From 6f26cf3dacf9d426bbd18cecb4d486da56c17351 Mon Sep 17 00:00:00 2001 From: Rich Felker Date: Sun, 18 Mar 2012 21:43:43 -0400 Subject: asm exponential functions for i386 --- src/math/i386/expm1.s | 36 ++++++++++++++++++++++++++++++++++++ 1 file changed, 36 insertions(+) create mode 100644 src/math/i386/expm1.s (limited to 'src/math/i386/expm1.s') diff --git a/src/math/i386/expm1.s b/src/math/i386/expm1.s new file mode 100644 index 00000000..d6d511ef --- /dev/null +++ b/src/math/i386/expm1.s @@ -0,0 +1,36 @@ +.global expm1f +.type expm1f,@function +expm1f: + flds 4(%esp) + jmp 1f + +.global expm1l +.type expm1l,@function +expm1l: + fldt 4(%esp) + jmp 1f + +.global expm1 +.type expm1,@function +expm1: + fldl 4(%esp) +1: fldl2e + fmulp + fld %st(0) + frndint + fldz + fcompi + jnz 1f + fstp %st(0) + f2xm1 + ret +1: fxch %st(1) + fsub %st(1) + f2xm1 + fld1 + faddp + fscale + fld1 + fsubrp + fstp %st(1) + ret -- cgit 1.4.1