[Numpy-discussion] NaN (Not a Number) occurs in calculation of complex number for Bessel functions

Happyman bahtiyor_zohidov@mail...
Fri Dec 21 06:56:08 CST 2012


 Thanks Pauli

But I have already very shortly built  for bessel function, but the code you gave me is in Fortran.. I also used f2py but I could not manage to read fortran codes..that is why I have asked in Python what is wrong??


Пятница, 21 декабря 2012, 12:46 UTC от Pauli Virtanen <pav@iki.fi>:
>Happyman <bahtiyor_zohidov <at> mail.ru> writes:
>[clip]
>> IF I GIVE ( it is necessary value for my program ): 
>> a , b = sph_jn ( 536 , 2513.2741228718346 + 201.0619298974676j )
>
>The implementation of the spherical Bessel functions is through
>this Fortran code:
>
>https://github.com/scipy/scipy/blob/master/scipy/special/specfun/specfun.f#L1091
>
>It does not have asymptotic expansions for dealing with parts of
>the complex plane where the computation via the recurrence does not
>work.
>
>-- 
>Pauli Virtanen
>
>
>_______________________________________________
>NumPy-Discussion mailing list
>NumPy-Discussion@scipy.org
>http://mail.scipy.org/mailman/listinfo/numpy-discussion

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.scipy.org/pipermail/numpy-discussion/attachments/20121221/4fdfc41d/attachment.html 


More information about the NumPy-Discussion mailing list