From seager@lll-zeus Thu Feb 12 13:27:32 1987
Received: from anl-mcs.ARPA by dasher.sun.com (3.2/SMI-3.2)
	id AA12212; Thu, 12 Feb 87 13:27:27 CST
Received: from lll-crg.ARPA (lll-crg.arpa.ARPA) by anl-mcs.ARPA (4.12/4.9)
	id AA29580; Wed, 11 Feb 87 13:06:05 cst
Received: Wed, 11 Feb 87 11:05:31 pst from lll-heracles.Arpa by lll-crg.ARPA (4.12/)
	id AA23988; Wed, 11 Feb 87 11:05:31 pst
Received: from lll-zeus.Arpa by lll-heracles.Arpa (1.1/SMI-3.0DEV3)
	id AA08304; Wed, 11 Feb 87 11:06:17 PST
Return-Path: <seager@lll-zeus>
Received: by lll-zeus.Arpa (1.1/SMI-3.0DEV3)
	id AA07597; Wed, 11 Feb 87 11:06:23 PST
Date: Wed, 11 Feb 87 11:06:23 PST
From: seager@lll-zeus (Mark Seager)
Message-Id: <8702111906.AA07597@lll-zeus.Arpa>
To: dongarra%dasher@anl-mcs.Arpa
In-Reply-To: Jack Dongarra's message of Wed, 11 Feb 87 09:48:43 CST <8702111548.AA09723@dasher.sun.com>
Subject: Mail addresses: Mine is seager@lll-crg.Arpa
Status: RO

Jack,
	Our system administrator has not put our sun workstation
addresses out on the arpanet.  
You can reach me at seager@lll-crg.Arpa.  Here are the machine
constants for the Suns in the form of i1mach.f, r1mach.f and d1mach.f
(in shar format).  Please let me know if you find any problems with
them !!! I will try to get the same numbers for the alliant... 
		Warmest Regards for you and your Family !!!
		Mark
======================================================================

# This is a shell archive.  Remove anything before this line,
# then unpack it by saving it in a file and typing "sh file".
#
# Wrapped by lll-zeus!seager on Wed Feb 11 11:03:43 PST 1987
# Contents:  d1mach.f i1mach.f r1mach.f
 
echo x - d1mach.f
sed 's/^@//' > "d1mach.f" <<'@//E*O*F d1mach.f//'
      double precision function d1mach(i)
c***begin prologue  d1mach
c***date written   750101   (yymmdd)
c***revision date  831014   (yymmdd)
c***category no.  r1
c***keywords  machine constants
c***author  fox, p. a., (bell labs)
c           hall, a. d., (bell labs)
c           schryer, n. l., (bell labs)
c***purpose  returns double precision machine dependent constants
c***description
c
c     d1mach can be used to obtain machine-dependent parameters
c     for the local machine environment.  it is a function
c     subprogram with one (input) argument, and can be called
c     as follows, for example
c
c          d = d1mach(i)
c
c     where i=1,...,5.  the (output) value of d above is
c     determined by the (input) value of i.  the results for
c     various values of i are discussed below.
c
c  double-precision machine constants
c  d1mach( 1) = b**(emin-1), the smallest positive magnitude.
c  d1mach( 2) = b**emax*(1 - b**(-t)), the largest magnitude.
c  d1mach( 3) = b**(-t), the smallest relative spacing.
c  d1mach( 4) = b**(1-t), the largest relative spacing.
c  d1mach( 5) = log10(b)
c***references  fox p.a., hall a.d., schryer n.l.,*framework for a
c                 portable library*, acm transactions on mathematical
c                 software, vol. 4, no. 2, june 1978, pp. 177-188.
c***routines called  xerror
c***end prologue  d1mach
c
      integer small(4)
      integer large(4)
      integer right(4)
      integer diver(4)
      integer log10(4)
c
      double precision dmach(5)
c
      equivalence (dmach(1),small(1))
      equivalence (dmach(2),large(1))
      equivalence (dmach(3),right(1))
      equivalence (dmach(4),diver(1))
      equivalence (dmach(5),log10(1))
c
c     machine constants for the sun
c
      data small / 2.22507385850720d-308 /
      data large / 1.79769313486232d+308 /
      data right / 1.1101827117665d-16 /
      data diver / 1.06021462628034d+292 /
      data log10 / 3.01029995663981e-1 /
c
c
c***first executable statement  d1mach
      if (i .lt. 1  .or.  i .gt. 5)
     1   call xerror( 'd1mach -- i out of bounds',25,1,2)
c
      d1mach = dmach(i)
      return
c
      end
@//E*O*F d1mach.f//
chmod u=rw,g=r,o=r d1mach.f
 
echo x - i1mach.f
sed 's/^@//' > "i1mach.f" <<'@//E*O*F i1mach.f//'
      integer function i1mach(i)
