#5175: pari error is thrown when computing Bessel function J0(0)
---------------------+------------------------------------------------------
Reporter: pdenapo | Owner: tbd
Type: defect | Status: new
Priority: major | Milestone: sage-3.4.1
Component: algebra | Resolution:
Keywords: |
---------------------+------------------------------------------------------
Changes (by mabshoff):
* summary: [bug] Bug computing Bessel function J0(0) => pari error is
thrown when computing Bessel function J0(0)
* milestone: => sage-3.4.1
Old description:
> A bug under sage 3.2.3 computing J_0(0)
> (which should be just "1")
>
> sage: g=Bessel(0)
> sage: g
> J_{0}
> sage: g(1)
> 0.765197686557967
> sage: g(0)
> ---------------------------------------------------------------------------
> PariError Traceback (most recent call
> last)
>
> /home/pablo/sage/sage-3.1/<ipython console> in <module>()
>
> /home/pablo/sage/sage-3.1/local/lib/python2.5/site-
> packages/sage/functions/special.pyc in __call__(self, z)
> 689 return bessel_I(nu,z,algorithm=s,prec=p)
> 690 if t == "J":
> --> 691 return bessel_J(nu,z,algorithm=s,prec=p)
> 692 if t == "K":
> 693 return bessel_K(nu,z,algorithm=s,prec=p)
>
> /home/pablo/sage/sage-3.1/local/lib/python2.5/site-
> packages/sage/functions/special.pyc in bessel_J(nu, z, algorithm, prec)
> 522 K = C
> 523 K = z.parent()
> --> 524 return K(pari(nu).besselj(z, precision=prec))
> 525 elif algorithm=="scipy":
> 526 if prec != 53:
>
> /home/pablo/sage/sage-3.1/local/lib/python2.5/site-
> packages/sage/libs/pari/gen.so in sage.libs.pari.gen._pari_trap
> (sage/libs/pari/gen.c:38645)()
>
> PariError: (8)
New description:
A bug under sage 3.2.3 computing J_0(0)
(which should be just "1")
{{{
sage: g=Bessel(0)
sage: g
J_{0}
sage: g(1)
0.765197686557967
sage: g(0)
---------------------------------------------------------------------------
PariError Traceback (most recent call
last)
/home/pablo/sage/sage-3.1/<ipython console> in <module>()
/home/pablo/sage/sage-3.1/local/lib/python2.5/site-
packages/sage/functions/special.pyc in __call__(self, z)
689 return bessel_I(nu,z,algorithm=s,prec=p)
690 if t == "J":
--> 691 return bessel_J(nu,z,algorithm=s,prec=p)
692 if t == "K":
693 return bessel_K(nu,z,algorithm=s,prec=p)
/home/pablo/sage/sage-3.1/local/lib/python2.5/site-
packages/sage/functions/special.pyc in bessel_J(nu, z, algorithm, prec)
522 K = C
523 K = z.parent()
--> 524 return K(pari(nu).besselj(z, precision=prec))
525 elif algorithm=="scipy":
526 if prec != 53:
/home/pablo/sage/sage-3.1/local/lib/python2.5/site-
packages/sage/libs/pari/gen.so in sage.libs.pari.gen._pari_trap
(sage/libs/pari/gen.c:38645)()
PariError: (8)
}}}}
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/5175#comment:1>
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
-~----------~----~----~----~------~----~------~--~---
