| #include <math.h> |
| #include <stdint.h> |
| |
| float hypotf(float x, float y) { |
| union { |
| float f; |
| uint32_t i; |
| } ux = {x}, uy = {y}, ut; |
| float_t z; |
| |
| ux.i &= -1U >> 1; |
| uy.i &= -1U >> 1; |
| if (ux.i < uy.i) { |
| ut = ux; |
| ux = uy; |
| uy = ut; |
| } |
| |
| x = ux.f; |
| y = uy.f; |
| if (uy.i == 0xff << 23) |
| return y; |
| if (ux.i >= 0xff << 23 || uy.i == 0 || ux.i - uy.i >= 25 << 23) |
| return x + y; |
| |
| z = 1; |
| if (ux.i >= (0x7f + 60) << 23) { |
| z = 0x1p90f; |
| x *= 0x1p-90f; |
| y *= 0x1p-90f; |
| } else if (uy.i < (0x7f - 60) << 23) { |
| z = 0x1p-90f; |
| x *= 0x1p90f; |
| y *= 0x1p90f; |
| } |
| return z * sqrtf((double)x * x + (double)y * y); |
| } |
