前些日子在使用Kinsol庫(求解非線性方程的數值庫)時，出現死循環。調試跟蹤，發現浮點型數據異常，異常的數據使得大小比較失效，無法退出循環。這個問題涉及到浮點型數據的具體構造，詳細介紹可參考下文：http://steve.hollasch.net/cgindex/coding/ieeefloat.html
本文討論怎樣的運算造成異常的浮點型數據，即異常數據的來源。浮點型數據分成三類，finite型，inf型，nan型。finite型好理解，是指有限數如0，1.1等。inf型是表示無窮大數，分為正無窮和負無窮，它由非0的有限數除0得到或者超過了浮點型能夠表示的最大最小數。浮點型不像整型那樣除零報錯。nan型即not a number， 表示這不是一個數。它是由一些無意義的運算引起的，比如0/0，inf-inf，inf*0等，這些運算不能確定結果是什么。
前兩種類型是可以比較大小，比如：1<inf，inf > 3，inf > -inf，inf == inf。這些比較是有意義的。但不能和nan型比較，如果比較則總返回0（返回什么值，大家最好自己做測試），我遇到的問題就是在比較表達式中出現了nan型數。
要判斷一個浮點數是否是inf或者nan，有函數int _isnan(double x)和int _finite(double x)，在float.h中聲明。