diff options
author | bde <bde@FreeBSD.org> | 2008-02-18 14:02:12 +0000 |
---|---|---|
committer | bde <bde@FreeBSD.org> | 2008-02-18 14:02:12 +0000 |
commit | ad78d666213a592bddae82c25e993dd47786dcef (patch) | |
tree | 390f2485d24a844d41ef64862adc9c01ce71db04 /lib/msun/src/s_sin.c | |
parent | 6c809a426ade64ce82f2f4480682bc1c520f6a2d (diff) | |
download | FreeBSD-src-ad78d666213a592bddae82c25e993dd47786dcef.zip FreeBSD-src-ad78d666213a592bddae82c25e993dd47786dcef.tar.gz |
Inline __ieee754__rem_pio2(). With gcc4-2, this gives an average
optimization of about 10% for cos(x), sin(x) and tan(x) on
|x| < 2**19*pi/2. We didn't do this before because __ieee754__rem_pio2()
is too large and complicated for gcc-3.3 to inline very well. We don't
do this for float precision because it interferes with optimization
of the usual (?) case (|x| < 9pi/4) which is manually inlined for float
precision only.
This has some rough edges:
- some static data is duplicated unnecessarily. There isn't much after
the recent move of large tables to k_rem_pio2.c, and some static data
is duplicated to good affect (all the data static const, so that the
compiler can evaluate expressions like 2*pio2 at compile time and
generate even more static data for the constant for this).
- extern inline is used (for the same reason as in previous inlining of
k_cosf.c etc.), but C99 apparently doesn't allow extern inline
functions with static data, and gcc will eventually warn about this.
Convert to __FBSDID().
Indent __ieee754_rem_pio2()'s declaration consistently (its style was
made inconsistent with fdlibm a while ago, so complete this).
Fix __ieee754_rem_pio2()'s return type to match its prototype. Someone
changed too many ints to int32_t's when fixing the assumption that all
ints are int32_t's.
Diffstat (limited to 'lib/msun/src/s_sin.c')
-rw-r--r-- | lib/msun/src/s_sin.c | 7 |
1 files changed, 4 insertions, 3 deletions
diff --git a/lib/msun/src/s_sin.c b/lib/msun/src/s_sin.c index 7cd3877..f6ed2b5 100644 --- a/lib/msun/src/s_sin.c +++ b/lib/msun/src/s_sin.c @@ -10,9 +10,8 @@ * ==================================================== */ -#ifndef lint -static char rcsid[] = "$FreeBSD$"; -#endif +#include <sys/cdefs.h> +__FBSDID("$FreeBSD$"); /* sin(x) * Return sine function of x. @@ -48,7 +47,9 @@ static char rcsid[] = "$FreeBSD$"; #include <float.h> #include "math.h" +#define INLINE_REM_PIO2 #include "math_private.h" +#include "e_rem_pio2.c" double sin(double x) |