Ticket #7007: math_changes.2.patch

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
@@ -56 +57 @@
 #define FP_ZERO         2
 #define FP_SUBNORMAL    3
 #define FP_NORMAL       4
+#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
 
+# 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
+
 #ifdef __cplusplus
 struct __exception;
 extern "C" int matherr(struct __exception *);
@@ -350 +370 @@
         : sizeof(value) == sizeof(double) ? __isinf(value)      \
         : __isinfl(value))
 
+/* Return nonzero value if X is greater than Y.  */
+#  define isgreater(x, y) \
+  (__extension__                                  \
+   ({ __typeof__(x) __x = (x); __typeof__(y) __y = (y);               \
+      !isunordered (__x, __y) && __x > __y; }))
+
+/* Return nonzero value if X is greater than or equal to Y.  */
+#  define isgreaterequal(x, y) \
+  (__extension__                                  \
+   ({ __typeof__(x) __x = (x); __typeof__(y) __y = (y);               \
+      !isunordered (__x, __y) && __x >= __y; }))
+
+/* Return nonzero value if X is less than Y.  */
+#  define isless(x, y) \
+  (__extension__                                  \
+   ({ __typeof__(x) __x = (x); __typeof__(y) __y = (y);               \
+      !isunordered (__x, __y) && __x < __y; }))
+
+/* Return nonzero value if X is less than or equal to Y.  */
+#  define islessequal(x, y) \
+  (__extension__                                  \
+   ({ __typeof__(x) __x = (x); __typeof__(y) __y = (y);               \
+      !isunordered (__x, __y) && __x <= __y; }))
+
+/* Return nonzero value if either X is less than Y or Y is less than X.  */
+#  define islessgreater(x, y) \
+  (__extension__                                  \
+   ({ __typeof__(x) __x = (x); __typeof__(y) __y = (y);               \
+      !isunordered (__x, __y) && (__x < __y || __y < __x); }))
+
+/* Return nonzero value if arguments are unordered.  */
+#  define isunordered(u, v) \
+  (__extension__                                  \
+   ({ __typeof__(u) __u = (u); __typeof__(v) __v = (v);               \
+      fpclassify (__u) == FP_NAN || fpclassify (__v) == FP_NAN; }))
+
 #ifdef __cplusplus
 }
 #endif