2011-11-15 68 views
1

我寫了一個小程序,它使用函數指針來做一些數值計算。C++代碼崩潰與「免費():無效的下一個大小」

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

該程序崩潰與一個不尋常的免費()錯誤。如果我將函數指針更改爲第一個函數(polynom),它可以正常工作,但它會與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

你能提供一個完整的小例子來證明這個問題嗎? –

+4

另外,'l'是可能的最差變量名! –

+1

也許你應該發佈你的代碼,使用函數指針... – stefan

回答

5

該代碼沒有錯誤,前提是您始終使用l_end >= 1調用該函數。

l_end == 0代替p_l[1] = xi;出現邊界外寫入操作時。

但是,請注意,您不能推斷這是具有問題的功能,因爲這是您遇到崩潰的位置,或者僅僅是因爲沒有調用此函數就沒有崩潰。

錯誤是錯誤,崩潰是崩潰。它們在C++中完全不同;越早意識到這個重要事實越好。其他地方可能有錯誤,這個功能可能只是受害者。

如果您看到崩潰,則會出現錯誤。如果你看不到崩潰,你什麼都不知道(錯誤可能仍然存在)。

+0

這使它。非常感謝。一旦我的限制下降,將會批准答案。 –

相關問題