Comment Fortran retourne-t-il les baies?

Question

La sous-routine Rule_Tn de la bibliothèque Fortran CUBPACK nécessite un paramètre Integrand décrivant la fonction vectorielle intégrée. Il est un

INTERFACE 
    FUNCTION Integrand(NF,X) RESULT(Value) 
     USE Precision_Model 
     INTEGER,      INTENT(IN) :: NF 
     REAL(KIND=STND), DIMENSION(:), INTENT(IN) :: X 
     REAL(KIND=STND), DIMENSION(NF)    :: Value 
    END FUNCTION Integrand 
END INTERFACE 

Depuis que je veux appeler Rule_Tn à partir du code C je dois définir un type de fonction en C correspondant exactement à cette interface ci-dessus. J'ai donc essayé de comprendre comment une fonction Fortran retourne des tableaux. Au début, j'ai supposé la signature C suivante

void Integrand(double* value, const int* nf, const int* x); 

correspond à l'interface ci-dessus. Mais loin faux! J'ai une erreur de segmentation. Et j'ai déjà testé que le double est le type correspondant à REAL(KIND=STND), le STND sort du module Precision_Model.

Maintenant, quelqu'un peut-il me dire quelle est la bonne signature? J'utilise les compilateurs GNU Fortran et C.

Answer 1

Voir GNU docs. Il semble que vous ayez fourni les arguments dans un ordre différent entre Fortran et C. Essayez de mettre value en dernier dans le prototype C. En outre, il vous manque bind(C) sur la ligne FUNCTION.