#3426: [with new patch, needs review] bessel_K function is broken
-----------------------------------------------+----------------------------
Reporter: bober | Owner: gfurnish
Type: defect | Status: new
Priority: major | Milestone: sage-3.2
Component: calculus | Resolution:
Keywords: bessel, bessel_K, editor_gfurnish |
-----------------------------------------------+----------------------------
Comment (by ddrake):
Uh oh:
{{{
sage: bessel_J(0,0)
---------------------------------------------------------------------------
PariError Traceback (most recent call
last)
/home/drake/.sage/temp/klee/32521/_home_drake__sage_init_sage_0.py in
<module>()
----> 1
2
3
4
5
/opt/sage-3.1.3.alpha2/local/lib/python2.5/site-
packages/sage/functions/special.pyc in bessel_J(nu, z, algorithm, prec)
570 z = C(z)
571 K = C
--> 572 pari_bes = pari(nu).besselj(z, precision=prec)
573 if K is R:
574 return fudge*K(pari_bes.real())
/opt/sage-3.1.3.alpha2/local/lib/python2.5/site-
packages/sage/libs/pari/gen.so in sage.libs.pari.gen._pari_trap
(sage/libs/pari/gen.c:34414)()
7865
7866
-> 7867
7868
7869
PariError: (8)
}}}
Doing `bessel_J(0, 0)` in 3.1.2 works fine. I get similar errors with this
patch for other Bessel functions, too.
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/3426#comment:28>
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
-~----------~----~----~----~------~----~------~--~---
