free():多次调用后无效的下一个尺寸(快速)

free(): invalid next size (fast) after several calls

我完全不确定这个函数怎么可能有内存分配错误。

template <typename T, typename U> // x -> ax+b
void transform_coords(double a, double b, int n, const T* oldc, U* newc) {
  int i, j, q0, q;
  if (n<2) return;

  int *p = new int[n]; // Pascal's triangle
  p[0] = 1;
  for (i=1; i<n; ++i) p[i] = 0;
  for (i=0; i<n; ++i) newc[i] = 0;

  for (j=0; j<n; ++j) {
    for (i=0; i<n; ++i) {
      if (p[i]==0) break;
      double m = p[i];
      if (i!=j) m *= std::pow(a,j-i);
      if (i!=0) m *= std::pow(b,i);
      newc[i] += m * oldc[j];
    }
    for (i=1, q=1; q; ++i) {
      q0 = p[i];
      p[i] += q;
      q = q0;
    }
  }

  delete[] p;
}

第三次 电话中,我得到了

*** glibc detected *** ./bin/bin_var: free(): invalid next size (fast)

如果我注释掉 delete[] p; 它会在所有调用中正常运行。

据我所知,关于内存分配,这个函数等同于

void transform_coords(int n) {
  if (n<2) return;
  int *p = new int[n];
  delete[] p;
}

问题出在你的一个循环中。正如 上所说的 in this post 错误:free():下一个尺寸无效(快速):

You may be overflowing a buffer or otherwise writing to memory to which you shouldn't be writing, causing heap corruption.

我最好的猜测是这个循环:

for (i=1, q=1; q; ++i) {
  q0 = p[i];
  p[i] += q;
  q = q0;
}

如果p数组没有0终止,这将覆盖堆中的内存。

也就是说,如果您使用 std::vector 而不是手动分配的数组,您的代码会更清晰。