[Fortran 2003] Erreurs d'exécution

Version imprimable

Bonjour,

Voici mon programme en Fortran 2003:
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 !####################################################################### program model real, dimension(:), allocatable :: arr1 real, dimension(8) :: arr2 = (/ 1., 2., 3., 4., 5., 6., 7., 8. /) real, dimension(3) :: arr3 = (/1., -2., 3. /) arr1 = 2. * arr2 print *, arr1 end program model !#######################################################################
C'est un exemple de la page 401 du livre Chapman, S. 2007. Fortran 95/2003 for Scientists And Engineers, McGraw-Hill.

Je le compile avec l'Intel Fortran Compiler:
ifort model.f90

La compilation s'effectue correctement, mais l'execution du programme affiche une ligne vide, et pas le contenu de arr1. Si je remplace par print *, arr2, l'affichage est correct.
Je ne comprends pas ce que je fais de faux.

J'ai essayé avec gfortran, en compilant avec
gfortran model.f90

et le programme donne une 'Segmentation fault'.
Je ne comprends pas non plus ce que je fais de faux.

Merci de votre aide, je débute en Fortran...:)

Guillaume

02/01/2009, 09h56
Ehouarn

Bonjour,

le tableau 'arr1' est déclaré comme 'allouable', mais n'est jamais alloué !!

D'où les erreurs que tu obtiens à l'exécution ...

Bonne continuation
02/01/2009, 11h34
Guillaume_C

Citation:

Envoyé par Ehouarn

le tableau 'arr1' est déclaré comme 'allouable', mais n'est jamais alloué !!

Je croyais qu'en Fortran 2003, on pouvait allouer un tableau a partir d'un autre (des lors qu'ils ont le same rank), sans explicitement ecrire l'etape d'allocation, qui est comprise dans arr1 = 2.*arr2

Ce n'est donc pas le cas ?
02/01/2009, 12h16
Ehouarn

Je ne sais pas (je ne suis pas très au fait des ajouts du F2003 par rapport au 90/95).
Si le F2003 permet effectivement cela, il semble bien que tes compilateurs ne le prennent pas en charge...
Une possible explication: ton programme se nomme 'model.f90' et l'extension 'f90' indique peut-être aux compilateurs qu'il s'agit de Fortran 90 et pas 2003.
Je te conseille en tous cas de compiler avec le plus d'options de débogage possible (par ex. 'gfortran -O2 -Wall -Wextra') pour voir ce qu'en disent tes compilateurs.
02/01/2009, 14h43
Guillaume_C

Citation:

Envoyé par Ehouarn

Je ne sais pas (je ne suis pas très au fait des ajouts du F2003 par rapport au 90/95).
Si le F2003 permet effectivement cela, il semble bien que tes compilateurs ne le prennent pas en charge...
Une possible explication: ton programme se nomme 'model.f90' et l'extension 'f90' indique peut-être aux compilateurs qu'il s'agit de Fortran 90 et pas 2003.
Je te conseille en tous cas de compiler avec le plus d'options de débogage possible (par ex. 'gfortran -O2 -Wall -Wextra') pour voir ce qu'en disent tes compilateurs.

J'ai trouvé :D Pour le compilateur intel, il faut compiler comme ceci, et ca marche:

ifort model.f90 -assume realloc_lhs