#4807: bug in exponential integral
----------------------+-----------------------------------------------------
Reporter: wdj | Owner: burcin
Type: defect | Status: new
Priority: major | Milestone: sage-3.4
Component: calculus | Resolution:
Keywords: |
----------------------+-----------------------------------------------------
Comment (by mabshoff):
The solution seems to be the following from the Scipy mailing list:
{{{
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:
> >
> > IMPLICIT DOUBLE PRECISION (A,D-H,O-Y)
> >
> > 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:
(25615628.4058-3.14159265359j)
(25615630.8316-3.14159265359j)
cheers,
David
}}}
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/4807#comment:2>
Sage <http://sagemath.org/>
Sage - Open Source Mathematical Software: Building the Car Instead of
Reinventing the Wheel
--~--~---------~--~----~------------~-------~--~----~
You received this message because you are subscribed to the Google Groups
"sage-trac" group.
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to
[email protected]
For more options, visit this group at
http://groups.google.com/group/sage-trac?hl=en
-~----------~----~----~----~------~----~------~--~---
