diff options
| author | Szabolcs Nagy <nsz@port70.net> | 2013-08-15 10:56:57 +0000 |
|---|---|---|
| committer | Szabolcs Nagy <nsz@port70.net> | 2013-08-15 10:56:57 +0000 |
| commit | 31c5fb80b9eae86f801be4f46025bc6532a554c5 (patch) | |
| tree | c5d0912699727ebd96bb5194334ee4dd33dc9c5b /src/math/i386/asin.s | |
| parent | 1b3973fb43fbef80dab1dfc9c788783e78ab5043 (diff) | |
| download | musl-31c5fb80b9eae86f801be4f46025bc6532a554c5.tar.gz | |
math: fix x86 asin, atan, exp, log1p to raise underflow
underflow is raised by an inexact subnormal float store,
since subnormal operations are slow, check the underflow
flag and skip the store if it's already raised
Diffstat (limited to 'src/math/i386/asin.s')
| -rw-r--r-- | src/math/i386/asin.s | 23 |
1 files changed, 22 insertions, 1 deletions
diff --git a/src/math/i386/asin.s b/src/math/i386/asin.s index 932c7542..a9f691bf 100644 --- a/src/math/i386/asin.s +++ b/src/math/i386/asin.s @@ -2,7 +2,18 @@ .type asinf,@function asinf: flds 4(%esp) - jmp 1f + mov 4(%esp),%eax + add %eax,%eax + cmp $0x01000000,%eax + jae 1f + # subnormal x, return x with underflow + fnstsw %ax + and $16,%ax + jnz 2f + fld %st(0) + fmul %st(1) + fstps 4(%esp) +2: ret .global asinl .type asinl,@function @@ -14,6 +25,16 @@ asinl: .type asin,@function asin: fldl 4(%esp) + mov 8(%esp),%eax + add %eax,%eax + cmp $0x00200000,%eax + jae 1f + # subnormal x, return x with underflow + fnstsw %ax + and $16,%ax + jnz 2f + fsts 4(%esp) +2: ret 1: fld %st(0) fld1 fsub %st(0),%st(1) |