March 29 2007   3:58:42.317 AM
 
MACHINE_PRB:
  FORTRAN90 version
  Test the routines in the MACHINE library,
  which returns the value of constants
  associated with computer arithmetic.
 
D1MACHINE_PRB
  D1MACH reports the value of constants associated
  with real double precision computer arithmetic.
 
  Assume that double precision numbers are stored 
  with a mantissa of T digits in base B, with an 
  exponent whose value must lie between EMIN and EMAX.
 
  For input arguments of 1 <= I <= 5,
  D1MACH will return the following values:
 
  D1MACH(1) = B**(EMIN-1), the smallest positive magnitude.
  2.225073858507201E-308
 
  D1MACH(2) = B**EMAX*(1-B**(-T)), the largest magnitude.
  1.797693134862316E+308
 
  D1MACH(3) = B**(-T), the smallest relative spacing.
  1.110223024625157E-016
 
  D1MACH(4) = B**(1-T), the largest relative spacing.
  2.220446049250313E-016
 
  D1MACH(5) = log10(B).
  0.301029995663981     
 
I1MACHINE_PRB
  I1MACH reports the value of constants associated
  with integer computer arithmetic.
 
  Numbers associated with input/output units:
 
  I1MACH(1) = the standard input unit.
           5
 
  I1MACH(2) = the standard output unit.
           6
 
  I1MACH(3) = the standard punch unit.
           6
 
  I1MACH(4) = the standard error message unit.
           6
 
  Numbers associated with words:
 
  I1MACH(5) = the number of bits per integer.
          32
 
  I1MACH(6) = the number of characters per integer.
           4
 
  Numbers associated with integer values:
 
  Assume integers are represented in the S digit 
  base A form:
 
    Sign * (X(S-1)*A**(S-1) + ... + X(1)*A + X(0))
 
  where the digits X satisfy 0 <= X(1:S-1) < A.
 
  I1MACH(7) = A, the base.
           2
 
  I1MACH(8) = S, the number of base A digits.
          31
 
  I1MACH(9) = A**S-1, the largest integer.
  2147483647
 
  Numbers associated with floating point values:
 
  Assume floating point numbers are represented 
  in the T digit base B form:
 
    Sign * (B**E) * ((X(1)/B) + ... + (X(T)/B**T) )
 
  where 
 
    0 <= X(1:T) < B,
    0 < X(1) (unless the value being represented is 0),
    EMIN <= E <= EMAX.
 
  I1MACH(10) = B, the base.
           2
 
  Numbers associated with single precision values:
 
  I1MACH(11) = T, the number of base B digits.
          24
 
  I1MACH(12) = EMIN, the smallest exponent E.
        -125
 
  I1MACH(13) = EMAX, the largest exponent E.
         128
 
  Numbers associated with double precision values:
 
  I1MACH(14) = T, the number of base B digits.
          53
 
  I1MACH(15) = EMIN, the smallest exponent E.
       -1021
 
  I1MACH(16) = EMAX, the largest exponent E.
        1024
 
R1MACHINE_PRB
  D1MACH reports the value of constants associated
  with real single precision computer arithmetic.
 
  Assume that single precision numbers are stored 
  with a mantissa of T digits in base B, with an 
  exponent whose value must lie between EMIN and EMAX.
 
  For input arguments of 1 <= I <= 5,
  R1MACH will return the following values:
 
  R1MACH(1) = B**(EMIN-1), the smallest positive magnitude.
  1.1754944E-38
 
  R1MACH(2) = B**EMAX*(1-B**(-T)), the largest magnitude.
  3.4028235E+38
 
  R1MACH(3) = B**(-T), the smallest relative spacing.
  5.9604645E-08
 
  R1MACH(4) = B**(1-T), the largest relative spacing.
  1.1920929E-07
 
  R1MACH(5) = log10(B).
  0.3010300    
 
MACHINE_PRB:
  Normal end of execution.
 
March 29 2007   3:58:42.322 AM