Index: fusl/src/math/fmaxl.c |
diff --git a/fusl/src/math/fmaxl.c b/fusl/src/math/fmaxl.c |
new file mode 100644 |
index 0000000000000000000000000000000000000000..4b03158e06292b28576615a2a1dd36bd6bf1d34a |
--- /dev/null |
+++ b/fusl/src/math/fmaxl.c |
@@ -0,0 +1,21 @@ |
+#include <math.h> |
+#include <float.h> |
+ |
+#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024 |
+long double fmaxl(long double x, long double y) |
+{ |
+ return fmax(x, y); |
+} |
+#else |
+long double fmaxl(long double x, long double y) |
+{ |
+ if (isnan(x)) |
+ return y; |
+ if (isnan(y)) |
+ return x; |
+ /* handle signed zeros, see C99 Annex F.9.9.2 */ |
+ if (signbit(x) != signbit(y)) |
+ return signbit(x) ? y : x; |
+ return x < y ? y : x; |
+} |
+#endif |