[Numpy-discussion] NaN (Not a Number) occurs in calculation of complex number for Bessel functions
bahtiyor_zohidov at mail.ru
Fri Dec 21 14:44:10 EST 2012
But I could find for Win64 bit windows????
Second question: Did you mean that I have to put lens limits of those number???
Пятница, 21 декабря 2012, 15:45 UTC от Pauli Virtanen <pav at iki.fi>:
>Your code tries to to evaluate
> z = 1263309.3633394379 + 101064.74910119522j
> jv(536, z)
> # -> (inf+inf*j)
>In reality, this number is not infinite, but
> jv(536, z) == -2.3955170861527422e+43888 + 9.6910119847300024e+43887
>These numbers (~ 10^43888) are too large for the floating point
>numbers that computers use (maximum ~ 10^308). This is why you get
>infinities and NaNs in the result. The same is true for the spherical
>You will not be able to do this calculation using any software
>that uses only floating point numbers (Scipy, Matlab, ...).
>You need to use analytical properties of your problem to
>get rid of such large numbers. Alternatively, you can use arbitrary
>precision numbers. Python has libraries for that:
>By the way, the proper place for this discussion is the following
>NumPy-Discussion mailing list
>NumPy-Discussion at scipy.org
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the NumPy-Discussion