#10239: Bessel method forces bessel_Y to use pari algorithm
-------------------------+--------------------------------------------------
Reporter: jvkersch | Owner: jason, jkantor
Type: defect | Status: new
Priority: major | Milestone: sage-4.6.1
Component: numerical | Keywords: functions, bessel, pari
Author: | Upstream: N/A
Reviewer: | Merged:
Work_issues: |
-------------------------+--------------------------------------------------
Comment(by jvkersch):
Too bad that there's no pari implementation: for Bessel-J it is just so
fast:
{{{
sage: timeit(r"bessel_J(1, 2.3, algorithm='pari')")
625 loops, best of 3: 110 µs per loop
sage: timeit(r"bessel_J(1, 2.3, algorithm='scipy')")
625 loops, best of 3: 480 µs per loop
sage: timeit(r"bessel_J(1, 2.3, algorithm='maxima')")
125 loops, best of 3: 5.27 ms per loop
}}}
Still, using the maxima implementation would be good for now, as the docs
mention that there are some accuracy issues with scipy. Also, thanks for
your quick response -- it feels like I never left Ghent and we are still
talking about Sage in the backyard of the department!
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/10239#comment:2>
Sage <http://www.sagemath.org>
Sage: Creating a Viable Open Source Alternative to Magma, Maple, Mathematica,
and MATLAB
--
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.
