#4626: [with fix, needs patch] error in bessel_J(0,0)
------------------------------+---------------------------------------------
Reporter: zimmerma | Owner: somebody
Type: defect | Status: new
Priority: major | Milestone: sage-3.4.1
Component: basic arithmetic | Resolution:
Keywords: |
------------------------------+---------------------------------------------
Changes (by rlm):
* summary: error in bessel_J(0,0) => [with fix, needs patch] error in
bessel_J(0,0)
Comment:
Aha!:
{{{
sage: pari('besselj(0,0)')
1.00000000000000
sage: pari('besselj(0.0,0.0)')
---------------------------------------------------------------------------
PariError Traceback (most recent call
last)
/Users/rlmill/<ipython console> in <module>()
/Users/rlmill/sage-3.2.2/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)
sage: pari('besselj(0,0.0)')
1.00000000000000
}}}
The problem is that Pari is expecting an integer:
{{{
sage: R = RealField(53)
sage: n = Integer(0)
sage: z = R(0)
sage: pari(n).besselj(z, precision=53)
1.00000000000000
}}}
I'll post a patch once my 3.3.alpha0 builds...
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/4626#comment:2>
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
-~----------~----~----~----~------~----~------~--~---
