This file is part of the GNU C Library.
Contributed by Jakub Jelinek <jj@ultra.linux.cz>
- This program is free software; you can redistribute it and/or modify
+ This program is free software: you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
- the Free Software Foundation; either version 2, or (at your option)
- any later version.
+ the Free Software Foundation; either version 3 of the License, or
+ (at your option) any later version.
This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License for more details.
- You should have received a copy of the GNU General Public License along
- with this program; if not, write to the Free Software Foundation,
- Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA. */
+ You should have received a copy of the GNU General Public License
+ along with this program. If not, see <http://www.gnu.org/licenses/>. */
#include <config.h>
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;
return -1;
}
- if (x + x == x || x != x) /* x is +=oo or NaN */
+ if (x + x == x) /* x is ±oo */
{
y[0] = x - x;
y[1] = y[0];