SpecialFunctions.jl icon indicating copy to clipboard operation
SpecialFunctions.jl copied to clipboard

Large negative real number arguments return NaN for expintx

Open weymouth opened this issue 2 years ago • 1 comments

expintx(complex(-x)) returns NaN for x more than around 700. It's almost like this is being handled by -exp(-x)*expinti(x) and we're getting an overflow error.

weymouth avatar Jan 12 '24 14:01 weymouth

It seems to be about -709.782712894 at which point the result goes from roughly -0.00141087244 + 1e-308*im to -Inf + NaN*im.

ararslan avatar Jan 18 '24 18:01 ararslan