diff options
| author | Alexander Monakov <amonakov@ispras.ru> | 2020-01-14 23:36:44 +0300 |
|---|---|---|
| committer | Rich Felker <dalias@aerifal.cx> | 2020-03-24 16:31:36 -0400 |
| commit | b173e4262f65800ea488b4494125284779a61547 (patch) | |
| tree | 46cca9389969b9c15bde136bcdd477c80a670c60 /src/math/i386/remainderl.c | |
| parent | 51f4f8c512d682fe0c1a7a891909e75f416f20f6 (diff) | |
| download | musl-b173e4262f65800ea488b4494125284779a61547.tar.gz | |
math: move x87-family remainder functions to C with inline asm
Diffstat (limited to 'src/math/i386/remainderl.c')
| -rw-r--r-- | src/math/i386/remainderl.c | 9 |
1 files changed, 9 insertions, 0 deletions
diff --git a/src/math/i386/remainderl.c b/src/math/i386/remainderl.c new file mode 100644 index 00000000..8cf75071 --- /dev/null +++ b/src/math/i386/remainderl.c @@ -0,0 +1,9 @@ +#include <math.h> + +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; +} |