[interfaçage C/Fortran] erreur de segmentation

Bonjour à tous,
J'ai hérité d'une librairie C écrite il y a des lustres et qui sert à allouer une plage mémoire.
Cette librairie est apellée par un programme fortran qui passe la variable à allouer en paramètre.
j'ai crée un petit programme de test qui reproduit le problème car malheureusement je ne vois pas ce qui cloche :/

Le programme fortran:

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12

program testlib   integer nelem,elsize,work,tbaddr,offset nelem=10 elsize=1 work=0 tbaddr=0 offset=0   call falloc(nelem,elsize,work,tbaddr,offset)   end program testlib

---

La librairie 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 26 27 28 29 30

#include <stdlib.h> #include <stdio.h>   falloc_(nelem,elsize,basevec,addr,offset) int *nelem,*elsize,*basevec,*addr,*offset; { int noct; int adbase; int *ptr; int rest;   if( *elsize > 8) noct=(*nelem + 1) * (*elsize); /* Nombre d'octets */ else noct=(*nelem)*(*elsize); /* Nombre d'octets */   if(noct<=0){ *addr=0; printf("noct < 0 bizarre \n"); return; }   ptr=(int *) malloc(noct); /* Allocation */ /*(int *)*addr=ptr; */ *addr=(int)ptr; adbase=(int)basevec; *offset=(*addr-adbase)/(*elsize); rest=(*addr-adbase)-(*offset) * (*elsize); if(rest!=0) *offset=*offset+1; }

---

Je compile avec sunstudio sans options particulières.
Merci à celui/celle qui pourrais jeter un oeil et me dire pourquoi j'ai des valeurs de sortie incohérentes ^^

par exemple: en entrée j'ai

Code:

---

1 2 3 4 5 6

nelem = 10 elsize=1 work=0 tbaddr=0 offset=0

---

en sortie:

Code:

---

1 2 3 4 5

nelem = 10 elsize = 1 work = 0 tbaddr = 415111376 offset = 408813608

---

Bonjour
Je ne connais pas fortran donc mes remarques ne peuvent se porter que sur la fonction falloc_(). Au passage on remarquera que la fonction se nomme falloc_ alors que le code fortran appelle juste falloc. De plus la fonction est sensée recevoir des adresses mais le code Fortran passe des valeurs mais je présume que c'est l'usage quand on interface Fortran/C.

Bon, ceci dit, voici ce qui me choque: déjà on ne caste plus un malloc depuis que celui-ci renvoie un void*. Au contraire, le cast d'un malloc peut en retour masquer un oubli #include. Mais je suis surtout étonné par cette ligne 25 *addr=(int)ptr;. ptr contient l'adresse allouée par malloc et tu vas copier cette adresse dans un int ??? Je ne sais pas comment on gère les adresse en Fortran mais ça ça me choque. De même cette copie de basevec (une adresse) dans adbase (un int)...

Au passage, bien que je n'en comprenne pas la raison, on peut optimiser un peu ton bloc if/else de cette façon

Code:

---

1 2 3

noct=(*nelem)*(*elsize); /* Nombre d'octets */ if( *elsize > 8) noct+=(*elsize);

---

Et la ligne 28 peut se raccourcir en rest=(*addr-adbase) % (*elsize); et la variable adbase ne sert à rien du tout vu que c'est le contenu de basevec

Maintenant, peut-être qu'en rajoutant quelques printf() dans ta fonction tu pourras affiner l'instruction qui pose soucis...

merci pour ton aide :)

En mettant des printf dans le code, on se rend compte que les valeurs

Code:

---

1 2	tbaddr = 415111376 offset = 408813608

---

que j'obtiens en sortie (et qui normalement sont les valeurs des variables) correspondent (si je comprend bien) à la valeur des adresses converties en entier dans la librairie...

