Discussion :

[Guillaume_C]

Bonjour,

Voici mon programme en Fortran 2003:

Code :

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

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

[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

[Guillaume_C]

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 ?

[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.

[Guillaume_C]

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é Pour le compilateur intel, il faut compiler comme ceci, et ca marche:

ifort model.f90 -assume realloc_lhs