*** glibc détecté *** gratuit(): invalide taille suivante (rapide) - devrait fonctionner?

Question

Je travaille sur un projet c existant (spglib sur sourceforge), et je suis en cours d'exécution dans le problème après le nettoyage des initialisations de tableau:

* glibc détectée * tests/spglibtest: free() : taille suivante invalide (rapide): 0x08ab46e0 ***

le est backtrace:

#0 0xb7fe1424 in __kernel_vsyscall() 
#1 0xb5cfdd61 in raise() from /lib/libc.so.6 
#2 0xb5cff5ee in abort() from /lib/libc.so.6 
#3 0xb5d397ed in ??() from /lib/libc.so.6 
#4 0xb5d3f7b1 in ??() from /lib/libc.so.6 
#5 0xb5d4052b in ??() from /lib/libc.so.6 
#6 0xb5d441cd in free() from /lib/libc.so.6 
#7 0xb6681484 in sym_get_multiplicity (cell=0xbfffe1f0, symprec=0.050000000000000003) at /git/xtalopt-public/src/spglib/symmetry.c:168 
#8 0xb6680550 in spg_find_primitive (lattice=0xbfffe2a8, position=0x813c6f0, types=0x813c700, num_atom=2, symprec=0.050000000000000003) 
    at /git/xtalopt-public/src/spglib/spglib.c:253 

l'erreur est dans la ligne "libre (trans)" ci-dessous:

int sym_get_multiplicity(const Cell *cell, const double symprec) 
{ 
    int i, rc; 
    double **trans; 
    trans = (double**)malloc(cell->size * sizeof(double*)); 
    for (i = 0; i < cell->size; i++) { 
    trans[i] = (double*)malloc(3 * sizeof(double)); 
    } 

    rc = get_translation(&trans[0][0], identity, cell, symprec); 

    for (i = 0; i < cell->size; i++) { 
    free(trans[i]); 
    } 
    free(trans); 

    return rc; 
} 

get_translation attribue des valeurs à trans comme ceci:

static int get_translation(double trans[][3], const int rot[3][3], const Cell *cell, 
          const double symprec) 
{ 
... 
    for (j = 0; j < 3; j++) { 
    trans[num_trans][j] = someDouble; 
    } 
... 
} 

Valgrind montre ce qui suit lors de l'écriture au tableau dans get_translation:

==17929== Invalid write of size 8 
==17929== at 0x56BE8A7: get_translation (symmetry.c:285) 
==17929== by 0x56BE44B: sym_get_multiplicity (symmetry.c:163) 
... 
==17929== Address 0x9cb5868 is 0 bytes after a block of size 8 alloc'd 
==17929== at 0x4024918: malloc (in /usr/lib/valgrind/vgpreload_memcheck-x86-linux.so) 
==17929== by 0x56BE3F7: sym_get_multiplicity (symmetry.c:158) 
.... 

Cela me donne à penser qu'il est en train d'écrire après la fin de la mémoire allouée pour trans, mais il écrit à trans [0] [0], et trans est de dimension [2] [3]. Cela devrait fonctionner, AFAIK, quelqu'un peut-il voir quelque chose qui me manque?

Answer 1

Vos types sont mauvais, vous ne pouvez pas passer un pointeur sur un tableau de pointeurs vers une fonction attendant un tableau de tableaux (par exemple un pointeur sur un tableau).

Pour la signature de get_translation que vous avez, vous avez besoin:

double (*trans)[3] = malloc(cell->size * sizeof(double[3])); 

Answer 2

est ici peut-être un SO problème. Votre fonction semble supposer (basé sur) trans [num_trans] [j] = someDouble; ce trans est en fait un tableau agencé séquentiellement, ce qui, comme je l'ai mentionné plus haut, n'est pas vrai dans ce cas. Vous allouez un tableau de pointeurs plutôt qu'un tableau à deux dimensions. Quelque chose de semblable à

double * trans = malloc (Cell-> Taille * 3); pourrait être mieux. En général, vous pouvez utiliser un tableau 1d au lieu d'un tableau 2d et l'utiliser comme tableau à deux dimensions.