林銳博士在《高質量C++編程指南》中說過:“不可將浮點變量用“==”或“!=”與數字比較,應該設法轉化成“>=”或“<=”此類形式”,因為那是“隱含錯誤的寫法”。應該設法轉化成下面這樣的寫法:
-------------------------------
const float EPSINON = 0.00001; //精度限制
if((x >= -EPSINON) && (x <= EPSINON)) //表示x為“零值”
-------------------------------
剛看到的時候沒太在意,因為感覺平時編程極少會遇到要比較兩個浮點值是否相等的情況。后來得空一想,又想不明白這樣寫為什么是錯的。正巧Wayne過來串門,問他一下,結果他也不太清楚。那好吧,不知道原因,那就先寫個程序試一下吧,看到底是不是錯誤寫法:
-------------------------------
#include <iostream>
using std::cout;
using std::endl;
int main()
{
double ff = 3;
double aa = 3.6;
double tmp = aa / 1.2 - ff;
if (tmp != 0)
{
cout << "1" << endl;
}
else
{
cout << "2" << endl;
}
return 0;
}
-------------------------------
在我們看來,上面程序的輸出結果應該是2的,但是很遺憾,結果是1。——我在Cygwin中用g++(gcc版本為3.4.4)測試時就是這個結果,但是用VC6.0測試時就是輸出2了。可見這種寫法不能保證一定正確,果然是隱含錯誤啊。再來試試林銳博士的辦法:
-------------------------------
#include <iostream>
using std::cout;
using std::endl;
int main()
{
double ff = 3;
double aa = 3.6;
double tmp = aa / 1.2 - ff;
if (tmp < 1e-9999999999999999 && tmp > -1e-9999999999999999)
{
cout << "1" << endl;
}
else
{
cout << "2" << endl;
}
return 0;
}
-------------------------------
這下正確了。但是精度最多可以到多少呢?根據我的測試,在上面的if語句中不論寫多少個“9”都不會有問題,只要不超出編譯器所約束的一行中最大字符數限制就可以了。我用Cygwin中的g++(gcc版本為3.4.4)測試的時候最多寫了30720個9也沒有問題,不知道后面到底還可以加多少;但把這個源文件用VC6.0打開時就不行了,提示說“這個文件包含最大2048個字符的長行。行將換行。”,編譯自然也是通不過的了。
現在再想想就明白前一種寫法為什么不對了,在計算機組成原理中我們學過計算機對浮點數的記錄和處理方法。在一系列地計算過程當中,總是少不了“對階”的操作,這就引入了計算精度的問題。所以就有很多我們人類看著是相等的東西,給計算機看就不等了。
posted on 2008-03-12 17:13
李陽 閱讀(653)
評論(0) 編輯 收藏 引用 所屬分類:
C++