[SciPy-user] exponential integral behaviour near -20

[David Cournapeau]

Robert Kern wrote:
>
> I cannot replicate this on OS X with a recent(ish) SVN scipy compiled
> with gfortran 4.2.1. This function is implemented in FORTRAN, so the
> FORTRAN compiler used may be relevant.

Isn't exp1 implemented from exp1n in cephes  ?

>  Can you give us more
> information about the platform, scipy version, and FORTRAN compiler?
>   

FWIW, I can replicate the reported behavior on Linux (x86 - gfortran
4.3.2 + recent SVN for both numpy and scipy as well).

cheers,

David