Thanks for your help. I tried your first and last suggestions, but
they yielded code that was slower than the original python
implementation. However, I'll take a look at sage.rings.real_mpfr and
try to use mpfr directly.
On Dec 22, 1:44 am, "William Stein" <wst...@gmail.com> wrote:
> On Sun, Dec 21, 2008 at 1:44 PM, M. Yurko <myu...@gmail.com> wrote:
>
> > I have recently been experimenting with converting some simple python
> > functions that I have made into cython. I have been quite impressed by
> > how simple it is for the massive speed increases that I have seen.
> > However, one thing that is mildly annoying at times is the limitation
> > to double precision computation. Is there any simple way to do
> > arbitrary precision floating-point arithmetic in cython? I tried using
> > importing mpfr, but I gave up after I saw the sheer number of
> > different things that seemed to need to be individually imported.
>
> I don't understand. Can't you just do, e.g.,
>
> from sage.all import RealNumber
> a = RealNumber('1.2939498029384028342983084203482093840283490823094829')
> # Now do stuff with a...
>
> Or do you want to completely avoid the Sage
> sage.rings.real_mpfr.RealNumber datatype and *directly* use the mpfr C
> library? If so, you might want to look at
> SAGE_ROOT/devel/sage/sage/rings/real_mpfr.pyx to see how
> sage.rings.real_mpfr.RealNumber is implemented. You could also do
>
> from sage.rings.real_mpfr cimport RealNumber
> import sage.all
> cdef RealNumber a = sage.all.RealNumber('1.29394980293840283429')
>
> Then in your code, a.value is a cdef's attribute of type mpfr_t.
>
> -- william
--~--~---------~--~----~------------~-------~--~----~
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to
sage-support-unsubscr...@googlegroups.com
For more options, visit this group at
http://groups.google.com/group/sage-support
URLs: http://www.sagemath.org
-~----------~----~----~----~------~----~------~--~---
