WebSVN – Kolibri OS – Diff – /contrib/sdk/sources/newlib/libc/include/math.h

 #ifndef  _MATH_H_
 #define  _MATH_H_
 #include
 #include
 #include
 extern double sqrt _PARAMS((double));
 extern double fmod _PARAMS((double, double));
 #endif /* ! defined (__math_68881) */
 #endif /* ! defined (_REENT_ONLY) */
-#if !defined(__STRICT_ANSI__) || defined(__cplusplus) || \
-  (defined(__STDC_VERSION__) && __STDC_VERSION__ >= 199901L)
 #if __ISO_C_VISIBLE >= 1999
 /* ISO C99 types and macros. */
 /* FIXME:  FLT_EVAL_METHOD should somehow be gotten from float.h (which is hard,
 extern int __fpclassifyf (float x);
 extern int __fpclassifyd (double x);
 extern int __signbitf (float x);
 extern int __signbitd (double x);
+/* Note: isinf and isnan were once functions in newlib that took double
+ *       arguments.  C99 specifies that these names are reserved for macros
+ *       supporting multiple floating point types.  Thus, they are
+ *       now defined as macros.  Implementations of the old functions
+ *       taking double arguments still exist for compatibility purposes
+ *       (prototypes for them are in ).  */
+#if __GNUC_PREREQ (4, 4)
+  #define fpclassify(__x) (__builtin_fpclassify (FP_NAN, FP_INFINITE, \
+						 FP_NORMAL, FP_SUBNORMAL, \
+						 FP_ZERO, __x))
+  #ifndef isfinite
+    #define isfinite(__x)	(__builtin_isfinite (__x))
+  #endif
+  #ifndef isinf
+    #define isinf(__x) (__builtin_isinf_sign (__x))
+  #endif
+  #ifndef isnan
+    #define isnan(__x) (__builtin_isnan (__x))
+  #endif
+  #define isnormal(__x) (__builtin_isnormal (__x))
 #else
-#define fpclassify(__x) \
+  #define fpclassify(__x) \
-	((sizeof(__x) == sizeof(float))  ? __fpclassifyf(__x) : \
+	  ((sizeof(__x) == sizeof(float))  ? __fpclassifyf(__x) : \
-	__fpclassifyd(__x))
 	  __fpclassifyd(__x))
-#ifndef isfinite
+  #ifndef isfinite
-  #define isfinite(__y) \
+    #define isfinite(__y) \
-          (__extension__ ({int __cy = fpclassify(__y); \
+	    (__extension__ ({int __cy = fpclassify(__y); \
-                           __cy != FP_INFINITE && __cy != FP_NAN;}))
+			     __cy != FP_INFINITE && __cy != FP_NAN;}))
-#endif
-/* Note: isinf and isnan were once functions in newlib that took double
- *       arguments.  C99 specifies that these names are reserved for macros
- *       supporting multiple floating point types.  Thus, they are
- *       now defined as macros.  Implementations of the old functions
- *       taking double arguments still exist for compatibility purposes
- *       (prototypes for them are in ).  */
+  #endif
-#ifndef isinf
+  #ifndef isinf
-  #define isinf(y) (fpclassify(y) == FP_INFINITE)
+    #define isinf(__x) (fpclassify(__x) == FP_INFINITE)
-#endif
   #endif
-#ifndef isnan
+  #ifndef isnan
+    #define isnan(__x) (fpclassify(__x) == FP_NAN)
+  #endif
-  #define isnan(y) (fpclassify(y) == FP_NAN)
+  #define isnormal(__x) (fpclassify(__x) == FP_NORMAL)
+#endif
+#if __GNUC_PREREQ (4, 0)
+  #if defined(_HAVE_LONG_DOUBLE)
+    #define signbit(__x) \
+	    ((sizeof(__x) == sizeof(float))  ? __builtin_signbitf(__x) : \
+	     (sizeof(__x) == sizeof(double)) ? __builtin_signbit (__x) : \
+					       __builtin_signbitl(__x))
+  #else
-#endif
+    #define signbit(__x) \
+	    ((sizeof(__x) == sizeof(float))  ? __builtin_signbitf(__x) : \
+					       __builtin_signbit (__x))
   #endif
