4.3BSD version dated 09/12/85.
authorGroup working with Kahn on math lib <elefunt@ucbvax.Berkeley.EDU>
Fri, 13 Sep 1985 11:22:16 +0000 (03:22 -0800)
committerGroup working with Kahn on math lib <elefunt@ucbvax.Berkeley.EDU>
Fri, 13 Sep 1985 11:22:16 +0000 (03:22 -0800)
SCCS-vsn: lib/libm/README 1.3
SCCS-vsn: lib/libm/vax/Makefile 1.2
SCCS-vsn: lib/libm/vax/argred.s 1.2
SCCS-vsn: lib/libm/vax/atan2.s 1.2
SCCS-vsn: lib/libm/vax/cabs.s 1.2
SCCS-vsn: lib/libm/vax/cbrt.s 1.2
SCCS-vsn: lib/libm/vax/infnan.s 1.2
SCCS-vsn: lib/libm/vax/sincos.s 1.2
SCCS-vsn: lib/libm/vax/sqrt.s 1.2
SCCS-vsn: lib/libm/vax/support.s 1.2
SCCS-vsn: lib/libm/vax/tan.s 1.2
SCCS-vsn: lib/libm/vax/mcount.sed 1.2
SCCS-vsn: lib/libm/ieee/Makefile 1.2
SCCS-vsn: lib/libm/common/atan2.c 1.2
SCCS-vsn: lib/libm/ieee/cabs.c 1.2
SCCS-vsn: lib/libm/ieee/cbrt.c 1.2
SCCS-vsn: lib/libm/ieee/support.c 1.2
SCCS-vsn: lib/libm/common/trig.c 1.2

17 files changed:
usr/src/lib/libm/README
usr/src/lib/libm/common/atan2.c
usr/src/lib/libm/common/trig.c
usr/src/lib/libm/ieee/Makefile
usr/src/lib/libm/ieee/cabs.c
usr/src/lib/libm/ieee/cbrt.c
usr/src/lib/libm/ieee/support.c
usr/src/lib/libm/vax/Makefile
usr/src/lib/libm/vax/argred.s
usr/src/lib/libm/vax/atan2.s
usr/src/lib/libm/vax/cabs.s
usr/src/lib/libm/vax/cbrt.s
usr/src/lib/libm/vax/infnan.s
usr/src/lib/libm/vax/sincos.s
usr/src/lib/libm/vax/sqrt.s
usr/src/lib/libm/vax/support.s
usr/src/lib/libm/vax/tan.s

index d06727d..b3e4cb9 100644 (file)
@@ -1,22 +1,38 @@
-# @(#)README   1.3 (ELEFUNT) %G%
--1.  The machine-independent Version 7 math library found in 4.2BSD
-     is now "/usr/lib/libom.a".  To compile with those routines use -lom.
+***************************************************************************
+*                                                                         * 
+* Copyright (c) 1985 Regents of the University of California.             *
+*                                                                         * 
+* Use and reproduction of this software are granted  in  accordance  with *
+* the terms and conditions specified in  the  Berkeley  Software  License *
+* Agreement (in particular, this entails acknowledgement of the programs' *
+* source, and inclusion of this notice) with the additional understanding *
+* that  all  recipients  should regard themselves as participants  in  an *
+* ongoing  research  project and hence should  feel  obligated  to report *
+* their  experiences (good or bad) with these elementary function  codes, *
+* using "sendbug 4bsd-bugs@BERKELEY", to the authors.                     *
+*                                                                         *
+* K.C. Ng, with Z-S. Alex Liu, S. McDonald, P. Tang, W. Kahan.            *
+* Revised on 5/10/85, 5/13/85, 6/14/85, 8/20/85, 8/27/85, 9/11/85.        *
+*                                                                         *
+***************************************************************************
 
 
-K.C. Ng, March 7, 1985, with Z-S. Alex Liu, S. McDonald, P. Tang, W. Kahan.
-Revised on 5/10/85, 5/13/85, 6/14/85, 8/20/85, 8/27/85.
+/*     @(#)README      1.6 (Berkeley) 9/12/85; 1.3 (ucb.elefunt) %G% */
+
+NB.  The machine-independent Version 7 math library found in 4.2BSD
+     is now /usr/lib/libom.a.  To compile with those routines use -lom.
 
 ******************************************************************************
 *  This is a description of the upgraded elementary functions (listed in 1). *
 *  Bessel functions (j0, j1, jn, y0, y1, yn), floor, and fabs passed over    *
 *  from 4.2BSD without change except perhaps for the way floating point      *
 
 ******************************************************************************
 *  This is a description of the upgraded elementary functions (listed in 1). *
 *  Bessel functions (j0, j1, jn, y0, y1, yn), floor, and fabs passed over    *
 *  from 4.2BSD without change except perhaps for the way floating point      *
