🚀 go-pugleaf

Some of you may be aware of a multiprecision computation package and a
multiprecision Fortran translator that I have developed during the
past year or two.  I am now pleased to announce that this software has
been approved by NASA Ames for unlimited international distribution,
at no cost.  A brief description of this software is as follows:

The first package contains approximately 85 Fortran subroutines that
perform a variety of arithmetic operations and transcendental
functions.  The basic arithmetic routines are very efficient, and for
extra high levels of precision (about about 1000 digits), some
additional routines are available that employ advanced algorithms,
such FFT-based multiplication and the Borwein transcendental
algorithms, for additional speedup.  On some problems these routines
are over 200 times faster than certain other multiprecision tools.

However, as with other such packages, manually converting a program to
use these routines is a tedious and error-prone process.  To
facilitate such conversions, I have now developed a multiprecision
translator for Fortran programs.  By means of source directives
(special comments) placed in a standard Fortran-77 program, the user
declares the precision level and specifies which variables in each
subprogram are to be treated as multiprecision.  The translator
program reads this source program and outputs a program with the
appropriate multiprecision subroutine calls, which can then be linked
with the multiprecision package and executed.

This translator program supports multiprecision integer, real and
complex datatypes.  The required array space for multiprecision data
types is automatically allocated.  In the evaluation of computational
expressions, all of the usual conventions for operator precedence and
mixed mode operations are upheld.  Furthermore, most of the Fortran-77
intrinsics, such as ABS, MOD, NINT, COS, EXP are supported and produce
true multiprecision values.

Both the multiprecision package and translator are written in highly
portable Fortran-77 and have been successfully implemented on a
variety of systems, ranging from parallel supercomputers to
workstations.

If anyone is interested in obtaining this software, plus the
accompanying documentation, send me a note, and I will respond by
e-mail.  I eventually intend to place this software on public media
such as netlib, but for the time being I am handling it myself, so
that any bug fixes or enhancements can be promptly distributed.  Since
it is a small but significant amount of trouble for me to respond,
please don't request unless you are somewhat serious about needing
this software.

David H. Bailey 
NASA Ames Research Center 
dbailey@nas.nasa.gov