Index: fusl/src/math/hypotf.c |
diff --git a/fusl/src/math/hypotf.c b/fusl/src/math/hypotf.c |
new file mode 100644 |
index 0000000000000000000000000000000000000000..2fc214b7232e09ab52d5f8b5bb1b10f024b9b8d3 |
--- /dev/null |
+++ b/fusl/src/math/hypotf.c |
@@ -0,0 +1,35 @@ |
+#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); |
+} |