LAPACK 3.3.1
Linear Algebra PACKage
|
00001 PROGRAM TEST3 00002 * 00003 * -- LAPACK test routine (version 3.2) -- 00004 * Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. 00005 * November 2006 00006 * 00007 * .. Local Scalars .. 00008 DOUBLE PRECISION BASE, EMAX, EMIN, EPS, PREC, RMAX, RMIN, RND, 00009 $ SFMIN, T 00010 * .. 00011 * .. External Functions .. 00012 DOUBLE PRECISION DLAMCH 00013 EXTERNAL DLAMCH 00014 * .. 00015 * .. Executable Statements .. 00016 * 00017 EPS = DLAMCH( 'Epsilon' ) 00018 SFMIN = DLAMCH( 'Safe minimum' ) 00019 BASE = DLAMCH( 'Base' ) 00020 PREC = DLAMCH( 'Precision' ) 00021 T = DLAMCH( 'Number of digits in mantissa' ) 00022 RND = DLAMCH( 'Rounding mode' ) 00023 EMIN = DLAMCH( 'Minimum exponent' ) 00024 RMIN = DLAMCH( 'Underflow threshold' ) 00025 EMAX = DLAMCH( 'Largest exponent' ) 00026 RMAX = DLAMCH( 'Overflow threshold' ) 00027 * 00028 WRITE( 6, * )' Epsilon = ', EPS 00029 WRITE( 6, * )' Safe minimum = ', SFMIN 00030 WRITE( 6, * )' Base = ', BASE 00031 WRITE( 6, * )' Precision = ', PREC 00032 WRITE( 6, * )' Number of digits in mantissa = ', T 00033 WRITE( 6, * )' Rounding mode = ', RND 00034 WRITE( 6, * )' Minimum exponent = ', EMIN 00035 WRITE( 6, * )' Underflow threshold = ', RMIN 00036 WRITE( 6, * )' Largest exponent = ', EMAX 00037 WRITE( 6, * )' Overflow threshold = ', RMAX 00038 WRITE( 6, * )' Reciprocal of safe minimum = ', 1 / SFMIN 00039 * 00040 END