-*  exception is signaled on a VAX.  Three lines that contain "errno in erf.c *
-*  (error function erf, erfc) have been deleted to prevent overriding the    *
+*  exception is signaled on a VAX.  Three lines that contain "errno" in erf.c*
+*  (error functions erf, erfc) have been deleted to prevent overriding the   *
 *  system "errno".                                                           *
 ******************************************************************************
 
 0. Total number of files: 40
 
 *  system "errno".                                                           *
 ******************************************************************************
 
 0. Total number of files: 40
 
-        IEEE/Makefile   VAX/Makefile    VAX/support.s   erf.c       lgama.c
+        IEEE/Makefile   VAX/Makefile    VAX/support.s   erf.c       lgamma.c
         IEEE/atan2.c    VAX/argred.s    VAX/tan.s       exp.c       log.c
         IEEE/cabs.c     VAX/atan2.s     acosh.c         exp__E.c    log10.c
         IEEE/cbrt.c     VAX/cabs.s      asincos.c       expm1.c     log1p.c
         IEEE/atan2.c    VAX/argred.s    VAX/tan.s       exp.c       log.c
         IEEE/cabs.c     VAX/atan2.s     acosh.c         exp__E.c    log10.c
         IEEE/cbrt.c     VAX/cabs.s      asincos.c       expm1.c     log1p.c
@@ -25,47 +41,44 @@ Revised on 5/10/85, 5/13/85, 6/14/85, 8/20/85, 8/27/85.
         Makefile        VAX/sincos.s    atanh.c         j1.c        sinh.c
         README          VAX/sqrt.s      cosh.c          jn.c        tanh.c
 
         Makefile        VAX/sincos.s    atanh.c         j1.c        sinh.c
         README          VAX/sqrt.s      cosh.c          jn.c        tanh.c
 
