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

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#L... 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