Totally rework the math library, this time based on the MacOs X

math library (which is itself based on the math lib from FreeBSD).
 -Erik

1 files changed, 0 insertions, 186 deletions

diff --git a/libm/ldouble/isnanl.c b/libm/ldouble/isnanl.c
deleted file mode 100644
index 44158ecc7..000000000
--- a/libm/ldouble/isnanl.c
+++ /dev/null
@@ -1,186 +0,0 @@
-/*							isnanl()
- *							isfinitel()
- *							signbitl()
- *
- *	Floating point IEEE special number tests
- *
- *
- *
- * SYNOPSIS:
- *
- * int signbitl(), isnanl(), isfinitel();
- * long double x, y;
- *
- * n = signbitl(x);
- * n = isnanl(x);
- * n = isfinitel(x);
- *
- *
- *
- * DESCRIPTION:
- *
- * These functions are part of the standard C run time library
- * for some but not all C compilers.  The ones supplied are
- * written in C for IEEE arithmetic.  They should
- * be used only if your compiler library does not already have
- * them.
- *
- */
-
-
-/*
-Cephes Math Library Release 2.7:  June, 1998
-Copyright 1992, 1998 by Stephen L. Moshier
-*/
-
-
-#include <math.h>
-
-/* This is defined in mconf.h. */
-/* #define DENORMAL 1 */
-
-#ifdef UNK
-/* Change UNK into something else.  */
-#undef UNK
-#if BIGENDIAN
-#define MIEEE 1
-#else
-#define IBMPC 1
-#endif
-#endif
-
-
-/* Return 1 if the sign bit of x is 1, else 0.  */
-
-int signbitl(x)
-long double x;
-{
-union
-	{
-	long double d;
-	short s[6];
-	int i[3];
-	} u;
-
-u.d = x;
-
-if( sizeof(int) == 4 )
-	{
-#ifdef IBMPC
-	return( u.s[4] < 0 );
-#endif
-#ifdef MIEEE
-	return( u.i[0] < 0 );
-#endif
-	}
-else
-	{
-#ifdef IBMPC
-	return( u.s[4] < 0 );
-#endif
-#ifdef MIEEE
-	return( u.s[0] < 0 );
-#endif
-	}
-}
-
-
-/* Return 1 if x is a number that is Not a Number, else return 0.  */
-
-int isnanl(x)
-long double x;
-{
-#ifdef NANS
-union
-	{
-	long double d;
-	unsigned short s[6];
-	unsigned int i[3];
-	} u;
-
-u.d = x;
-
-if( sizeof(int) == 4 )
-	{
-#ifdef IBMPC
-	if( ((u.s[4] & 0x7fff) == 0x7fff)
-	    && (((u.i[1] & 0x7fffffff)!= 0) || (u.i[0] != 0)))
-		return 1;
-#endif
-#ifdef MIEEE
-	if( ((u.i[0] & 0x7fff0000) == 0x7fff0000)
-	    && (((u.i[1] & 0x7fffffff) != 0) || (u.i[2] != 0)))
-		return 1;
-#endif
-	return(0);
-	}
-else
-	{ /* size int not 4 */
-#ifdef IBMPC
-	if( (u.s[4] & 0x7fff) == 0x7fff)
-		{
-		if((u.s[3] & 0x7fff) | u.s[2] | u.s[1] | u.s[0])
-			return(1);
-		}
-#endif
-#ifdef MIEEE
-	if( (u.s[0] & 0x7fff) == 0x7fff)
-		{
-		if((u.s[2] & 0x7fff) | u.s[3] | u.s[4] | u.s[5])
-			return(1);
-		}
-#endif
-	return(0);
-	} /* size int not 4 */
-
-#else
-/* No NANS.  */
-return(0);
-#endif
-}
-
-
-/* Return 1 if x is not infinite and is not a NaN.  */
-
-int isfinitel(x)
-long double x;
-{
-#ifdef INFINITIES
-union
-	{
-	long double d;
-	unsigned short s[6];
-	unsigned int i[3];
-	} u;
-
-u.d = x;
-
-if( sizeof(int) == 4 )
-	{
-#ifdef IBMPC
-	if( (u.s[4] & 0x7fff) != 0x7fff)
-		return 1;
-#endif
-#ifdef MIEEE
-	if( (u.i[0] & 0x7fff0000) != 0x7fff0000)
-		return 1;
-#endif
-	return(0);
-	}
-else
-	{
-#ifdef IBMPC
-	if( (u.s[4] & 0x7fff) != 0x7fff)
-		return 1;
-#endif
-#ifdef MIEEE
-	if( (u.s[0] & 0x7fff) != 0x7fff)
-		return 1;
-#endif
-	return(0);
-	}
-#else
-/* No INFINITY.  */
-return(1);
-#endif
-}