Question bête : tu n'as aucune option pour débugger ou mettre le numéro de la ligne qui segfault ?
Si tu as sunstudio, je suppose que tu es sur une station Sun (thx cpt. obvious ! pour moi... :P), n'as-tu aucun port de gcc ou gdb dessus ? Ni l'option pour faire du pas à pas ? (si sunstudio possède une IDE, ce que j'espère !)

Ça fait bizarre, tous tes paramètres déclarés en int*, sans aucune notion de type (tbaddr devrait être au minimum void**) ou de constante (pour les paramètres que la fonction ne modifie pas).

Et qu'est censé signifier basevec?

je suis sous Linux avec un debogueur (totalview) qui m'affiche les valeurs des variables au petit bonheur la chance (apparemment il aime pas du tout les pointeurs :/).

Déjà, le même code avec un minimum de rigueur (et de modernité):

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

#include <stdlib.h> #include <stdio.h>   void falloc_impl(int nelem, int elsize, int basevec, void **pAddr, int *pOffset) { size_t noct; intptr_t adbase, adalloc; void *ptr; int rest;     if( elsize > 8) noct = (nelem + 1) * elsize; /* Nombre d'octets */ else noct = nelem*elsize; /* Nombre d'octets */   if(noct<=0){ *pAddr=NULL; printf("noct < 0 bizarre \n"); return; }   ptr = malloc(noct); /* Allocation */ *pAddr = ptr; adbase = (intptr_t)basevec; adalloc = (intptr_t)ptr; *pOffset = (adalloc-adbase)/elsize; rest=(adalloc-adbase)-(*pOffset)*elsize; if(rest!=0) (*pOffset)++; }   void falloc_(int *nelem, int *elsize, int *basevec, void **addr, int *offset) { falloc_impl(*nelem, *elsize, *basevec, addr, offset); }

---

Je me suis renseigné, il parait que le gars qui a écrit la librairie est mort :whistle:

Merci pour votre aide, je vais essayer de me débrouiller avec les infos que vous m'avez données :)

Citation:

---

Envoyé par tiresias54

Je me suis renseigné, il parait que le gars qui a écrit la librairie est mort :whistle:

---

Ouais, de honte...:aie:

@tiresias54 : tu fait du FORTRAN ou du Fortran ? :P
Dis autrement, tu fais du Fortran 77 (ou inférieur) ou du Fortran 90 (ou supérieur) ?

Citation:

---

Envoyé par Sve@r

Au passage on remarquera que la fonction se nomme falloc_ alors que le code fortran appelle juste falloc.

---

Les compilateurs Fortran ajoute un underscore aux noms de fonctions. Ajouter cet underscore à la main dans le code C est une des façons d’interfacer les deux.

Citation:

---

Envoyé par Sve@r

De plus la fonction est sensée recevoir des adresses mais le code Fortran passe des valeurs

---

Fortran passe les paramètres par référence en général.

Citation:

---

Envoyé par Sve@r

mais je présume que c'est l'usage quand on interface Fortran/C.

---

C’est tout à fait ça.

Citation:

---

Envoyé par Sve@r

Bon, ceci dit, voici ce qui me choque: déjà on ne caste plus un malloc depuis que celui-ci renvoie un void*.

---

Vu la tête du code, c’est du vieux C (pré-ANSI), donc malloc devait renvoyer char* à cette époque.

Citation:

---

Envoyé par Médinoc

Ça fait bizarre, tous tes paramètres déclarés en int*, sans aucune notion de type

---

Déclaration à la syntaxe K&R (pré-ANSI donc), toujours valide mais dépréciée (Cf. la normedu C99 par exemple, 6.9.1 Function definitions, page 155) car il y a pas mal d’inconvénient (dont pas de vérification de types entre autres).

Bon, je viens de faire un test et, à part le fait que le code soit écrit dans un style très vieux, il semble fonctionner correctement (à condition que la taille d’un pointeur soit inférieure ou égale à la taille d’un int).

Soit le code suivant (je n’ai rien modifié par rapport au code d’origine, juste ajouter des printf/print :
testlib.f

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13

program testlib   integer nelem,elsize,work,tbaddr,offset nelem=10 elsize=1 work=0 tbaddr=0 offset=0   print *, nelem, elsize, work, tbaddr, offset call falloc(nelem,elsize,work,tbaddr,offset) print *, nelem, elsize, work, tbaddr, offset end program testlib

---

falloc.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 26 27 28 29 30 31

#include <stdlib.h> #include <stdio.h>   falloc_(nelem,elsize,basevec,addr,offset) int *nelem,*elsize,*basevec,*addr,*offset; { int noct; int adbase; int *ptr; int rest;   if( *elsize > 8) noct=(*nelem + 1) * (*elsize); /* Nombre d'octets */ else noct=(*nelem)*(*elsize); /* Nombre d'octets */   if(noct<=0){ *addr=0; printf("noct < 0 bizarre \n"); return; }   ptr=(int *) malloc(noct); /* Allocation */ /*(int *)*addr=ptr; */ *addr=(int)ptr; adbase=(int)basevec; *offset=(*addr-adbase)/(*elsize); rest=(*addr-adbase)-(*offset) * (*elsize); if(rest!=0) *offset=*offset+1; printf("%d %d\n", *addr, *offset); }

---

Résultat :

Code:

---

1 2 3 4 5 6

% gcc -m32 -c falloc.c % gfortran --std=legacy -ffree-form -m32 -o tiresias54 testlib.f falloc.o % ./tiresias54 10 1 0 0 0 135003072 135968036 10 1 0 135003072 135968036

---

Les valeurs obtenues dans le codeFotran sont bien celles produites dans le code C donc l’interfaçage semble OK.
Du coup, est ce que tu aurais un exemple de code qui reproduit une erreur de segmentation ?

Bon j’ai dû compiler en 32-bit, sinon la conversion pointeur->int aurait été foireuse.

(Désolé pour le double post).

merci grim, ton idée m'a mis sur la piste: en remplaçant les integer par des integer*4 et les real par des real*4 dans mon code appellant, j'ai résolu le probème :)

Pour info le code a plus de 20 ans.

Je pense que tu devrais plutôt faire l'inverse: Utiliser des size_t et intptr_t côté C au lieu des int (plus un void** pour le retour du pointeur, au lieu d'un int*).

