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/float/tandgf.c | 206 ---------------------------------------------------- 1 file changed, 206 deletions(-) delete mode 100644 libm/float/tandgf.c (limited to 'libm/float/tandgf.c') 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 - -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) ); -} - -- cgit v1.2.3