-1. Functions implemented:
-    (A). Standard elementary functions (total 22):
-        acos(x)                 ... in file  "asincos.c"
-        asin(x)                 ... in file  "asincos.c"
-        atan(x)                 ... in file  "atan.c"
-        atan2(x,y)              ... in files "IEEE/atan2.c", "VAX/atan2.s"
-        sin(x)                  ... in files "IEEE/trig.c",  "VAX/sincos.s"
-        cos(x)                  ... in files "IEEE/trig.c",  "VAX/sincos.s"
-        tan(x)                  ... in files "IEEE/trig.c",  "VAX/tan.s"
-        cabs(x,y)               ... in files "IEEE/cabs.c",  "VAX/cabs.s"
-        hypot(x,y)              ... in files "IEEE/cabs.c",  "VAX/cabs.s"
-        cbrt(x)                 ... in files "IEEE/cbrt.c",  "VAX/cbrt.s"
-        exp(x)                  ... in file  "exp.c"
-        expm1(x):=exp(x)-1      ... in file  "expm1.c"
-        log(x)                  ... in file  "log.c"
-        log10(x)                ... in file  "log10.c"
-        log1p(x):=log(1+x)      ... in file  "log1p.c"
-        pow(x,y)                ... in file  "pow.c"
-        sinh(x)                 ... in file  "sinh.c"
-        cosh(x)                 ... in file  "cosh.c"
-        tanh(x)                 ... in file  "tanh.c"
-        asinh(x)                ... in file  "asinh.c"
-        acosh(x)                ... in file  "acosh.c"
-        atanh(x)                ... in file  "atanh.c"
-
+1. Functions implemented :
+    (A). Standard elementary functions (total 22) :
+        acos(x)                 ...in file  asincos.c 
+        asin(x)                 ...in file  asincos.c
+        atan(x)                 ...in file  atan.c
+        atan2(x,y)              ...in files IEEE/atan2.c, VAX/atan2.s
+        sin(x)                  ...in files IEEE/trig.c,  VAX/sincos.s
+        cos(x)                  ...in files IEEE/trig.c,  VAX/sincos.s
+        tan(x)                  ...in files IEEE/trig.c,  VAX/tan.s
+        cabs(x,y)               ...in files IEEE/cabs.c,  VAX/cabs.s
+        hypot(x,y)              ...in files IEEE/cabs.c,  VAX/cabs.s
+        cbrt(x)                 ...in files IEEE/cbrt.c,  VAX/cbrt.s
+        exp(x)                  ...in file  exp.c
+        expm1(x):=exp(x)-1      ...in file  expm1.c
+        log(x)                  ...in file  log.c
+        log10(x)                ...in file  log10.c
+        log1p(x):=log(1+x)      ...in file  log1p.c
+        pow(x,y)                ...in file  pow.c
+        sinh(x)                 ...in file  sinh.c
+        cosh(x)                 ...in file  cosh.c
+        tanh(x)                 ...in file  tanh.c
+        asinh(x)                ...in file  asinh.c
+        acosh(x)                ...in file  acosh.c
+        atanh(x)                ...in file  atanh.c
+                        
     (B). Kernel functions :
     (B). Kernel functions :
-        exp__E(x,c) ... in file "exp__E.c", used by
-                       expm1(), exp(), pow() and cosh()
-        log__L(s)   ... in file "log__L.c", used by
-                       log1p(), log() and pow()
-        libm$argred ... in file "VAX/argred.s", used by VAX version of
-                       sin(), cos() and tan()
+        exp__E(x,c) ...in file exp__E.c, used by expm1/exp/pow/cosh
+        log__L(s)   ...in file log__L.c, used by log1p/log/pow
+        libm$argred ...in file VAX/argred.s, used by VAX version of sin/cos/tan
 
     (C). System supported functions :
 
     (C). System supported functions :
-        sqrt()      ... in files "IEEE/support.c", "VAX/sqrt.s"
-        drem()      ... in files "IEEE/support.c", "VAX/support.s"
-        finite()    ... in files "IEEE/support.c", "VAX/support.s"
-        logb()      ... in files "IEEE/support.c", "VAX/support.s"
-        scalb()     ... in files "IEEE/support.c", "VAX/support.s"
-        copysign()  ... in files "IEEE/support.c", "VAX/support.s"
-        rint()      ... in file  "floor.c"
+        sqrt()      ...in files IEEE/support.c, VAX/sqrt.s
+        drem()      ...in files IEEE/support.c, VAX/support.s
+        finite()    ...in files IEEE/support.c, VAX/support.s
+        logb()      ...in files IEEE/support.c, VAX/support.s
+        scalb()     ...in files IEEE/support.c, VAX/support.s
+        copysign()  ...in files IEEE/support.c, VAX/support.s
+        rint()      ...in file  floor.c
 
 
    Notes: 
 
 
    Notes: 
@@ -77,33 +90,33 @@ Revised on 5/10/85, 5/13/85, 6/14/85, 8/20/85, 8/27/85.
           standard 754 for double precision floating-point arithmetic.
 
       ii. On other than VAX or IEEE machines, run the original math 
           standard 754 for double precision floating-point arithmetic.
 
       ii. On other than VAX or IEEE machines, run the original math 
