Mélange Fortran et C

Version imprimable

13/04/2007, 11h08
vinzzzz

Mélange Fortran et C

Hello,
Je précise davance, je connais pas le fortran (enfin je découvre depuis hiers)...
J'ai donc récupéré un programme fortran relativement long. Dedan il y a une fonction qui m'interesse et que je souhaiterai pouvoir appeler a partir d'un code C.
Première question:
Sur ce lien:
http://www.nersc.gov/nusers/resource...m/c_and_f.html
On donne des exemples d'appel de C a partir du fortran et vice versa. L'édition des liens se fait avec soit un compilateur C, soit un fortran, selon le cas. Or moi je souhaite pouvoir a la fois appeler des routines fortran a partir du C, et des fonctions C a partir du fortran le tout dans un même programme. Donc première question, est-ce possible, et puis-je utiliser indifférement l'un des deux compilateurs? (ce qui me parait fort peu probable) Et si non, comment faire?

Ensuite, la routine principale fortran que je souhaite appeler prends en paramètres des tableaux de double, matrices, des structures (ou l'équivalent fortran...), des fonctions (enfin ds pointeurs de fonction je suppose) et des arguments optionnels...
Là dessu, sur les structures, pointeurs de fonction et stt sur les arguments optionnels, je souhaiterai savoir si c'est possible d'appeler ce genre de fonction à partir d'un programme C (pour les détails techniques je chercherai plus tard), histoire que je ne perde pas mon temps a essayer de faire quelque chose d'impossible...

Voila (pour l'instant ;) )
Merci d'avance ;)
13/04/2007, 13h51
poulacou

1. Il est tout a fait possible de mélanger le C et le fortran !

2. génralement dans le source C il faut rajouter un _ (underscore) dans le nom des subroutine fortran que tu appeles.

3. Je te conseille de compiler chaque fichier avec le compilateur associé au language et d'utiliser pour le link l'editeur de lien fortran en rajoutant un option pour indiquer que le main est en C ( exemple pour ifort => -norformain). sinon essayer avec le linker C.
explication :
pour l'instant un main en C et fonction fortran est possible mais plus pour longtemps.
pour s'execute le C n'a pas besoin de runtime, le fortran non plus mais le C++ oui.
Le Fortran aura bientot besoin d'un runtime pour cause des évolutions objets et pour lancer ce runtime le main devra obligatoirement etre dans un fichier fortran.
Donc creation d'un main fortran qui appelle l'ancien main du C.

Tiens moi au courant de tes tests
@+

Pour être plus précis, du programme C j'appèlerai normalement uniquement une fonction qui retourne une structure contenant pas mal de chose.
A cette fonction fortran, je donnerais un pointeur sur fonction qui devrait appeler une fonction qui elle seraié écrite en C... Voila les deux intéractions principales.

En me fiant aux exemple que j'ai vu dans le lien indiqué, je fais un test bidon déja pr voir si j'arrive a faire les appels. Je souligne que je dois utiliser des bibliotheques de GSL et quelques fichiers minimaux, dc je les inclus dirèctement pr mes premiers tests. J'en suis donc là:
C:

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
 
extern void __dde_NMOD_essai_(double*, double*, double *, int *, int *, double (*fct)(double*)) ;
extern void essai_(double*, double*, double *, int *, int *, double (*fct)(double*)) ;
extern void essai(double*, double*, double *, int *, int *, double (*fct)(double*)) ;
 
double fct(double *x) {
	return (*x)*(*x) ;
}
 
int main(int argc, char **argv){
	printf("\n\n- Program %s begins- \n\n", argv[0]);
	double t ;
	double y[] = {1.0, 2.0, 3.0} ;
	double dy[] = {0., 0., 0.};
	int neqn = 3. ;
	int nlag = 3. ;
 
	__dde_NMOD_essai_(&t, y, dy, &neqn, &nlag, &fct) ;	
	essai(&t, y, dy, &neqn, &nlag, &fct) ;
	essai_(&t, y, dy, &neqn, &nlag, &fct) ;
 
	printf("\n\n- End of %s - \n\n", argv[0]);
 
	return 0 ;
}

Fortran:

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
 
MODULE dde
  IMPLICIT NONE
  CONTAINS
 
  SUBROUTINE ESSAI(T,Y,DY,NEQN,NLAGS,FCT)
      INTERFACE
	  DOUBLE PRECISION FUNCTION FCT(T)  
              DOUBLE PRECISION, INTENT(IN) :: T
          END FUNCTION FCT
      END INTERFACE
 
      ! DOUBLE PRECISION, EXTERNAL :: FCT
      INTEGER :: I, NEQN, NLAGS
      DOUBLE PRECISION :: T
      DOUBLE PRECISION, DIMENSION(NEQN) :: Y,DY
 
      INTENT(IN)  :: T,Y,NEQN,NLAGS
      INTENT(OUT) :: DY
 
      DY(1) = FCT(T+1.0D0)
 
      RETURN
  END SUBROUTINE ESSAI
END MODULE dde

Makefile:

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
44
 
#------------------------------------------------------------
# Makefile. Auteurs: Vince
#------------------------------------------------------------
 
PATH_GSL    = /home/vince/gsl-1.9/
PATH_OBJ    = obj/
PATH_SCR    = src/
PATH_HEADER = headers/
PATH_BIN    = bin/
 
PROG        = dde
CC          = gcc
FCC         = gfortran
 
CFLAGS      = -Wall -Wextra -Wwrite-strings -Wstrict-prototypes -D MD_OS_LINUX -I$(PATH_GSL)include
LFLAGS	    = -L$(PATH_GSL)lib -lm -lgsl -lgslcblas 
 
OBJ = $(PATH_OBJ)main.o $(PATH_OBJ)utils.o $(PATH_OBJ)dde.o $(PATH_OBJ)essaI.o
 
#------------------------------------------------------------
# REGLE GENERALE DE COMPILATION
#------------------------------------------------------------
 
$(PATH_OBJ)%.o: $(PATH_SCR)%.c $(PATH_HEADER)%.h
	$(CC) $(CFLAGS) -c $< -o $@
 
$(PATH_OBJ)%.o: $(PATH_SCR)%.f90
	$(FCC) -c $< -o $@
 
#-----------------------------------------------------------
# REGLE POUR L'EXECUTABLE
#-----------------------------------------------------------
 
$(PATH_BIN)$(PROG): $(OBJ)
	$(CC) $(LFLAGS) $^ -o $@
 
#-----------------------------------------------------------
# DEPENDANCES
#-----------------------------------------------------------
 
$(PATH_OBJ)main.o: $(PATH_HEADER)main.h $(PATH_HEADER)utils.h $(PATH_HEADER)dde.h
$(PATH_OBJ)dde.o: $(PATH_HEADER)dde.h $(PATH_HEADER)utils.h
$(PATH_OBJ)utils.o: $(PATH_HEADER)utils.h

A l'édition des liens:

Citation:

gcc -L/home/vince/gsl-1.9/lib -lm -lgsl -lgslcblas obj/main.o obj/utils.o obj/dde.o obj/essaI.o -o bin/dde
obj/main.o: In function `main':main.c:(.text+0x99): undefined reference to `__dde_NMOD_essai_'
:main.c:(.text+0xc8): undefined reference to `essai'
:main.c:(.text+0xf7): undefined reference to `essai_'
collect2: ld returned 1 exit status

Bref ca marche pas ^^ aucun des prototypes ne passe, je continue mes recherches, tout commentaire est le bienvenu..

13/04/2007, 15h15
poulacou

j'ai réutilisé a compilé ton exemple avec quelques modifications :

1. j'ai utilisé ifort a la place de gfortran
2. pour gcc j'ai forcé la compilation en 32 bits -m32 ; pas besoin pour le fortran
3. dans dde.f90 j'ai supprimé le module et est declaré directement une subroutine dde
4. dans main j'appelle dde_
5. les pointeurs sur fonctions pas sur que ca marche.

Je peux pas plus t'aider
mais tiens moi au courant de tes investigations et conclusion sur l'utilisations des modules et des pointeurs sur fonctions.

A+
Fabrice
13/04/2007, 15h32
vinzzzz
Ok ca marche de mon coté, en revanche j'ai juste retiré les histoires de modules... Ce qui est assez embettant car le programme que j'ai récupéré est quelque peu complexe et la fonction que je veux appeler est déclarée ainsi:
Code:

1 2 3 4 5 ! .. Generic Interface Blocks .. INTERFACE DDE_SOLVER MODULE PROCEDURE DKL_1, DKL_2, DKL_3, DKL_4, ODEAVG END INTERFACE
Ce qui me trouble c'est que j'ai pas trouvé de déclaration de type FUNCTION DDE_SOLVER ou SUBROUTINE DDE_SOLVER dans le fichier...

En tout cas pparamaent l'option -m32 n'est pas necessaire (du moins pr l'exemple présent) et le copilateur que j'utilise a la base fonctionne également.

Pour le pointeur de fonction, ca semble fonctionner (je met un 'coucou' dans la fonction fct et il apparait bien). Je vas continuer histoire de voir que les valeurs collent bien etc... et je passerait aux structure & aux arguments optionnels je sens que ca va etre prise de tête...

En tt cas merci a toi ;)
13/04/2007, 15h41
poulacou

pour l'histoire des modules , tu as la réponse ici : http://www.developpez.net/forums/sho...d.php?t=311098
13/04/2007, 15h49
vinzzzz

Ca passe pas, c'est peut être spécifique au c++...
EDIT: Enfin je dit que ca passe pas ms j'ai pas encore testé la première méthode ;)
18/04/2007, 16h49
Sylvain Bergeron
Citation:
Envoyé par vinzzzz

Code:

1 2 3 4 5 ! .. Generic Interface Blocks .. INTERFACE DDE_SOLVER MODULE PROCEDURE DKL_1, DKL_2, DKL_3, DKL_4, ODEAVG END INTERFACE

Ce qui me trouble c'est que j'ai pas trouvé de déclaration de type FUNCTION DDE_SOLVER ou SUBROUTINE DDE_SOLVER dans le fichier...
Tu n'en trouveras pas. Ce que dit le bout de code, c'est que DDE_SOLVER est une fonction ou routine générique. Les fonctions spécifiques étant fournies par la ligne MODULE PROCEDURE. La selection de la fonction ou routine spécifique étant résolue par le nombre et le type des arguments.

Par exemple, quand tu fais CALL DDE_SOLVER(real), le compilateur utilisera la routine de la liste dont l'argument est real.
18/04/2007, 17h54
vinzzzz

Merci pour cette précision ;)