Discussion :

[bored]

Bonjour !

Je poste ici car je ne comprends pas bien un programme Fortran que l'on m'a fait (et qui ne fonctionne pas ). Je précise que je suis complètement novice en la matière, je n'ai eu que 2-3 séances d'initiation à Fortran pour le moment.

Il s'agit de la programmation de la méthode de Gauss pour les matrices tridiagonales.
La ligne qui me pose problème est la première du corps du programme, après les définitions :

IF(DABS(B(1)).LT.1.D-7) WRITE(*,1000)
1000 FORMAT(' LE SYSTEME EST NON INVERSIBLE')

Où B est le terme de la diagonale principale de la matrice.
Je ne comprends pas 'DABS', ni '1.D-7', je sais que LT est Lower Than mais je ne comprends pas les termes qui l'entourent, je n'ai aucun D de défini nulle part. Par contre au début du programme j'ai 'IMPLICIT DOUBLE PRECISION (A-H,O-Z)' : est-ce que cela définit implicitement un D ?

Voilà si quelqu'un pouvait m'éclairer, ça serait chouette

Merci d'avance !

[FR119492]

Salut!
Commence par lire le cours de P. Corde et A. Fouilloux qui est disponible sur ce site.
Jean-Marc Blanc

[ZeNoob]

C'est juste :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

|B(1)| < 1.0E-7 = 0.0000001

dabs est la valeur absolue pour les variables de type double precision
et le d-7 c'est la façon de faire les exposants en fortran.

[Sylvain Bergeron]

...
IF(DABS(B(1)).LT.1.D-7) WRITE(*,1000)
1000 FORMAT(' LE SYSTEME EST NON INVERSIBLE')

...
Je ne comprends pas 'DABS', ni '1.D-7', ... j'ai 'IMPLICIT DOUBLE PRECISION (A-H,O-Z)' : est-ce que cela définit implicitement un D ?

Je seconde Jean-Marc.

C'est un très ancien programme, ou à tout le moins, un programme écrit dans un très ancien style. Ça fait au moins 30 ans qu'on sait que le typage implicite (implicit double ...), c'est de la mauvaise programmation.

Comme ZeNoob a dit, DABS est la fonction valeur absolue pour un double. Techniquement, c'est la fonction explicite substituée à la fonction générique abs() lorsqu'appliquée à un double. Sauf dans le cas très précis du passage de fonction à un sous-programme (argument dummy external), on devrait toujours utiliser les fonctions génériques, et jamais les fonctions explicites. Le DABS laisse supposer que ton programme est écrit en Fortran IV ou en Fortran 1966 (les fonctions génériques existent depuis Fortran 77). Ton if devrait donc se lire IF(ABS(B(...

1D-7 : C'est de la notation scientifique, comme sur une calculette, mais pour de la double précision. Ça correspond à 1e-7 en simple précision, soit 1 x 10^-7