林銳博士在《高質(zhì)量C++編程指南》中說(shuō)過(guò):“不可將浮點(diǎn)變量用“==”或“!=”與數(shù)字比較,應(yīng)該設(shè)法轉(zhuǎn)化成“>=”或“<=”此類形式”,因?yàn)槟鞘?#8220;隱含錯(cuò)誤的寫法”。應(yīng)該設(shè)法轉(zhuǎn)化成下面這樣的寫法:
-------------------------------
const float EPSINON = 0.00001; //精度限制
if((x >= -EPSINON) && (x <= EPSINON)) //表示x為“零值”
-------------------------------
剛看到的時(shí)候沒(méi)太在意,因?yàn)楦杏X(jué)平時(shí)編程極少會(huì)遇到要比較兩個(gè)浮點(diǎn)值是否相等的情況。后來(lái)得空一想,又想不明白這樣寫為什么是錯(cuò)的。正巧Wayne過(guò)來(lái)串門,問(wèn)他一下,結(jié)果他也不太清楚。那好吧,不知道原因,那就先寫個(gè)程序試一下吧,看到底是不是錯(cuò)誤寫法:
-------------------------------
#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;
}
-------------------------------
在我們看來(lái),上面程序的輸出結(jié)果應(yīng)該是2的,但是很遺憾,結(jié)果是1。——我在Cygwin中用g++(gcc版本為3.4.4)測(cè)試時(shí)就是這個(gè)結(jié)果,但是用VC6.0測(cè)試時(shí)就是輸出2了。可見這種寫法不能保證一定正確,果然是隱含錯(cuò)誤啊。再來(lái)試試林銳博士的辦法:
-------------------------------
#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;
}
-------------------------------
這下正確了。但是精度最多可以到多少呢?根據(jù)我的測(cè)試,在上面的if語(yǔ)句中不論寫多少個(gè)“9”都不會(huì)有問(wèn)題,只要不超出編譯器所約束的一行中最大字符數(shù)限制就可以了。我用Cygwin中的g++(gcc版本為3.4.4)測(cè)試的時(shí)候最多寫了30720個(gè)9也沒(méi)有問(wèn)題,不知道后面到底還可以加多少;但把這個(gè)源文件用VC6.0打開時(shí)就不行了,提示說(shuō)“這個(gè)文件包含最大2048個(gè)字符的長(zhǎng)行。行將換行。”,編譯自然也是通不過(guò)的了。
現(xiàn)在再想想就明白前一種寫法為什么不對(duì)了,在計(jì)算機(jī)組成原理中我們學(xué)過(guò)計(jì)算機(jī)對(duì)浮點(diǎn)數(shù)的記錄和處理方法。在一系列地計(jì)算過(guò)程當(dāng)中,總是少不了“對(duì)階”的操作,這就引入了計(jì)算精度的問(wèn)題。所以就有很多我們?nèi)祟惪粗窍嗟鹊臇|西,給計(jì)算機(jī)看就不等了。
posted on 2008-03-12 17:13
李陽(yáng) 閱讀(653)
評(píng)論(0) 編輯 收藏 引用 所屬分類:
C++