Ticket #7007: math_changes.patch

src/system/libroot/posix/glibc/math/math.h

@@ -129 +129 @@
 /* ISO C99 defines some generic macros which work on any data type.  */
 #if __USE_ISOC99
 
-/* Get the architecture specific values describing the floating-point
-   evaluation.  The following symbols will get defined:
+# if FLT_EVAL_METHOD == 1
+typedef double float_t;
+typedef double double_t;
+# elif FLT_EVAL_METHOD == 2
+typedef long double float_t;
+typedef long double double_t;
+# else
+typedef float float_t;
+typedef double double_t;
+# endif          
 
-    float_t floating-point type at least as wide as `float' used
-        to evaluate `float' expressions
-    double_t    floating-point type at least as wide as `double' used
-        to evaluate `double' expressions
+#define FP_FAST_FMA         1
+#define FP_FAST_FMAL        1
+#define FP_FAST_FMAF        1
+#define FP_ILOGB0           2147483648
+#define FP_ILOGBNAN         2147483647
+#define MATH_ERRNO          1
+#define MATH_ERREXCEPT      2
+#define math_errhandling    3
 
-    FLT_EVAL_METHOD
-        Defined to
-          0 if `float_t' is `float' and `double_t' is `double'
-          1 if `float_t' and `double_t' are `double'
-          2 if `float_t' and `double_t' are `long double'
-          else  `float_t' and `double_t' are unspecified
-
-    INFINITY    representation of the infinity value of type `float'
-
-    FP_FAST_FMA
-    FP_FAST_FMAF
-    FP_FAST_FMAL
-        If defined it indicates that the `fma' function
-        generally executes about as fast as a multiply and an add.
-        This macro is defined only iff the `fma' function is
-        implemented directly with a hardware multiply-add instructions.
-
-    FP_ILOGB0   Expands to a value returned by `ilogb (0.0)'.
-    FP_ILOGBNAN Expands to a value returned by `ilogb (NAN)'.
-
-    DECIMAL_DIG Number of decimal digits supported by conversion between
-        decimal and all internal floating-point formats.
-
-*/
-
 /* All floating-point numbers can be put in one of these categories.  */
 enum
   {

src/system/libroot/posix/glibc/arch/generic/w_gammaf_r.c

@@ +1 @@
+/* w_gammaf_r.c -- float version of w_gamma_r.c.
+ * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com.
+ */
+
+/*
+ * ====================================================
+ * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
+ *
+ * Developed at SunPro, a Sun Microsystems, Inc. business.
+ * Permission to use, copy, modify, and distribute this
+ * software is freely granted, provided that this notice
+ * is preserved.
+ * ====================================================
+ */
+
+#if defined(LIBM_SCCS) && !defined(lint)
+static char rcsid[] = "$NetBSD: w_gammaf_r.c,v 1.8 2002/05/26 22:02:01 wiz Exp $";
+#endif
+
+/*
+ * wrapper float gammaf_r(float x, int *signgamp)
+ */
+
+#include "math.h"
+#include "math_private.h"
+
+float gammaf_r(float x, int *signgamp) /* wrapper lgammaf_r */
+{
+#ifdef _IEEE_LIBM
+    return __ieee754_lgammaf_r(x,signgamp);
+#else
+        float y;
+        y = __ieee754_lgammaf_r(x,signgamp);
+        if(_LIB_VERSION == _IEEE_) return y;
+        if(!finitef(y)&&finitef(x)) {
+            if(floorf(x)==x&&x<=(float)0.0)
+            /* gammaf pole */
+                return (float)__kernel_standard((double)x,(double)x,141);
+            else
+            /* gamma overflow */
+                return (float)__kernel_standard((double)x,(double)x,140);
+        } else
+            return y;
+#endif
+}
+weak_alias (__gammaf_r, gammaf_r)

src/system/libroot/posix/glibc/arch/generic/w_gamma_r.c

@@ +1 @@
+/* @(#)w_gamma_r.c 1.3 95/01/18 */
+/*
+ * ====================================================
+ * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
+ *
+ * Developed at SunSoft, a Sun Microsystems, Inc. business.
+ * Permission to use, copy, modify, and distribute this
+ * software is freely granted, provided that this notice 
+ * is preserved.
+ * ====================================================
+ */
+
+#include "math.h"
+#include "math_private.h"
+
+
+#ifdef __STDC__
+    double gamma_r(double x, int *signgamp) /* wrapper lgamma_r */
+#else
+    double gamma_r(x,signgamp)              /* wrapper lgamma_r */
+        double x; int *signgamp;
+#endif
+{
+#ifdef _IEEE_LIBM
+    return __ieee754_gamma_r(x,signgamp);
+#else
+        double y;
+        y = __ieee754_gamma_r(x,signgamp);
+        if(_LIB_VERSION == _IEEE_) return y;
+        if(!finite(y)&&finite(x)) {
+            if(floor(x)==x&&x<=0.0)
+                return __kernel_standard(x,x,41); /* gamma pole */
+            else
+                return __kernel_standard(x,x,40); /* gamma overflow */
+        } else
+            return y;
+#endif
+}
+weak_alias (__gamma_r, gamma_r)
+#ifdef NO_LONG_DOUBLE
+strong_alias (__gamma_r, __gammal_r)
+weak_alias (__gamma_r, gammal_r)
+#endif

headers/posix/math.h

@@ -19 +19 @@
 #define M_2_SQRTPI      1.12837916709551257390  /* 2/sqrt(pi) */
 #define M_SQRT2         1.41421356237309504880  /* sqrt(2) */
 #define M_SQRT1_2       0.70710678118654752440  /* 1/sqrt(2) */
+#define MAXFLOAT            3.40282346638528860e+38 /*max float value */
 
 /* platform independent IEEE floating point special values */
 #if __GNUC__ >= 4