-          library, formerly "/usr/lib/libm.a", now "/usr/lib/libom.a",
-         if nothing better is available.
+          library, formerly "/usr/lib/libm.a", now "/usr/lib/libom.a", if
+         nothing better is available.
 
 
-     iii. The trigonometric functions sin(), cos(), tan() and atan2() in files
-         "VAX/sincos.s", "VAX/tan.s" and "VAX/atan2.s" are different from
-         those in "IEEE/trig.c" and "IEEE/atan2.c".  The VAX assembler code
-         uses the true value of pi to perform argument reduction, while the
-         C code uses the machine's value of PI rounded (see "IEEE/trig.c").
+     iii. The trigonometric functions sin/cos/tan/atan2 in files "VAX/sincos.s",
+          "VAX/tan.s" and "VAX/atan2.s" are different from those in
+          "IEEE/trig.c" and "IEEE/atan2.c".  The VAX assembler code uses the
+          true value of pi to perform argument reduction, while the C code uses
+          a machine value of PI (see "IEEE/trig.c").
 
 
 2. A computer system that conforms to IEEE standard 754 should provide 
 
 
 2. A computer system that conforms to IEEE standard 754 should provide 
-         sqrt(x),
-         drem(x,p), (double precision remainder function)
-         copysign(x,y),
-         finite(x),
-         scalb(x,N),
-         logb(x) and
-         rint(x).
+                sqrt(x),
+                drem(x,p), (double precision remainder function)
+                copysign(x,y),
+                finite(x),
+                scalb(x,N),
+                logb(x) and
+                rint(x).
    These functions are either required or recommended by the standard.
    For convenience, a (slow) C implementation of these functions is 
    provided in the file "IEEE/support.c".
 
    These functions are either required or recommended by the standard.
    For convenience, a (slow) C implementation of these functions is 
    provided in the file "IEEE/support.c".
 
-   Warning: The functions in "IEEE/support.c" are somewhat machine dependent.
+   Warning: The functions in IEEE/support.c are somewhat machine dependent.
    Some modifications may be necessary to run them on a different machine.
    Some modifications may be necessary to run them on a different machine.
-   Currently, if compiled with a suitable flag, "IEEE/support.c" will work on a
-   National 32000, a Zilog 8000, a VAX, and a SUN (cf. the "Makefile" in
-   this directory). Invoke the C compiler thus:
+   Currently, if compiled with a suitable flag, "IEEE/support.c" will work
+   on a National 32000, a Zilog 8000, a VAX, and a SUN (cf. the "Makefile"
+   in this directory). Invoke the C compiler thus:
 
         cc -c -DVAX IEEE/support.c              ... on a VAX, D-format
         cc -c -DNATIONAL IEEE/support.c         ... on a National 32000
 
         cc -c -DVAX IEEE/support.c              ... on a VAX, D-format
         cc -c -DNATIONAL IEEE/support.c         ... on a National 32000
@@ -125,8 +138,8 @@ Revised on 5/10/85, 5/13/85, 6/14/85, 8/20/85, 8/27/85.
    the VAX D-format (56-bit precision), and the IEEE double format 
    (53-bit precision).  The cbrt() in "IEEE/cbrt.c" is for IEEE machines 
    only. The functions in files that end with ".s" are for VAX computers 
    the VAX D-format (56-bit precision), and the IEEE double format 
    (53-bit precision).  The cbrt() in "IEEE/cbrt.c" is for IEEE machines 
    only. The functions in files that end with ".s" are for VAX computers 
-   only. The functions in files that end with ".c" (except "IEEE/cbrt.c") are
-   for VAX and IEEE machines. To use the VAX D-format, compile the code 
+   only. The functions in files that end with ".c" (except "IEEE/cbrt.c")
+   are for VAX and IEEE machines. To use the VAX D-format, compile the code 
    with -DVAX; to use IEEE double format on various IEEE machines, see 
    "Makefile" in this directory). 
 
    with -DVAX; to use IEEE double format on various IEEE machines, see 
    "Makefile" in this directory). 
 
