Discussion :

[Yu Cloud]

Bonjours à tous,

Je suis en train d'élaborer un programme qui contient deux variables principales X et Y. J'utilise des équations qui relie X à Y et vis versa, ces variables sont calculées à chaque instant t. Dans ce programme je résous plusieurs équations dont des équations du second degré.

A chaque résolution, le résultat change à 10^-6 (ex : 10.000001 devient 10.0000015). Le problème c'est que je fais des itérations par rapport au temps, donc ce petit changement répété influence au bout d'un certain temps le résultat que j'obtiens.

J'aimerai donc savoir quels sont les moyens de réduire ces fluctuation pour stabiliser mes résultats ?

En vous remerciant par avance,

Yu Cloud

[rogue-spectre]

Bonjour,

Est-ce que tu travailles en simple ou en double précision ? Tu peux forcer la précision au niveau de la compilation avec les bons flags.

Avec gfortran l'option fdefault-double-8, plus d'infos :
http://linux.die.net/man/1/gfortran
Avec ifort l'option -r8, plus d'infos :
http://docs.notur.no/uit/files-uit/i...rt-9-x-options

Il y a également des options qui "améliorent" les calculs semble-t-il. L'option -mp1 pour ifort notamment. Après si tu fais des mélanges de simple et de double précision ça peut également jouer sur la précision des tes calculs. Une autre idée est de regarder si certaines opérations peuvent être réalisée de façon différentes

( 10E42 + 1 ) - 10E42 = 0

Bonne soirée.

[Yu Cloud]

En effet j'utilise des mélanges entre simple et double précision. Certain de ces mélanges me permet d'éliminer des erreurs en 19^-19 que j'obtenais.

[rogue-spectre]

Quand je parlais de mélange double-simple pour jouer sur la précision, c'était plutôt dans le sens négatif du terme, je te conseilles donc d'utiliser soit la simple, soit la double mais pas les deux en même temps.

[Yu Cloud]

Quand je parlais de mélange double-simple pour jouer sur la précision, c'était plutôt dans le sens négatif du terme

Je sais mais je n'ai pas trouvé d'autres solutions afin d'enlever les erreurs numériques en 10^-19 que j'obtenais. J'ai passé le reste de mes variables qui n'étaient pas touchées par l'erreur numérique en double précision. Par contre si tu as d'autres solutions à me proposer, je suis preneur !

De plus, je fais appelle à des fichiers extérieurs dont certains points varie légèrement (de l'ordre 10^-5). Et tout ceci rajoute des légère erreurs numériques.

D'où l'importance de pouvoir stabiliser le programme lorsque j'ai des légères variations. Donc tous conseils est le bienvenue

Bonne journée,

Yu Cloud

[Nanzilla]

Chaque opération de calcul introduit une erreur et les erreurs se multiplient entre elles. Plus tu as d'opérations, plus l'erreur augmente.
A partir de cela tu peux:

• Soit minimiser le nombre d'opérations pour diminuer l'erreur (optimisation méthodologique)
• Soit considérer qu'à chaque étape une partie de ta précision est superflue. Par exemple pour le nombre 10.00015000004 on pourrait ne garder que 10.00015 ou 10.0 (tout dépend de la précision voulue). Par cette approche on corrige manuellement les valeurs, mais il faut faire très attention aux données que l'on manipulent pour ne pas introduire un effet inverse.

En tout état de cause, il faudrait en savoir un peu plus sur ton programme pour pouvoire t'aider un peu plus.