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