From bc87299ce72a52f4debf9fc19d859abe34dbdf43 Mon Sep 17 00:00:00 2001 From: Alexander Monakov Date: Wed, 15 Jan 2020 18:42:46 +0300 Subject: math: move x87-family fmod functions to C with inline asm --- src/math/i386/fmod.c | 10 ++++++++++ 1 file changed, 10 insertions(+) create mode 100644 src/math/i386/fmod.c (limited to 'src/math/i386/fmod.c') diff --git a/src/math/i386/fmod.c b/src/math/i386/fmod.c new file mode 100644 index 00000000..ea0c58d9 --- /dev/null +++ b/src/math/i386/fmod.c @@ -0,0 +1,10 @@ +#include + +double fmod(double x, double y) +{ + unsigned short fpsr; + // fprem does not introduce excess precision into x + do __asm__ ("fprem; fnstsw %%ax" : "+t"(x), "=a"(fpsr) : "u"(y)); + while (fpsr & 0x400); + return x; +} -- cgit v1.2.1