This library offers a Fortran binding to the GNU MPFR Library multi-precision library. The aim of the library is to make the use of MPFR convenient for a Fortran user.

This is an initial work, expect bugs and (especially) missing features and possibly some instability, at least in binary compatibility. Source code compatibility should be fine.

As of now, you need a compiler which supports `ERROR STOP for the testsuite.

• Arithmetic expressions can be used the usual way; expressions can be mixed with real and integer expressions.

• Mathematical intrinsics, like sin, cos etc.

• Array expressions

• Setting and changing a default precision of variables

• List-directed output

• Reading values from strings

• Input (user-defined I/O is too unevenly implemented in existing compilers)

• Array intrinsics like MAXLOC and MAXVAL

Here is a short exampe, if your compiler supports user-defined I/O:

program memain
  use fmpfr_oper
  implicit none
  type (fmpfr), dimension(:), allocatable :: a
  call set_default_prec (128)
  allocate (a(2))
  a(1) = fmpfr("1.3")
  a(2) = a(1) + 2
  print *,a
  print *,sin(a)
  print *,get_str(cos(a(1)),10)
end program memain

whose output is

 1.299999999999999999999999999999999999998E0  3.299999999999999999999999999999999999998E0 
 9.635581854171929647013486300395548153418E-1  -1.577456941432483820116542776024823708430E-1 
 2.674988286E-1

Type ./configure and make. If everything was successful, type make check.

If you want to modify the code, it is best if you make your changes to the src/generate.pl file and run

$ grep ^@deftypefun mfpr.texi | perl generate.pl

where mpfr.texi is the MPFR documentation.

This is covered under the MIT license.