#4807: bug in exponential integral
----------------------+-----------------------------------------------------
Reporter: wdj | Owner: burcin
Type: defect | Status: new
Priority: major | Milestone: sage-3.4.1
Component: calculus | Resolution:
Keywords: |
----------------------+-----------------------------------------------------
Comment (by was):
Mike Hansen observes that one could implement Ei in *general* and with
arbitrary precision by using the library "mpmath", which just happens to
already be in Sage. That library can be made fast if we included gmpy or
hack it to use our integers.
{{{
sage: import sympy.mpmath as a
sage: a.ei(complex(2r,3r))
mpc(real='-0.3615519445996403', imag='5.2705484358136943')
sage: a.mp.prec = 1000
sage:
a.ei('20.1000000000000000000000000000000000000000000000000000000000000000000')
mpf('28160453.3833994153950012661507902048168259435687244439573563929054452528467189574883919398422855919158391624626109323992762879680063235714851415880668470266487381859473426294489489391883476549497496312702303650499027252495570681212146125439331860229395375385300595969122490770779143241578042047266845835')
}}}
So this is an alternative approach to fixing this bug and greatly
improving
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/4807#comment:3>
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
-~----------~----~----~----~------~----~------~--~---
