summaryrefslogtreecommitdiff
path: root/libm/ldouble/fltestl.c
diff options
context:
space:
mode:
authorEric Andersen <andersen@codepoet.org>2001-05-10 00:40:28 +0000
committerEric Andersen <andersen@codepoet.org>2001-05-10 00:40:28 +0000
commit1077fa4d772832f77a677ce7fb7c2d513b959e3f (patch)
tree579bee13fb0b58d2800206366ec2caecbb15f3fc /libm/ldouble/fltestl.c
parent22358dd7ce7bb49792204b698f01a6f69b9c8e08 (diff)
uClibc now has a math library. muahahahaha!
-Erik
Diffstat (limited to 'libm/ldouble/fltestl.c')
-rw-r--r--libm/ldouble/fltestl.c265
1 files changed, 265 insertions, 0 deletions
diff --git a/libm/ldouble/fltestl.c b/libm/ldouble/fltestl.c
new file mode 100644
index 000000000..963e92467
--- /dev/null
+++ b/libm/ldouble/fltestl.c
@@ -0,0 +1,265 @@
+/* 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 );
+}