Discussion :

[info!]

bonsoir,
j'utilise en ce moment le fortran 77 de cygwin sous windows pour calculer des matrices de triangles numériques d'entiers positifs (exple : triangle de pascal) avec des lois de récurrence du genre :
A(i+1,j+1)=A(i,j) + (j+1)A(i,j+1)
Tout semble marcher mais je suis étonné lorsque j'affiche ces triangles de voir que seules les 10 ou 12 premières lignes sont correctes et qu'ensuite apparaissent des entiers négatifs, ce qui est mathématiquement impossible, et les lignes suivantes sont complètement fausses.
Il y a t-il déjà une instabilité numérique dans fortran même après un calcul de triangles aussi simple ??
merci beaucoup !

[info!]

bonsoir,
dernièrement, je crois avoir compris qu'il s'agissait d'un problème de précision en termes de nombre d'octets utilisés :
en effet quand je déclare mes variables entieres comme integer*8, j'obtiens plus de lignes exactes (et positives!) qu'avec integer*1.
Donc il semble que la taille des entiers dans mes matrices jouent beaucoup, or j'ai vraiment besoin de calculer de très grand entiers, comment faire alors pour augmenter la précision avec g77 ??
merci beaucoup

[Ladgalen]

Ce n'est pas un problème de fortran, c'est un problème général de mémoire. Un nombre (qu'il soit entier ou réel) est stocké sur un nombre fini de d'octet (8bits).
Il en découle obligatoirement un limite du nombre maximum que tu peux enregistrer.

Avec integer*1 tu enregistre tes entiers sur 1 octet soit 8bit. Le premier étant un bit de signe il te reste 7 bit pour enregistrer ton entier, le nombre max que tu peux enregistrer est donc 256 ! D'où le fait que tu ne peux pas aller très loin dans le tableau de pascal !

regarde ce post pour plus de détails :

http://www.developpez.net/forums/d62...ble-precision/

[info!]

merci Lagdalen,
j'ai aussi vu ton lien, ma question désespérée est alors : peut-on augmenter de beaucoup la capacité de mémoire des entiers de fortran 77 ( genre integer*32) ou alors... changer carrément de langage ??
merci !

[Ladgalen]

Je ne suis pas le plus qualifié pour répondre à cette question. Cependant il me semble qu'en fortran 90 (ou peut être 95) tu peux donner le nombre d'octet sur lequel tu enregistres tes nombres sans contraintes.

Sinon je crois que c'est un des atouts de C ... mais là j'y connais encore moins.

Sinon au passage si tu te mets en réel double précision, tu as un éventail de valeurs entre 10^-300 et 10^+300 ce qui te laisse de la marge !!!!

[info!]

merci encore Lagdalen,
je ne suis pas sûr que des fonctions comme mod(a,b) affichent parfaitement zéro en double précision : je vais réessayer.
Par contre, puis-je trouver le fortran 90 ou 95 dans les mêmes packages cygwin où j'ai extrait le f77 ??
merci, si d'autres ont plus de précisions...