C ++ код падает с «free (): неверный следующий размер» - PullRequest
Новая
       6

C ++ код падает с «free (): неверный следующий размер»

1 голос
/ 16 ноября 2011

Я написал небольшую программу, которая использует функциональные указатели для выполнения некоторых численных расчетов. 

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]

Итак, мой вопрос, что здесь не так?

1 Ответ

5 голосов
/ 16 ноября 2011

В этом коде нет ошибки, при условии, что вы всегда вызываете эту функцию с помощью l_end >= 1.

Когда вместо l_end == 0 есть операция записи за границу в p_l[1] = xi;.

Обратите внимание, однако, что вы не можете сделать вывод, что эта функция имеет проблему только потому, что именно здесь происходит сбой, или просто потому, что не вызывая эту функцию, у вас нет сбоя.

Ошибка - это ошибка, икрушение крушение.Они полностью различны в C ++;чем раньше вы поймете этот важный факт, тем лучше.Может быть ошибка где-то еще, и эта функция может быть просто жертвой.

Если вы видите сбой, то есть ошибка.Если вы не видите сбоя, вы ничего не знаете (возможно, ошибка все еще присутствует).

...