free():多次调用后下一个大小(快速)无效

Question

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

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

Answer 1

问题出在您的一个循环中。正如这篇文章中所说错误：free((：无效的下一个大小(快速(：

您可能会溢出缓冲区或以其他方式写入不应写入的内存，从而导致堆损坏。

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

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

如果数组p未0终止，这将覆盖堆中的内存。

话虽如此，如果您使用std::vector而不是手动分配的数组，您的代码会更清晰。