12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273 |
- #include "math.h"
- #include "math_private.h"
- double sin(double x)
- {
- double y[2],z=0.0;
- int32_t n, ix;
-
- GET_HIGH_WORD(ix,x);
-
- ix &= 0x7fffffff;
- if(ix <= 0x3fe921fb) return __kernel_sin(x,z,0);
-
- else if (ix>=0x7ff00000) return x-x;
-
- else {
- n = __ieee754_rem_pio2(x,y);
- switch(n&3) {
- case 0: return __kernel_sin(y[0],y[1],1);
- case 1: return __kernel_cos(y[0],y[1]);
- case 2: return -__kernel_sin(y[0],y[1],1);
- default:
- return -__kernel_cos(y[0],y[1]);
- }
- }
- }
- libm_hidden_def(sin)
|