C++コードが "free（）：次のサイズが無効"でクラッシュします。

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 == 0の代わりに、p_l[1] = xi;に境界外書き込み操作があります。

ただし、これはクラッシュする場所であるという理由だけで問題が発生している関数であると推測することはできません。

エラーはエラーであり、クラッシュはクラッシュです。それらはC++では完全に異なっています。この重要な事実をより早く理解することができれば早くなります。他の場所でエラーが発生している可能性があり、この機能はただの被害者である可能性があります。

クラッシュが発生した場合は、エラーが発生します。クラッシュがない場合は、何も知らない（エラーがまだ存在する可能性があります）。