[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,