@@ -140,7 +153,7 @@ Revised on 5/10/85, 5/13/85, 6/14/85, 8/20/85, 8/27/85.
                 enough to produce the hexadecimal values shown. If the
                 conversion is inaccurate, then one must know the exact machine 
                 representation of the constants and alter the assembly
                 enough to produce the hexadecimal values shown. If the
                 conversion is inaccurate, then one must know the exact machine 
                 representation of the constants and alter the assembly
-                language output from the compiler, or play tricks like 
+                language output from the compiler, or play tricks like
                 the following in a C program.
 
                         Example: to store the floating-point constant 
                 the following in a C program.
 
                         Example: to store the floating-point constant 
@@ -151,11 +164,11 @@ Revised on 5/10/85, 5/13/85, 6/14/85, 8/20/85, 8/27/85.
                         machine value and define p1 to be the double constant 
                         at the location of these two longwords:
 
                         machine value and define p1 to be the double constant 
                         at the location of these two longwords:
 
-                        static long  p1x[] = {0x3abe3d78, 0x066a67e1};
+                        static long  p1x[] = { 0x3abe3d78, 0x066a67e1};
                         #define      p1      (*(double*)p1x)
 
     Note:  On a VAX, some functions have two codes. For example, cabs() has
                         #define      p1      (*(double*)p1x)
 
     Note:  On a VAX, some functions have two codes. For example, cabs() has
-          one implementation in "IEEE/cabs.c", and another in "VAX/cabs.s".
+          one implementation in "IEEE/cabs.c", and another in "VAX/cabs.s". 
            In this case, the assembly language version is preferred.
 
 
            In this case, the assembly language version is preferred.
 
 
@@ -169,15 +182,15 @@ Revised on 5/10/85, 5/13/85, 6/14/85, 8/20/85, 8/27/85.
             on all tested machines (VAX, SUN, NATIONAL, ZILOG), provided that 
             x to the power of y is representable exactly.
 
             on all tested machines (VAX, SUN, NATIONAL, ZILOG), provided that 
             x to the power of y is representable exactly.
 
-            cosh(), sinh(), acosh(), asinh(), tanh(), atanh() and log10() have
-           errors below about 3 ULPs. 
+            cosh, sinh, acosh, asinh, tanh, atanh and log10 have errors below
+            about 3 ULPs. 
 
 
-            For trigonometric and inverse trigonometric functions, let
-           [trig(x)] denote the value actually computed for trig(x).
+            For trigonometric and inverse trigonometric functions: 
+
+                Let [trig(x)] denote the value actually computed for trig(x),
 
                 1) Those codes using the machine's value PI (true pi rounded):
 
                 1) Those codes using the machine's value PI (true pi rounded):
-                   (in files "IEEE/trig.c", "IEEE/atan2.c", "asincos.c" and
-                  "atan.c".)
+                   (source codes: IEEE/{trig.c,atan2.c}, asincos.c and atan.c)
 
                    The errors in [sin(x)], [cos(x)], and [atan(x)] are below 
                    1 ULP compared with sin(x*pi/PI), cos(x*pi/PI), and 
 
                    The errors in [sin(x)], [cos(x)], and [atan(x)] are below 
                    1 ULP compared with sin(x*pi/PI), cos(x*pi/PI), and 
@@ -187,14 +200,16 @@ Revised on 5/10/85, 5/13/85, 6/14/85, 8/20/85, 8/27/85.
                    return acos(x)*PI/pi, asin(x)*PI/pi, and atan2(y,x)*PI/pi
                    respectively to similar accuracy.
 
                    return acos(x)*PI/pi, asin(x)*PI/pi, and atan2(y,x)*PI/pi
                    respectively to similar accuracy.
 
+
                 2) Those using true pi (for VAX D-format only):
                 2) Those using true pi (for VAX D-format only):
-                   (in files "VAX/sincos.s", "VAX/tan.s", "VAX/atan2.s",
-                  "asincos.c" and "atan.c".)
+                   (source codes: VAX/{sincos.s,tan.s,atan2.s}, asincos.c and
+                   atan.c)
 
                    The errors in [sin(x)], [cos(x)], and [atan(x)] are below
 
                    The errors in [sin(x)], [cos(x)], and [atan(x)] are below
