Discussion :

[tiresias54]

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 :

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

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 :

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

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 :

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

1
2
3
4
5
6
 
nelem = 10
elsize=1
work=0
tbaddr=0
offset=0

en sortie:

Code :

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

1
2
3
4
5
 nelem =  10
 elsize =  1
 work =  0
 tbaddr =  415111376
 offset =  408813608

[Sve@r]

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 c :

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

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...

[tiresias54]

merci pour ton aide

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

Code :

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

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...

[Metalman]

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... ), 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 !)

[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).

Et qu'est censé signifier basevec?

[tiresias54]

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 :/).

[Médinoc]

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

Code C :

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

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);
}

[tiresias54]

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

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

[Sve@r]

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

Ouais, de honte...

[grim7reaper]

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

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.

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.

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

C’est tout à fait ça.

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.

Ç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).

[grim7reaper]

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 fortran :

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

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 C :

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

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 :

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

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).

[tiresias54]

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.

[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*).

[tiresias54]

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 :

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

#include <stdint.h>

en début de fichier *.c

Code :

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

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 :

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

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 :

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

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

[Médinoc]

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

[souviron34]

Ç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...