#10239: Bessel method forces bessel_Y to use pari algorithm
-------------------------+--------------------------------------------------
Reporter: jvkersch | Owner: jason, jkantor
Type: defect | Status: needs_review
Priority: major | Milestone: sage-4.6.1
Component: numerical | Keywords: functions, bessel, pari
Author: | Upstream: N/A
Reviewer: | Merged:
Work_issues: |
-------------------------+--------------------------------------------------
Changes (by jvkersch):
* status: needs_work => needs_review
Comment:
The problem with scipy is the way in which the argument is presented to
the Bessel function: when the argument is real, scipy returns a real value
for the Bessel function:
{{{
sage: import scipy.special
sage: scipy.special.jv(6, complex(float(pi)))
(0.014545966982505558+5.3441015716885544e-18j)
sage: scipy.special.jv(6, float(pi))
0.014545966982505555
}}}
For now I would suggest adding an ellipsis to the doctest to obscure the
complex part (e.g write the output as `0.0145459669825000...`). I've
modified the patch to take this into account.
As everything else works fine, I would change the status to positive
review but I want to hear your opinion about the modification first.
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/10239#comment:4>
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.
