Thank you for this bug report. I wonder if you would be kind enough to include some commented code in an email, so someone can fill out a trac report for the issue?
On Mon, Dec 15, 2008 at 7:13 PM, M. Yurko <myu...@gmail.com> wrote: > > I have noticed recently that when evaluating the EI function at any > number over 20, it adds 2pi i, which shouldn't be there. After looking > at the code it appears that it uses scipy, so the error is probably > there. However, when looking for alternatives, I found the > exponential_integral_1 from PARI function which doesn't seem to suffer > from the same issue as scipy. When I did a test using the identity Ei > (x) = -exponential_integral_1(-x), PARI's implementation was about > twice as fast as scipy and had a much greater precision. > > > --~--~---------~--~----~------------~-------~--~----~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URLs: http://www.sagemath.org -~----------~----~----~----~------~----~------~--~---
