1
DOUBLE PRECISION FUNCTION spmpar(i)
2
C-----------------------------------------------------------------------
3
C SPMPAR PROVIDES THE SINGLE PRECISION MACHINE CONSTANTS FOR
4
C THE COMPUTER BEING USED. IT IS ASSUMED THAT THE ARGUMENT
5
C I IS AN INTEGER HAVING ONE OF THE VALUES 1, 2, OR 3. IF THE
6
C SINGLE PRECISION ARITHMETIC BEING USED HAS M BASE B DIGITS AND
7
C ITS SMALLEST AND LARGEST EXPONENTS ARE EMIN AND EMAX, THEN
9
C SPMPAR(1) = B**(1 - M), THE MACHINE PRECISION,
11
C SPMPAR(2) = B**(EMIN - 1), THE SMALLEST MAGNITUDE,
13
C SPMPAR(3) = B**EMAX*(1 - B**(-M)), THE LARGEST MAGNITUDE.
15
C-----------------------------------------------------------------------
16
C RWRITTEN BY JPC to use lapack dlamch
17
C-----------------------------------------------------------------------
18
DOUBLE PRECISION DLAMCH