Discussion :

[Matthieu Brucher]

Bonjour,

Ma question est aussi simple que ça : comment transformer un pointeur récupéré par un appel extérieur comme un tableau à 2 dimensions en Fortran 90. Aussi simple que :

Code :

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

real, intent(in)    :: tab(dim1,dim2)

?

[Sylvain Bergeron]

Solution #1 - Cray Pointer (non standard mais assez disponible)

Je l'utilise pour récupérer le data sous forme de tableau 1D. Je n'ai jamais essayé avec un tableau 2D.

Par exemple : real rData(*); pointer(pData, rData)

--> Cherche le web pour "Cray Pointer Fortran"

Solution #2 - Cascade en Fortran 77

Il faut passer le pointeur en "value" à un sous-programme qui s'attend à recevoir un tableau Fortran (il s'attend en fait à recevoir un pointeur sur l'origine du tableau). Pour que ça fonctionne, il ne faut évidemment pas que l'appelant voit l'interface du sous-programme; on est donc en Fortran 77...

Solution #3 - Fortran 2003 - Interface to C

Une solution très standard existe très certainement, ton besoin est trop commun. Je ne sais malheureusement pas comment !

[Sylvain Bergeron]

Ajout à la solution #3

Vas voir le document ftp://ftp-developpez.com/fortran/cours/Fortran-2003.pdf, page 21, premier point en haut.

[Matthieu Brucher]

Merci, je vais regarder tout ça !
La solution 1, j'évite maintenant car on a eu pas mal de soucis avec les pointeurs Cray, et notamment tout la partie check bound qui ne fonctionne pas dans ce cas.
Mon gros problème, c'est que le code d'interface est en Fortran 90 (mais sans module). Je sais que c'est très courant et un des gars avec qui je travaille est halluciné que ce n'est plus vraiment possible en Fortran 90 simplement. Avant, il y avait la possibilité de changer le type avec je ne sais plus quelle instruction, mais elle est dépréciée, et dans les bouquins de Fortran, je ne trouve plus rien !

[Sylvain Bergeron]

Mon gros problème, c'est que le code d'interface est en Fortran 90 (mais sans module). Je sais que c'est très courant et un des gars avec qui je travaille est halluciné que ce n'est plus vraiment possible en Fortran 90 simplement.

Qu'est-ce qui n'est plus vraiment possible ?

Avant, il y avait la possibilité de changer le type avec je ne sais plus quelle instruction, mais elle est dépréciée, et dans les bouquins de Fortran, je ne trouve plus rien !

Equivalence ?

[Matthieu Brucher]

Faire simplement des changements de types, genre tu alloues un gros tableau de flottants de travail, puis tu utilises un morceaux pour en faire des complexes, un autre pour des entiers finalement...
On peut faire ça ?

Oui, c'était equivalence.

J'ai finalement utilisé ISO_C_BINDING, le C_F_POINTER devrait faire plus que l'affaire, d'autant que j'ai vu qu'on pouvait exporter la fonction en C de manière transparente, c'est vraiment bien (OK, c'est vieux, mais tout de même !)

En tout cas, encore merci pour tes pistes, elles ont été précieuses.

[Sylvain Bergeron]

Faire simplement des changements de types, genre tu alloues un gros tableau de flottants de travail, puis tu utilises un morceaux pour en faire des complexes, un autre pour des entiers finalement...
On peut faire ça ?

Traditionnellement, on faisait ça par des allocations statiques et des equivalences. Je ne pense pas qu'on puisse le faire avec des allocations dynamiques de façon propre.

[Matthieu Brucher]

Bien ce que je me disais