From b173e4262f65800ea488b4494125284779a61547 Mon Sep 17 00:00:00 2001 From: Alexander Monakov Date: Tue, 14 Jan 2020 23:36:44 +0300 Subject: math: move x87-family remainder functions to C with inline asm --- src/math/x86_64/remainderl.c | 9 +++++++++ src/math/x86_64/remainderl.s | 11 ----------- 2 files changed, 9 insertions(+), 11 deletions(-) create mode 100644 src/math/x86_64/remainderl.c delete mode 100644 src/math/x86_64/remainderl.s (limited to 'src/math/x86_64') diff --git a/src/math/x86_64/remainderl.c b/src/math/x86_64/remainderl.c new file mode 100644 index 00000000..8cf75071 --- /dev/null +++ b/src/math/x86_64/remainderl.c @@ -0,0 +1,9 @@ +#include + +long double remainderl(long double x, long double y) +{ + unsigned short fpsr; + do __asm__ ("fprem1; fnstsw %%ax" : "+t"(x), "=a"(fpsr) : "u"(y)); + while (fpsr & 0x400); + return x; +} diff --git a/src/math/x86_64/remainderl.s b/src/math/x86_64/remainderl.s deleted file mode 100644 index cb3857b4..00000000 --- a/src/math/x86_64/remainderl.s +++ /dev/null @@ -1,11 +0,0 @@ -.global remainderl -.type remainderl,@function -remainderl: - fldt 24(%rsp) - fldt 8(%rsp) -1: fprem1 - fnstsw %ax - testb $4,%ah - jnz 1b - fstp %st(1) - ret -- cgit 1.4.1