mirrors/ecl
1
Fork 0
mirror of https://gitlab.com/embeddable-common-lisp/ecl.git synced 2026-01-25 14:01:07 -08:00
Code Issues Projects Releases Packages Wiki Activity

atan: fix if-clause order to work on signed zero (2)

Browse source 
Fixes #329. This time for real :-)
This commit is contained in:
Daniel Kochmański 2016-12-20 17:56:20 +01:00
parent a302a2e36c
commit 4ba367531e
1 changed files with 20 additions and 18 deletions
Download patch file Download diff file Expand all files Collapse all files

38
src/c/numbers/atan.d
View file

@ -23,28 +23,28 @@ static double
ecl_atan2_double(double y, double x)
{
if (x == 0) {
if (signbit(y)) {
return -ECL_PI2_D;
} else if (y == 0) {
if (y == 0) {
/* Signals floating-point-invalid-operation. If not trapped,
produces NaN. */
return x / y;
} else if (signbit(y)) {
return -ECL_PI2_D;
} else {
return ECL_PI2_D;
}
} else if (signbit(x)) {
if (signbit(y)) {
return -ECL_PI_D + atan(-y / -x);
} else if (y == 0) {
if (y == 0) {
return ECL_PI_D;
} else if (signbit(y)) {
return -ECL_PI_D + atan(-y / -x);
} else {
return ECL_PI_D - atan(y / -x);
}
} else {
if (signbit(y)) {
return -atan(-y / x);
} else if (y == 0) {
if (y == 0) {
return (double)0;
} else if (signbit(y)) {
return -atan(-y / x);
} else {
return atan(y / x);
}
@ -56,29 +56,31 @@ static long double
ecl_atan2_long_double(long double y, long double x)
{
if (x == 0) {
if (signbit(y)) {
return -ECL_PI2_L;
} else if (y == 0) {
if (y == 0) {
/* Signals floating-point-invalid-operation. If not trapped,
produces NaN. */
return x / y;
}
else if (signbit(y)) {
return -ECL_PI2_L;
} else {
return ECL_PI2_L;
}
}
else if (signbit(x)) {
if (signbit(y)) {
return -ECL_PI_L + atanl(-y / -x);
} else if (y == 0) {
if (y == 0) {
return ECL_PI_L;
}
else if (signbit(y)) {
return -ECL_PI_L + atanl(-y / -x);
} else {
return ECL_PI_L - atanl(y / -x);
}
} else {
if (signbit(y)) {
return -atanl(-y / x);
} else if (y == 0) {
if (y == 0) {
return (long double)0;
} else if (signbit(y)) {
return -atanl(-y / x);
} else {
return atanl(y / x);
}