c***begin prologue  i1mach
c***date written   750101   (yymmdd)
c***revision date  860106   (yymmdd)
c***category no.  r1
c***keywords  machine constants
c***author  fox, p. a., (bell labs)
c           hall, a. d., (bell labs)
c           schryer, n. l., (bell labs)
c***purpose  return integer machine dependent constants.
c***description
c
c * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
c   these machine constant routines must be activated for
c   a particular environment.
c * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
c
c     i1mach can be used to obtain machine-dependent parameters
c     for the local machine environment.  it is a function
c     subroutine with one (input) argument, and can be called
c     as follows, for example
c
c          k = i1mach(i)
c
c     where i=1,...,16.  the (output) value of k above is
c     determined by the (input) value of i.  the results for
c     various values of i are discussed below.
c
c  i/o unit numbers.
c    i1mach( 1) = the standard input unit.
c    i1mach( 2) = the standard output unit.
c    i1mach( 3) = the standard punch unit.
c    i1mach( 4) = the standard error message unit.
c
c  words.
c    i1mach( 5) = the number of bits per integer storage unit.
c    i1mach( 6) = the number of characters per integer storage unit.
c
c  integers.
c    assume integers are represented in the s-digit, base-a form
c
c               sign ( x(s-1)*a**(s-1) + ... + x(1)*a + x(0) )
c
c               where 0 .le. x(i) .lt. a for i=0,...,s-1.
c    i1mach( 7) = a, the base.
c    i1mach( 8) = s, the number of base-a digits.
c    i1mach( 9) = a**s - 1, the largest magnitude.
c
c  floating-point numbers.
c    assume floating-point numbers are represented in the t-digit,
c    base-b form
c               sign (b**e)*( (x(1)/b) + ... + (x(t)/b**t) )
c
c               where 0 .le. x(i) .lt. b for i=1,...,t,
c               0 .lt. x(1), and emin .le. e .le. emax.
c    i1mach(10) = b, the base.
c
c  single-precision
c    i1mach(11) = t, the number of base-b digits.
c    i1mach(12) = emin, the smallest exponent e.
c    i1mach(13) = emax, the largest exponent e.
c
c  double-precision
c    i1mach(14) = t, the number of base-b digits.
c    i1mach(15) = emin, the smallest exponent e.
c    i1mach(16) = emax, the largest exponent e.
c
c  to alter this function for a particular environment,
c  the desired set of data statements should be activated by
c  removing the c from column 1.  also, the values of
c  i1mach(1) - i1mach(4) should be checked for consistency
c  with the local operating system.
c***references  fox p.a., hall a.d., schryer n.l.,*framework for a
c                 portable library*, acm transactions on mathematical
c                 software, vol. 4, no. 2, june 1978, pp. 177-188.
c***routines called  (none)
c***end prologue  i1mach
c
      integer imach(16),output
      equivalence (imach(4),output)
c
c
c     machine constants for the sun.
c
      data imach( 1) /     5 /
      data imach( 2) /     6 /
      data imach( 3) /     6 /
      data imach( 4) /     0 /
      data imach( 5) /    32 /
      data imach( 6) /     4 /
      data imach( 7) /     2 /
      data imach( 8) /    32 /
      data imach( 9) /2147483647/
      data imach(10) /     2 /
      data imach(11) /    23 /
      data imach(12) /   -38 /
      data imach(13) /    38 /
      data imach(14) /    52 /
      data imach(15) /  -308 /
      data imach(16) /   308 /
c
c
c***first executable statement  i1mach
      if (i .lt. 1  .or.  i .gt. 16) go to 10
c
      i1mach=imach(i)
      return
c
   10 continue
      write(output,9000)
9000  format('1error    1 in i1mach - i out of bounds ')
c
c     call fdump
c
c
      stop
      end
c      the routine i1mach has been edited.  the routine fdump
c      has been deleted from the original slatec release.
@//E*O*F i1mach.f//
chmod u=rw,g=r,o=r i1mach.f
 
echo x - r1mach.f
sed 's/^@//' > "r1mach.f" <<'@//E*O*F r1mach.f//'
      real function r1mach(i)
c***begin prologue  r1mach
c***date written   790101   (yymmdd)
c***revision date  831014   (yymmdd)
c***category no.  r1
c***keywords  machine constants
c***author  fox, p. a., (bell labs)
c           hall, a. d., (bell labs)
c           schryer, n. l., (bell labs)
c***purpose  returns single precision machine dependent constants
c***description
c
c     r1mach can be used to obtain machine-dependent parameters
c     for the local machine environment.  it is a function
c     subroutine with one (input) argument, and can be called
c     as follows, for example
c
c          a = r1mach(i)
c
c     where i=1,...,5.  the (output) value of a above is
c     determined by the (input) value of i.  the results for
c     various values of i are discussed below.
c
c  single-precision machine constants
c  r1mach(1) = b**(emin-1), the smallest positive magnitude.
c  r1mach(2) = b**emax*(1 - b**(-t)), the largest magnitude.
c  r1mach(3) = b**(-t), the smallest relative spacing.
c  r1mach(4) = b**(1-t), the largest relative spacing.
c  r1mach(5) = log10(b)
c***references  fox, p.a., hall, a.d., schryer, n.l, *framework for
c                 a portable library*, acm transactions on mathe-
c                 matical software, vol. 4, no. 2, june 1978,
c                 pp. 177-188.
c***routines called  xerror
c***end prologue  r1mach
c
      integer small(2)
      integer large(2)
      integer right(2)
      integer diver(2)
      integer log10(2)
c
      real rmach(5)
c
      equivalence (rmach(1),small(1))
      equivalence (rmach(2),large(1))
      equivalence (rmach(3),right(1))
      equivalence (rmach(4),diver(1))
      equivalence (rmach(5),log10(1))
c
c
c     machine constants for the sun
c
      data rmach(1) / 1.17549435e-38 /
      data rmach(2) / 3.40282347e+38 /
      data rmach(3) / 5.96016605e-8  /
      data rmach(4) / 1.10795077e+31 /
      data rmach(5) / 3.0103001e-1   /
c
c
c***first executable statement  r1mach
      if (i .lt. 1  .or.  i .gt. 5)
     1   call xerror ( 'r1mach -- i out of bounds',25,1,2)
c
      r1mach = rmach(i)
      return
c
      end
@//E*O*F r1mach.f//
chmod u=rw,g=r,o=r r1mach.f
 
exit 0


.