From 7ce331c01ce6eb7b3f5c715a38a24359da9c6ee2 Mon Sep 17 00:00:00 2001 From: Eric Andersen 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/double/eltst.c | 37 ------------------------------------- 1 file changed, 37 deletions(-) delete mode 100644 libm/double/eltst.c (limited to 'libm/double/eltst.c') diff --git a/libm/double/eltst.c b/libm/double/eltst.c deleted file mode 100644 index cef249eaf..000000000 --- a/libm/double/eltst.c +++ /dev/null @@ -1,37 +0,0 @@ -extern double MACHEP, PIO2, PI; -double ellie(), ellpe(), floor(), fabs(); -double ellie2(); - -main() -{ -double y, m, phi, e, E, phipi, y1; -int i, j, npi; - -/* dprec(); */ -m = 0.9; -E = ellpe(0.1); -for( j=-10; j<=10; j++ ) - { - printf( "%d * PIO2\n", j ); - for( i=-2; i<=2; i++ ) - { - phi = PIO2 * j + 50 * MACHEP * i; - npi = floor(phi/PIO2); - if( npi & 1 ) - npi += 1; - phipi = phi - npi * PIO2; - npi = floor(phi/PIO2); - if( npi & 1 ) - npi += 1; - phipi = phi - npi * PIO2; - printf( "phi %.9e npi %d ", phi, npi ); - y1 = E * npi + ellie(phipi,m); - y = ellie2( phi, m ); - printf( "y %.9e ", y ); - e = fabs(y - y1); - if( y1 != 0.0 ) - e /= y1; - printf( "e %.4e\n", e ); - } - } -} -- cgit v1.2.3