summaryrefslogtreecommitdiff
path: root/libm/float/tandgf.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/float/tandgf.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/float/tandgf.c')
-rw-r--r--libm/float/tandgf.c206
1 files changed, 0 insertions, 206 deletions
diff --git a/libm/float/tandgf.c b/libm/float/tandgf.c
deleted file mode 100644
index dc55ad5e4..000000000
--- a/libm/float/tandgf.c
+++ /dev/null
@@ -1,206 +0,0 @@
-/* tandgf.c
- *
- * Circular tangent of angle in degrees
- *
- *
- *
- * SYNOPSIS:
- *
- * float x, y, tandgf();
- *
- * y = tandgf( x );
- *
- *
- *
- * DESCRIPTION:
- *
- * Returns the circular tangent of the radian argument x.
- *
- * Range reduction is into intervals of 45 degrees.
- *
- *
- *
- *
- * ACCURACY:
- *
- * Relative error:
- * arithmetic domain # trials peak rms
- * IEEE +-2^24 50000 2.4e-7 4.8e-8
- *
- * ERROR MESSAGES:
- *
- * message condition value returned
- * tanf total loss x > 2^24 0.0
- *
- */
- /* cotdgf.c
- *
- * Circular cotangent of angle in degrees
- *
- *
- *
- * SYNOPSIS:
- *
- * float x, y, cotdgf();
- *
- * y = cotdgf( x );
- *
- *
- *
- * DESCRIPTION:
- *
- * Range reduction is into intervals of 45 degrees.
- * A common routine computes either the tangent or cotangent.
- *
- *
- *
- * ACCURACY:
- *
- * Relative error:
- * arithmetic domain # trials peak rms
- * IEEE +-2^24 50000 2.4e-7 4.8e-8
- *
- *
- * ERROR MESSAGES:
- *
- * message condition value returned
- * cot total loss x > 2^24 0.0
- * cot singularity x = 0 MAXNUMF
- *
- */
-
-/*
-Cephes Math Library Release 2.2: June, 1992
-Copyright 1984, 1987, 1989, 1992 by Stephen L. Moshier
-Direct inquiries to 30 Frost Street, Cambridge, MA 02140
-*/
-
-/* Single precision circular tangent
- * test interval: [-pi/4, +pi/4]
- * trials: 10000
- * peak relative error: 8.7e-8
- * rms relative error: 2.8e-8
- */
-#include <math.h>
-
-extern float MAXNUMF;
-
-static float T24M1 = 16777215.;
-static float PI180 = 0.0174532925199432957692; /* pi/180 */
-
-static float tancotf( float xx, int cotflg )
-{
-float x, y, z, zz;
-long j;
-int sign;
-
-
-/* make argument positive but save the sign */
-if( xx < 0.0 )
- {
- x = -xx;
- sign = -1;
- }
-else
- {
- x = xx;
- sign = 1;
- }
-
-if( x > T24M1 )
- {
- if( cotflg )
- mtherr( "cotdgf", TLOSS );
- else
- mtherr( "tandgf", TLOSS );
- return(0.0);
- }
-
-/* compute x mod PIO4 */
-j = 0.022222222222222222222 * x; /* integer part of x/45 */
-y = j;
-
-/* map zeros and singularities to origin */
-if( j & 1 )
- {
- j += 1;
- y += 1.0;
- }
-
-z = x - y * 45.0;
-z *= PI180; /* multiply by pi/180 to convert to radians */
-
-zz = z * z;
-
-if( x > 1.0e-4 )
- {
-/* 1.7e-8 relative error in [-pi/4, +pi/4] */
- y =
- ((((( 9.38540185543E-3 * zz
- + 3.11992232697E-3) * zz
- + 2.44301354525E-2) * zz
- + 5.34112807005E-2) * zz
- + 1.33387994085E-1) * zz
- + 3.33331568548E-1) * zz * z
- + z;
- }
-else
- {
- y = z;
- }
-
-if( j & 2 )
- {
- if( cotflg )
- y = -y;
- else
- {
- if( y != 0.0 )
- {
- y = -1.0/y;
- }
- else
- {
- mtherr( "tandgf", SING );
- y = MAXNUMF;
- }
- }
- }
-else
- {
- if( cotflg )
- {
- if( y != 0.0 )
- y = 1.0/y;
- else
- {
- mtherr( "cotdgf", SING );
- y = MAXNUMF;
- }
- }
- }
-
-if( sign < 0 )
- y = -y;
-
-return( y );
-}
-
-
-float tandgf( float x )
-{
-
-return( tancotf(x,0) );
-}
-
-float cotdgf( float x )
-{
-
-if( x == 0.0 )
- {
- mtherr( "cotdgf", SING );
- return( MAXNUMF );
- }
-return( tancotf(x,1) );
-}
-