-                   1 ULP.  [tan(x)], [atan2(y,x)], [acos(x)], and [asin(x)] 
+                   1 ULP. [tan(x)], [atan2(y,x)], [acos(x)], and [asin(x)] 
                    have errors below about 2 ULPs. 
 
                    have errors below about 2 ULPs. 
 
+
             Here are the results of some test runs to find worst errors on 
             the VAX :
 
             Here are the results of some test runs to find worst errors on 
             the VAX :
 
@@ -243,25 +258,7 @@ Revised on 5/10/85, 5/13/85, 6/14/85, 8/20/85, 8/27/85.
        be called via short subroutine calls.
 
 
        be called via short subroutine calls.
 
 
-6. j0,j1,jn.
+6. j0, j1, jn.
 
         The modifications to these routines were only in how an invalid
 
         The modifications to these routines were only in how an invalid
-        floating point operation is signaled on a VAX.
-
-
-7. Copyright notice, and Disclaimer:
-
-***************************************************************************
-*                                                                         * 
-* Copyright (c) 1985 Regents of the University of California.             *
-*                                                                         * 
-* Use and reproduction of this software are granted  in  accordance  with *
-* the terms and conditions specified in  the  Berkeley  Software  License *
-* Agreement (in particular, this entails acknowledgement of the programs' *
-* source, and inclusion of this notice) with the additional understanding *
-* that  all  recipients  should regard themselves as participants  in  an *
-* ongoing  research  project and hence should  feel  obligated  to report *
-* their  experiences (good or bad) with these elementary function  codes, *
-* using "sendbug 4bsd-bugs@BERKELEY", to the authors.                     *
-*                                                                         *
-***************************************************************************
+        floating point operations is signaled.
index 5d637e5..8b0233b 100644 (file)
@@ -12,7 +12,8 @@
  */
 
 #ifndef lint
  */
 
 #ifndef lint
-static char sccsid[] = "@(#)atan2.c    1.1 (ELEFUNT) %G%";
+static char sccsid[] =
+"@(#)atan2.c   1.3 (Berkeley) 8/21/85; 1.2 (ucb.elefunt) %G%";
 #endif not lint
 
 /* ATAN2(Y,X)
 #endif not lint
 
 /* ATAN2(Y,X)
index bd1676b..f684a3f 100644 (file)
@@ -12,7 +12,8 @@
  */
 
 #ifndef lint
  */
 
 #ifndef lint
-static char sccsid[] = "@(#)trig.c     1.1 (ELEFUNT) %G%";
+static char sccsid[] =
+"@(#)trig.c    1.2 (Berkeley) 8/22/85; 1.2 (ucb.elefunt) %G%";
 #endif not lint
 
 /* SIN(X), COS(X), TAN(X)
 #endif not lint
 
 /* SIN(X), COS(X), TAN(X)
index 3cf94b6..3af6327 100644 (file)
@@ -1,5 +1,5 @@
-# @(#)Makefile 1.1 (ELEFUNT) %G%
-SCCSID = "@(#)Makefile 1.1 (ELEFUNT) %G%"
+# @(#)Makefile 1.2 (ucb.elefunt) %G%
+SCCSID = "@(#)Makefile 1.2 (ucb.elefunt) %G%"
 MORE=atan2.o cbrt.o trig.o cabs.o support.o
 
 .c.o:
 MORE=atan2.o cbrt.o trig.o cabs.o support.o
 
 .c.o:
index 79fcabe..9eea151 100644 (file)
@@ -12,7 +12,8 @@
  */
 
 #ifndef lint
  */
 
 #ifndef lint
-static char sccsid[] = "@(#)cabs.c     1.1 (ELEFUNT) %G%";
+static char sccsid[] =
+"@(#)cabs.c    1.2 (Berkeley) 8/21/85; 1.2 (ucb.elefunt) %G%";
 #endif not lint
 
 /* CABS(Z)
 #endif not lint
 
 /* CABS(Z)
index fb2b248..d5908a7 100644 (file)
@@ -12,7 +12,8 @@
  */
 
 #ifndef lint
  */
 
 #ifndef lint
