|
1
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 |
PROGRAM TEST3
* * -- LAPACK test routine (version 3.2) -- * Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. * November 2006 * * .. Local Scalars .. DOUBLE PRECISION BASE, EMAX, EMIN, EPS, PREC, RMAX, RMIN, RND, $ SFMIN, T * .. * .. External Functions .. DOUBLE PRECISION DLAMCH EXTERNAL DLAMCH * .. * .. Executable Statements .. * EPS = DLAMCH( 'Epsilon' ) SFMIN = DLAMCH( 'Safe minimum' ) BASE = DLAMCH( 'Base' ) PREC = DLAMCH( 'Precision' ) T = DLAMCH( 'Number of digits in mantissa' ) RND = DLAMCH( 'Rounding mode' ) EMIN = DLAMCH( 'Minimum exponent' ) RMIN = DLAMCH( 'Underflow threshold' ) EMAX = DLAMCH( 'Largest exponent' ) RMAX = DLAMCH( 'Overflow threshold' ) * WRITE( 6, * )' Epsilon = ', EPS WRITE( 6, * )' Safe minimum = ', SFMIN WRITE( 6, * )' Base = ', BASE WRITE( 6, * )' Precision = ', PREC WRITE( 6, * )' Number of digits in mantissa = ', T WRITE( 6, * )' Rounding mode = ', RND WRITE( 6, * )' Minimum exponent = ', EMIN WRITE( 6, * )' Underflow threshold = ', RMIN WRITE( 6, * )' Largest exponent = ', EMAX WRITE( 6, * )' Overflow threshold = ', RMAX WRITE( 6, * )' Reciprocal of safe minimum = ', 1 / SFMIN * END |