int exp, n;
if (x >= -0.78539816339744830961566084581987572104929234984377
- && x < 0.78539816339744830961566084581987572104929234984377)
+ && x <= 0.78539816339744830961566084581987572104929234984377)
/* x in <-pi/4, pi/4> */
{
y[0] = x;
return 0;
}
- if (x >= 2.35619449019234492884698253745962716314787704953131
- && x < 2.35619449019234492884698253745962716314787704953131)
- if (x > 0)
+ if (x > 0 && x < 2.35619449019234492884698253745962716314787704953131)
{
/* 113 + 93 bit PI is ok */
z = x - PI_2_1;
y[1] = (z - y[0]) - PI_2_1t;
return 1;
}
- else
+
+ if (x < 0 && x > -2.35619449019234492884698253745962716314787704953131)
{
/* 113 + 93 bit PI is ok */
z = x + PI_2_1;