-static char sccsid[] = "@(#)cbrt.c     1.1 (ELEFUNT) %G%";
+static char sccsid[] =
+"@(#)cbrt.c    1.1 (Berkeley) 5/23/85; 1.2 (ucb.elefunt) %G%";
 #endif not lint
 
 /* kahan's cube root (53 bits IEEE double precision)
 #endif not lint
 
 /* kahan's cube root (53 bits IEEE double precision)
index 350c9a4..6f8dd62 100644 (file)
@@ -12,7 +12,8 @@
  */
 
 #ifndef lint
  */
 
 #ifndef lint
-static char sccsid[] = "@(#)support.c  1.1 (ELEFUNT) %G%";
+static char sccsid[] =
+"@(#)support.c 1.1 (Berkeley) 5/23/85; 1.2 (ucb.elefunt) %G%";
 #endif not lint
 
 /* 
 #endif not lint
 
 /* 
index 1c4db2d..f17ed6a 100644 (file)
@@ -1,5 +1,5 @@
-# @(#)Makefile 1.1 (ELEFUNT) %G%
-SCCSID = "@(#)Makefile 1.1 (ELEFUNT) %G%"
+# @(#)Makefile 1.2 (ucb.elefunt) %G%
+# SCCSID = "@(#)Makefile       1.2 (ucb.elefunt) %G%"
 MORE=atan2.o cabs.o cbrt.o sqrt.o sincos.o tan.o argred.o support.o infnan.o
 
 .s.o:
 MORE=atan2.o cabs.o cbrt.o sqrt.o sincos.o tan.o argred.o support.o infnan.o
 
 .s.o:
index 1e07fdf..d771a76 100644 (file)
@@ -11,7 +11,7 @@
 # using "sendbug 4bsd-bugs@BERKELEY", to the authors.
 #
 
 # using "sendbug 4bsd-bugs@BERKELEY", to the authors.
 #
 
-# @(#)argred.s 1.1 (ELEFUNT) %G%
+# @(#)argred.s 1.1 (Berkeley) 8/21/85; 1.2 (ucb.elefunt) %G%
 
 #  libm$argred implements Bob Corbett's argument reduction and
 #  libm$sincos implements Peter Tang's double precision sin/cos.
 
 #  libm$argred implements Bob Corbett's argument reduction and
 #  libm$sincos implements Peter Tang's double precision sin/cos.
index 3ff0285..430adcd 100644 (file)
@@ -11,7 +11,7 @@
 # using "sendbug 4bsd-bugs@BERKELEY", to the authors.
 #
 
 # using "sendbug 4bsd-bugs@BERKELEY", to the authors.
 #
 
-# @(#)atan2.s  1.1 (ELEFUNT) %G%
+# @(#)atan2.s  1.2 (Berkeley) 8/21/85; 1.2 (ucb.elefunt) %G%
 
 # ATAN2(Y,X)
 # RETURN ARG (X+iY)
 
 # ATAN2(Y,X)
 # RETURN ARG (X+iY)
index e7935c0..3088fc6 100644 (file)
@@ -11,7 +11,7 @@
 # using "sendbug 4bsd-bugs@BERKELEY", to the authors.
 #
 
 # using "sendbug 4bsd-bugs@BERKELEY", to the authors.
 #
 
-# @(#)cabs.s   1.1 (ELEFUNT) %G%
+# @(#)cabs.s   1.2 (Berkeley) 8/21/85; 1.2 (ucb.elefunt) %G%
 
 # double precision complex absolute value
 # CABS by W. Kahan, 9/7/80.
 
 # double precision complex absolute value
 # CABS by W. Kahan, 9/7/80.
index e2921a8..953ebc1 100644 (file)
@@ -11,7 +11,7 @@
 # using "sendbug 4bsd-bugs@BERKELEY", to the authors.
 #
 
 # using "sendbug 4bsd-bugs@BERKELEY", to the authors.
 #
 
-# @(#)cbrt.s   1.1 (ELEFUNT) %G%
+# @(#)cbrt.s   1.1 (Berkeley) 5/23/85; 1.2 (ucb.elefunt) %G%
 
 # double cbrt(double arg)
 # W. Kahan, 10/13/80. revised 1/13/84 for keeping sign symmetry
 
 # double cbrt(double arg)
 # W. Kahan, 10/13/80. revised 1/13/84 for keeping sign symmetry
index f28c16c..4b6f02c 100644 (file)
@@ -10,7 +10,7 @@
  * their  experiences (good or bad) with these elementary function  codes,
  * using "sendbug 4bsd-bugs@BERKELEY", to the authors.
  *
  * their  experiences (good or bad) with these elementary function  codes,
  * using "sendbug 4bsd-bugs@BERKELEY", to the authors.
  *
- * @(#)infnan.s        1.1 (ELEFUNT) %G%
+ * @(#)infnan.s        1.1 (Berkeley) 8/21/85; 1.2 (ucb.elefunt) %G%
  *
  * infnan(arg) int arg;
  * where arg :=    EDOM        if result is  NaN
  *
  * infnan(arg) int arg;
  * where arg :=    EDOM        if result is  NaN
index a2c1721..d8e4f57 100644 (file)
@@ -11,7 +11,7 @@
 # using "sendbug 4bsd-bugs@BERKELEY", to the authors.
 #
 
 # using "sendbug 4bsd-bugs@BERKELEY", to the authors.
 #
 
-# @(#)sincos.s 1.1 (ELEFUNT) %G%
+# @(#)sincos.s 1.1 (Berkeley) 8/21/85; 1.2 (ucb.elefunt) %G%
 
 #  This is the implementation of Peter Tang's double precision  
 #  sine and cosine for the VAX using Bob Corbett's argument reduction.
 
 #  This is the implementation of Peter Tang's double precision  
 #  sine and cosine for the VAX using Bob Corbett's argument reduction.
index 8befae7..764c7f4 100644 (file)
@@ -11,7 +11,7 @@
  * using "sendbug 4bsd-bugs@BERKELEY", to the authors.
  *
  *
  * using "sendbug 4bsd-bugs@BERKELEY", to the authors.
  *
  *
- * @(#)sqrt.s  1.1 (ELEFUNT) %G%
+ * @(#)sqrt.s  1.1 (Berkeley) 8/21/85; 1.2 (ucb.elefunt) %G%
  *
  * double sqrt(arg)   revised August 15,1982
  * double arg;
  *
  * double sqrt(arg)   revised August 15,1982
  * double arg;
index 7190229..b4b444c 100644 (file)
@@ -10,7 +10,7 @@
  * their  experiences (good or bad) with these elementary function  codes,
  * using "sendbug 4bsd-bugs@BERKELEY", to the authors.
  *
  * their  experiences (good or bad) with these elementary function  codes,
  * using "sendbug 4bsd-bugs@BERKELEY", to the authors.
  *
- * @(#)support.s       1.1 (ELEFUNT) %G%
+ * @(#)support.s       1.3 (Berkeley) 8/21/85; 1.2 (ucb.elefunt) %G%
  *
  * copysign(x,y),
  * logb(x),
  *
  * copysign(x,y),
  * logb(x),
index c31c0e5..028e13e 100644 (file)
@@ -11,7 +11,7 @@
 # using "sendbug 4bsd-bugs@BERKELEY", to the authors.
 #
 
 # using "sendbug 4bsd-bugs@BERKELEY", to the authors.
 #
 
-# @(#)tan.s    1.1 (ELEFUNT) %G%
+# @(#)tan.s    1.1 (Berkeley) 8/21/85; 1.2 (ucb.elefunt) %G%
 
 #  This is the implementation of Peter Tang's double precision  
 #  tangent for the VAX using Bob Corbett's argument reduction.
 
 #  This is the implementation of Peter Tang's double precision  
 #  tangent for the VAX using Bob Corbett's argument reduction.