Test d’arrêt numérique

[invite6236f8da]

Bonjour,
J'ai élaboré un code numérique ( en fortran90) pour approcher la valeur des températures dans une plaque 2D discrétisée par la méthode des volumes finis. Et vu que j'adopte la méthode ADI, j'ai utilisé un test d'arrêt ( l'erreur absolue ).
Le souci est que le test d'arrêt a l'air de ne pas bien marcher, en effet, quand je fais tourner le code, le test marche bien au départ, mais à un certain temps, il se stabilise à une valeur ( barrière que l'erreur ne peut pas franchir)
J'ai augmenté le nombre d'itérations mais en vain, je ne peux jamais atteindre la précision souhaitable.
sauriez vous quel est le problème ici ?
vous pouvez me recommander d'autres tests d'arrêt ?

Merci à vous.
-----

[SULREN]

Bonjour,
Assurez vous que toutes vos variables et paramètres ont bien été déclarés en flottant double précision, (voir mieux; je ne connais pas le Fortran90) et qu'il n'y en ait pas une déclarée par mégarde en variable logique, par exemple.

[invite6236f8da]

Bonjour,

Toutes les variables sont déclarées en real/int.
Je vais les déclarer en double et voir ce que ça va donner

[SULREN]

Re,
Dans le langage que j’utilise :
Integer est codé :
- sur 4 octets sur compilateur 32 bits et donc va de -2147483648 à +2147483647
- sur 8 octets sur compilateur 64 bits et donc va de 9223372036854775808 à +9223372036854775807
C’est insuffisant pour certains calculs où on cherche de la précision.

Le flottant double est aussi codé sur 8 octets, mais avec mantisse et exposant, c’est autre chose.

[A voir en vidéo sur Futura]

[invite6236f8da]

merci pour l'info, je viens de changer la nature des variables mais j'ai toujours le même problème.
apparemment, je dois changer le test d'arrêt

[SULREN]

Re,
Laissez quand même en Flottant Double ce qui mérite de l'être et en Integer ce qui doit l'être: indices des boucles, etc.

Je ne connais pas la méthode ADI: ne présente t'elle pas intrinsèquement, quelque part, une troncature qui fait qu'on ne peut pas descendre en dessous d'une certaine limite?

[invite6236f8da]

D'accord.
La méthode ADI ou encore la méthode des directions alternées implicites est un schéma implicite à deux pas.
L’idée derrière cette méthode est de scinder les équations aux différences finies en deux étapes pour trouver des systèmes linéaires dont les matrices sont tridiagonales ( pratiquement simples à résoudre )