Discussion :

[panda31]

Bonjour à tous,

je ne suis pas sûr que mon post est bien placé puisqu'il concerne deux langages : C & Fortran.

Bref, je possède un logiciel dont je possède aussi la source (cf AMBER). C'est intégralement écrit en Fortran. Dans un des modules, une fonction fait le calcul d'une matrice de dérivées secondes pour ensuite s'en servir dans d'autres calculs. Cette matrice n'est jamais renvoyée en sortie de fonction.

Jusque là, tou va bien dans mon explication ?

Mon problème est que je veux écrire cette matrice de dérivées secondes dans un fichier texte pour qu'un autre programme puisse s'en servir... Si AMBER avait été écrit en C, aucun problème... Mais c'est du fortran et, comment dire pour être convenable, je suis un buse en Fortran. J'ai fait des essais d'écriture externe dans un fichier mais cela ne veut même pas compiler !!!

Du coup, je viens au fait de suite, j'ai pensé (oui cela m'arrive des fois ) à faire un module en C et d'utiliser ses fonctions dans le programme Fortran. Je sais qu'il y a des pré-requis. Par exemple, il faut définir les types en fonction du langage Fortran comme suit :

Code :

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

#define int INTEGER

par exemple...

Encore, j'ai appris dernièrement que toutes les paramètres en Fortran se passaient par référence. Donc il faut se servir des pointeurs...Ca, ça va.

Mais y a-t-il ici des personnes ayant déjà eu ce genre d'interfaçage à faire et qui pourrait me conseiller SVP ?

Un immense merci à tous !

[DaZumba]

Cette matrice n'est jamais renvoyée en sortie de fonction.

Tu peux modifier la liste d'argument de la fonction pour qu'elle renvoie la matrice calculee.

J'ai fait des essais d'écriture externe dans un fichier mais cela ne veut même pas compiler !

Tu n'es pas tres perseverant. Ce n'est pas tellement dur de faire (Fortran 77)

Code :

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

1
2
3
4
5
6
7
8
 
      OPEN(unit=17, file='matrix.txt')
      DO i = 1, imax
         DO j = 1, jmax
           WRITE(17,*) matrix(i,j)
         ENDDO
      ENDDO
      CLOSE(17)

en supposant que c'est une matrice a deux dimensions.

Mais y a-t-il ici des personnes ayant déjà eu ce genre d'interfaçage à faire et qui pourrait me conseiller SVP ?

Je fais cela assez souvent - en supposant que l'on utilise des compilateurs C et Fortran faisant partie de la meme famille (genre g77 et gcc). Il y a pas mal d'info sur Google. Recherche 'mixing C and Fortran'. Remarque que ce n'est pas trivial, notamment lors du passage de tableau a plusieurs dimensions.

[panda31]

Tu peux modifier la liste d'argument de la fonction pour qu'elle renvoie la matrice calculee.


Tu n'es pas tres perseverant. Ce n'est pas tellement dur de faire (Fortran 77)

Code :

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

1
2
3
4
5
6
7
8
 
      OPEN(unit=17, file='matrix.txt')
      DO i = 1, imax
         DO j = 1, jmax
           WRITE(17,*) matrix(i,j)
         ENDDO
      ENDDO
      CLOSE(17)

en supposant que c'est une matrice a deux dimensions.


Je fais cela assez souvent - en supposant que l'on utilise des compilateurs C et Fortran faisant partie de la meme famille (genre g77 et gcc). Il y a pas mal d'info sur Google. Recherche 'mixing C and Fortran'. Remarque que ce n'est pas trivial, notamment lors du passage de tableau a plusieurs dimensions.

Si tu as vu mon code, tu as bien vu que ma syntaxe est la même que celle que tu me conseilles.


Pour ma persévérance, je travaille sur trois projets en ce moment ce qui me force à faire des choix !!!

En tout cas, ça compile toujours pas...J'enlève mon code, ça compile...

ou alors, je compile mal... Mais vu que j'ai rien changé d'autre, j'ai laissé les options telles quelles...

[DaZumba]

Si tu as vu mon code, tu as bien vu que ma syntaxe est la même que celle que tu me conseilles.

J'ai vu ton post dans le forum Fortran, mais comme tu ne te plains de rien, je n'ai rien a repondre. Comme d'habitude, 'ca compile pas' est une declaration vague qui ne permet pas de faire grand'chose. Il faudrait au moins me dire ce qui ne compile pas, quels sont les messages donnes par le compilateur, et si c'est le code C ou Fortran qui pose probleme...

[panda31]

J'ai vu ton post dans le forum Fortran, mais comme tu ne te plains de rien, je n'ai rien a repondre. Comme d'habitude, 'ca compile pas' est une declaration vague qui ne permet pas de faire grand'chose. Il faudrait au moins me dire ce qui ne compile pas, quels sont les messages donnes par le compilateur, et si c'est le code C ou Fortran qui pose probleme...

1°) Pour le moment, je n'ai pas interfacé C et Fortran
2°) Je ne serais au boulot que mardi ! Du coup, je ne pourrais te donner le message d'erreur précis que mardi soir !!!