Wilco Dijkstra ea5c662c62 Improve performance of sincosf ... This patch is a complete rewrite of sincosf. The new version is significantly faster, as well as simple and accurate. The worst-case ULP is 0.5607, maximum relative error is 0.5303 * 2^-23 over all 4 billion inputs. In non-nearest rounding modes the error is 1ULP. The algorithm uses 3 main cases: small inputs which don't need argument reduction, small inputs which need a simple range reduction and large inputs requiring complex range reduction. The code uses approximate integer comparisons to quickly decide between these cases. The small range reducer uses a single reduction step to handle values up to 120.0. It is fastest on targets which support inlined round instructions. The large range reducer uses integer arithmetic for simplicity. It does a 32x96 bit multiply to compute a 64-bit modulo result. This is more than accurate enough to handle the worst-case cancellation for values close to an integer multiple of PI/4. It could be further optimized, however it is already much faster than necessary. sincosf throughput gains on Cortex-A72: * |x| < 0x1p-12 : 1.6x * |x| < M_PI_4 : 1.7x * |x| < 2 * M_PI: 1.5x * |x| < 120.0 : 1.8x * |x| < Inf : 2.3x * math/Makefile: Add s_sincosf_data.c. * sysdeps/ia64/fpu/s_sincosf_data.c: New file. * sysdeps/ieee754/flt-32/s_sincosf.h (abstop12): Add new function. (sincosf_poly): Likewise. (reduce_small): Likewise. (reduce_large): Likewise. * sysdeps/ieee754/flt-32/s_sincosf.c (sincosf): Rewrite. * sysdeps/ieee754/flt-32/s_sincosf_data.c: New file with sincosf data. * sysdeps/m68k/m680x0/fpu/s_sincosf_data.c: New file. * sysdeps/x86_64/fpu/s_sincosf_data.c: New file.		2018-08-10 17:34:39 +01:00
..
e_acosf.c	Rename all __ieee754_sqrt(f/l) calls to sqrt(f/l)	2018-03-15 19:21:36 +00:00
e_acoshf.c	Rename all __ieee754_sqrt(f/l) calls to sqrt(f/l)	2018-03-15 19:21:36 +00:00
e_asinf.c	Move math_check_force_underflow macros to separate math-underflow.h.	2018-05-10 00:53:04 +00:00
e_atan2f.c	Consistently use uintN_t not u_intN_t in libm.	2017-08-03 19:55:04 +00:00
e_atanhf.c	Do not include math-barriers.h in math_private.h.	2018-05-11 15:11:38 +00:00
e_coshf.c	Move math_narrow_eval to separate math-narrow-eval.h.	2018-05-09 00:15:10 +00:00
e_exp2f_data.c	Update copyright dates with scripts/update-copyrights.	2018-01-01 00:32:25 +00:00
e_exp2f.c	Move math_narrow_eval to separate math-narrow-eval.h.	2018-05-09 00:15:10 +00:00
e_expf.c	Clean up converttoint handling and document the semantics	2018-08-10 17:23:16 +01:00
e_fmodf.c	Consistently use uintN_t not u_intN_t in libm.	2017-08-03 19:55:04 +00:00
e_gammaf_r.c	Move math_check_force_underflow macros to separate math-underflow.h.	2018-05-10 00:53:04 +00:00
e_hypotf.c	Rename all __ieee754_sqrt(f/l) calls to sqrt(f/l)	2018-03-15 19:21:36 +00:00
e_ilogbf.c	Remove useless __ilogb*_finite aliases	2012-04-18 00:40:13 +02:00
e_j0f.c	Do not include math-barriers.h in math_private.h.	2018-05-11 15:11:38 +00:00
e_j1f.c	Move math_check_force_underflow macros to separate math-underflow.h.	2018-05-10 00:53:04 +00:00
e_jnf.c	Move math_check_force_underflow macros to separate math-underflow.h.	2018-05-10 00:53:04 +00:00
e_lgammaf_r.c	Move math_narrow_eval to separate math-narrow-eval.h.	2018-05-09 00:15:10 +00:00
e_log2f_data.c	Update copyright dates with scripts/update-copyrights.	2018-01-01 00:32:25 +00:00
e_log2f.c	Update copyright dates with scripts/update-copyrights.	2018-01-01 00:32:25 +00:00
e_log10f.c	Use fabs(f/l) rather than __fabs	2017-09-29 18:54:24 +01:00
e_logf_data.c	Update copyright dates with scripts/update-copyrights.	2018-01-01 00:32:25 +00:00
e_logf.c	Update copyright dates with scripts/update-copyrights.	2018-01-01 00:32:25 +00:00
e_powf_log2_data.c	Update copyright dates with scripts/update-copyrights.	2018-01-01 00:32:25 +00:00
e_powf.c	Use uint32_t sign in single precision math error handling functions	2018-07-02 09:29:04 +01:00
e_rem_pio2f.c	Consistently use uintN_t not u_intN_t in libm.	2017-08-03 19:55:04 +00:00
e_remainderf.c	Consistently use uintN_t not u_intN_t in libm.	2017-08-03 19:55:04 +00:00
e_sinhf.c	Move math_check_force_underflow macros to separate math-underflow.h.	2018-05-10 00:53:04 +00:00
e_sqrtf.c	Consistently use uintN_t not u_intN_t in libm.	2017-08-03 19:55:04 +00:00
k_cosf.c	Remove trailing whitespace.	2013-06-05 20:44:03 +00:00
k_rem_pio2f.c	Fix i686-linux-gnu build with GCC mainline.	2018-05-22 16:55:04 +00:00
k_sinf.c	Move math_check_force_underflow macros to separate math-underflow.h.	2018-05-10 00:53:04 +00:00
k_tanf.c	Move math_check_force_underflow macros to separate math-underflow.h.	2018-05-10 00:53:04 +00:00
lgamma_negf.c	Move math_narrow_eval to separate math-narrow-eval.h.	2018-05-09 00:15:10 +00:00
lgamma_productf.c	Fix lgamma (negative) inaccuracy (bug 2542, bug 2543, bug 2558).	2015-09-10 22:27:58 +00:00
math_config.h	Clean up converttoint handling and document the semantics	2018-08-10 17:23:16 +01:00
math_errf.c	Use uint32_t sign in single precision math error handling functions	2018-07-02 09:29:04 +01:00
mpn2flt.c	Update copyright dates with scripts/update-copyrights.	2018-01-01 00:32:25 +00:00
s_asinhf.c	Move math_check_force_underflow macros to separate math-underflow.h.	2018-05-10 00:53:04 +00:00
s_atanf.c	Move math_check_force_underflow macros to separate math-underflow.h.	2018-05-10 00:53:04 +00:00
s_cbrtf.c	Update copyright dates with scripts/update-copyrights.	2018-01-01 00:32:25 +00:00
s_ceilf.c	Use libm_alias_float in flt-32.	2017-09-22 20:24:12 +00:00
s_copysignf.c	Use libm_alias_float in flt-32.	2017-09-22 20:24:12 +00:00
s_cosf.c	Update copyright dates with scripts/update-copyrights.	2018-01-01 00:32:25 +00:00
s_erff.c	Move math_check_force_underflow macros to separate math-underflow.h.	2018-05-10 00:53:04 +00:00
s_expm1f.c	Do not include math-barriers.h in math_private.h.	2018-05-11 15:11:38 +00:00
s_fabsf.c	Use libm_alias_float in flt-32.	2017-09-22 20:24:12 +00:00
s_finitef.c	Consistently use uintN_t not u_intN_t in libm.	2017-08-03 19:55:04 +00:00
s_floorf.c	Use libm_alias_float in flt-32.	2017-09-22 20:24:12 +00:00
s_fpclassifyf.c	Update copyright dates with scripts/update-copyrights.	2018-01-01 00:32:25 +00:00
s_frexpf.c	Use libm_alias_float in flt-32.	2017-09-22 20:24:12 +00:00
s_fromfpf_main.c	Update copyright dates with scripts/update-copyrights.	2018-01-01 00:32:25 +00:00
s_fromfpf.c	Use libm_alias_float in flt-32.	2017-09-22 20:24:12 +00:00
s_fromfpxf.c	Use libm_alias_float in flt-32.	2017-09-22 20:24:12 +00:00
s_getpayloadf.c	Update copyright dates with scripts/update-copyrights.	2018-01-01 00:32:25 +00:00
s_isinff.c	Use <> for math.h and math_private.h everywhere.	2012-03-09 16:09:10 -08:00
s_isnanf.c	Consistently use uintN_t not u_intN_t in libm.	2017-08-03 19:55:04 +00:00
s_issignalingf.c	Update copyright dates with scripts/update-copyrights.	2018-01-01 00:32:25 +00:00
s_llrintf.c	Move math_narrow_eval to separate math-narrow-eval.h.	2018-05-09 00:15:10 +00:00
s_llroundf.c	Update copyright dates with scripts/update-copyrights.	2018-01-01 00:32:25 +00:00
s_log1pf.c	Do not include math-barriers.h in math_private.h.	2018-05-11 15:11:38 +00:00
s_logbf.c	Use libm_alias_float in flt-32.	2017-09-22 20:24:12 +00:00
s_lrintf.c	Move math_narrow_eval to separate math-narrow-eval.h.	2018-05-09 00:15:10 +00:00
s_lroundf.c	Update copyright dates with scripts/update-copyrights.	2018-01-01 00:32:25 +00:00
s_modff.c	Use libm_alias_float in flt-32.	2017-09-22 20:24:12 +00:00
s_nearbyintf.c	Do not include math-barriers.h in math_private.h.	2018-05-11 15:11:38 +00:00
s_nextafterf.c	Do not include math-barriers.h in math_private.h.	2018-05-11 15:11:38 +00:00
s_nextupf.c	Move math_check_force_underflow macros to separate math-underflow.h.	2018-05-10 00:53:04 +00:00
s_remquof.c	Update copyright dates with scripts/update-copyrights.	2018-01-01 00:32:25 +00:00
s_rintf.c	Use libm_alias_float in flt-32.	2017-09-22 20:24:12 +00:00
s_roundevenf.c	Update copyright dates with scripts/update-copyrights.	2018-01-01 00:32:25 +00:00
s_roundf.c	Update copyright dates with scripts/update-copyrights.	2018-01-01 00:32:25 +00:00
s_scalblnf.c	[BZ #6803] Set errno for scalbln, scalbn	2014-06-20 07:48:20 +05:30
s_scalbnf.c	Make scalbn set errno (bug 6803).	2015-09-16 21:11:00 +00:00
s_setpayloadf_main.c	Update copyright dates with scripts/update-copyrights.	2018-01-01 00:32:25 +00:00
s_setpayloadf.c	Use libm_alias_float in flt-32.	2017-09-22 20:24:12 +00:00
s_setpayloadsigf.c	Use libm_alias_float in flt-32.	2017-09-22 20:24:12 +00:00
s_signbitf.c	Update copyright dates with scripts/update-copyrights.	2018-01-01 00:32:25 +00:00
s_sincosf_data.c	Improve performance of sincosf	2018-08-10 17:34:39 +01:00
s_sincosf.c	Improve performance of sincosf	2018-08-10 17:34:39 +01:00
s_sincosf.h	Improve performance of sincosf	2018-08-10 17:34:39 +01:00
s_sinf.c	Update copyright dates with scripts/update-copyrights.	2018-01-01 00:32:25 +00:00
s_tanf.c	Use libm_alias_float in flt-32.	2017-09-22 20:24:12 +00:00
s_tanhf.c	Move math_check_force_underflow macros to separate math-underflow.h.	2018-05-10 00:53:04 +00:00
s_totalorderf.c	Update copyright dates with scripts/update-copyrights.	2018-01-01 00:32:25 +00:00
s_totalordermagf.c	Update copyright dates with scripts/update-copyrights.	2018-01-01 00:32:25 +00:00
s_truncf.c	Update copyright dates with scripts/update-copyrights.	2018-01-01 00:32:25 +00:00
s_ufromfpf.c	Use libm_alias_float in flt-32.	2017-09-22 20:24:12 +00:00
s_ufromfpxf.c	Use libm_alias_float in flt-32.	2017-09-22 20:24:12 +00:00
w_exp2f.c	Do not wrap expf and exp2f	2017-10-02 14:38:54 +01:00
w_expf.c	Do not wrap expf and exp2f	2017-10-02 14:38:54 +01:00
w_log2f.c	Do not wrap logf, log2f and powf	2017-10-02 14:39:38 +01:00
w_logf.c	Do not wrap logf, log2f and powf	2017-10-02 14:39:38 +01:00
w_powf.c	Do not wrap logf, log2f and powf	2017-10-02 14:39:38 +01:00