Discussion :

[PadawanDuDelphi]

Bonjour à tous,

Voici mon problème, j'ai repris un ancien programme que je suis censé débogger et à l'intérieur se trouvent un tas de structures du genre:

Code :

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

1
2
3
4
5
6
7
8
9
10
TYPE PARALLELE
    integer::indiceParallele
    complex(8)::Zparallele
    complex(8),dimension(:,:),allocatable::YP
END TYPE PARALLELE
 
TYPE PARALLELE_LIST
    integer::nbParallele
    type(PARALLELE),dimension(:),allocatable::TParallele
END TYPE PARALLELE_LIST

que j'alloue dans mon programme initiale en fonction de paramètres issus d'un base de donnée.

Et j'aimerais simplement savoir si un deallocate(TParallele) suffit à tous supprimer ou si je suis obligé de désallouer d'abord tous les YP puis les TParallele.

A noter que dans le programme c'est la première solution qui est utilisée, alors que j'opterais pour la seconde (mais avant de modifier un grand nombre de lignes de codes, je préfère demander l'avis d'experts)

Merci beaucoup pour votre aide,

@+.

[DaZumba]

Et j'aimerais simplement savoir si un deallocate(TParallele) suffit à tous supprimer ou si je suis obligé de désallouer d'abord tous les YP puis les TParallele.

C'est la seconde solution. Typiquement, on désalloue dans l'ordre inverse des allocations.

EDIT: Il semble en fait que j'ai tort. Voir la suite du fil.

[PadawanDuDelphi]

Ok,

merci beaucoup pour ton aide, il me semblait aussi qu'on avait vu ça en cours en C. Mais je voulais être sûr avant de me lancer dans la modification des 4000 lignes de codes...

Encore merci,

@+.

[Sylvain Bergeron]

Tu n'as pas à désallouer les membres allocatables. C'est une différence importante et un grand avantage des membres allocatables sur les membres pointers. Ta première solution est donc suffisante.

[DaZumba]

Tu n'as pas à désallouer les membres allocatables. C'est une différence importante et un grand avantage des membres allocatables sur les membres pointers. Ta première solution est donc suffisante.

Interessant. Tu as une source pour cela ?

[Sylvain Bergeron]

J'ai malheureusement une journée de fou. J'ai cherché très rapidement sous google group comp.lang.fortran. Je n'ai pas trouvé explicitement ce que tu demande, mais une phrase qui reflète l'étendu de la cascade...

It's also worth mentioning that the TR support in the particular case of the Lahey compiler is an "unsupported" feature. While most of it works very well, there is one particular area which is broken: the automatic deallocation, on subroutine entry, of an INTENT(out) allocatable dummy argument. I submitted this bug to Lahey, but since the TR capabilities are unsupported, they filed my bug report as a product suggestion (d'oh!).

La discussion complète : http://groups.google.com/group/comp....7fbd71d24113f4