-#define isnormal(y) (fpclassify(y) == FP_NORMAL)
+#else
-#define signbit(__x) \
+  #define signbit(__x) \
+	((sizeof(__x) == sizeof(float))  ?  __signbitf(__x) : \
+		__signbitd(__x))
+#endif
+#if __GNUC_PREREQ (2, 97)
+#define isgreater(__x,__y)	(__builtin_isgreater (__x, __y))
+#define isgreaterequal(__x,__y)	(__builtin_isgreaterequal (__x, __y))
+#define isless(__x,__y)		(__builtin_isless (__x, __y))
+#define islessequal(__x,__y)	(__builtin_islessequal (__x, __y))
-	((sizeof(__x) == sizeof(float))  ?  __signbitf(__x) : \
+#define islessgreater(__x,__y)	(__builtin_islessgreater (__x, __y))
-		__signbitd(__x))
+#define isunordered(__x,__y)	(__builtin_isunordered (__x, __y))
 #else
 #define isgreater(x,y) \
           (__extension__ ({__typeof__(x) __x = (x); __typeof__(y) __y = (y); \
                            !isunordered(__x,__y) && (__x < __y || __x > __y);}))
 #define isunordered(a,b) \
           (__extension__ ({__typeof__(a) __a = (a); __typeof__(b) __b = (b); \
+                           fpclassify(__a) == FP_NAN || fpclassify(__b) == FP_NAN;}))
-                           fpclassify(__a) == FP_NAN || fpclassify(__b) == FP_NAN;}))
+#endif
 extern float erfcf _PARAMS((float));
 extern float log2f _PARAMS((float));
 extern float hypotf _PARAMS((float, float));
 #endif /* ! defined (_REENT_ONLY) */
+/* Newlib doesn't fully support long double math functions so far.
+   On platforms where long double equals double the long double functions
+   simply call the double functions.  On Cygwin the long double functions
+   are implemented independently from newlib to be able to use optimized
-/* On platforms where long double equals double.  */
+   assembler functions despite using the Microsoft x86_64 ABI. */
-#ifdef _LDBL_EQ_DBL
+#if defined (_LDBL_EQ_DBL) || defined (__CYGWIN__)
 /* Reentrant ANSI C functions.  */
 #ifndef __math_68881
 extern long double atanl _PARAMS((long double));
 extern long double cosl _PARAMS((long double));
 extern long double remainderl _PARAMS((long double, long double));
 extern long double lgammal _PARAMS((long double));
 extern long double erfl _PARAMS((long double));
 extern long double erfcl _PARAMS((long double));
 #endif /* ! defined (_REENT_ONLY) */
-#else /* !_LDBL_EQ_DBL */
+#else /* !_LDBL_EQ_DBL && !__CYGWIN__ */
 extern long double hypotl _PARAMS((long double, long double));
 extern long double sqrtl _PARAMS((long double));
 #ifdef __i386__
 /* Other long double precision functions.  */
 extern _LONG_DOUBLE rintl _PARAMS((_LONG_DOUBLE));
 extern long int lrintl _PARAMS((_LONG_DOUBLE));
 extern long long int llrintl _PARAMS((_LONG_DOUBLE));
 #endif /* __i386__ */
-#endif /* !_LDBL_EQ_DBL */
+#endif /* !_LDBL_EQ_DBL && !__CYGWIN__ */
-#endif /* !defined (__STRICT_ANSI__) || defined(__cplusplus) || (defined(__STDC_VERSION__) && __STDC_VERSION__ >= 199901L) */
+#endif /* __ISO_C_VISIBLE >= 1999 */
 #if !defined (__STRICT_ANSI__) || defined(__cplusplus)
+#if __MISC_VISIBLE
+extern double drem _PARAMS((double, double));
+extern float dremf _PARAMS((float, float));
+#ifdef __CYGWIN__
-extern double drem _PARAMS((double, double));
+extern float dreml _PARAMS((long double, long double));
-extern void sincos _PARAMS((double, double *, double *));
+#endif /* __CYGWIN__ */
+extern double gamma_r _PARAMS((double, int *));
+extern double lgamma_r _PARAMS((double, int *));
+extern float gammaf_r _PARAMS((float, int *));
+extern float lgammaf_r _PARAMS((float, int *));
-extern double gamma_r _PARAMS((double, int *));
+#endif
 extern double lgamma_r _PARAMS((double, int *));
 #if __MISC_VISIBLE || __XSI_VISIBLE
 extern double y0 _PARAMS((double));
 extern double y1 _PARAMS((double));
 extern double yn _PARAMS((int, double));
