newlib-cygwin/newlib/libm/machine/i386/f_pow.c

48 lines
1.4 KiB
C
Raw Normal View History

* Makefile.shared: New file. * libc/sys/linux/shared.ld: New file. * libm/machine/*: New files. * libm/machine/i386/*: New files. * Makefile.am: Add libtool support. Change math and mathfp references to variables. * configure.host: Add variables for libtool support. Add libm_machine_dir variable. * configure.in: Add objectlist variables, for libtool support. Add CC_FOR_BUILD tests. * libc/Makefile.am: Add libtool support. Change crt0.o reference to be a variable reference. * libc/configure.in: Add libtool support. Change sublib names to be lib${subdir}.la when using libtool. * libc/ctype/Makefile.am: Add libtool support. * libc/errno/Makefile.am: Likewise. * libc/locale/Makefile.am: Likewise. * libc/machine/Makefile.am: Likewise. * libc/machine/configure.in: Likewise. * libc/machine/i386/Makefile.am: Likewise. * libc/machine/i386/configure.in: Likewise. * libc/misc/Makefile.am: Likewise. * libc/posix/Makefile.am: Likewise. * libc/reent/Makefile.am: Likewise. * libc/signal/Makefile.am: Likewise. * libc/stdio/Makefile.am: Likewise. * libc/stdlib/Makefile.am: Likewise. * libc/string/Makefile.am: Likewise. * libc/sys/Makefile.am: Likewise. * libc/sys/configure.in: Likewise. * libc/sys/linux/Makefile.am: Add libtool support. Change awk reference to a variable reference. Change signal.h reference to a variable reference. * libc/sys/linux/configure.in: Add libtool support. * libc/syscalls/Makefile.am: Likewise. * libc/time/Makefile.am: Likewise. * libc/unix/Makefile.am: Likewise. * libm/Makefile.am: Add libtool support. Change math and mathfp references to variables. * libm/configure.in: Add libtool support. Add LIBM_MACHINE_LIB variable. * libm/common/Makefile.am: Add libtool support. * libm/math/Makefile.am: Likewise. * libm/mathfp/Makefile.am: Likewise. Regenerate all Makefile.in, aclocal.m4, and configure.
2001-12-14 07:50:11 +08:00
/*
* ====================================================
2002-12-20 Jeff Johnston <jjohnstn@redhat.com> * libc/machine/i386/f_atan2.S: Change copyright from Cygnus Solutions to Red Hat Inc. * libc/machine/i386/f_atan2f.S: Ditto. * libc/machine/i386/f_exp.c: Ditto. * libc/machine/i386/f_expf.c: Ditto. * libc/machine/i386/f_frexp.S: Ditto. * libc/machine/i386/f_frexpf.S: Ditto. * libc/machine/i386/f_ldexp.S: Ditto. * libc/machine/i386/f_ldexpf.S: Ditto. * libc/machine/i386/f_log.S: Ditto. * libc/machine/i386/f_log10.S: Ditto. * libc/machine/i386/f_log10f.S: Ditto. * libc/machine/i386/f_logf.S: Ditto. * libc/machine/i386/f_pow.c: Ditto. * libc/machine/i386/f_powf.c: Ditto. * libc/machine/i386/f_tan.S: Ditto. * libc/machine/i386/f_tanf.S: Ditto. * libc/machine/i386/memchr.S: Ditto. * libc/machine/i386/memcmp.S: Ditto. * libc/machine/i386/memcpy.S: Ditto. * libc/machine/i386/memmove.S: Ditto. * libc/machine/i386/memset.S: Ditto. * libc/machine/i386/strchr.S: Ditto. * libc/machine/i386/strlen.S: Ditto. * libm/machine/i386/f_atan2.S: Ditto. * libm/machine/i386/f_atan2f.S: Ditto. * libm/machine/i386/f_exp.c: Ditto. * libm/machine/i386/f_expf.c: Ditto. * libm/machine/i386/f_frexp.S: Ditto. * libm/machine/i386/f_frexpf.S: Ditto. * libm/machine/i386/f_ldexp.S: Ditto. * libm/machine/i386/f_ldexpf.S: Ditto. * libm/machine/i386/f_log.S: Ditto. * libm/machine/i386/f_log10.S: Ditto. * libm/machine/i386/f_log10f.S: Ditto. * libm/machine/i386/f_logf.S: Ditto. * libm/machine/i386/f_pow.c: Ditto. * libm/machine/i386/f_powf.c: Ditto. * libm/machine/i386/f_tan.S: Ditto. * libm/machine/i386/f_tanf.S: Ditto.
2002-12-21 05:31:20 +08:00
* Copyright (C) 1998, 2002 by Red Hat Inc. All rights reserved.
* Makefile.shared: New file. * libc/sys/linux/shared.ld: New file. * libm/machine/*: New files. * libm/machine/i386/*: New files. * Makefile.am: Add libtool support. Change math and mathfp references to variables. * configure.host: Add variables for libtool support. Add libm_machine_dir variable. * configure.in: Add objectlist variables, for libtool support. Add CC_FOR_BUILD tests. * libc/Makefile.am: Add libtool support. Change crt0.o reference to be a variable reference. * libc/configure.in: Add libtool support. Change sublib names to be lib${subdir}.la when using libtool. * libc/ctype/Makefile.am: Add libtool support. * libc/errno/Makefile.am: Likewise. * libc/locale/Makefile.am: Likewise. * libc/machine/Makefile.am: Likewise. * libc/machine/configure.in: Likewise. * libc/machine/i386/Makefile.am: Likewise. * libc/machine/i386/configure.in: Likewise. * libc/misc/Makefile.am: Likewise. * libc/posix/Makefile.am: Likewise. * libc/reent/Makefile.am: Likewise. * libc/signal/Makefile.am: Likewise. * libc/stdio/Makefile.am: Likewise. * libc/stdlib/Makefile.am: Likewise. * libc/string/Makefile.am: Likewise. * libc/sys/Makefile.am: Likewise. * libc/sys/configure.in: Likewise. * libc/sys/linux/Makefile.am: Add libtool support. Change awk reference to a variable reference. Change signal.h reference to a variable reference. * libc/sys/linux/configure.in: Add libtool support. * libc/syscalls/Makefile.am: Likewise. * libc/time/Makefile.am: Likewise. * libc/unix/Makefile.am: Likewise. * libm/Makefile.am: Add libtool support. Change math and mathfp references to variables. * libm/configure.in: Add libtool support. Add LIBM_MACHINE_LIB variable. * libm/common/Makefile.am: Add libtool support. * libm/math/Makefile.am: Likewise. * libm/mathfp/Makefile.am: Likewise. Regenerate all Makefile.in, aclocal.m4, and configure.
2001-12-14 07:50:11 +08:00
*
* Permission to use, copy, modify, and distribute this
* software is freely granted, provided that this notice
* is preserved.
* ====================================================
*/
#if !defined(_SOFT_FLOAT)
* Makefile.shared: New file. * libc/sys/linux/shared.ld: New file. * libm/machine/*: New files. * libm/machine/i386/*: New files. * Makefile.am: Add libtool support. Change math and mathfp references to variables. * configure.host: Add variables for libtool support. Add libm_machine_dir variable. * configure.in: Add objectlist variables, for libtool support. Add CC_FOR_BUILD tests. * libc/Makefile.am: Add libtool support. Change crt0.o reference to be a variable reference. * libc/configure.in: Add libtool support. Change sublib names to be lib${subdir}.la when using libtool. * libc/ctype/Makefile.am: Add libtool support. * libc/errno/Makefile.am: Likewise. * libc/locale/Makefile.am: Likewise. * libc/machine/Makefile.am: Likewise. * libc/machine/configure.in: Likewise. * libc/machine/i386/Makefile.am: Likewise. * libc/machine/i386/configure.in: Likewise. * libc/misc/Makefile.am: Likewise. * libc/posix/Makefile.am: Likewise. * libc/reent/Makefile.am: Likewise. * libc/signal/Makefile.am: Likewise. * libc/stdio/Makefile.am: Likewise. * libc/stdlib/Makefile.am: Likewise. * libc/string/Makefile.am: Likewise. * libc/sys/Makefile.am: Likewise. * libc/sys/configure.in: Likewise. * libc/sys/linux/Makefile.am: Add libtool support. Change awk reference to a variable reference. Change signal.h reference to a variable reference. * libc/sys/linux/configure.in: Add libtool support. * libc/syscalls/Makefile.am: Likewise. * libc/time/Makefile.am: Likewise. * libc/unix/Makefile.am: Likewise. * libm/Makefile.am: Add libtool support. Change math and mathfp references to variables. * libm/configure.in: Add libtool support. Add LIBM_MACHINE_LIB variable. * libm/common/Makefile.am: Add libtool support. * libm/math/Makefile.am: Likewise. * libm/mathfp/Makefile.am: Likewise. Regenerate all Makefile.in, aclocal.m4, and configure.
2001-12-14 07:50:11 +08:00
/*
Fast version of pow using Intel float instructions.
double _f_pow (double x, double y);
Function calculates x to power of y.
The function optimizes the case where x is >0.0 and y is finite.
In such a case, there is no error checking or setting of errno.
All other cases defer to normal pow() function which will
set errno as normal.
*/
#include <math.h>
#include <ieeefp.h>
#include "f_math.h"
double _f_pow (double x, double y)
{
/* following sequence handles the majority of cases for pow() */
if (x > 0.0 && check_finite(y))
{
double result;
/* calculate x ** y as 2 ** (y log2(x)). On Intel, can only
raise 2 to an integer or a small fraction, thus, we have
to perform two steps 2**integer portion * 2**fraction. */
asm ("fldl 8(%%ebp); fyl2x; fld %%st; frndint; fsub %%st,%%st(1);" \
"fxch; fchs; f2xm1; fld1; faddp; fxch; fld1; fscale; fstp %%st(1);"\
"fmulp" : "=t" (result) : "0" (y));
return result;
}
else /* all other strange cases, defer to normal pow() */
return pow (x,y);
}
#endif