diff options
| author | Vladimir Azarov <avm@intermediate-node.net> | 2024-10-01 15:47:05 +0200 |
|---|---|---|
| committer | Vladimir Azarov <avm@intermediate-node.net> | 2024-10-01 15:47:05 +0200 |
| commit | 4abab5ad6c8465a7528ccdd5f49367da05f78bbd (patch) | |
| tree | ebf009bf1376a5a223a915bc27cbbd791a1316bc /src/math/i386/log1p.s | |
Initial version
Diffstat (limited to 'src/math/i386/log1p.s')
| -rw-r--r-- | src/math/i386/log1p.s | 25 |
1 files changed, 25 insertions, 0 deletions
diff --git a/src/math/i386/log1p.s b/src/math/i386/log1p.s new file mode 100644 index 0000000..f3c95f8 --- /dev/null +++ b/src/math/i386/log1p.s @@ -0,0 +1,25 @@ +.global log1p +.type log1p,@function +log1p: + mov 8(%esp),%eax + fldln2 + and $0x7fffffff,%eax + fldl 4(%esp) + cmp $0x3fd28f00,%eax + ja 1f + 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: fsts 4(%esp) + fstp %st(1) + ret |