#4229: special functions should use mpfr when available
-------------------------+--------------------------------------------------
Reporter: AlexGhitza | Owner: burcin
Type: enhancement | Status: new
Priority: major | Milestone: sage-3.2
Component: calculus | Keywords:
-------------------------+--------------------------------------------------
MPFR has fast implementations for restricted types of arguments in some
special functions, e.g. Bessel J and Y with integer order and positive
real argument. We should be using these instead of Pari or Maxima or
Scipy whenever that is feasible.
Example:
{{{
sage: a = RR(2)
sage: timeit("bessel_J(1, a)")
625 loops, best of 3: 370 µs per loop
sage: timeit("a.j1()")
625 loops, best of 3: 13.9 µs per loop
}}}
That's 26 times faster than Pari.
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/4229>
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
-~----------~----~----~----~------~----~------~--~---
