From 845e4f6692336ea7f7672045dee05fa39c4022ec Mon Sep 17 00:00:00 2001
From: Rich Felker <dalias@aerifal.cx>
Date: Sun, 2 Aug 2020 23:31:51 -0400
Subject: add m68k sqrtl using native instruction
this is actually a functional fix at present, since the C sqrtl does
not support ld80 and just wraps double sqrt. once that's fixed it will
just be an optimization.
---
src/math/m68k/sqrtl.c | 15 +++++++++++++++
1 file changed, 15 insertions(+)
create mode 100644 src/math/m68k/sqrtl.c
(limited to 'src/math')
diff --git a/src/math/m68k/sqrtl.c b/src/math/m68k/sqrtl.c
new file mode 100644
index 00000000..b1c303c7
--- /dev/null
+++ b/src/math/m68k/sqrtl.c
@@ -0,0 +1,15 @@
+#include <math.h>
+
+#if __HAVE_68881__
+
+long double sqrtl(long double x)
+{
+ __asm__ ("fsqrt.x %1,%0" : "=f"(x) : "fm"(x));
+ return x;
+}
+
+#else
+
+#include "../sqrtl.c"
+
+#endif
--
cgit v1.2.3-70-g09d2