[Matrix-SIG] NaN (almost)

Yoon, Hoon (CICG - NY Program Trading) HYoon@exchange.ml.com
Wed, 7 Oct 1998 18:16:10 -0400


Thanks Paul, but not quite....

>>> excP
array([  9.12500000e+000,   9.27519777e+000,   9.25000000e+000,
9.06250000e+000,
              8.93750000e+000,               NaN,   9.00000000e+000,
              9.00000000e+000,   9.00000000e+000,               NaN,
              9.00000000e+000,               NaN,               NaN,
              9.00000000e+000,               NaN,   9.00000000e+000,
                          NaN,   9.00000000e+000,               NaN,
              8.93750000e+000,   8.93750000e+000,               NaN,
              9.02269591e+000,               NaN,               NaN,
                          NaN,               NaN,   0.00000000e+000])
>>> infini = (array([1.0])/array([0.0]))/[0]
>>> x = array([1.,2.,3.])
>>> y = array([1.,0.,1.])
>>> equal(x/y, infini)
array([0, 1, 0])
>>> NaNv = 1e1000000 - 1e10000000
>>> NaNv = array((1e1000000 - 1e10000000))
>>> NaNv
NaN
>>> NaNv = infini - infini
>>> NaNv
array([              NaN])
>>> equal(excP, NaNv)
array([0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0,
            0, 0])

No wonder

equal(NaNv, NaNv)  is 0 too:
Appreciate it neverthless,

Hoon,