summaryrefslogtreecommitdiff
path: root/sysdeps/m68k/fpu/e_atan2.c
diff options
context:
space:
mode:
Diffstat (limited to 'sysdeps/m68k/fpu/e_atan2.c')
-rw-r--r--sysdeps/m68k/fpu/e_atan2.c52
1 files changed, 28 insertions, 24 deletions
diff --git a/sysdeps/m68k/fpu/e_atan2.c b/sysdeps/m68k/fpu/e_atan2.c
index c012070a93..59bc990f5b 100644
--- a/sysdeps/m68k/fpu/e_atan2.c
+++ b/sysdeps/m68k/fpu/e_atan2.c
@@ -35,63 +35,67 @@ float_type
s(__ieee754_atan2) (float_type y, float_type x)
{
float_type pi, pi_2, z;
+ unsigned long y_cond, x_cond;
__asm ("fmovecr%.x %#0, %0" : "=f" (pi));
__asm ("fscale%.w %#-1, %0" : "=f" (pi_2) : "0" (pi));
- if (x != x || y != y)
+ y_cond = __m81_test (y);
+ x_cond = __m81_test (x);
+
+ if ((x_cond | y_cond) & __M81_COND_NAN)
z = x + y;
- else if (y == 0)
+ else if (y_cond & __M81_COND_ZERO)
{
- if (m81(__signbit) (x))
- z = m81(__signbit) (y) ? -pi : pi;
+ if (x_cond & __M81_COND_NEG)
+ z = y_cond & __M81_COND_NEG ? -pi : pi;
else
z = y;
}
- else if (m81(__isinf) (x))
+ else if (x_cond & __M81_COND_INF)
{
- if (m81(__isinf) (y))
+ if (y_cond & __M81_COND_INF)
{
float_type pi_4;
__asm ("fscale%.w %#-2, %0" : "=f" (pi_4) : "0" (pi));
- z = x > 0 ? pi_4 : 3 * pi_4;
+ z = x_cond & __M81_COND_NEG ? 3 * pi_4 : pi_4;
}
else
- z = x > 0 ? 0 : pi;
- if (m81(__signbit) (y))
+ z = x_cond & __M81_COND_NEG ? pi : 0;
+ if (y_cond & __M81_COND_NEG)
z = -z;
}
- else if (m81(__isinf) (y))
- z = y > 0 ? pi_2 : -pi_2;
- else if (x > 0)
+ else if (y_cond & __M81_COND_INF)
+ z = y_cond & __M81_COND_NEG ? -pi_2 : pi_2;
+ else if (x_cond & __M81_COND_NEG)
{
- if (y > 0)
+ if (y_cond & __M81_COND_NEG)
{
- if (x > y)
- z = m81(__atan) (y / x);
+ if (-x > -y)
+ z = -pi + m81(__atan) (y / x);
else
- z = pi_2 - m81(__atan) (x / y);
+ z = -pi_2 - m81(__atan) (x / y);
}
else
{
- if (x > -y)
- z = m81(__atan) (y / x);
+ if (-x > y)
+ z = pi + m81(__atan) (y / x);
else
- z = -pi_2 - m81(__atan) (x / y);
+ z = pi_2 - m81(__atan) (x / y);
}
}
else
{
- if (y < 0)
+ if (y_cond & __M81_COND_NEG)
{
- if (-x > y)
- z = -pi + m81(__atan) (y / x);
+ if (x > -y)
+ z = m81(__atan) (y / x);
else
z = -pi_2 - m81(__atan) (x / y);
}
else
{
- if (-x > y)
- z = pi + m81(__atan) (y / x);
+ if (x > y)
+ z = m81(__atan) (y / x);
else
z = pi_2 - m81(__atan) (x / y);
}