Question

我编写了一个小程序，它使用函数指针进行一些数值计算。

double polynom(const int j, const double xi) {
  return pow(xi, j);
}

/**
 * Calculate the legendre_polynom l_end on a certain position xi.
 */
double legendre_polynom(const int l_end, const double xi) {
  vector <double> p_l(l_end+1);
  p_l[0] = 1.0;
  p_l[1] = xi;

  for (int x = 2; x <= l_end; x++) {
    // p_l = ((2x-1) * p_{x-1} - (x-1) * p_{x-2}) / l
    p_l[x] = ((2 * x - 1) * p_l[x - 1] - (x - 1) * p_l[x - 2]) / x;
  }

  double result = p_l[l_end];
  return result;
}

程序因异常的free（）错误而崩溃。如果我将函数指针更改为第一个函数（多项式），它可以正常工作，但它会因legendre_polynom而失败。

我已经调试了那么远，它在退出该函数之后和其他代码继续之前就会中断。

*** glibc detected *** blub: free(): invalid next size (fast): 0x0804f248 ***
======= Backtrace: ========= /lib/i386-linux-gnu/libc.so.6(+0x6ebc2)[0xb7d70bc2]
/lib/i386-linux-gnu/libc.so.6(+0x6f862)[0xb7d71862]
/lib/i386-linux-gnu/libc.so.6(cfree+0x6d)[0xb7d7494d]

...

number2(_ZN9__gnu_cxx13new_allocatorIdE10deallocateEPdj+0x11)[0x804bc8b]
number2(_ZNSt12_Vector_baseIdSaIdEE13_M_deallocateEPdj+0x25)[0x804bbc3]
number2(_ZNSt12_Vector_baseIdSaIdEED1Ev+0x37)[0x804ba33]
number2(_ZNSt6vectorIdSaIdEED1Ev+0x38)[0x804b8a0]
number2(_Z16legendre_polynomid+0x13f)[0x804af9b]

所以我的问题是这里有什么问题？

Answer 1

该代码中没有错误，前提是您始终使用l_end >= 1调用该函数。

当l_end == 0改为p_l[1] = xi;时出现边界外写操作。

但请注意，您无法推断这是出现此问题的函数，因为这是您遇到崩溃的原因，或者仅因为没有调用此函数而导致崩溃。

错误是错误，崩溃是崩溃。它们在C ++中完全不同;你越早意识到这个重要的事实就越好。其他地方可能存在错误，此功能可能只是受害者。

如果发现崩溃，则会出现错误。如果你没有看到任何崩溃，你什么都不知道（错误可能仍然存在）。

C ++代码崩溃与“free（）：无效的下一个大小”

1 个答案: