Module et tableau dynamique

Version imprimable

Salut à tous,

Avant de me lancer dans l'écriture d'une longue subroutine, et n'ayant jamais touché aux modules avec fortran95, j'ai fait un petit programme test afin de me familiariser avec ce nouvel outil.
Mais, au malheur, l'allocation dynamique d'un tableau dans un module, n'aboutit qu'à recevoir les insultes du compilateur...

Je laisse le module en copie, si quelqu'un a une idée, je le remercie par avance de sauver l'amitié que j'avais lié avec mon compilateur.
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 MODULE conf_ini IMPLICIT NONE REAL (KIND=8), DIMENSION(:), ALLOCATABLE :: atom INTEGER (KIND=4) :: ntot INTEGER (KIND=4) :: x,compt CONTAINS SUBROUTINE CPA(atom,ntot) REAL (KIND=8), DIMENSION(:), ALLOCATABLE :: atom INTEGER (KIND=4) :: ntot INTEGER (KIND=4) :: x,compt !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ntot=9 allocate(atom(ntot)) compt=0 DO x=1,ntot compt=compt+1 atom(x)=real(compt) ENDDO deallocate(atom) END SUBROUTINE CPA END MODULE conf_ini

02/09/2007, 20h23
FR119492

Module et tableau dynamique

Salut.

Parfois, dans la vie, il faut faire des choix: si tu veux "faire moderne", utilise l'allocation dynamique, mais ne t'étonne pas si tes programmes trainent les pieds.

Tu peux aussi imiter le vieux croulant que je suis et te restreindre à une allocation strictement statique. C'est souvent beaucoup plus efficace. C'est une règle tout à fait générale, comme tu le verras par exemple dans la documentation de MatLab.

A toi de choisir.
Jean-Marc Blanc
04/09/2007, 12h47
yann_m

Salut,

Merci pour l'info, c'est un problème auquel je dois effectivement me préocuper puisque le code en question est un Monte Carlo et autant le faire converger le plus rapidement possible.

Cependant mon problème initial demeure, comment remplir un tableau dans un module2, dont la taille est fixé dans un premier (module1) ?

Est-ce ainsi qu'il faut s'y prendre ?

MODULE1 : Je détermine la taille de mon tableau (par exemple l'entier ntot)

MODULE2 : Je déclare mon tableau
SUBROUTINE(montableau)
REAL KIND(8), DIMENSION(ntot) :: montableau

PROGRAM : je déclare un tableau que j'alloue en dynamique ou autre ? j'appel le module1 puis le module2 ?
Il y a quelque chose à cette étape qui coince et dans lequel je me fourvoie...

Merci encore ton message.
04/09/2007, 16h32
FR119492

Module et tableau dynamique

Re-salut!

Juste une petite suggestion: va sur le site www.netlib.org et cherche la librairie LinPack. Ensuite, tu ouvres un sous-programme, par exemple DGEFA et tu regardes comment ils font. C'est peut-être vieillot, mais ça marche, et plutôt bien.

Bon succès
Jean-Marc Blanc

Bonjour,

Pour info, ton module compile très bien chez moi (avec g95) et le code test suivant:
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 MODULE conf_ini IMPLICIT NONE REAL (KIND=8), DIMENSION(:), ALLOCATABLE :: atom INTEGER (KIND=4) :: ntot INTEGER (KIND=4) :: x,compt CONTAINS SUBROUTINE CPA(atom,ntot) REAL (KIND=8), DIMENSION(:), ALLOCATABLE :: atom INTEGER (KIND=4) :: ntot INTEGER (KIND=4) :: x,compt !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ntot=9 allocate(atom(ntot)) compt=0 DO x=1,ntot compt=compt+1 atom(x)=real(compt) write(*,*)'x=',x,'atom(x)=',atom(x) ENDDO deallocate(atom) END SUBROUTINE CPA END MODULE conf_ini !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! program test_module use conf_ini call CPA(atom,ntot) end program test_module
affiche bien le contenu du tableau comme on s'y attend.

Par contre, vu que le tableau est désalloué dans la foulée, s'est assez contre-productif...

Vu tes messages, le but que tu recherche est de partager l'info entre plusieurs routines et le programme principal; les modules sont effectivement ce qu'il te faut mais il faut bien voir que c'est le fait d'inclure un module ici et là qui mène au partage de l'information.
Peut-être que l'exemple suivant sera plus parlant qu'un long discours:
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 module common_stuff implicit none ! Variables du module: integer :: taille real(kind=8),dimension(:),allocatable :: tableau contains subroutine ini_tableau implicit none ! variable locale integer :: i write(*,*) 'Taille du tableau?' read(*,*) taille ! Allocation du tableau + initialisation des elements allocate(tableau(taille)) do i=1,taille tableau(i)=i enddo end subroutine ini_tableau end module common_stuff !!!!!!!!!!!!!!!!!!!!!!!! program share_data use common_stuff implicit none ! variables locales integer :: i ! initialisation call ini_tableau do i=1,taille write(*,'(a9,i2,a2,g15.5)')'tableau(',i,')=',tableau(i) enddo ! NB: 'taille' et 'tableau sont connu grace au 'use common_stuff' end program share_data
La subroutine 'ini_tableau' connait 'taille' et 'tableau' à cause du contains; de même pour le programme principal grâce au use common_stuff.

Bonne continuation.

05/09/2007, 16h54
Sylvain Bergeron

Je seconde Ehouarn pour les modifications proposées. La routine ini_tableau pourrait cependant avoir un paramètre de taille en paramètre. La saisie de la taille requise pourrait alors être déplacée ailleurs (dans ton Module1).

De l'externe, tu dois uniquement fournir la taille du tableau. Il me semble illogique de passer en paramètre le tableau à dimensionner puisqu'il fait partie du module.
07/09/2007, 00h13
yann_m

Bonjour,

merci à vous, le programme proposé fonctionne bien. L'erreur venait du fait que le tableau était alloué et désalloué dans un même module, puis appelé par le programme principal. L'erreur alors retourné était un message de segmentation... Logique. Bon, vos messages et ces quelques jours de galères m'auront au moins permi de me familiariser avec les modules et d'en arriver à la conclusion que ce tableau, je dois le déclarer dans le programme principal.

Merci beaucoup pour vos messages, à un de ces jours dans les méandres du forum.