(Si l'anglais cause problème) Ça dit que si tu passe une variable (type) contenant un membre allocatable à une fonction dont l'argument est intent out, le membre sera désalloué automatiquement à l'entrée.

Si tu veux chercher, la piste est probablement Fortran 2003 ou TR 15581 (extension standard ISO à F95 intégré dans F2003)

[Mathlus]

Je ne sais pas si ça peut faire avancer les choses mais je travaille en ce moment avec des tableaux (allocatables) de structures contenant des variables "allocatable" et je suis plutôt de l'avis de Sylvain Bergeron. Je viens même de rencontrer un cas où la désallocation d'un tableau contenu dans un structure (type) qui est elle-même l'élèment d'un tableau "allocatable" causait des erreurs à l'execution. Alors qu'en désallouant uniquement le tableau de structures contenant tout le reste, aucun problème ! Toutes les variables se réinitialisent comme il faut à la réallocation. (Mon programme tourne en boucle sur plusieurs fichiers et je réalloue plusieurs fois d'affilée les tableaux de structures).

Voilou pour mon expérience perso... ça vaut ce que ça vaut.

@ +,

Mathlus.

[DaZumba]

Je ne sais pas si ça peut faire avancer les choses mais je travaille en ce moment avec des tableaux (allocatables) de structures contenant des variables "allocatable" et je suis plutôt de l'avis de Sylvain Bergeron.

Ok, j'ai regarde les quelques rares informations disponibles la-dessus sur l'Internet, et il semble que Sylvain ait raison, en effet. My mistake! J'edite mon premier post. Par contre, je suis surpris que desallouer les allocatable a la main induise un comportement different, comme tu l'ecris, Mathlus. La gestion de la memoire en Fortran90 me semble assez etrange, tout d'un coup...

[Sylvain Bergeron]

Ça s'est un peu calmé, alors j'ai repris ma recherche.

Voici le document TR 15581 : ftp://ftp.nag.co.uk/sc22wg5/N1351-N1400/N1379.pdf

Première phrase en haut de la page 5 :

Deallocation of a variable containing an ultimate allocatable array component automatically deallocates all such components of the variable that are currently allocated.

[PadawanDuDelphi]

OK,

Merci beaucoup pour votre aide (surtout que ca va m'eviter de modifier ce long code)...

Dernière petite question: Si on effectue une desallocation des elements internes avant de supprimer l'element principal (solution 1), il y a aucune raison pour que ça plante, non?

@+.

[Sylvain Bergeron]

Je ne vois pas pourquoi ça ne marcherait pas, mais Mathlus semble avoir une expérience contraire. Faudrait voir son code pour juger...

[PadawanDuDelphi]

De toute façon ça parait plus logique et plus facile d'opter pour la désallocation d'un élément plutôt que de désallouer l'intérieur des stucture auparavant...

J'avais juste continuer le programme avec la methode numérp 2 et je voulais savoir si je devais le modifier ou pas, mais je pense que je vais le faire pour uniformiser.

Merci pour tout,

@+.

[Mathlus]

En fait, la déallocation qui m'a posé problème est :

Code :

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

Deallocate(MyMesureAcous(i)%Response%LocalColumnX%Values)

où MyMesureAcous est un tableau allocatable et Values aussi.
Pas d'erreur à la compilation mais un "access violation" à l'execution. J'ai donc opté pour la solution "je déalloue juste MyMesureAcous" et là pas de problème...

J'ai aussi eu une erreur (à la compilation cette fois, du style "vous ne pouvez pas déallouer un tableau qui n'est pas allocatable") avec :

Code :

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

Deallocate(MyChampAcous%Response%LocalColumnX%Values)

où comme dans l'exemple précédent, Values est un tableau allocatable (alloué et utilisé pour des calculs pour tôt dans le code) et MyChampAcous un type simple (pas de tableau).
Du coup, j'ai laissé tomber le deallocate et j'ai réussi à réallouer le tableau (Allocate(MyChampAcous%...%Values)) sans problème à l'itération suivante, avec des valeurs tout à fait cohérentes et un tableau de taille diffénrete à l'itération précédente.

Si vous pouvez m'expliquer tout ça, c'est pas de refus !!

@ +,

Mathlus.

[PadawanDuDelphi]

Salut,

J'ai déjà un message de ce genre, mais c'était parcque mon tableau n'était pas alloué précédemment...
Sinon je vois pas...

@+.

[Mathlus]

Je te confirme que mon tableau était bel et bien alloué au moment où j'essayais de le désallouer. Je ne comprends pas non plus...

[Sylvain Bergeron]

Mathlus :

Il est difficile de comprendre ce qui ne va pas avec aussi peu d'information. Si tu as besoin de la réponse au problème, tu devras nous fournir soit le code source complet, soit une version simplifiée qui reproduit le problème.

As-tu compilé en mode « debug » ?

Il y a aussi une autre possibilité pour que tu trouve par toi même : essaie avec un autre compilateur. Il est assez facile d'en avoir un gratuit pour une période d'essai, ce qui est suffisant pour un tel besoin. Un compilateur différent produit souvent un diagnostic différent, ce qui permet d'isoler le problème.