On Sat, May 16, 2009 at 11:56 AM, Ondrej Certik <ond...@certik.cz> wrote: > That'd be cool -- but it will still be slower than scipy -- sometimes > I just need it fast, but it's ok that maybe the 10th digit is wrong. > So I want to have both.
Some mpmath functions have turned out to be faster than their scipy counterparts when I reimplemented them for float. I don't know what the case will be for Bessel functions. > How will you implement jn()? One can use the same thing as in scipy, > but I find it pretty cool that sympy can now (I mean after the patches > are accepted) find an exact formula for it. Numerically, there's no point implementing spherical Bessel functions using any other formula than sqrt(pi/2/z)*besselj(v+0.5,z), except adding some additional code to handle special cases. By the way, jn() denotes an ordinary Bessel function in C99 math.h and scipy. You might want to change the names. Fredrik --~--~---------~--~----~------------~-------~--~----~ You received this message because you are subscribed to the Google Groups "sympy-patches" group. To post to this group, send email to sympy-patches@googlegroups.com To unsubscribe from this group, send email to sympy-patches+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sympy-patches?hl=en -~----------~----~----~----~------~----~------~--~---