Re: arbitrary precision

On 09/26/2011 03:52 PM, Lurkos wrote:
I'm right now investigating the possibility of using arbitrary precision
arithmetic with Fortran.

I've found three libraries:
+ ARPREC (http://crd.lbl.gov/~dhbailey/mpdist/)
+ MPFUN90 (http://crd.lbl.gov/~dhbailey/mpdist/)
+ FMLIB (http://myweb.lmu.edu/dmsmith/FMLIB.html)

MPFUN90 seems the easiest to use, while ARPREC might supersede it
(although it needs g++ post compilation and redefinition of the main).
FMLIB supports more intrinsic function, but seems more complicated.

Do you have any experience with this libraries and give me an thoughts
about them?


Don't forget the GNU projects: GMP/MPFR/MPC

These are C libraries, but since they are requirements of the gcc compiler suite, they may actually be automatically available when you are using gfortran (some gfortran expert may correct me).

By the way, what is the standard return value of HUGE, EPSILON etc. if the KIND happens to be arbitrary-digits/arbitrary-precision ...? (I guess that's just theory right now)

-- Wolfgang


--
E-mail: firstnameinitial.lastname@xxxxxxxxx
Domain: yahoo
.