ALGORITHM 719 - MULTIPRECISION TRANSLATION AND EXECUTION OF FORTRAN PROGRAMS

This paper describes two Fortran utilities for multiprecision computat ion. The first is a package of Fortran subroutines that perform a vari ety of arithmetic operations and transcendental functions on floating point numbers of arbitrarily high precision. This package is in some c ases over 200 times faster than that of certain other packages that ha ve been developed for this purpose. The second utility is a translator program, which facilitates the conversion of ordinary Fortran program s to use this package. By means of source directives (special comments ) in the original Fortran program, the user declares the precision lev el and specifies which variables in each subprogram are to be treated as multiprecision. The translator program reads this source program an d outputs a program with the appropriate multiprecision subroutine cal ls. This translator supports multiprecision integer, real, and complex datatypes. The required array space for multiprecision data types is automatically allocated. In the evaluation of computational expression s, 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 mult iprecision values.