Citation:

---

Envoyé par Médinoc

Je pense que tu devrais plutôt faire l'inverse: Utiliser des size_t et intptr_t côté C au lieu des int (plus un void** pour le retour du pointeur, au lieu d'un int*).

---

du coup ce code te parait-il pertinent (et sûr)? (les lignes que j'ai modifiées sont en gras). Et j'ai du ajouter un

Code:

---

#include <stdint.h>

---

en début de fichier *.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 26

int *nelem,*elsize,*basevec,*addr,*offset; { size_t noct; size_t adbase; size_t *ptr; size_t rest; if( *elsize > 8) noct=(*nelem + 1) * (*elsize); /* Nombre d'octets */ else noct=(*nelem)*(*elsize); /* Nombre d'octets */ if(noct<=0){ *addr=0; printf("noct < 0 bizarre \n"); return; } ptr=(void **) malloc(noct); /* Allocation */ *addr=(intptr_t)ptr; adbase=(intptr_t)basevec; *offset=(*addr-adbase)/(*elsize); rest=(*addr-adbase)-(*offset) * (*elsize); if(rest!=0) *offset=*offset+1; printf("%d %d\n", *addr, *offset); }

---

quoiqu'il en soit lorsque je compile j'ai des warnings du genre:

Code:

---

1 2	"falloc.c", line 19: warning: assignment type mismatch: pointer to unsigned long "=" pointer to pointer to void

---

avec sunstudio cc

ou bien égalemert (si je compile avec gnu gcc):

Code:

---

1 2	falloc.c: In function ‘falloc_’: falloc.c:19: warning: assignment from incompatible pointer type

---

Retirer le cast de malloc() devrait corriger le warning.

Citation:

---

Envoyé par Médinoc

Ça fait bizarre, tous tes paramètres déclarés en int*, sans aucune notion de type (tbaddr devrait être au minimum void**) ou de constante (pour les paramètres que la fonction ne modifie pas).

---

manière de déclaration pré C90....


Quant au problème, outre le fait que Fortran passe par référence, il y a un cas particulier pour les tableaux (et donc les chaines aussi) : Fortran attend à la suite la taille..

Mais le mieux est de poster sur le forum Fortran si des questions subsistent.. Eux sont habitués, ça arrive tout le temps les conversions dans un sesn et dans l'autre...