d/freebsd-skq
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

* Use ENTERI() and RETURNI() to toggle the rounding precision if

Browse Source 
necessary, so that cosl(), sinl() and tanl() work on i386 even
  for naive callers.

Suggested by:	bde
Reviewed by:	bde
Approved by: 	das (mentor)
This commit is contained in:
Steve Kargl 2012-09-22 15:38:29 +00:00
parent e6f9129aa0
commit 340076f0f6
3 changed files with 21 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
lib/msun/src/s_cosl.c
View File

@ -33,6 +33,9 @@ __FBSDID("$FreeBSD$");
*/
#include <float.h>
#ifdef __i386__
#include <ieeefp.h>
#endif
#include "math.h"
#include "math_private.h"
@ -63,9 +66,11 @@ cosl(long double x)
if (z.bits.exp == 32767)
return ((x - x) / (x - x));
ENTERI();
/* Optimize the case where x is already within range. */
if (z.e < M_PI_4)
return (__kernel_cosl(z.e, 0));
RETURNI(__kernel_cosl(z.e, 0));
e0 = __ieee754_rem_pio2l(x, y);
hi = y[0];
@ -86,5 +91,5 @@ cosl(long double x)
break;
}
return (hi);
RETURNI(hi);
}

9
lib/msun/src/s_sinl.c
View File

@ -28,6 +28,9 @@
__FBSDID("$FreeBSD$");
#include <float.h>
#ifdef __i386__
#include <ieeefp.h>
#endif
#include "math.h"
#include "math_private.h"
@ -59,10 +62,12 @@ sinl(long double x)
if (z.bits.exp == 32767)
return ((x - x) / (x - x));
ENTERI();
/* Optimize the case where x is already within range. */
if (z.e < M_PI_4) {
hi = __kernel_sinl(z.e, 0, 0);
return (s ? -hi : hi);
RETURNI(s ? -hi : hi);
}
e0 = __ieee754_rem_pio2l(x, y);
@ -84,5 +89,5 @@ sinl(long double x)
break;
}
return (hi);
RETURNI(hi);
}

9
lib/msun/src/s_tanl.c
View File

@ -34,6 +34,9 @@ __FBSDID("$FreeBSD$");
*/
#include <float.h>
#ifdef __i386__
#include <ieeefp.h>
#endif
#include "math.h"
#include "math_private.h"
@ -65,10 +68,12 @@ tanl(long double x)
if (z.bits.exp == 32767)
return ((x - x) / (x - x));
ENTERI();
/* Optimize the case where x is already within range. */
if (z.e < M_PI_4) {
hi = __kernel_tanl(z.e, 0, 0);
return (s ? -hi : hi);
RETURNI(s ? -hi : hi);
}
e0 = __ieee754_rem_pio2l(x, y);
@ -86,5 +91,5 @@ tanl(long double x)
break;
}
return (hi);
RETURNI(hi);
}