http://www.netlib.org/bmp/

Multiple precision arithmetic

This is Richard Brent's classic MP package, which was published as TOMS algorithm 524. Contains some features such as the Bernoulli numbers, which are not ...
www.netlib.org/bmp/ - 2k - Cached - Similar pages

            This was found by google using a search for “multiple precision”

The exponent diffs

            In sminv, a check is made that diff s greater than 10-12smt to ensure that there are 4 digits left for the matrix inverse

16384**x=10**60

x ln 16384 = 60 ln 10

x = 60 ln(10)/ln(16384) = 14

            Item 2 in the multi-precision vector is the power of 16384.  To make the equivalent chek we need to see thaf diff(2) is > smt(2)-14 so that the subtraction loses only 50 digits.

            The original as downloaded is in ..\brent0.zip

Digits

            Real*8 uses 52 bits for the mantissa or 252=4503599627370496

                                                                              123456789abcdefgh ~ 17 digits

            Real*16 uses 52+64=116 bits for the mantissa or 2116=8.3076749736557242056487941267522e+34 ~ 34 digits

            MP uses 16 214 integers for the mantissa or 214*16=224=2.695994666715063979466701508702e+67 ~67 digits – mp setup asked for 64 digits.

Documents

07/10/2007  09:15a             328,192 MP USER Guide.doc This is an unchanged transformation of Brent’s buide into word format.

 

07/26/2007  12:28p              47,104 ..\Bob\BitsInFortran.doc .htm  This is a set of quotes from various helps about bit routines in Fortran extensions.

07/19/2007  06:24p             369,664 Example2a.doc htm  Contains example2a.zip, a working verision of a modification of the example given in Brent’s user guide. example.wpj – extends into Example 2b.doc .htm.

 

07/27/2007  01:01p              30,208 MPCDB.doc .htm This converts a double precision word to multiple precision with a series of left and right shifts.   The needed shifts are arithmetic statement functions.  The watfor version of this is tested in WFTMPCDBW.zip.

 

07/29/2007  02:57p              25,600 MPADD.doc .htm  Converts and adds in the same routine.  Does not normalize.

 

 

Wpj’s

07/12/2007  10:10p                 246 example.wpj

07/06/2007  12:38a                 238 test.wpj

07/18/2007  03:37p                 236 bif.wpj tests the make up of a double precision word

07/27/2007  04:30p                 245 example2.wpj

07/29/2007  02:44p                 239 mpadd.wpj

07/29/2007  05:54p                 238 mpmul.wpj

07/26/2007  01:56p                 239 titob.wpj

 

Zips

07/24/2007  07:57a               3,681 BITSINFORTRAN.zip     See BitsInFortran.doc  .htm

07/18/2007  03:38p               1,385 BIF.zip                             See BitsInFortran.doc  .htm

07/18/2007  10:14a                 771 OBIF.zip                           See BitsInFortran.doc  .htm

07/24/2007  08:00a                 669 itob.zip                               See BitsInFortran.doc  .htm

07/26/2007  10:28a               3,467 INLINE.zip                      See BitsInFortran.doc  .htm

07/26/2007  10:56a               3,194 WFTITOB.zip                  See BitsInFortran.doc  .htm

07/13/2007  02:57p              49,254 example2a.zip                  See Example2a.doc htm

07/19/2007  09:44a              55,084 example2b1.zip                See Example 2b.doc htm

07/19/2007  01:49p              55,600 example2b2.zip                See Example 2b.doc htm

07/28/2007  09:54a              12,753 EXAMPLE2c.zip             See MPCDB.doc

07/28/2007  09:54a              56,498 EXAMPLE2B3.zip          See Example 2b.doc htm

07/28/2007  09:53a              56,021 EXAMPLE2B4.zip

07/28/2007  09:54a              53,654 EXAMPLE2d.zip

07/31/2007  09:26a               8,146 wtitob.zip                          See MPCDB.doc htm

07/29/2007  02:56p              54,127 mpadd.zip                        See MPADD.doc htm

ISHFT locations

1.                  bif.for 07/18/2007  04:31p

2.                  ishft.for 07/23/2007  03:46p

3.                  ispan.for 07/22/2007  09:09a    

4.                  ispan2.for 07/22/2007  04:29p

5.                  mpaddb.for 07/29/2007  02:40p

6.                  mpaddbe.for 07/29/2007  02:22p ß

7.                  mpaddbw.for 08/01/2007  04:04p

8.                  mpcdbe.for 07/29/2007  01:36p

9.                  mpcdbe1.for  07/27/2007  09:10a

10.              mpcdbw.for 08/01/2007  03:04p

11.              tibset.for 07/23/2007  03:46p

12.              wfshfts.for 08/04/2007  02:17p