On Sat, May 16, 2009 at 11:53 AM, Fredrik Johansson <fredrik.johans...@gmail.com> wrote: > >>> Do you think some of the functions should be rather implemented in mpmath? >> >> which functions do you mean exactly? The jn_zeros? I don't mind. But >> mpmath is sometimes very slow, e.g. I want to be able to use scipy >> optionally. It can be in mpmath, I have no problems with that. But >> then you would need the jn() function too in there, but that one uses >> sympy a lot, so I don't think it's that easy. > > Bessel function zeros, along with a couple of other inverse functions, > will benefit greatly from the more generic high-level interface > being added to mpmath. > > Basically, the first 16 digits can be computed entirely using machine > precision arithmetic, and then in the high precision case, mpmath can > use the machine precision estimate as the initial argument to > findroot.
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. 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. Ondrej --~--~---------~--~----~------------~-------~--~----~ 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 -~----------~----~----~----~------~----~------~--~---