BSD 4_4 release

[unix-history] / usr / src / lib / libm / common_source / gamma.c

diff --git a/usr/src/lib/libm/common_source/gamma.c b/usr/src/lib/libm/common_source/gamma.c
index cc52aa8..5d270f0 100644 (file)
--- a/usr/src/lib/libm/common_source/gamma.c
+++ b/usr/src/lib/libm/common_source/gamma.c
@@ -1,12 +1,38 @@
 /*-
 /*-

- * Copyright (c) 1992 The Regents of the University of California.
- * All rights reserved.
+ * Copyright (c) 1992, 1993
+ *     The Regents of the University of California.  All rights reserved.

  *
  *

- * %sccs.include.redist.c%
+ * Redistribution and use in source and binary forms, with or without
+ * modification, are permitted provided that the following conditions
+ * are met:
+ * 1. Redistributions of source code must retain the above copyright
+ *    notice, this list of conditions and the following disclaimer.
+ * 2. Redistributions in binary form must reproduce the above copyright
+ *    notice, this list of conditions and the following disclaimer in the
+ *    documentation and/or other materials provided with the distribution.
+ * 3. All advertising materials mentioning features or use of this software
+ *    must display the following acknowledgement:
+ *     This product includes software developed by the University of
+ *     California, Berkeley and its contributors.
+ * 4. Neither the name of the University nor the names of its contributors
+ *    may be used to endorse or promote products derived from this software
+ *    without specific prior written permission.
+ *
+ * THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND
+ * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+ * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+ * ARE DISCLAIMED.  IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE
+ * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
+ * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
+ * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
+ * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
+ * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
+ * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
+ * SUCH DAMAGE.

  */
 
 #ifndef lint
  */
 
 #ifndef lint

-static char sccsid[] = "@(#)gamma.c    5.3 (Berkeley) %G%";
+static char sccsid[] = "@(#)gamma.c    8.1 (Berkeley) 6/4/93";

 #endif /* not lint */
 
 /*
 #endif /* not lint */
 
 /*


@@ -120,7 +146,7 @@ gamma(x)
                if(x > 171.63)
                        return(one/zero);
                u = large_gam(x);
                if(x > 171.63)
                        return(one/zero);
                u = large_gam(x);

-               return(exp__D(u.a, u.b));
+               return(__exp__D(u.a, u.b));

        } else if (x >= 1.0 + LEFT + x0)
                return (small_gam(x));
        else if (x > 1.e-17)
        } else if (x >= 1.0 + LEFT + x0)
                return (small_gam(x));
        else if (x > 1.e-17)


@@ -154,7 +180,7 @@ large_gam(x)
        p = Pa0+z*(Pa1+z*(Pa2+z*(Pa3+z*(Pa4+z*(Pa5+z*(Pa6+z*Pa7))))));
        p = p/x;
 
        p = Pa0+z*(Pa1+z*(Pa2+z*(Pa3+z*(Pa4+z*(Pa5+z*(Pa6+z*Pa7))))));
        p = p/x;
 

-       u = log__D(x);
+       u = __log__D(x);

        u.a -= one;
        v.a = (x -= .5);
        TRUNC(v.a);
        u.a -= one;
        v.a = (x -= .5);
        TRUNC(v.a);


@@ -291,12 +317,12 @@ neg_gam(x)
                        return ((double)sgn*tiny*tiny);
                y = one - x;            /* exact: 128 < |x| < 255 */
                lg = large_gam(y);
                        return ((double)sgn*tiny*tiny);
                y = one - x;            /* exact: 128 < |x| < 255 */
                lg = large_gam(y);

-               lsine = log__D(M_PI/z); /* = TRUNC(log(u)) + small */
-               lg.a -= lsine.a;        /* exact (opposite signs) */
+               lsine = __log__D(M_PI/z);       /* = TRUNC(log(u)) + small */
+               lg.a -= lsine.a;                /* exact (opposite signs) */

                lg.b -= lsine.b;
                y = -(lg.a + lg.b);
                z = (y + lg.a) + lg.b;
                lg.b -= lsine.b;
                y = -(lg.a + lg.b);
                z = (y + lg.a) + lg.b;

-               y = exp__D(y, z);
+               y = __exp__D(y, z);

                if (sgn < 0) y = -y;
                return (y);
        }
                if (sgn < 0) y = -y;
                return (y);
        }