Page: 22 prev next

Extended Precision Calculations

D.H. Bailey's multiple-precision FORTRAN preprocessor and subroutine library (available from NETLIB) was used for all calculations.

Standard FORTRAN code was modified by defining certain variables to be multiple-precision using specially-formatted FORTRAN comments. The number of digits of precision desired was also specified.

The preprocessor program then translated the "ordinary" program to a multiple-precision one by generating subroutine calls and other FORTRAN code from the multiple precision directives.

The source code output by the translation program was compiled with a standard FORTRAN compiler, and linked with the subroutine library.