#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):
Now `bessel_J(0,0)` works but I'm seeing other problems. I'm concentrating
here on the "K" functions since that's what this ticket is about. This is
all with the current patch applied to 3.1.4.
* `bessel_K(0, -1)` doesn't work at all; it just soaks up all the CPU and
doesn't return. The correct answer is about `0.421024438240708 -
3.97746326050642*I`.
* `bessel_K(-1*I - 1, 0)` gives an uninformative Pari error; the function
isn't defined there. Mathematica says "ComplexInfinity"; can we give a
better error message? Even just "function not defined at 0" would be fine.
* `bessel_K(-1, -1)` says it can't convert
`-0.601907230197235-1.77549968921218*I` to a real number, which I agree
with. :) It looks like Pari has the right answer but we're trying to
convert it to a real when we shouldn't be.
* `bessel_K(0, -1 - I)` says "incorrect type", but Mathematica and Maple
evaluate it just fine.
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/3426#comment:31>
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
-~----------~----~----~----~------~----~------~--~---