+extern double j0 _PARAMS((double));
-extern double j0 _PARAMS((double));
-extern double j1 _PARAMS((double));
-extern double jn _PARAMS((int, double));
+extern double j1 _PARAMS((double));
-extern float dremf _PARAMS((float, float));
-extern void sincosf _PARAMS((float, float *, float *));
+extern double jn _PARAMS((int, double));
-extern float gammaf_r _PARAMS((float, int *));
+#endif
 extern float lgammaf_r _PARAMS((float, int *));
 #if __MISC_VISIBLE || __XSI_VISIBLE >= 600
 extern float y0f _PARAMS((float));
 extern float y1f _PARAMS((float));
+extern float ynf _PARAMS((int, float));
-extern float ynf _PARAMS((int, float));
+extern float j0f _PARAMS((float));
+extern float j1f _PARAMS((float));
+extern float jnf _PARAMS((int, float));
+#endif
+/* GNU extensions */
+#if __GNU_VISIBLE
-extern float j0f _PARAMS((float));
+extern void sincos _PARAMS((double, double *, double *));
-extern float j1f _PARAMS((float));
+extern void sincosf _PARAMS((float, float *, float *));
-extern float jnf _PARAMS((int, float));
+#ifdef __CYGWIN__
 extern void sincosl _PARAMS((long double, long double *, long double *));
-/* GNU extensions */
+#endif /* __CYGWIN__ */
 extern float exp10f _PARAMS((float));
 # endif
 # ifndef pow10f
 extern float pow10f _PARAMS((float));
 # endif
+#ifdef __CYGWIN__
+# ifndef exp10l
+extern float exp10l _PARAMS((float));
+# endif
+# ifndef pow10l
+extern float pow10l _PARAMS((float));
+# endif
+#endif /* __CYGWIN__ */
+#endif /* __GNU_VISIBLE */
-#endif /* !defined (__STRICT_ANSI__) || defined(__cplusplus) */
-#ifndef __STRICT_ANSI__
 #if __MISC_VISIBLE || __XSI_VISIBLE
 /* The gamma functions use a global variable, signgam.  */
 #ifndef _REENT_ONLY
 #define signgam (*__signgam())
 extern int *__signgam _PARAMS((void));
 #endif /* ! defined (_REENT_ONLY) */
+#define __signgam_r(ptr) _REENT_SIGNGAM(ptr)
 #endif /* __MISC_VISIBLE || __XSI_VISIBLE */
 #define __signgam_r(ptr) _REENT_SIGNGAM(ptr)
 #if __SVID_VISIBLE
 /* The exception structure passed to the matherr routine.  */
 /* We have a problem when using C++ since `exception' is a reserved
 #define OVERFLOW 3
 #define UNDERFLOW 4
 #define TLOSS 5
 #define PLOSS 6
-#endif /* ! defined (__STRICT_ANSI__) */
+#endif /* __SVID_VISIBLE */
 /* Useful constants.  */
-#if !defined(__STRICT_ANSI__) || ((_XOPEN_SOURCE - 0) >= 500)
+#if __BSD_VISIBLE || __XSI_VISIBLE >= 500
 #define M_SQRT2		1.41421356237309504880
 #define M_SQRT1_2	0.70710678118654752440
 #endif
-#ifndef __STRICT_ANSI__
+#if __BSD_VISIBLE
 #define M_TWOPI         (M_PI * 2.0)
 #define _IEEE_  __fdlibm_ieee
 #define _SVID_  __fdlibm_svid
 #define _XOPEN_ __fdlibm_xopen
 #define _POSIX_ __fdlibm_posix
-#endif /* ! defined (__STRICT_ANSI__) */
+#endif /* __BSD_VISIBLE */
 _END_STD_C

Subversion Repositories Kolibri OS

(root)/contrib/sdk/sources/newlib/libc/include/math.h @ 6897 – Rev 6099 → 6536