Index: fusl/src/complex/catanl.c |
diff --git a/fusl/src/complex/catanl.c b/fusl/src/complex/catanl.c |
index 5ace7704fdc086461216389aa83b980c161e2c3f..4a24ff578ff05a4565a797e40ab2d1ad9388ced9 100644 |
--- a/fusl/src/complex/catanl.c |
+++ b/fusl/src/complex/catanl.c |
@@ -62,9 +62,8 @@ |
#include "libm.h" |
#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024 |
-long double complex catanl(long double complex z) |
-{ |
- return catan(z); |
+long double complex catanl(long double complex z) { |
+ return catan(z); |
} |
#else |
static const long double PIL = 3.141592653589793238462643383279502884197169L; |
@@ -72,55 +71,53 @@ static const long double DP1 = 3.14159265358979323829596852490908531763125L; |
static const long double DP2 = 1.6667485837041756656403424829301998703007e-19L; |
static const long double DP3 = 1.8830410776607851167459095484560349402753e-39L; |
-static long double redupil(long double x) |
-{ |
- long double t; |
- long i; |
+static long double redupil(long double x) { |
+ long double t; |
+ long i; |
- t = x / PIL; |
- if (t >= 0.0L) |
- t += 0.5L; |
- else |
- t -= 0.5L; |
+ t = x / PIL; |
+ if (t >= 0.0L) |
+ t += 0.5L; |
+ else |
+ t -= 0.5L; |
- i = t; /* the multiple */ |
- t = i; |
- t = ((x - t * DP1) - t * DP2) - t * DP3; |
- return t; |
+ i = t; /* the multiple */ |
+ t = i; |
+ t = ((x - t * DP1) - t * DP2) - t * DP3; |
+ return t; |
} |
-long double complex catanl(long double complex z) |
-{ |
- long double complex w; |
- long double a, t, x, x2, y; |
+long double complex catanl(long double complex z) { |
+ long double complex w; |
+ long double a, t, x, x2, y; |
- x = creall(z); |
- y = cimagl(z); |
+ x = creall(z); |
+ y = cimagl(z); |
- if ((x == 0.0L) && (y > 1.0L)) |
- goto ovrf; |
+ if ((x == 0.0L) && (y > 1.0L)) |
+ goto ovrf; |
- x2 = x * x; |
- a = 1.0L - x2 - (y * y); |
- if (a == 0.0L) |
- goto ovrf; |
+ x2 = x * x; |
+ a = 1.0L - x2 - (y * y); |
+ if (a == 0.0L) |
+ goto ovrf; |
- t = atan2l(2.0L * x, a) * 0.5L; |
- w = redupil(t); |
+ t = atan2l(2.0L * x, a) * 0.5L; |
+ w = redupil(t); |
- t = y - 1.0L; |
- a = x2 + (t * t); |
- if (a == 0.0L) |
- goto ovrf; |
+ t = y - 1.0L; |
+ a = x2 + (t * t); |
+ if (a == 0.0L) |
+ goto ovrf; |
- t = y + 1.0L; |
- a = (x2 + (t * t)) / a; |
- w = w + (0.25L * logl(a)) * I; |
- return w; |
+ t = y + 1.0L; |
+ a = (x2 + (t * t)) / a; |
+ w = w + (0.25L * logl(a)) * I; |
+ return w; |
ovrf: |
- // FIXME |
- w = LDBL_MAX + LDBL_MAX * I; |
- return w; |
+ // FIXME |
+ w = LDBL_MAX + LDBL_MAX * I; |
+ return w; |
} |
#endif |