*** glibc检测到***免费（）：无效的下一个大小（快） - 应该工作？

Question

我工作的一个现有的C项目（sourceforge上spglib），和我清理一些阵列初始化后运行到以下问题：

* glibc的检测*测试/ spglibtest：免费（） ：无效的下一个尺寸（快速）：0x08ab46e0 ***

回溯是：

#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 

的错误是在 “自由（反）” 下面一行：

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值分配给反像这样：

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

在get_translation写入到阵列时，Valgrind是显示以下内容：

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

这表明，我认为它试图写过去的结束为trans分配的内存，但它写入trans [0] [0]，trans是dimension [2] [3]。这应该工作，AFAIK，任何人都可以看到我失踪的东西？

Answer 1

您的类型错误，您无法将指针传递给期望数组数组（即指向数组的指针）的函数的指针数组。

对于你有get_translation签名，您需要：

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

Answer 2

所以这里可能是一个问题。你的函数似乎假设（基于）trans [num_trans] [j] = someDouble;这个trans实际上是一个按顺序排列的数组，正如我上面提到的在这种情况下不是真的。您正在分配一个指针数组而不是2维数组。类似于

double * trans = malloc（cell-> size * 3）;可能会更好。一般而言，您可能想使用1d数组而不是2d数组，并将其用作2维数组。