爲什麼gdb將sqrt（3）評估爲0？

Question

平方根的3中，作爲估計由鎢阿爾法：

1.7320508075688772935274463415058723669428052538103806280558... 

當我用C做sqrt(3)，它的計算結果爲0。爲什麼呢？

EDIT4：以下介紹如何在GDB中重現此問題。創建test.c如下：

#include <stdio.h>                                          
#include <math.h> 

int main() 
{ 
    printf("sqrt(3): %f\n", sqrt(3)); 
    return 0; 
} 

編譯：

gcc -O0 -g -Wall -pedantic -ansi -lm -o test test.c 

運行調試器：

gdb test 

在控制檯輸入以下內容：

(gdb) break test.c:6 
Breakpoint 1 at 0x400578: file test.c, line 6. 
(gdb) r 
Starting program: /home/pdedecker/Desktop/test 
Breakpoint 1, main() at test.c:6 
6   printf("sqrt(3): %f\n", sqrt(3)); 
(gdb) print sqrt(3) 
$1 = 0 
(gdb) s 
sqrt(3): 1.732051 

我的GDB版本爲GNU gdb (GDB) SUSE (7.1-3.12)。

Answer 1

的問題不是缺少函數聲明（這ISN不要錯過，因爲你確實包括了<math.h>）。

的問題是缺少對sqrt你實際使用的是調試信息。如果沒有這些調試信息，GDB不知道要傳遞給sqrt()的參數類型以及返回的內容。

您可以通過安裝的libc-debuginfo軟包得到很多Linux發行版所需要的調試信息。以下是我看到這樣一個系統：

gdb -q ./a.out 
Reading symbols from /tmp/a.out...done. 
(gdb) b main 
Breakpoint 1 at 0x400558: file t.c, line 6. 
(gdb) r 

Breakpoint 1, main() at t.c:6 
6  printf("sqrt(3): %f\n", sqrt(3)); 
(gdb) p sqrt 
$1 = {<text variable, no debug info>} 0x7ffff7b7fb50 <__sqrt> 

注： 「沒有調試信息」

(gdb) p sqrt(3) 
$2 = 0 
(gdb) p sqrt(3.0) 
$3 = 0 

注意：符合你的行爲。 什麼sqrt函數做有調試信息？

(gdb) info func sqrt 
All functions matching regular expression "sqrt": 

File ../sysdeps/x86_64/fpu/e_sqrt.c: 
double __ieee754_sqrt(double); 

File s_csqrt.c: 
complex double __csqrt(complex double); 

File ../sysdeps/x86_64/fpu/e_sqrtf.c: 
float __ieee754_sqrtf(float); 

File w_sqrtf.c: 
float __sqrtf(float); 

File s_csqrtf.c: 
complex float __csqrtf(complex float); 

File ../sysdeps/i386/fpu/e_sqrtl.c: 
long double __ieee754_sqrtl(long double); 

File w_sqrtl.c: 
long double __sqrtl(long double); 

File s_csqrtl.c: 
complex long double __csqrtl(complex long double); 

File ../sysdeps/ieee754/dbl-64/mpsqrt.c: 
void __mpsqrt(mp_no *, mp_no *, int); 

File w_sqrt.c: 
double __sqrt(double); 

(gdb) p __sqrt 
$4 = {double (double)} 0x7ffff7b7fb50 <__sqrt> 

注：__sqrt是在同一地址爲sqrt，但GDB知道它的類型！

(gdb) p __sqrt(3) 
$5 = 1.7320508075688772 
(gdb) p __sqrt(3.0) 
$6 = 1.7320508075688772 

可以合理地認爲這是GDB中的一個bug。隨意創建一個在GDB bugzilla。

Answer 2

#include <stdio.h>                                          
#include <math.h> 

int main() 
{ 
    printf("sqrt(3): %f\n", sqrt(3)); 

    return 0; 
} 

輸出：

josh@josh-ubuntu:~/scratch$ ./a.out 
sqrt(3): 1.732051 

Answer 3

我預測，你沒有做#include <math.h>

沒有一個函數聲明，C將一個函數的返回值默認爲int。取決於int的大小，浮點數可能會回到0。 C也不知道如何轉換函數參數。它會默認將參數傳遞爲它所發生的類型。如果您傳遞一個整數到sqrt()它不會被轉換爲雙精度，但sqrt()函數將解釋位模式爲double。

Answer 4

也許不支持調用sqrt！也許是因爲它是一個libc函數。我不知道的深層原因，但下面的測試顯示了一個有趣的現象：

double mysqrt(double x) { return sqrt(x) }; 

然後在GDB會話：

(gdb) p mysqrt(3) 
$1 = 1.7320508075688772 
(gdb) p sqrt(3) 
$2 = -1209775368 

Answer 5

要調用沒有調試信息的函數，必須使用函數指針cast來顯式地告訴gdb返回值和參數的類型。因此，例如：

(gdb) print ((double (*) (double)) sqrt) (3) 
$1 = 1.7320508075688772