summaryrefslogtreecommitdiff
path: root/libm/double/eltst.c
diff options
context:
space:
mode:
authorEric Andersen <andersen@codepoet.org>2001-11-22 14:04:29 +0000
committerEric Andersen <andersen@codepoet.org>2001-11-22 14:04:29 +0000
commit7ce331c01ce6eb7b3f5c715a38a24359da9c6ee2 (patch)
tree3a7e8476e868ae15f4da1b7ce26b2db6f434468c /libm/double/eltst.c
parentc117dd5fb183afb1a4790a6f6110d88704be6bf8 (diff)
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
Diffstat (limited to 'libm/double/eltst.c')
-rw-r--r--libm/double/eltst.c37
1 files changed, 0 insertions, 37 deletions
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 );
- }
- }
-}