Article View: sci.math.num-analysis
Article #159Bailey's \*new\* M.P. Package (was: Brent's M.P.)
From: Frederick W. Chapman
Date: Fri, 31 Jul 1992 18:36
Date: Fri, 31 Jul 1992 18:36
120 lines
5480 bytes
5480 bytes
In article <Bs5p9K.76L@undergrad.math.waterloo.edu>,
drpatton@undergrad.math.waterloo.edu (Dave Patton) writes:
>I am attempting to use Brent's Multiple-precision library with a MIPS
>Fortran 77
>compiler. I would like to run a program which is fairly simple in nature but
>which requires very high precision for certain choices of parameters.
>
>Does anyone know how to run a simple program, which makes use of the AUGMENT
>preprocessor?
I have never had any luck even *finding* a copy of the AUGMENT
preprocessor much less using it!
>Also, I seem to be encountering numerous compile errors with the library
>routines. Is this as a result of machine-dependent code, or am I missing
>something more obvious?
...........................................................................
For anyone who may be interested, there is a *new* FORTRAN package for
doing multiple-precision floating-point and integer arithmetic; the
package includes an expression translator (program preprocessor). The
package was written by David H. Bailey, of the NASA Ames Research
Center. I have attached an announcement (below) that originally
appeared in the "NA Digest" electronic newsletter.
Fred Chapman
...........................................................................
NA Digest Sunday, April 26, 1992 Volume 92 : Issue 17
Today's Editor:
Cleve Moler
The MathWorks, Inc.
moler@mathworks.com
Today's Topics:
Multiprecision Computation Translator
Multivariate Optimization
SIAM Meeting Advance Registration Deadlines
Position at University of Vienna
Contents: Journal of Supercomputer Applications
Submissions for NA Digest:
Mail to na.digest@na-net.ornl.gov.
Information about NA-NET:
Mail to na.help@na-net.ornl.gov.
-------------------------------------------------------
From: David Bailey <dbailey@nas.nasa.gov>
Date: Tue, 21 Apr 92 10:46:38 -0700
Subject: Multiprecision Computation Translator
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, please send me a note. I eventually
intend to place this software on public media such as netlib, but for
the time being I prefer to keep a list of users, so that any bug fixes
or enhancements can be promptly distributed.
David H. Bailey
NASA Ames Research Center
dbailey@nas.nasa.gov
-------------------------------------------------------
--
o ------------------------------------------------------------------------- o
| Frederick W. Chapman, User Services, Computing Center, Lehigh University |
| Campus Phone: 8-3218 Preferred E-mail Address: fc03@Lehigh.Edu |
| "I do comedy and magic; what you don't find funny -- that's the magic." |
o ------------------------------------------------------------------------- o
Technical Details
Message-ID:
<1992Jul31.183619.39396@ns1.cc.lehigh.edu>Path:
arc1-ca.pugleaf.net!sparky!uunet!dtix!darwin.sura.net!jvnc.net!netnews.upenn.edu!netnews.cc.lehigh.edu!ns1.cc.lehigh.edu!fc03Show Full Headers
Complete Article Headers
Path: arc1-ca.pugleaf.net!sparky!uunet!dtix!darwin.sura.net!jvnc.net!netnews.upenn.edu!netnews.cc.lehigh.edu!ns1.cc.lehigh.edu!fc03 From: fc03@ns1.cc.lehigh.edu (Frederick W. Chapman) Newsgroups: sci.math.num-analysis Subject: Bailey's \*new\* M.P. Package (was: Brent's M.P.) Message-ID: <1992Jul31.183619.39396@ns1.cc.lehigh.edu> Date: 31 Jul 92 18:36:19 GMT Organization: Lehigh University Lines: 120