From 7ce331c01ce6eb7b3f5c715a38a24359da9c6ee2 Mon Sep 17 00:00:00 2001
From: Eric Andersen <andersen@codepoet.org>
Date: Thu, 22 Nov 2001 14:04:29 +0000
Subject: 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

---
 libm/ldouble/fltestl.c | 265 -------------------------------------------------
 1 file changed, 265 deletions(-)
 delete mode 100644 libm/ldouble/fltestl.c

(limited to 'libm/ldouble/fltestl.c')

diff --git a/libm/ldouble/fltestl.c b/libm/ldouble/fltestl.c
deleted file mode 100644
index 963e92467..000000000
--- a/libm/ldouble/fltestl.c
+++ /dev/null
@@ -1,265 +0,0 @@
-/* fltest.c
- * Test program for floor(), frexp(), ldexp()
- */
-
-/*
-Cephes Math Library Release 2.1:  December, 1988
-Copyright 1984, 1987, 1988 by Stephen L. Moshier (moshier@world.std.com)
-*/
-
-
-
-/*#include <math.h>*/
-#define MACHEPL  5.42101086242752217003726400434970855712890625E-20L
-#define N 16300
-
-void flierr();
-int printf();
-void exit();
-
-int
-main()
-{
-long double x, y, y0, z, f, x00, y00;
-int i, j, e, e0;
-int errfr, errld, errfl, underexp, err, errth, e00;
-long double frexpl(), ldexpl(), floorl();
-
-
-/*
-if( 1 )
-	goto flrtst;
-*/
-
-printf( "Testing frexpl() and ldexpl().\n" );
-errth = 0.0L;
-errfr = 0;
-errld = 0;
-underexp = 0;
-f = 1.0L;
-x00 = 2.0L;
-y00 = 0.5L;
-e00 = 2;
-
-for( j=0; j<20; j++ )
-{
-if( j == 10 )
-	{
-	f = 1.0L;
-	x00 = 2.0L;
-	e00 = 1;
-/* Find 2**(2**14) / 2 */
-	for( i=0; i<13; i++ )
-		{
-		x00 *= x00;
-		e00 += e00;
-		}
-	y00 = x00/2.0L;
-	x00 = x00 * y00;
-	e00 += e00;
-	y00 = 0.5L;
-	}
-x = x00 * f;
-y0 = y00 * f;
-e0 = e00;
-
-#if 1
-/* If ldexp, frexp support denormal numbers, this should work.  */
-for( i=0; i<16448; i++ )
-#else
-for( i=0; i<16383; i++ )
-#endif
-	{
-	x /= 2.0L;
-	e0 -= 1;
-	if( x == 0.0L )
-		{
-		if( f == 1.0L )
-			underexp = e0;
-		y0 = 0.0L;
-		e0 = 0;
-		}
-	y = frexpl( x, &e );
-	if( (e0 < -16383) && (e != e0) )
-		{
-		if( e == (e0 - 1) )
-			{
-			e += 1;
-			y /= 2.0L;
-			}
-		if( e == (e0 + 1) )
-			{
-			e -= 1;
-			y *= 2.0L;
-			}
-		}
-	err = y - y0;
-	if( y0 != 0.0L )
-		err /= y0;
-	if( err < 0.0L )
-		err = -err;
-	if( e0 > -1023 )
-		errth = 0.0L;
-	else
-		{/* Denormal numbers may have rounding errors */
-		if( e0 == -16383 )
-			{
-			errth = 2.0L * MACHEPL;
-			}
-		else
-			{
-			errth *= 2.0L;
-			}
-		}
-
-	if( (x != 0.0L) && ((err > errth) || (e != e0)) )
-		{
-		printf( "Test %d: ", j+1 );
-		printf( " frexpl( %.20Le) =?= %.20Le * 2**%d;", x, y, e );
-		printf( " should be %.20Le * 2**%d\n", y0, e0 );
-		errfr += 1;
-		}
-	y = ldexpl( x, 1-e0 );
-	err = y - 1.0L;
-	if( err < 0.0L )
-		err = -err;
-	if( (err > errth) && ((x == 0.0L) && (y != 0.0L)) )
-		{
-		printf( "Test %d: ", j+1 );
-		printf( "ldexpl( %.15Le, %d ) =?= %.15Le;", x, 1-e0, y );
-		if( x != 0.0L )
-			printf( " should be %.15Le\n", f );
-		else
-			printf( " should be %.15Le\n", 0.0L );
-		errld += 1;
-		}
-	if( x == 0.0L )
-		{
-		break;
-		}
-	}
-f = f * 1.08005973889L;
-}
-
-if( (errld == 0) && (errfr == 0) )
-	{
-	printf( "No errors found.\n" );
-	}
-
-/*flrtst:*/
-
-printf( "Testing floorl().\n" );
-errfl = 0;
-
-f = 1.0L/MACHEPL;
-x00 = 1.0L;
-for( j=0; j<57; j++ )
-{
-x = x00 - 1.0L;
-for( i=0; i<128; i++ )
-	{
-	y = floorl(x);
-	if( y != x )
-		{
-		flierr( x, y, j );
-		errfl += 1;
-		}
-/* Warning! the if() statement is compiler dependent,
- * since x-0.49 may be held in extra precision accumulator
- * so would never compare equal to x!  The subroutine call
- * y = floor() forces z to be stored as a double and reloaded
- * for the if() statement.
- */
-	z = x - 0.49L;
-	y = floorl(z);
-	if( z == x )
-		break;
-	if( y != (x - 1.0L) )
-		{
-		flierr( z, y, j );
-		errfl += 1;
-		}
-
-	z = x + 0.49L;
-	y = floorl(z);
-	if( z != x )
-		{
-		if( y != x )
-			{
-			flierr( z, y, j );
-			errfl += 1;
-			}
-		}
-	x = -x;
-	y = floorl(x);
-	if( z != x )
-		{
-		if( y != x )
-			{
-			flierr( x, y, j );
-			errfl += 1;
-			}
-		}
-	z = x + 0.49L;
-	y = floorl(z);
-	if( z != x )
-		{
-		if( y != x )
-			{
-			flierr( z, y, j );
-			errfl += 1;
-			}
-		}
-	z = x - 0.49L;
-	y = floorl(z);
-	if( z != x )
-		{
-		if( y != (x - 1.0L) )
-			{
-			flierr( z, y, j );
-			errfl += 1;
-			}
-		}
-	x = -x;
-	x += 1.0L;
-	}
-x00 = x00 + x00;
-}
-y = floorl(0.0L);
-if( y != 0.0L )
-	{
-	flierr( 0.0L, y, 57 );
-	errfl += 1;
-	}
-y = floorl(-0.0L);
-if( y != 0.0L )
-	{
-	flierr( -0.0L, y, 58 );
-	errfl += 1;
-	}
-y = floorl(-1.0L);
-if( y != -1.0L )
-	{
-	flierr( -1.0L, y, 59 );
-	errfl += 1;
-	}
-y = floorl(-0.1L);
-if( y != -1.0l )
-	{
-	flierr( -0.1L, y, 60 );
-	errfl += 1;
-	}
-
-if( errfl == 0 )
-	printf( "No errors found in floorl().\n" );
-exit(0);
-return 0;
-}
-
-void flierr( x, y, k )
-long double x, y;
-int k;
-{
-printf( "Test %d: ", k+1 );
-printf( "floorl(%.15Le) =?= %.15Le\n", x, y );
-}
-- 
cgit v1.2.3