URL
https://opencores.org/ocsvn/openrisc_2011-10-31/openrisc_2011-10-31/trunk
Subversion Repositories openrisc_2011-10-31
[/] [openrisc/] [trunk/] [gnu-src/] [newlib-1.17.0/] [newlib/] [libc/] [machine/] [i386/] [f_exp.c] - Rev 194
Go to most recent revision | Compare with Previous | Blame | View Log
/* * ==================================================== * Copyright (C) 1998,2002 by Red Hat Inc. All rights reserved. * * Permission to use, copy, modify, and distribute this * software is freely granted, provided that this notice * is preserved. * ==================================================== */ #if !defined(_SOFT_FLOAT) /* Fast version of exp using Intel float instructions. double _f_exp (double x); Function computes e ** x. The following special cases exist: 1. if x is 0.0 ==> return 1.0 2. if x is infinity ==> return infinity 3. if x is -infinity ==> return 0.0 4. if x is NaN ==> return x There is no error checking or setting of errno. */ #include <math.h> #include <ieeefp.h> #include "f_math.h" double _f_exp (double x) { if (check_finite(x)) { double result; asm ("fldl2e; fmulp; fld %%st; frndint; fsub %%st,%%st(1); fxch;" \ "fchs; f2xm1; fld1; faddp; fxch; fld1; fscale; fstp %%st(1); fmulp" : "=t"(result) : "0"(x)); return result; } else if (x == -infinity()) return 0.0; return x; } #endif
Go to most recent revision | Compare with Previous | Blame | View Log