|
Lines 97-108
Link Here
|
| 97 |
} else if (ix<0x3fe00000) { /* |x|<0.5 */ |
97 |
} else if (ix<0x3fe00000) { /* |x|<0.5 */ |
| 98 |
if(ix<0x3e400000) { /* if |x| < 2**-27 */ |
98 |
if(ix<0x3e400000) { /* if |x| < 2**-27 */ |
| 99 |
if(huge+x>one) return x;/* return x with inexact if x!=0*/ |
99 |
if(huge+x>one) return x;/* return x with inexact if x!=0*/ |
| 100 |
} else |
100 |
} else { |
| 101 |
t = x*x; |
101 |
t = x*x; |
| 102 |
p = t*(pS0+t*(pS1+t*(pS2+t*(pS3+t*(pS4+t*pS5))))); |
102 |
p = t*(pS0+t*(pS1+t*(pS2+t*(pS3+t*(pS4+t*pS5))))); |
| 103 |
q = one+t*(qS1+t*(qS2+t*(qS3+t*qS4))); |
103 |
q = one+t*(qS1+t*(qS2+t*(qS3+t*qS4))); |
| 104 |
w = p/q; |
104 |
w = p/q; |
| 105 |
return x+x*w; |
105 |
return x+x*w; |
|
|
106 |
} |
| 106 |
} |
107 |
} |
| 107 |
/* 1> |x|>= 0.5 */ |
108 |
/* 1> |x|>= 0.5 */ |
| 108 |
w = one-fabs(x); |
109 |
w = one-fabs(x); |