math: move x87-family rint functions to C with inline asm

[musl] / src / math / i386 / log1p.s

diff --git a/src/math/i386/log1p.s b/src/math/i386/log1p.s
index 6b6929c..f3c95f8 100644 (file)
--- a/src/math/i386/log1p.s
+++ b/src/math/i386/log1p.s
@@ -10,15 +10,16 @@ log1p:
        cmp $0x00100000,%eax
        jb 2f
        fyl2xp1
+       fstpl 4(%esp)
+       fldl 4(%esp)
        ret
 1:     fld1
        faddp
        fyl2x
+       fstpl 4(%esp)
+       fldl 4(%esp)
        ret
                # subnormal x, return x with underflow
-2:     fnstsw %ax
-       and $16,%ax
-       jnz 1f
-       fsts 4(%esp)
+2:     fsts 4(%esp)
        fstp %st(1)
-1:     ret
+       ret