[SciPy-user] exponential integral behaviour near -20

David Cournapeau david@ar.media.kyoto-u.ac...
Tue Dec 16 02:34:53 CST 2008

Robert Kern wrote:
> Ah, I think found it using this clue. It's a bug in SPECFUN. The
> "IMPLICIT DOUBLE PRECISION" statement is missing "A" so A0 is REAL
> rather than DOUBLE. Fixing that makes both of them go through the same
> code path. Can you change the line to this:
> in your specfun.f file, and rebuild scipy?

Sorry for the delay: you're right, this seems to fix the problem, at
least for me. The example now gives me:




More information about the SciPy-user mailing list