[comp.unix.aix] Mathematical constants

AER7101%TECHNION@TAUNIVM.TAU.AC.IL (Zvika Bar-Deroma) (06/18/91)

Hi,
I'm installing  a subroutine library  and need to  know the values  of a
number of mathematical  constants. I would be grateful  if someone could
provide me with the constants for the  RS/6000 (AIX 3.1 w/ xlf) and also
for the PS/2 (AIX 1.2 w/ vs fortran).

Here are the definitions of the required constants:


     ASSUME INTEGERS ARE REPRESENTED IN THE N-DIGIT, BASE-A FORM
               SIGN ( X(N-1)*A**(N-1) + ... + X(1)*A + X(0) )
               WHERE 0 .LE. X(I) .LT. A FOR I=0,...,N-1.
     IPMPAR(1) = A, THE BASE.
     IPMPAR(2) = N, THE NUMBER OF BASE-A DIGITS.
     IPMPAR(3) = A**N - 1, THE LARGEST MAGNITUDE.

  FLOATING-POINT NUMBERS.

     IT IS ASSUMED THAT THE SINGLE AND DOUBLE PRECISION FLOATING
     POINT ARITHMETICS HAVE THE SAME BASE, SAY B, AND THAT THE
     NONZERO NUMBERS ARE REPRESENTED IN THE FORM

               SIGN (B**E) * (X(1)/B + ... + X(M)/B**M)

               WHERE X(I) = 0,1,...,B-1 FOR I=1,...,M,
               X(1) .GE. 1, AND EMIN .LE. E .LE. EMAX.

     IPMPAR(4) = B, THE BASE.

  SINGLE-PRECISION

     IPMPAR(5) = M, THE NUMBER OF BASE-B DIGITS.
     IPMPAR(6) = EMIN, THE SMALLEST EXPONENT E.
     IPMPAR(7) = EMAX, THE LARGEST EXPONENT E.

  DOUBLE-PRECISION

     IPMPAR(8) = M, THE NUMBER OF BASE-B DIGITS.
     IPMPAR(9) = EMIN, THE SMALLEST EXPONENT E.
     IPMPAR(10) = EMAX, THE LARGEST EXPONENT E.

Thanks,
/Zvika


Zvika Bar-Deroma                                  Phone: (+972)-4-292706
Faculty of Aerospace Engineering,                 Fax  : (+972)-4-231848
Technion
Haifa 32000
Israel

BITNET        :   AER7101@TECHNION
Internet      :   AER7101@TECHNION.TECHNION.AC.IL
UUCP          :   ...!uunet!pucc.princeton.edu!technion!aer7101