| Index: fusl/src/math/tanhf.c
|
| diff --git a/fusl/src/math/tanhf.c b/fusl/src/math/tanhf.c
|
| index 10636fbd7be6c21d1c58326eba2bf7179b0dea8f..7c57ccf37b98691196355889bb7c07d42ef3e994 100644
|
| --- a/fusl/src/math/tanhf.c
|
| +++ b/fusl/src/math/tanhf.c
|
| @@ -1,39 +1,41 @@
|
| #include "libm.h"
|
|
|
| -float tanhf(float x)
|
| -{
|
| - union {float f; uint32_t i;} u = {.f = x};
|
| - uint32_t w;
|
| - int sign;
|
| - float t;
|
| +float tanhf(float x) {
|
| + union {
|
| + float f;
|
| + uint32_t i;
|
| + } u = {.f = x};
|
| + uint32_t w;
|
| + int sign;
|
| + float t;
|
|
|
| - /* x = |x| */
|
| - sign = u.i >> 31;
|
| - u.i &= 0x7fffffff;
|
| - x = u.f;
|
| - w = u.i;
|
| + /* x = |x| */
|
| + sign = u.i >> 31;
|
| + u.i &= 0x7fffffff;
|
| + x = u.f;
|
| + w = u.i;
|
|
|
| - if (w > 0x3f0c9f54) {
|
| - /* |x| > log(3)/2 ~= 0.5493 or nan */
|
| - if (w > 0x41200000) {
|
| - /* |x| > 10 */
|
| - t = 1 + 0/x;
|
| - } else {
|
| - t = expm1f(2*x);
|
| - t = 1 - 2/(t+2);
|
| - }
|
| - } else if (w > 0x3e82c578) {
|
| - /* |x| > log(5/3)/2 ~= 0.2554 */
|
| - t = expm1f(2*x);
|
| - t = t/(t+2);
|
| - } else if (w >= 0x00800000) {
|
| - /* |x| >= 0x1p-126 */
|
| - t = expm1f(-2*x);
|
| - t = -t/(t+2);
|
| - } else {
|
| - /* |x| is subnormal */
|
| - FORCE_EVAL(x*x);
|
| - t = x;
|
| - }
|
| - return sign ? -t : t;
|
| + if (w > 0x3f0c9f54) {
|
| + /* |x| > log(3)/2 ~= 0.5493 or nan */
|
| + if (w > 0x41200000) {
|
| + /* |x| > 10 */
|
| + t = 1 + 0 / x;
|
| + } else {
|
| + t = expm1f(2 * x);
|
| + t = 1 - 2 / (t + 2);
|
| + }
|
| + } else if (w > 0x3e82c578) {
|
| + /* |x| > log(5/3)/2 ~= 0.2554 */
|
| + t = expm1f(2 * x);
|
| + t = t / (t + 2);
|
| + } else if (w >= 0x00800000) {
|
| + /* |x| >= 0x1p-126 */
|
| + t = expm1f(-2 * x);
|
| + t = -t / (t + 2);
|
| + } else {
|
| + /* |x| is subnormal */
|
| + FORCE_EVAL(x * x);
|
| + t = x;
|
| + }
|
| + return sign ? -t : t;
|
| }
|
|
|