2

我正在开发一个现有的 c 项目(sourceforge 上的 spglib),在清理了一些数组初始化后我遇到了以下问题:

* 检测到 glibc *测试/spglibtest:free():下一个大小无效(快速):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

错误在下面的“free(trans)”行中:

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 将值分配给 trans ,如下所示:

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 的维度为 [2][3]。这应该可行,AFAIK,任何人都可以看到我缺少的东西吗?

4

2 回答 2

6

您的类型错误,您不能将指向指针数组的指针传递给期望数组数组的函数(即指向数组的指针)。

对于get_translation您拥有的签名,您需要:

double (*trans)[3] = malloc(cell->size * sizeof(double[3]));
于 2010-09-15T22:18:40.227 回答
2

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

double* trans = malloc(cell->size * 3); 可能会更好。通常,您可能希望使用一维数组而不是二维数组,并将其用作二维数组。

于 2010-09-15T22:25:59.680 回答