e_expf.S revision 1.4
1 1.1 jtc /* 2 1.1 jtc * Written by J.T. Conklin <jtc (at) netbsd.org>. 3 1.1 jtc * Public domain. 4 1.1 jtc */ 5 1.1 jtc 6 1.1 jtc #include <machine/asm.h> 7 1.1 jtc 8 1.4 fvdl #include "abi.h" 9 1.4 fvdl 10 1.4 fvdl 11 1.4 fvdl RCSID("$NetBSD: e_expf.S,v 1.4 2001/06/19 00:26:30 fvdl Exp $") 12 1.1 jtc 13 1.1 jtc /* e^x = 2^(x * log2(e)) */ 14 1.1 jtc ENTRY(__ieee754_expf) 15 1.4 fvdl XMM_ONE_ARG_FLOAT_PROLOGUE 16 1.4 fvdl flds ARG_FLOAT_ONE 17 1.1 jtc fldl2e 18 1.1 jtc fmulp /* x * log2(e) */ 19 1.4 fvdl fld %st(0) 20 1.1 jtc frndint /* int(x * log2(e)) */ 21 1.3 kleink fsubr %st(0),%st(1) /* fract(x * log2(e)) */ 22 1.3 kleink fxch 23 1.2 simonb f2xm1 /* 2^(fract(x * log2(e))) - 1 */ 24 1.1 jtc fld1 25 1.1 jtc faddp /* 2^(fract(x * log2(e))) */ 26 1.1 jtc fscale /* e^x */ 27 1.4 fvdl fstp %st(1) 28 1.4 fvdl XMM_FLOAT_EPILOGUE 29 1.1 jtc ret 30
Indexes created Tue Sep 30 17:09:57 GMT 2025