lib/glibc
1
0
Fork 0
mirror of https://sourceware.org/git/glibc.git synced 2025-07-28 00:21:52 +03:00
Code Activity

Refactor code forcing underflow exceptions.

Browse Source 
Various floating-point functions have code to force underflow
exceptions if a tiny result was computed in a way that might not have
resulted in such exceptions even though the result is inexact.  This
typically uses math_force_eval to ensure that the underflowing
expression is evaluated, but sometimes uses volatile.

This patch refactors such code to use three new macros
math_check_force_underflow, math_check_force_underflow_nonneg and
math_check_force_underflow_complex (which in turn use
math_force_eval).  In the limited number of cases not suited to a
simple conversion to these macros, existing uses of volatile are
changed to use math_force_eval instead.  The converted code does not
always execute exactly the same sequence of operations as the original
code, but the overall effects should be the same.

Tested for x86_64, x86, mips64 and powerpc.

	* sysdeps/generic/math_private.h (fabs_tg): New macro.
	(min_of_type): Likewise.
	(math_check_force_underflow): Likewise.
	(math_check_force_underflow_nonneg): Likewise.
	(math_check_force_underflow_complex): Likewise.
	* math/e_exp2l.c (__ieee754_exp2l): Use
	math_check_force_underflow_nonneg.
	* math/k_casinh.c (__kernel_casinh): Likewise.
	* math/k_casinhf.c (__kernel_casinhf): Likewise.
	* math/k_casinhl.c (__kernel_casinhl): Likewise.
	* math/s_catan.c (__catan): Use
	math_check_force_underflow_complex.
	* math/s_catanf.c (__catanf): Likewise.
	* math/s_catanh.c (__catanh): Likewise.
	* math/s_catanhf.c (__catanhf): Likewise.
	* math/s_catanhl.c (__catanhl): Likewise.
	* math/s_catanl.c (__catanl): Likewise.
	* math/s_ccosh.c (__ccosh): Likewise.
	* math/s_ccoshf.c (__ccoshf): Likewise.
	* math/s_ccoshl.c (__ccoshl): Likewise.
	* math/s_cexp.c (__cexp): Likewise.
	* math/s_cexpf.c (__cexpf): Likewise.
	* math/s_cexpl.c (__cexpl): Likewise.
	* math/s_clog.c (__clog): Use math_check_force_underflow_nonneg.
	* math/s_clog10.c (__clog10): Likewise.
	* math/s_clog10f.c (__clog10f): Likewise.
	* math/s_clog10l.c (__clog10l): Likewise.
	* math/s_clogf.c (__clogf): Likewise.
	* math/s_clogl.c (__clogl): Likewise.
	* math/s_csin.c (__csin): Use math_check_force_underflow_complex.
	* math/s_csinf.c (__csinf): Likewise.
	* math/s_csinh.c (__csinh): Likewise.
	* math/s_csinhf.c (__csinhf): Likewise.
	* math/s_csinhl.c (__csinhl): Likewise.
	* math/s_csinl.c (__csinl): Likewise.
	* math/s_csqrt.c (__csqrt): Use math_check_force_underflow.
	* math/s_csqrtf.c (__csqrtf): Likewise.
	* math/s_csqrtl.c (__csqrtl): Likewise.
	* math/s_ctan.c (__ctan): Use math_check_force_underflow_complex.
	* math/s_ctanf.c (__ctanf): Likewise.
	* math/s_ctanh.c (__ctanh): Likewise.
	* math/s_ctanhf.c (__ctanhf): Likewise.
	* math/s_ctanhl.c (__ctanhl): Likewise.
	* math/s_ctanl.c (__ctanl): Likewise.
	* stdlib/strtod_l.c (round_and_return): Use math_force_eval
	instead of volatile.
	* sysdeps/ieee754/dbl-64/e_asin.c (__ieee754_asin): Use
	math_check_force_underflow.
	* sysdeps/ieee754/dbl-64/e_atanh.c (__ieee754_atanh): Likewise.
	* sysdeps/ieee754/dbl-64/e_exp.c (__ieee754_exp): Do not use
	volatile when forcing underflow.
	* sysdeps/ieee754/dbl-64/e_exp2.c (__ieee754_exp2): Use
	math_check_force_underflow_nonneg.
	* sysdeps/ieee754/dbl-64/e_gamma_r.c (__ieee754_gamma_r):
	Likewise.
	* sysdeps/ieee754/dbl-64/e_j1.c (__ieee754_j1): Use
	math_check_force_underflow.
	* sysdeps/ieee754/dbl-64/e_jn.c (__ieee754_jn): Likewise.
	* sysdeps/ieee754/dbl-64/e_sinh.c (__ieee754_sinh): Likewise.
	* sysdeps/ieee754/dbl-64/s_asinh.c (__asinh): Likewise.
	* sysdeps/ieee754/dbl-64/s_atan.c (atan): Use
	math_check_force_underflow_nonneg.
	* sysdeps/ieee754/dbl-64/s_erf.c (__erf): Use
	math_check_force_underflow.
	* sysdeps/ieee754/dbl-64/s_expm1.c (__expm1): Likewise.
	* sysdeps/ieee754/dbl-64/s_fma.c (__fma): Use math_force_eval
	instead of volatile.
	* sysdeps/ieee754/dbl-64/s_log1p.c (__log1p): Use
	math_check_force_underflow.
	* sysdeps/ieee754/dbl-64/s_sin.c (__sin): Likewise.
	* sysdeps/ieee754/dbl-64/s_tan.c (tan): Use
	math_check_force_underflow_nonneg.
	* sysdeps/ieee754/dbl-64/s_tanh.c (__tanh): Use
	math_check_force_underflow.
	* sysdeps/ieee754/flt-32/e_asinf.c (__ieee754_asinf): Likewise.
	* sysdeps/ieee754/flt-32/e_atanhf.c (__ieee754_atanhf): Likewise.
	* sysdeps/ieee754/flt-32/e_exp2f.c (__ieee754_exp2f): Use
	math_check_force_underflow_nonneg.
	* sysdeps/ieee754/flt-32/e_gammaf_r.c (__ieee754_gammaf_r):
	Likewise.
	* sysdeps/ieee754/flt-32/e_j1f.c (__ieee754_j1f): Use
	math_check_force_underflow.
	* sysdeps/ieee754/flt-32/e_jnf.c (__ieee754_jnf): Likewise.
	* sysdeps/ieee754/flt-32/e_sinhf.c (__ieee754_sinhf): Likewise.
	* sysdeps/ieee754/flt-32/k_sinf.c (__kernel_sinf): Likewise.
	* sysdeps/ieee754/flt-32/k_tanf.c (__kernel_tanf): Likewise.
	* sysdeps/ieee754/flt-32/s_asinhf.c (__asinhf): Likewise.
	* sysdeps/ieee754/flt-32/s_atanf.c (__atanf): Likewise.
	* sysdeps/ieee754/flt-32/s_erff.c (__erff): Likewise.
	* sysdeps/ieee754/flt-32/s_expm1f.c (__expm1f): Likewise.
	* sysdeps/ieee754/flt-32/s_log1pf.c (__log1pf): Likewise.
	* sysdeps/ieee754/flt-32/s_tanhf.c (__tanhf): Likewise.
	* sysdeps/ieee754/ldbl-128/e_asinl.c (__ieee754_asinl): Likewise.
	* sysdeps/ieee754/ldbl-128/e_atanhl.c (__ieee754_atanhl):
	Likewise.
	* sysdeps/ieee754/ldbl-128/e_expl.c (__ieee754_expl): Use
	math_check_force_underflow_nonneg.
	* sysdeps/ieee754/ldbl-128/e_gammal_r.c (__ieee754_gammal_r):
	Likewise.
	* sysdeps/ieee754/ldbl-128/e_j1l.c (__ieee754_j1l): Use
	math_check_force_underflow.
	* sysdeps/ieee754/ldbl-128/e_jnl.c (__ieee754_jnl): Likewise.
	* sysdeps/ieee754/ldbl-128/e_sinhl.c (__ieee754_sinhl): Likewise.
	* sysdeps/ieee754/ldbl-128/k_sincosl.c (__kernel_sincosl):
	Likewise.
	* sysdeps/ieee754/ldbl-128/k_sinl.c (__kernel_sinl): Likewise.
	* sysdeps/ieee754/ldbl-128/k_tanl.c (__kernel_tanl): Likewise.
	* sysdeps/ieee754/ldbl-128/s_asinhl.c (__asinhl): Likewise.
	* sysdeps/ieee754/ldbl-128/s_atanl.c (__atanl): Likewise.
	* sysdeps/ieee754/ldbl-128/s_erfl.c (__erfl): Likewise.
	* sysdeps/ieee754/ldbl-128/s_expm1l.c (__expm1l): Likewise.
	* sysdeps/ieee754/ldbl-128/s_fmal.c (__fmal): Use math_force_eval
	instead of volatile.
	* sysdeps/ieee754/ldbl-128/s_log1pl.c (__log1pl): Use
	math_check_force_underflow.
	* sysdeps/ieee754/ldbl-128/s_tanhl.c (__tanhl): Likewise.
	* sysdeps/ieee754/ldbl-128ibm/e_asinl.c (__ieee754_asinl): Use
	math_check_force_underflow.
	* sysdeps/ieee754/ldbl-128ibm/e_atanhl.c (__ieee754_atanhl):
	Likewise.
	* sysdeps/ieee754/ldbl-128ibm/e_gammal_r.c (__ieee754_gammal_r):
	Use math_check_force_underflow_nonneg.
	* sysdeps/ieee754/ldbl-128ibm/e_jnl.c (__ieee754_jnl): Use
	math_check_force_underflow.
	* sysdeps/ieee754/ldbl-128ibm/e_sinhl.c (__ieee754_sinhl):
	Likewise.
	* sysdeps/ieee754/ldbl-128ibm/k_sincosl.c (__kernel_sincosl):
	Likewise.
	* sysdeps/ieee754/ldbl-128ibm/k_sinl.c (__kernel_sinl): Likewise.
	* sysdeps/ieee754/ldbl-128ibm/k_tanl.c (__kernel_tanl): Likewise.
	* sysdeps/ieee754/ldbl-128ibm/s_asinhl.c (__asinhl): Likewise.
	* sysdeps/ieee754/ldbl-128ibm/s_atanl.c (__atanl): Likewise.
	* sysdeps/ieee754/ldbl-128ibm/s_erfl.c (__erfl): Likewise.
	* sysdeps/ieee754/ldbl-128ibm/s_tanhl.c (__tanhl): Likewise.
	* sysdeps/ieee754/ldbl-96/e_asinl.c (__ieee754_asinl): Likewise.
	* sysdeps/ieee754/ldbl-96/e_atanhl.c (__ieee754_atanhl): Likewise.
	* sysdeps/ieee754/ldbl-96/e_gammal_r.c (__ieee754_gammal_r): Use
	math_check_force_underflow_nonneg.
	* sysdeps/ieee754/ldbl-96/e_j1l.c (__ieee754_j1l): Use
	math_check_force_underflow.
	* sysdeps/ieee754/ldbl-96/e_jnl.c (__ieee754_jnl): Likewise.
	* sysdeps/ieee754/ldbl-96/e_sinhl.c (__ieee754_sinhl): Likewise.
	* sysdeps/ieee754/ldbl-96/k_sinl.c (__kernel_sinl): Likewise.
	* sysdeps/ieee754/ldbl-96/k_tanl.c (__kernel_tanl): Use
	math_check_force_underflow_nonneg.
	* sysdeps/ieee754/ldbl-96/s_asinhl.c (__asinhl): Use
	math_check_force_underflow.
	* sysdeps/ieee754/ldbl-96/s_erfl.c (__erfl): Likewise.
	* sysdeps/ieee754/ldbl-96/s_fmal.c (__fmal): Use math_force_eval
	instead of volatile.
	* sysdeps/ieee754/ldbl-96/s_tanhl.c (__tanhl): Use
	math_check_force_underflow.
This commit is contained in:
Joseph Myers
2015-09-23 22:42:30 +00:00
parent 54142c44e9
commit d96164c330
113 changed files with 334 additions and 723 deletions
Download Patch File Download Diff File Expand all files Collapse all files

153
ChangeLog
View File

@ -1,5 +1,158 @@
2015-09-23 Joseph Myers <joseph@codesourcery.com> 2015-09-23 Joseph Myers <joseph@codesourcery.com>
* sysdeps/generic/math_private.h (fabs_tg): New macro.
(min_of_type): Likewise.
(math_check_force_underflow): Likewise.
(math_check_force_underflow_nonneg): Likewise.
(math_check_force_underflow_complex): Likewise.
* math/e_exp2l.c (__ieee754_exp2l): Use
math_check_force_underflow_nonneg.
* math/k_casinh.c (__kernel_casinh): Likewise.
* math/k_casinhf.c (__kernel_casinhf): Likewise.
* math/k_casinhl.c (__kernel_casinhl): Likewise.
* math/s_catan.c (__catan): Use
math_check_force_underflow_complex.
* math/s_catanf.c (__catanf): Likewise.
* math/s_catanh.c (__catanh): Likewise.
* math/s_catanhf.c (__catanhf): Likewise.
* math/s_catanhl.c (__catanhl): Likewise.
* math/s_catanl.c (__catanl): Likewise.
* math/s_ccosh.c (__ccosh): Likewise.
* math/s_ccoshf.c (__ccoshf): Likewise.
* math/s_ccoshl.c (__ccoshl): Likewise.
* math/s_cexp.c (__cexp): Likewise.
* math/s_cexpf.c (__cexpf): Likewise.
* math/s_cexpl.c (__cexpl): Likewise.
* math/s_clog.c (__clog): Use math_check_force_underflow_nonneg.
* math/s_clog10.c (__clog10): Likewise.
* math/s_clog10f.c (__clog10f): Likewise.
* math/s_clog10l.c (__clog10l): Likewise.
* math/s_clogf.c (__clogf): Likewise.
* math/s_clogl.c (__clogl): Likewise.
* math/s_csin.c (__csin): Use math_check_force_underflow_complex.
* math/s_csinf.c (__csinf): Likewise.
* math/s_csinh.c (__csinh): Likewise.
* math/s_csinhf.c (__csinhf): Likewise.
* math/s_csinhl.c (__csinhl): Likewise.
* math/s_csinl.c (__csinl): Likewise.
* math/s_csqrt.c (__csqrt): Use math_check_force_underflow.
* math/s_csqrtf.c (__csqrtf): Likewise.
* math/s_csqrtl.c (__csqrtl): Likewise.
* math/s_ctan.c (__ctan): Use math_check_force_underflow_complex.
* math/s_ctanf.c (__ctanf): Likewise.
* math/s_ctanh.c (__ctanh): Likewise.
* math/s_ctanhf.c (__ctanhf): Likewise.
* math/s_ctanhl.c (__ctanhl): Likewise.
* math/s_ctanl.c (__ctanl): Likewise.
* stdlib/strtod_l.c (round_and_return): Use math_force_eval
instead of volatile.
* sysdeps/ieee754/dbl-64/e_asin.c (__ieee754_asin): Use
math_check_force_underflow.
* sysdeps/ieee754/dbl-64/e_atanh.c (__ieee754_atanh): Likewise.
* sysdeps/ieee754/dbl-64/e_exp.c (__ieee754_exp): Do not use
volatile when forcing underflow.
* sysdeps/ieee754/dbl-64/e_exp2.c (__ieee754_exp2): Use
math_check_force_underflow_nonneg.
* sysdeps/ieee754/dbl-64/e_gamma_r.c (__ieee754_gamma_r):
Likewise.
* sysdeps/ieee754/dbl-64/e_j1.c (__ieee754_j1): Use
math_check_force_underflow.
* sysdeps/ieee754/dbl-64/e_jn.c (__ieee754_jn): Likewise.
* sysdeps/ieee754/dbl-64/e_sinh.c (__ieee754_sinh): Likewise.
* sysdeps/ieee754/dbl-64/s_asinh.c (__asinh): Likewise.
* sysdeps/ieee754/dbl-64/s_atan.c (atan): Use
math_check_force_underflow_nonneg.
* sysdeps/ieee754/dbl-64/s_erf.c (__erf): Use
math_check_force_underflow.
* sysdeps/ieee754/dbl-64/s_expm1.c (__expm1): Likewise.
* sysdeps/ieee754/dbl-64/s_fma.c (__fma): Use math_force_eval
instead of volatile.
* sysdeps/ieee754/dbl-64/s_log1p.c (__log1p): Use
math_check_force_underflow.
* sysdeps/ieee754/dbl-64/s_sin.c (__sin): Likewise.
* sysdeps/ieee754/dbl-64/s_tan.c (tan): Use
math_check_force_underflow_nonneg.
* sysdeps/ieee754/dbl-64/s_tanh.c (__tanh): Use
math_check_force_underflow.
* sysdeps/ieee754/flt-32/e_asinf.c (__ieee754_asinf): Likewise.
* sysdeps/ieee754/flt-32/e_atanhf.c (__ieee754_atanhf): Likewise.
* sysdeps/ieee754/flt-32/e_exp2f.c (__ieee754_exp2f): Use
math_check_force_underflow_nonneg.
* sysdeps/ieee754/flt-32/e_gammaf_r.c (__ieee754_gammaf_r):
Likewise.
* sysdeps/ieee754/flt-32/e_j1f.c (__ieee754_j1f): Use
math_check_force_underflow.
* sysdeps/ieee754/flt-32/e_jnf.c (__ieee754_jnf): Likewise.
* sysdeps/ieee754/flt-32/e_sinhf.c (__ieee754_sinhf): Likewise.
* sysdeps/ieee754/flt-32/k_sinf.c (__kernel_sinf): Likewise.
* sysdeps/ieee754/flt-32/k_tanf.c (__kernel_tanf): Likewise.
* sysdeps/ieee754/flt-32/s_asinhf.c (__asinhf): Likewise.
* sysdeps/ieee754/flt-32/s_atanf.c (__atanf): Likewise.
* sysdeps/ieee754/flt-32/s_erff.c (__erff): Likewise.
* sysdeps/ieee754/flt-32/s_expm1f.c (__expm1f): Likewise.
* sysdeps/ieee754/flt-32/s_log1pf.c (__log1pf): Likewise.
* sysdeps/ieee754/flt-32/s_tanhf.c (__tanhf): Likewise.
* sysdeps/ieee754/ldbl-128/e_asinl.c (__ieee754_asinl): Likewise.
* sysdeps/ieee754/ldbl-128/e_atanhl.c (__ieee754_atanhl):
Likewise.
* sysdeps/ieee754/ldbl-128/e_expl.c (__ieee754_expl): Use
math_check_force_underflow_nonneg.
* sysdeps/ieee754/ldbl-128/e_gammal_r.c (__ieee754_gammal_r):
Likewise.
* sysdeps/ieee754/ldbl-128/e_j1l.c (__ieee754_j1l): Use
math_check_force_underflow.
* sysdeps/ieee754/ldbl-128/e_jnl.c (__ieee754_jnl): Likewise.
* sysdeps/ieee754/ldbl-128/e_sinhl.c (__ieee754_sinhl): Likewise.
* sysdeps/ieee754/ldbl-128/k_sincosl.c (__kernel_sincosl):
Likewise.
* sysdeps/ieee754/ldbl-128/k_sinl.c (__kernel_sinl): Likewise.
* sysdeps/ieee754/ldbl-128/k_tanl.c (__kernel_tanl): Likewise.
* sysdeps/ieee754/ldbl-128/s_asinhl.c (__asinhl): Likewise.
* sysdeps/ieee754/ldbl-128/s_atanl.c (__atanl): Likewise.
* sysdeps/ieee754/ldbl-128/s_erfl.c (__erfl): Likewise.
* sysdeps/ieee754/ldbl-128/s_expm1l.c (__expm1l): Likewise.
* sysdeps/ieee754/ldbl-128/s_fmal.c (__fmal): Use math_force_eval
instead of volatile.
* sysdeps/ieee754/ldbl-128/s_log1pl.c (__log1pl): Use
math_check_force_underflow.
* sysdeps/ieee754/ldbl-128/s_tanhl.c (__tanhl): Likewise.
* sysdeps/ieee754/ldbl-128ibm/e_asinl.c (__ieee754_asinl): Use
math_check_force_underflow.
* sysdeps/ieee754/ldbl-128ibm/e_atanhl.c (__ieee754_atanhl):
Likewise.
* sysdeps/ieee754/ldbl-128ibm/e_gammal_r.c (__ieee754_gammal_r):
Use math_check_force_underflow_nonneg.
* sysdeps/ieee754/ldbl-128ibm/e_jnl.c (__ieee754_jnl): Use
math_check_force_underflow.
* sysdeps/ieee754/ldbl-128ibm/e_sinhl.c (__ieee754_sinhl):
Likewise.
* sysdeps/ieee754/ldbl-128ibm/k_sincosl.c (__kernel_sincosl):
Likewise.
* sysdeps/ieee754/ldbl-128ibm/k_sinl.c (__kernel_sinl): Likewise.
* sysdeps/ieee754/ldbl-128ibm/k_tanl.c (__kernel_tanl): Likewise.
* sysdeps/ieee754/ldbl-128ibm/s_asinhl.c (__asinhl): Likewise.
* sysdeps/ieee754/ldbl-128ibm/s_atanl.c (__atanl): Likewise.
* sysdeps/ieee754/ldbl-128ibm/s_erfl.c (__erfl): Likewise.
* sysdeps/ieee754/ldbl-128ibm/s_tanhl.c (__tanhl): Likewise.
* sysdeps/ieee754/ldbl-96/e_asinl.c (__ieee754_asinl): Likewise.
* sysdeps/ieee754/ldbl-96/e_atanhl.c (__ieee754_atanhl): Likewise.
* sysdeps/ieee754/ldbl-96/e_gammal_r.c (__ieee754_gammal_r): Use
math_check_force_underflow_nonneg.
* sysdeps/ieee754/ldbl-96/e_j1l.c (__ieee754_j1l): Use
math_check_force_underflow.
* sysdeps/ieee754/ldbl-96/e_jnl.c (__ieee754_jnl): Likewise.
* sysdeps/ieee754/ldbl-96/e_sinhl.c (__ieee754_sinhl): Likewise.
* sysdeps/ieee754/ldbl-96/k_sinl.c (__kernel_sinl): Likewise.
* sysdeps/ieee754/ldbl-96/k_tanl.c (__kernel_tanl): Use
math_check_force_underflow_nonneg.
* sysdeps/ieee754/ldbl-96/s_asinhl.c (__asinhl): Use
math_check_force_underflow.
* sysdeps/ieee754/ldbl-96/s_erfl.c (__erfl): Likewise.
* sysdeps/ieee754/ldbl-96/s_fmal.c (__fmal): Use math_force_eval
instead of volatile.
* sysdeps/ieee754/ldbl-96/s_tanhl.c (__tanhl): Use
math_check_force_underflow.
* math/s_nexttowardf.c (__nexttowardf): Use math_narrow_eval. * math/s_nexttowardf.c (__nexttowardf): Use math_narrow_eval.
* stdlib/strtod_l.c: Include <math_private.h>. * stdlib/strtod_l.c: Include <math_private.h>.
(overflow_value): Use math_narrow_eval. (overflow_value): Use math_narrow_eval.

6
math/e_exp2l.c
View File

@ -43,11 +43,7 @@ __ieee754_exp2l (long double x)
result = __scalbnl (1.0L + fractx, intx); result = __scalbnl (1.0L + fractx, intx);
else else
result = __scalbnl (__ieee754_expl (M_LN2l * fractx), intx); result = __scalbnl (__ieee754_expl (M_LN2l * fractx), intx);
if (result < LDBL_MIN) math_check_force_underflow_nonneg (result);
{
long double force_underflow = result * result;
math_force_eval (force_underflow);
}
return result; return result;
} }
else else

6
math/k_casinh.c
View File

@ -180,11 +180,7 @@ __kernel_casinh (__complex__ double x, int adj)
else else
__imag__ res = __ieee754_atan2 (ix, s); __imag__ res = __ieee754_atan2 (ix, s);
} }
if (__real__ res < DBL_MIN) math_check_force_underflow_nonneg (__real__ res);
{
volatile double force_underflow = __real__ res * __real__ res;
(void) force_underflow;
}
} }
else else
{ {

6
math/k_casinhf.c
View File

@ -182,11 +182,7 @@ __kernel_casinhf (__complex__ float x, int adj)
else else
__imag__ res = __ieee754_atan2f (ix, s); __imag__ res = __ieee754_atan2f (ix, s);
} }
if (__real__ res < FLT_MIN) math_check_force_underflow_nonneg (__real__ res);
{
volatile float force_underflow = __real__ res * __real__ res;
(void) force_underflow;
}
} }
else else
{ {

6
math/k_casinhl.c
View File

@ -189,11 +189,7 @@ __kernel_casinhl (__complex__ long double x, int adj)
else else
__imag__ res = __ieee754_atan2l (ix, s); __imag__ res = __ieee754_atan2l (ix, s);
} }
if (__real__ res < LDBL_MIN) math_check_force_underflow_nonneg (__real__ res);
{
volatile long double force_underflow = __real__ res * __real__ res;
(void) force_underflow;
}
} }
else else
{ {

11
math/s_catan.c
View File

@ -131,16 +131,7 @@ __catan (__complex__ double x)
} }
} }
if (fabs (__real__ res) < DBL_MIN) math_check_force_underflow_complex (res);
{
volatile double force_underflow = __real__ res * __real__ res;
(void) force_underflow;
}
if (fabs (__imag__ res) < DBL_MIN)
{
volatile double force_underflow = __imag__ res * __imag__ res;
(void) force_underflow;
}
} }
return res; return res;

11
math/s_catanf.c
View File

@ -133,16 +133,7 @@ __catanf (__complex__ float x)
} }
} }
if (fabsf (__real__ res) < FLT_MIN) math_check_force_underflow_complex (res);
{
volatile float force_underflow = __real__ res * __real__ res;
(void) force_underflow;
}
if (fabsf (__imag__ res) < FLT_MIN)
{
volatile float force_underflow = __imag__ res * __imag__ res;
(void) force_underflow;
}
} }
return res; return res;

11
math/s_catanh.c
View File

@ -125,16 +125,7 @@ __catanh (__complex__ double x)
__imag__ res = 0.5 * __ieee754_atan2 (2.0 * __imag__ x, den); __imag__ res = 0.5 * __ieee754_atan2 (2.0 * __imag__ x, den);
} }
if (fabs (__real__ res) < DBL_MIN) math_check_force_underflow_complex (res);
{
volatile double force_underflow = __real__ res * __real__ res;
(void) force_underflow;
}
if (fabs (__imag__ res) < DBL_MIN)
{
volatile double force_underflow = __imag__ res * __imag__ res;
(void) force_underflow;
}
} }
return res; return res;

11
math/s_catanhf.c
View File

@ -127,16 +127,7 @@ __catanhf (__complex__ float x)
__imag__ res = 0.5f * __ieee754_atan2f (2.0f * __imag__ x, den); __imag__ res = 0.5f * __ieee754_atan2f (2.0f * __imag__ x, den);
} }
if (fabsf (__real__ res) < FLT_MIN) math_check_force_underflow_complex (res);
{
volatile float force_underflow = __real__ res * __real__ res;
(void) force_underflow;
}
if (fabsf (__imag__ res) < FLT_MIN)
{
volatile float force_underflow = __imag__ res * __imag__ res;
(void) force_underflow;
}
} }
return res; return res;

11
math/s_catanhl.c
View File

@ -133,16 +133,7 @@ __catanhl (__complex__ long double x)
__imag__ res = 0.5L * __ieee754_atan2l (2.0L * __imag__ x, den); __imag__ res = 0.5L * __ieee754_atan2l (2.0L * __imag__ x, den);
} }
if (fabsl (__real__ res) < LDBL_MIN) math_check_force_underflow_complex (res);
{
volatile long double force_underflow = __real__ res * __real__ res;
(void) force_underflow;
}
if (fabsl (__imag__ res) < LDBL_MIN)
{
volatile long double force_underflow = __imag__ res * __imag__ res;
(void) force_underflow;
}
} }
return res; return res;

11
math/s_catanl.c
View File

@ -139,16 +139,7 @@ __catanl (__complex__ long double x)
} }
} }
if (fabsl (__real__ res) < LDBL_MIN) math_check_force_underflow_complex (res);
{
volatile long double force_underflow = __real__ res * __real__ res;
(void) force_underflow;
}
if (fabsl (__imag__ res) < LDBL_MIN)
{
volatile long double force_underflow = __imag__ res * __imag__ res;
(void) force_underflow;
}
} }
return res; return res;

13
math/s_ccosh.c
View File

@ -83,18 +83,7 @@ __ccosh (__complex__ double x)
__imag__ retval = __ieee754_sinh (__real__ x) * sinix; __imag__ retval = __ieee754_sinh (__real__ x) * sinix;
} }
if (fabs (__real__ retval) < DBL_MIN) math_check_force_underflow_complex (retval);
{
volatile double force_underflow
= __real__ retval * __real__ retval;
(void) force_underflow;
}
if (fabs (__imag__ retval) < DBL_MIN)
{
volatile double force_underflow
= __imag__ retval * __imag__ retval;
(void) force_underflow;
}
} }
else else
{ {

13
math/s_ccoshf.c
View File

@ -83,18 +83,7 @@ __ccoshf (__complex__ float x)
__imag__ retval = __ieee754_sinhf (__real__ x) * sinix; __imag__ retval = __ieee754_sinhf (__real__ x) * sinix;
} }
if (fabsf (__real__ retval) < FLT_MIN) math_check_force_underflow_complex (retval);
{
volatile float force_underflow
= __real__ retval * __real__ retval;
(void) force_underflow;
}
if (fabsf (__imag__ retval) < FLT_MIN)
{
volatile float force_underflow
= __imag__ retval * __imag__ retval;
(void) force_underflow;
}
} }
else else
{ {

13
math/s_ccoshl.c
View File

@ -83,18 +83,7 @@ __ccoshl (__complex__ long double x)
__imag__ retval = __ieee754_sinhl (__real__ x) * sinix; __imag__ retval = __ieee754_sinhl (__real__ x) * sinix;
} }
if (fabsl (__real__ retval) < LDBL_MIN) math_check_force_underflow_complex (retval);
{
volatile long double force_underflow
= __real__ retval * __real__ retval;
(void) force_underflow;
}
if (fabsl (__imag__ retval) < LDBL_MIN)
{
volatile long double force_underflow
= __imag__ retval * __imag__ retval;
(void) force_underflow;
}
} }
else else
{ {

13
math/s_cexp.c
View File

@ -74,18 +74,7 @@ __cexp (__complex__ double x)
__real__ retval = exp_val * cosix; __real__ retval = exp_val * cosix;
__imag__ retval = exp_val * sinix; __imag__ retval = exp_val * sinix;
} }
if (fabs (__real__ retval) < DBL_MIN) math_check_force_underflow_complex (retval);
{
volatile double force_underflow
= __real__ retval * __real__ retval;
(void) force_underflow;
}
if (fabs (__imag__ retval) < DBL_MIN)
{
volatile double force_underflow
= __imag__ retval * __imag__ retval;
(void) force_underflow;
}
} }
else else
{ {

13
math/s_cexpf.c
View File

@ -74,18 +74,7 @@ __cexpf (__complex__ float x)
__real__ retval = exp_val * cosix; __real__ retval = exp_val * cosix;
__imag__ retval = exp_val * sinix; __imag__ retval = exp_val * sinix;
} }
if (fabsf (__real__ retval) < FLT_MIN) math_check_force_underflow_complex (retval);
{
volatile float force_underflow
= __real__ retval * __real__ retval;
(void) force_underflow;
}
if (fabsf (__imag__ retval) < FLT_MIN)
{
volatile float force_underflow
= __imag__ retval * __imag__ retval;
(void) force_underflow;
}
} }
else else
{ {

13
math/s_cexpl.c
View File

@ -74,18 +74,7 @@ __cexpl (__complex__ long double x)
__real__ retval = exp_val * cosix; __real__ retval = exp_val * cosix;
__imag__ retval = exp_val * sinix; __imag__ retval = exp_val * sinix;
} }
if (fabsl (__real__ retval) < LDBL_MIN) math_check_force_underflow_complex (retval);
{
volatile long double force_underflow
= __real__ retval * __real__ retval;
(void) force_underflow;
}
if (fabsl (__imag__ retval) < LDBL_MIN)
{
volatile long double force_underflow
= __imag__ retval * __imag__ retval;
(void) force_underflow;
}
} }
else else
{ {

11
math/s_clog.c
View File

@ -65,15 +65,8 @@ __clog (__complex__ double x)
if (absx == 1.0 && scale == 0) if (absx == 1.0 && scale == 0)
{ {
double absy2 = absy * absy; __real__ result = __log1p (absy * absy) / 2.0;
if (absy2 <= DBL_MIN * 2.0) math_check_force_underflow_nonneg (__real__ result);
{
double force_underflow = absy2 * absy2;
__real__ result = absy2 / 2.0;
math_force_eval (force_underflow);
}
else
__real__ result = __log1p (absy2) / 2.0;
} }
else if (absx > 1.0 && absx < 2.0 && absy < 1.0 && scale == 0) else if (absx > 1.0 && absx < 2.0 && absy < 1.0 && scale == 0)
{ {

11
math/s_clog10.c
View File

@ -71,15 +71,8 @@ __clog10 (__complex__ double x)
if (absx == 1.0 && scale == 0) if (absx == 1.0 && scale == 0)
{ {
double absy2 = absy * absy; __real__ result = __log1p (absy * absy) * (M_LOG10E / 2.0);
if (absy2 <= DBL_MIN * 2.0 * M_LN10) math_check_force_underflow_nonneg (__real__ result);
{
double force_underflow = absy2 * absy2;
__real__ result = absy2 * (M_LOG10E / 2.0);
math_force_eval (force_underflow);
}
else
__real__ result = __log1p (absy2) * (M_LOG10E / 2.0);
} }
else if (absx > 1.0 && absx < 2.0 && absy < 1.0 && scale == 0) else if (absx > 1.0 && absx < 2.0 && absy < 1.0 && scale == 0)
{ {

11
math/s_clog10f.c
View File

@ -71,15 +71,8 @@ __clog10f (__complex__ float x)
if (absx == 1.0f && scale == 0) if (absx == 1.0f && scale == 0)
{ {
float absy2 = absy * absy; __real__ result = __log1pf (absy * absy) * ((float) M_LOG10E / 2.0f);
if (absy2 <= FLT_MIN * 2.0f * (float) M_LN10) math_check_force_underflow_nonneg (__real__ result);
{
float force_underflow = absy2 * absy2;
__real__ result = absy2 * ((float) M_LOG10E / 2.0f);
math_force_eval (force_underflow);
}
else
__real__ result = __log1pf (absy2) * ((float) M_LOG10E / 2.0f);
} }
else if (absx > 1.0f && absx < 2.0f && absy < 1.0f && scale == 0) else if (absx > 1.0f && absx < 2.0f && absy < 1.0f && scale == 0)
{ {

11
math/s_clog10l.c
View File

@ -78,15 +78,8 @@ __clog10l (__complex__ long double x)
if (absx == 1.0L && scale == 0) if (absx == 1.0L && scale == 0)
{ {
long double absy2 = absy * absy; __real__ result = __log1pl (absy * absy) * (M_LOG10El / 2.0L);
if (absy2 <= LDBL_MIN * 2.0L * M_LN10l) math_check_force_underflow_nonneg (__real__ result);
{
long double force_underflow = absy2 * absy2;
__real__ result = absy2 * (M_LOG10El / 2.0);
math_force_eval (force_underflow);
}
else
__real__ result = __log1pl (absy2) * (M_LOG10El / 2.0L);
} }
else if (absx > 1.0L && absx < 2.0L && absy < 1.0L && scale == 0) else if (absx > 1.0L && absx < 2.0L && absy < 1.0L && scale == 0)
{ {

11
math/s_clogf.c
View File

@ -65,15 +65,8 @@ __clogf (__complex__ float x)
if (absx == 1.0f && scale == 0) if (absx == 1.0f && scale == 0)
{ {
float absy2 = absy * absy; __real__ result = __log1pf (absy * absy) / 2.0f;
if (absy2 <= FLT_MIN * 2.0f) math_check_force_underflow_nonneg (__real__ result);
{
float force_underflow = absy2 * absy2;
__real__ result = absy2 / 2.0f;
math_force_eval (force_underflow);
}
else
__real__ result = __log1pf (absy2) / 2.0f;
} }
else if (absx > 1.0f && absx < 2.0f && absy < 1.0f && scale == 0) else if (absx > 1.0f && absx < 2.0f && absy < 1.0f && scale == 0)
{ {

11
math/s_clogl.c
View File

@ -72,15 +72,8 @@ __clogl (__complex__ long double x)
if (absx == 1.0L && scale == 0) if (absx == 1.0L && scale == 0)
{ {
long double absy2 = absy * absy; __real__ result = __log1pl (absy * absy) / 2.0L;
if (absy2 <= LDBL_MIN * 2.0L) math_check_force_underflow_nonneg (__real__ result);
{
long double force_underflow = absy2 * absy2;
__real__ result = absy2 / 2.0L;
math_force_eval (force_underflow);
}
else
__real__ result = __log1pl (absy2) / 2.0L;
} }
else if (absx > 1.0L && absx < 2.0L && absy < 1.0L && scale == 0) else if (absx > 1.0L && absx < 2.0L && absy < 1.0L && scale == 0)
{ {

13
math/s_csin.c
View File

@ -89,18 +89,7 @@ __csin (__complex__ double x)
__imag__ retval = __ieee754_sinh (__imag__ x) * cosix; __imag__ retval = __ieee754_sinh (__imag__ x) * cosix;
} }
if (fabs (__real__ retval) < DBL_MIN) math_check_force_underflow_complex (retval);
{
volatile double force_underflow
= __real__ retval * __real__ retval;
(void) force_underflow;
}
if (fabs (__imag__ retval) < DBL_MIN)
{
volatile double force_underflow
= __imag__ retval * __imag__ retval;
(void) force_underflow;
}
} }
else else
{ {

13
math/s_csinf.c
View File

@ -89,18 +89,7 @@ __csinf (__complex__ float x)
__imag__ retval = __ieee754_sinhf (__imag__ x) * cosix; __imag__ retval = __ieee754_sinhf (__imag__ x) * cosix;
} }
if (fabsf (__real__ retval) < FLT_MIN) math_check_force_underflow_complex (retval);
{
volatile float force_underflow
= __real__ retval * __real__ retval;
(void) force_underflow;
}
if (fabsf (__imag__ retval) < FLT_MIN)
{
volatile float force_underflow
= __imag__ retval * __imag__ retval;
(void) force_underflow;
}
} }
else else
{ {

13
math/s_csinh.c
View File

@ -89,18 +89,7 @@ __csinh (__complex__ double x)
__imag__ retval = __ieee754_cosh (__real__ x) * sinix; __imag__ retval = __ieee754_cosh (__real__ x) * sinix;
} }
if (fabs (__real__ retval) < DBL_MIN) math_check_force_underflow_complex (retval);
{
volatile double force_underflow
= __real__ retval * __real__ retval;
(void) force_underflow;
}
if (fabs (__imag__ retval) < DBL_MIN)
{
volatile double force_underflow
= __imag__ retval * __imag__ retval;
(void) force_underflow;
}
} }
else else
{ {

13
math/s_csinhf.c
View File

@ -89,18 +89,7 @@ __csinhf (__complex__ float x)
__imag__ retval = __ieee754_coshf (__real__ x) * sinix; __imag__ retval = __ieee754_coshf (__real__ x) * sinix;
} }
if (fabsf (__real__ retval) < FLT_MIN) math_check_force_underflow_complex (retval);
{
volatile float force_underflow
= __real__ retval * __real__ retval;
(void) force_underflow;
}
if (fabsf (__imag__ retval) < FLT_MIN)
{
volatile float force_underflow
= __imag__ retval * __imag__ retval;
(void) force_underflow;
}
} }
else else
{ {

13
math/s_csinhl.c
View File

@ -89,18 +89,7 @@ __csinhl (__complex__ long double x)
__imag__ retval = __ieee754_coshl (__real__ x) * sinix; __imag__ retval = __ieee754_coshl (__real__ x) * sinix;
} }
if (fabsl (__real__ retval) < LDBL_MIN) math_check_force_underflow_complex (retval);
{
volatile long double force_underflow
= __real__ retval * __real__ retval;
(void) force_underflow;
}
if (fabsl (__imag__ retval) < LDBL_MIN)
{
volatile long double force_underflow
= __imag__ retval * __imag__ retval;
(void) force_underflow;
}
} }
else else
{ {

13
math/s_csinl.c
View File

@ -89,18 +89,7 @@ __csinl (__complex__ long double x)
__imag__ retval = __ieee754_sinhl (__imag__ x) * cosix; __imag__ retval = __ieee754_sinhl (__imag__ x) * cosix;
} }
if (fabsl (__real__ retval) < LDBL_MIN) math_check_force_underflow_complex (retval);
{
volatile long double force_underflow
= __real__ retval * __real__ retval;
(void) force_underflow;
}
if (fabsl (__imag__ retval) < LDBL_MIN)
{
volatile long double force_underflow
= __imag__ retval * __imag__ retval;
(void) force_underflow;
}
} }
else else
{ {

12
math/s_csqrt.c
View File

@ -148,16 +148,8 @@ __csqrt (__complex__ double x)
s = __scalbn (s, scale); s = __scalbn (s, scale);
} }
if (fabs (r) < DBL_MIN) math_check_force_underflow (r);
{ math_check_force_underflow (s);
double force_underflow = r * r;
math_force_eval (force_underflow);
}
if (fabs (s) < DBL_MIN)
{
double force_underflow = s * s;
math_force_eval (force_underflow);
}
__real__ res = r; __real__ res = r;
__imag__ res = __copysign (s, __imag__ x); __imag__ res = __copysign (s, __imag__ x);

12
math/s_csqrtf.c
View File

@ -148,16 +148,8 @@ __csqrtf (__complex__ float x)
s = __scalbnf (s, scale); s = __scalbnf (s, scale);
} }
if (fabsf (r) < FLT_MIN) math_check_force_underflow (r);
{ math_check_force_underflow (s);
float force_underflow = r * r;
math_force_eval (force_underflow);
}
if (fabsf (s) < FLT_MIN)
{
float force_underflow = s * s;
math_force_eval (force_underflow);
}
__real__ res = r; __real__ res = r;
__imag__ res = __copysignf (s, __imag__ x); __imag__ res = __copysignf (s, __imag__ x);

12
math/s_csqrtl.c
View File

@ -148,16 +148,8 @@ __csqrtl (__complex__ long double x)
s = __scalbnl (s, scale); s = __scalbnl (s, scale);
} }
if (fabsl (r) < LDBL_MIN) math_check_force_underflow (r);
{ math_check_force_underflow (s);
long double force_underflow = r * r;
math_force_eval (force_underflow);
}
if (fabsl (s) < LDBL_MIN)
{
long double force_underflow = s * s;
math_force_eval (force_underflow);
}
__real__ res = r; __real__ res = r;
__imag__ res = __copysignl (s, __imag__ x); __imag__ res = __copysignl (s, __imag__ x);

11
math/s_ctan.c
View File

@ -117,16 +117,7 @@ __ctan (__complex__ double x)
__real__ res = sinrx * cosrx / den; __real__ res = sinrx * cosrx / den;
__imag__ res = sinhix * coshix / den; __imag__ res = sinhix * coshix / den;
} }
if (fabs (__real__ res) < DBL_MIN) math_check_force_underflow_complex (res);
{
double force_underflow = __real__ res * __real__ res;
math_force_eval (force_underflow);
}
if (fabs (__imag__ res) < DBL_MIN)
{
double force_underflow = __imag__ res * __imag__ res;
math_force_eval (force_underflow);
}
} }
return res; return res;

11
math/s_ctanf.c
View File

@ -117,16 +117,7 @@ __ctanf (__complex__ float x)
__real__ res = sinrx * cosrx / den; __real__ res = sinrx * cosrx / den;
__imag__ res = sinhix * coshix / den; __imag__ res = sinhix * coshix / den;
} }
if (fabsf (__real__ res) < FLT_MIN) math_check_force_underflow_complex (res);
{
float force_underflow = __real__ res * __real__ res;
math_force_eval (force_underflow);
}
if (fabsf (__imag__ res) < FLT_MIN)
{
float force_underflow = __imag__ res * __imag__ res;
math_force_eval (force_underflow);
}
} }
return res; return res;

11
math/s_ctanh.c
View File

@ -117,16 +117,7 @@ __ctanh (__complex__ double x)
__real__ res = sinhrx * coshrx / den; __real__ res = sinhrx * coshrx / den;
__imag__ res = sinix * cosix / den; __imag__ res = sinix * cosix / den;
} }
if (fabs (__real__ res) < DBL_MIN) math_check_force_underflow_complex (res);
{
double force_underflow = __real__ res * __real__ res;
math_force_eval (force_underflow);
}
if (fabs (__imag__ res) < DBL_MIN)
{
double force_underflow = __imag__ res * __imag__ res;
math_force_eval (force_underflow);
}
} }
return res; return res;

11
math/s_ctanhf.c
View File

@ -117,16 +117,7 @@ __ctanhf (__complex__ float x)
__real__ res = sinhrx * coshrx / den; __real__ res = sinhrx * coshrx / den;
__imag__ res = sinix * cosix / den; __imag__ res = sinix * cosix / den;
} }
if (fabsf (__real__ res) < FLT_MIN) math_check_force_underflow_complex (res);
{
float force_underflow = __real__ res * __real__ res;
math_force_eval (force_underflow);
}
if (fabsf (__imag__ res) < FLT_MIN)
{
float force_underflow = __imag__ res * __imag__ res;
math_force_eval (force_underflow);
}
} }
return res; return res;

11
math/s_ctanhl.c
View File

@ -124,16 +124,7 @@ __ctanhl (__complex__ long double x)
__real__ res = sinhrx * coshrx / den; __real__ res = sinhrx * coshrx / den;
__imag__ res = sinix * cosix / den; __imag__ res = sinix * cosix / den;
} }
if (fabsl (__real__ res) < LDBL_MIN) math_check_force_underflow_complex (res);
{
long double force_underflow = __real__ res * __real__ res;
math_force_eval (force_underflow);
}
if (fabsl (__imag__ res) < LDBL_MIN)
{
long double force_underflow = __imag__ res * __imag__ res;
math_force_eval (force_underflow);
}
} }
return res; return res;

11
math/s_ctanl.c
View File

@ -124,16 +124,7 @@ __ctanl (__complex__ long double x)
__real__ res = sinrx * cosrx / den; __real__ res = sinrx * cosrx / den;
__imag__ res = sinhix * coshix / den; __imag__ res = sinhix * coshix / den;
} }
if (fabsl (__real__ res) < LDBL_MIN) math_check_force_underflow_complex (res);
{
long double force_underflow = __real__ res * __real__ res;
math_force_eval (force_underflow);
}
if (fabsl (__imag__ res) < LDBL_MIN)
{
long double force_underflow = __imag__ res * __imag__ res;
math_force_eval (force_underflow);
}
} }
return res; return res;

4
stdlib/strtod_l.c
View File

@ -297,8 +297,8 @@ round_and_return (mp_limb_t *retval, intmax_t exponent, int negative,
|| (round_limb & ((((mp_limb_t) 1) << round_bit) - 1)) != 0)) || (round_limb & ((((mp_limb_t) 1) << round_bit) - 1)) != 0))
{ {
__set_errno (ERANGE); __set_errno (ERANGE);
volatile FLOAT force_underflow_exception = MIN_VALUE * MIN_VALUE; FLOAT force_underflow = MIN_VALUE * MIN_VALUE;
(void) force_underflow_exception; math_force_eval (force_underflow);
} }
} }

52
sysdeps/generic/math_private.h
View File

@ -429,6 +429,58 @@ extern long double __lgamma_productl (long double t, long double x,
}) })
#endif #endif
#define fabs_tg(x) __builtin_choose_expr \
(__builtin_types_compatible_p (__typeof (x), float), \
__builtin_fabsf (x), \
__builtin_choose_expr \
(__builtin_types_compatible_p (__typeof (x), double), \
__builtin_fabs (x), __builtin_fabsl (x)))
#define min_of_type(type) __builtin_choose_expr \
(__builtin_types_compatible_p (type, float), \
FLT_MIN, \
__builtin_choose_expr \
(__builtin_types_compatible_p (type, double), \
DBL_MIN, LDBL_MIN))
/* If X (which is not a NaN) is subnormal, force an underflow
exception. */
#define math_check_force_underflow(x) \
do \
{ \
__typeof (x) force_underflow_tmp = (x); \
if (fabs_tg (force_underflow_tmp) \
< min_of_type (__typeof (force_underflow_tmp))) \
{ \
__typeof (force_underflow_tmp) force_underflow_tmp2 \
= force_underflow_tmp * force_underflow_tmp; \
math_force_eval (force_underflow_tmp2); \
} \
} \
while (0)
/* Likewise, but X is also known to be nonnegative. */
#define math_check_force_underflow_nonneg(x) \
do \
{ \
__typeof (x) force_underflow_tmp = (x); \
if (force_underflow_tmp \
< min_of_type (__typeof (force_underflow_tmp))) \
{ \
__typeof (force_underflow_tmp) force_underflow_tmp2 \
= force_underflow_tmp * force_underflow_tmp; \
math_force_eval (force_underflow_tmp2); \
} \
} \
while (0)
/* Likewise, for both real and imaginary parts of a complex
result. */
#define math_check_force_underflow_complex(x) \
do \
{ \
__typeof (x) force_underflow_complex_tmp = (x); \
math_check_force_underflow (__real__ force_underflow_complex_tmp); \
math_check_force_underflow (__imag__ force_underflow_complex_tmp); \
} \
while (0)
/* The standards only specify one variant of the fenv.h interfaces. /* The standards only specify one variant of the fenv.h interfaces.
But at least for some architectures we can be more efficient if we But at least for some architectures we can be more efficient if we

6
sysdeps/ieee754/dbl-64/e_asin.c
View File

@ -71,11 +71,7 @@ __ieee754_asin(double x){
if (k < 0x3e500000) if (k < 0x3e500000)
{ {
if (fabs (x) < DBL_MIN) math_check_force_underflow (x);
{
double force_underflow = x * x;
math_force_eval (force_underflow);
}
return x; /* for x->0 => sin(x)=x */ return x; /* for x->0 => sin(x)=x */
} }
/*----------------------2^-26 <= |x| < 2^ -3 -----------------*/ /*----------------------2^-26 <= |x| < 2^ -3 -----------------*/

6
sysdeps/ieee754/dbl-64/e_atanh.c
View File

@ -52,11 +52,7 @@ __ieee754_atanh (double x)
if (__glibc_unlikely (xa < 0x1.0p-28)) if (__glibc_unlikely (xa < 0x1.0p-28))
{ {
math_force_eval (huge + x); math_force_eval (huge + x);
if (fabs (x) < DBL_MIN) math_check_force_underflow (x);
{
double force_underflow = x * x;
math_force_eval (force_underflow);
}
return x; return x;
} }

5
sysdeps/ieee754/dbl-64/e_exp.c
View File

@ -200,10 +200,7 @@ __ieee754_exp (double x)
check_uflow_ret: check_uflow_ret:
if (retval < DBL_MIN) if (retval < DBL_MIN)
{ {
#if FLT_EVAL_METHOD != 0 double force_underflow = tiny * tiny;
volatile
#endif
double force_underflow = tiny * tiny;
math_force_eval (force_underflow); math_force_eval (force_underflow);
} }
if (retval == 0) if (retval == 0)

6
sysdeps/ieee754/dbl-64/e_exp2.c
View File

@ -122,11 +122,7 @@ __ieee754_exp2 (double x)
else else
{ {
result *= scale_u.d; result *= scale_u.d;
if (result < DBL_MIN) math_check_force_underflow_nonneg (result);
{
double force_underflow = result * result;
math_force_eval (force_underflow);
}
return result; return result;
} }
} }

6
sysdeps/ieee754/dbl-64/e_gamma_r.c
View File

@ -187,11 +187,7 @@ __ieee754_gamma_r (double x, int *signgamp)
double tret = M_PI / (-x * sinpix double tret = M_PI / (-x * sinpix
* gamma_positive (-x, &exp2_adj)); * gamma_positive (-x, &exp2_adj));
ret = __scalbn (tret, -exp2_adj); ret = __scalbn (tret, -exp2_adj);
if (ret < DBL_MIN) math_check_force_underflow_nonneg (ret);
{
double force_underflow = ret * ret;
math_force_eval (force_underflow);
}
} }
} }
ret = math_narrow_eval (ret); ret = math_narrow_eval (ret);

6
sysdeps/ieee754/dbl-64/e_j1.c
View File

@ -128,11 +128,7 @@ __ieee754_j1 (double x)
if (huge + x > one) /* inexact if x!=0 necessary */ if (huge + x > one) /* inexact if x!=0 necessary */
{ {
double ret = 0.5 * x; double ret = 0.5 * x;
if (fabs (ret) < DBL_MIN) math_check_force_underflow (ret);
{
double force_underflow = ret * ret;
math_force_eval (force_underflow);
}
return ret; return ret;
} }
} }

7
sysdeps/ieee754/dbl-64/e_jn.c
View File

@ -246,11 +246,8 @@ __ieee754_jn (int n, double x)
} }
if (ret == 0) if (ret == 0)
ret = __copysign (DBL_MIN, ret) * DBL_MIN; ret = __copysign (DBL_MIN, ret) * DBL_MIN;
else if (fabs (ret) < DBL_MIN) else
{ math_check_force_underflow (ret);
double force_underflow = ret * ret;
math_force_eval (force_underflow);
}
return ret; return ret;
} }
strong_alias (__ieee754_jn, __jn_finite) strong_alias (__ieee754_jn, __jn_finite)

6
sysdeps/ieee754/dbl-64/e_sinh.c
View File

@ -60,11 +60,7 @@ __ieee754_sinh (double x)
if (ix < 0x40360000) /* |x|<22 */ if (ix < 0x40360000) /* |x|<22 */
{ {
if (__glibc_unlikely (ix < 0x3e300000)) { /* |x|<2**-28 */ if (__glibc_unlikely (ix < 0x3e300000)) { /* |x|<2**-28 */
if (fabs (x) < DBL_MIN) math_check_force_underflow (x);
{
double force_underflow = x * x;
math_force_eval (force_underflow);
}
if (shuge + x > one) if (shuge + x > one)
return x; return x;
/* sinh(tiny) = tiny with inexact */ /* sinh(tiny) = tiny with inexact */

6
sysdeps/ieee754/dbl-64/s_asinh.c
View File

@ -39,11 +39,7 @@ __asinh (double x)
ix = hx & 0x7fffffff; ix = hx & 0x7fffffff;
if (__glibc_unlikely (ix < 0x3e300000)) /* |x|<2**-28 */ if (__glibc_unlikely (ix < 0x3e300000)) /* |x|<2**-28 */
{ {
if (fabs (x) < DBL_MIN) math_check_force_underflow (x);
{
double force_underflow = x * x;
math_force_eval (force_underflow);
}
if (huge + x > one) if (huge + x > one)
return x; /* return x inexact except 0 */ return x; /* return x inexact except 0 */
} }

6
sysdeps/ieee754/dbl-64/s_atan.c
View File

@ -90,11 +90,7 @@ atan (double x)
{ {
if (u < A) if (u < A)
{ {
if (u < DBL_MIN) math_check_force_underflow_nonneg (u);
{
double force_underflow = x * x;
math_force_eval (force_underflow);
}
return x; return x;
} }
else else

6
sysdeps/ieee754/dbl-64/s_erf.c
View File

@ -213,11 +213,7 @@ __erf (double x)
{ {
/* Avoid spurious underflow. */ /* Avoid spurious underflow. */
double ret = 0.0625 * (16.0 * x + (16.0 * efx) * x); double ret = 0.0625 * (16.0 * x + (16.0 * efx) * x);
if (fabs (ret) < DBL_MIN) math_check_force_underflow (ret);
{
double force_underflow = ret * ret;
math_force_eval (force_underflow);
}
return ret; return ret;
} }
return x + efx * x; return x + efx * x;

6
sysdeps/ieee754/dbl-64/s_expm1.c
View File

@ -195,11 +195,7 @@ __expm1 (double x)
} }
else if (hx < 0x3c900000) /* when |x|<2**-54, return x */ else if (hx < 0x3c900000) /* when |x|<2**-54, return x */
{ {
if (fabs (x) < DBL_MIN) math_check_force_underflow (x);
{
double force_underflow = x * x;
math_force_eval (force_underflow);
}
t = huge + x; /* return x with inexact flags when x!=0 */ t = huge + x; /* return x with inexact flags when x!=0 */
return x - (t - (huge + x)); return x - (t - (huge + x));
} }

4
sysdeps/ieee754/dbl-64/s_fma.c
View File

@ -90,8 +90,8 @@ __fma (double x, double y, double z)
&& w.ieee.mantissa1 == 0 && w.ieee.mantissa1 == 0
&& w.ieee.mantissa0 == 0))) && w.ieee.mantissa0 == 0)))
{ {
volatile double force_underflow = x * y; double force_underflow = x * y;
(void) force_underflow; math_force_eval (force_underflow);
} }
return v.d * 0x1p-54; return v.d * 0x1p-54;
} }

6
sysdeps/ieee754/dbl-64/s_log1p.c
View File

@ -120,11 +120,7 @@ __log1p (double x)
math_force_eval (two54 + x); /* raise inexact */ math_force_eval (two54 + x); /* raise inexact */
if (ax < 0x3c900000) /* |x| < 2**-54 */ if (ax < 0x3c900000) /* |x| < 2**-54 */
{ {
if (fabs (x) < DBL_MIN) math_check_force_underflow (x);
{
double force_underflow = x * x;
math_force_eval (force_underflow);
}
return x; return x;
} }
else else

6
sysdeps/ieee754/dbl-64/s_sin.c
View File

@ -297,11 +297,7 @@ __sin (double x)
k = 0x7fffffff & m; /* no sign */ k = 0x7fffffff & m; /* no sign */
if (k < 0x3e500000) /* if x->0 =>sin(x)=x */ if (k < 0x3e500000) /* if x->0 =>sin(x)=x */
{ {
if (fabs (x) < DBL_MIN) math_check_force_underflow (x);
{
double force_underflow = x * x;
math_force_eval (force_underflow);
}
retval = x; retval = x;
} }
/*---------------------------- 2^-26 < |x|< 0.25 ----------------------*/ /*---------------------------- 2^-26 < |x|< 0.25 ----------------------*/

6
sysdeps/ieee754/dbl-64/s_tan.c
View File

@ -92,11 +92,7 @@ tan (double x)
/* (I) The case abs(x) <= 1.259e-8 */ /* (I) The case abs(x) <= 1.259e-8 */
if (w <= g1.d) if (w <= g1.d)
{ {
if (w < DBL_MIN) math_check_force_underflow_nonneg (w);
{
double force_underflow = x * x;
math_force_eval (force_underflow);
}
retval = x; retval = x;
goto ret; goto ret;
} }

6
sysdeps/ieee754/dbl-64/s_tanh.c
View File

@ -70,11 +70,7 @@ __tanh (double x)
return x; /* x == +-0 */ return x; /* x == +-0 */
if (ix < 0x3c800000) /* |x|<2**-55 */ if (ix < 0x3c800000) /* |x|<2**-55 */
{ {
if (fabs (x) < DBL_MIN) math_check_force_underflow (x);
{
double force_underflow = x * x;
math_force_eval (force_underflow);
}
return x * (one + x); /* tanh(small) = small */ return x * (one + x); /* tanh(small) = small */
} }
if (ix >= 0x3ff00000) /* |x|>=1 */ if (ix >= 0x3ff00000) /* |x|>=1 */

6
sysdeps/ieee754/flt-32/e_asinf.c
View File

@ -73,11 +73,7 @@ float __ieee754_asinf(float x)
return (x-x)/(x-x); /* asin(|x|>1) is NaN */ return (x-x)/(x-x); /* asin(|x|>1) is NaN */
} else if (ix<0x3f000000) { /* |x|<0.5 */ } else if (ix<0x3f000000) { /* |x|<0.5 */
if(ix<0x32000000) { /* if |x| < 2**-27 */ if(ix<0x32000000) { /* if |x| < 2**-27 */
if (fabsf (x) < FLT_MIN) math_check_force_underflow (x);
{
float force_underflow = x * x;
math_force_eval (force_underflow);
}
if(huge+x>one) return x;/* return x with inexact if x!=0*/ if(huge+x>one) return x;/* return x with inexact if x!=0*/
} else { } else {
t = x*x; t = x*x;

6
sysdeps/ieee754/flt-32/e_atanhf.c
View File

@ -52,11 +52,7 @@ __ieee754_atanhf (float x)
if (__glibc_unlikely (xa < 0x1.0p-28f)) if (__glibc_unlikely (xa < 0x1.0p-28f))
{ {
math_force_eval (huge + x); math_force_eval (huge + x);
if (fabsf (x) < FLT_MIN) math_check_force_underflow (x);
{
float force_underflow = x * x;
math_force_eval (force_underflow);
}
return x; return x;
} }

6
sysdeps/ieee754/flt-32/e_exp2f.c
View File

@ -111,11 +111,7 @@ __ieee754_exp2f (float x)
else else
{ {
result *= scale_u.f; result *= scale_u.f;
if (result < FLT_MIN) math_check_force_underflow_nonneg (result);
{
float force_underflow = result * result;
math_force_eval (force_underflow);
}
return result; return result;
} }
} }

6
sysdeps/ieee754/flt-32/e_gammaf_r.c
View File

@ -179,11 +179,7 @@ __ieee754_gammaf_r (float x, int *signgamp)
float tret = (float) M_PI / (-x * sinpix float tret = (float) M_PI / (-x * sinpix
* gammaf_positive (-x, &exp2_adj)); * gammaf_positive (-x, &exp2_adj));
ret = __scalbnf (tret, -exp2_adj); ret = __scalbnf (tret, -exp2_adj);
if (ret < FLT_MIN) math_check_force_underflow_nonneg (ret);
{
float force_underflow = ret * ret;
math_force_eval (force_underflow);
}
} }
} }
ret = math_narrow_eval (ret); ret = math_narrow_eval (ret);

5
sysdeps/ieee754/flt-32/e_j1f.c
View File

@ -72,10 +72,7 @@ __ieee754_j1f(float x)
if(__builtin_expect(ix<0x32000000, 0)) { /* |x|<2**-27 */ if(__builtin_expect(ix<0x32000000, 0)) { /* |x|<2**-27 */
if(huge+x>one) { /* inexact if x!=0 necessary */ if(huge+x>one) { /* inexact if x!=0 necessary */
float ret = (float) 0.5 * x; float ret = (float) 0.5 * x;
if (fabsf (ret) < FLT_MIN) { math_check_force_underflow (ret);
float force_underflow = ret * ret;
math_force_eval (force_underflow);
}
return ret; return ret;
} }
} }

6
sysdeps/ieee754/flt-32/e_jnf.c
View File

@ -170,10 +170,8 @@ __ieee754_jnf(int n, float x)
} }
if (ret == 0) if (ret == 0)
ret = __copysignf (FLT_MIN, ret) * FLT_MIN; ret = __copysignf (FLT_MIN, ret) * FLT_MIN;
else if (fabsf (ret) < FLT_MIN) { else
float force_underflow = ret * ret; math_check_force_underflow (ret);
math_force_eval (force_underflow);
}
return ret; return ret;
} }
strong_alias (__ieee754_jnf, __jnf_finite) strong_alias (__ieee754_jnf, __jnf_finite)

6
sysdeps/ieee754/flt-32/e_sinhf.c
View File

@ -36,11 +36,7 @@ __ieee754_sinhf(float x)
/* |x| in [0,22], return sign(x)*0.5*(E+E/(E+1))) */ /* |x| in [0,22], return sign(x)*0.5*(E+E/(E+1))) */
if (ix < 0x41b00000) { /* |x|<22 */ if (ix < 0x41b00000) { /* |x|<22 */
if (__builtin_expect(ix<0x31800000, 0)) { /* |x|<2**-28 */ if (__builtin_expect(ix<0x31800000, 0)) { /* |x|<2**-28 */
if (fabsf (x) < FLT_MIN) math_check_force_underflow (x);
{
float force_underflow = x * x;
math_force_eval (force_underflow);
}
if(shuge+x>one) return x;/* sinh(tiny) = tiny with inexact */ if(shuge+x>one) return x;/* sinh(tiny) = tiny with inexact */
} }
t = __expm1f(fabsf(x)); t = __expm1f(fabsf(x));

6
sysdeps/ieee754/flt-32/k_sinf.c
View File

@ -38,11 +38,7 @@ float __kernel_sinf(float x, float y, int iy)
ix &= 0x7fffffff; /* high word of x */ ix &= 0x7fffffff; /* high word of x */
if(ix<0x32000000) /* |x| < 2**-27 */ if(ix<0x32000000) /* |x| < 2**-27 */
{ {
if (fabsf (x) < FLT_MIN) math_check_force_underflow (x);
{
float force_underflow = x * x;
math_force_eval (force_underflow);
}
if ((int) x == 0) if ((int) x == 0)
return x; /* generate inexact */ return x; /* generate inexact */
} }

6
sysdeps/ieee754/flt-32/k_tanf.c
View File

@ -51,11 +51,7 @@ float __kernel_tanf(float x, float y, int iy)
if((ix|(iy+1))==0) return one/fabsf(x); if((ix|(iy+1))==0) return one/fabsf(x);
else if (iy == 1) else if (iy == 1)
{ {
if (fabsf (x) < FLT_MIN) math_check_force_underflow (x);
{
float force_underflow = x * x;
math_force_eval (force_underflow);
}
return x; return x;
} }
else else

6
sysdeps/ieee754/flt-32/s_asinhf.c
View File

@ -30,11 +30,7 @@ __asinhf(float x)
GET_FLOAT_WORD(hx,x); GET_FLOAT_WORD(hx,x);
ix = hx&0x7fffffff; ix = hx&0x7fffffff;
if(__builtin_expect(ix< 0x38000000, 0)) { /* |x|<2**-14 */ if(__builtin_expect(ix< 0x38000000, 0)) { /* |x|<2**-14 */
if (fabsf (x) < FLT_MIN) math_check_force_underflow (x);
{
float force_underflow = x * x;
math_force_eval (force_underflow);
}
if(huge+x>one) return x; /* return x inexact except 0 */ if(huge+x>one) return x; /* return x inexact except 0 */
} }
if(__builtin_expect(ix>0x47000000, 0)) { /* |x| > 2**14 */ if(__builtin_expect(ix>0x47000000, 0)) { /* |x| > 2**14 */

6
sysdeps/ieee754/flt-32/s_atanf.c
View File

@ -67,11 +67,7 @@ float __atanf(float x)
else return -atanhi[3]-atanlo[3]; else return -atanhi[3]-atanlo[3];
} if (ix < 0x3ee00000) { /* |x| < 0.4375 */ } if (ix < 0x3ee00000) { /* |x| < 0.4375 */
if (ix < 0x31000000) { /* |x| < 2^-29 */ if (ix < 0x31000000) { /* |x| < 2^-29 */
if (fabsf (x) < FLT_MIN) math_check_force_underflow (x);
{
float force_underflow = x * x;
math_force_eval (force_underflow);
}
if(huge+x>one) return x; /* raise inexact */ if(huge+x>one) return x; /* raise inexact */
} }
id = -1; id = -1;

6
sysdeps/ieee754/flt-32/s_erff.c
View File

@ -113,11 +113,7 @@ float __erff(float x)
{ {
/* Avoid spurious underflow. */ /* Avoid spurious underflow. */
float ret = 0.0625f * (16.0f * x + (16.0f * efx) * x); float ret = 0.0625f * (16.0f * x + (16.0f * efx) * x);
if (fabsf (ret) < FLT_MIN) math_check_force_underflow (ret);
{
float force_underflow = ret * ret;
math_force_eval (force_underflow);
}
return ret; return ret;
} }
return x + efx*x; return x + efx*x;

6
sysdeps/ieee754/flt-32/s_expm1f.c
View File

@ -81,11 +81,7 @@ __expm1f(float x)
c = (hi-x)-lo; c = (hi-x)-lo;
} }
else if(hx < 0x33000000) { /* when |x|<2**-25, return x */ else if(hx < 0x33000000) { /* when |x|<2**-25, return x */
if (fabsf (x) < FLT_MIN) math_check_force_underflow (x);
{
float force_underflow = x * x;
math_force_eval (force_underflow);
}
t = huge+x; /* return x with inexact flags when x!=0 */ t = huge+x; /* return x with inexact flags when x!=0 */
return x - (t-(huge+x)); return x - (t-(huge+x));
} }

6
sysdeps/ieee754/flt-32/s_log1pf.c
View File

@ -50,11 +50,7 @@ __log1pf(float x)
math_force_eval(two25+x); /* raise inexact */ math_force_eval(two25+x); /* raise inexact */
if (ax<0x24800000) /* |x| < 2**-54 */ if (ax<0x24800000) /* |x| < 2**-54 */
{ {
if (fabsf (x) < FLT_MIN) math_check_force_underflow (x);
{
float force_underflow = x * x;
math_force_eval (force_underflow);
}
return x; return x;
} }
else else

6
sysdeps/ieee754/flt-32/s_tanhf.c
View File

@ -43,11 +43,7 @@ float __tanhf(float x)
return x; /* x == +-0 */ return x; /* x == +-0 */
if (ix<0x24000000) /* |x|<2**-55 */ if (ix<0x24000000) /* |x|<2**-55 */
{ {
if (fabsf (x) < FLT_MIN) math_check_force_underflow (x);
{
float force_underflow = x * x;
math_force_eval (force_underflow);
}
return x*(one+x); /* tanh(small) = small */ return x*(one+x); /* tanh(small) = small */
} }
if (ix>=0x3f800000) { /* |x|>=1 */ if (ix>=0x3f800000) { /* |x|>=1 */

6
sysdeps/ieee754/ldbl-128/e_asinl.c
View File

@ -153,11 +153,7 @@ __ieee754_asinl (long double x)
{ {
if (ix < 0x3fc60000) /* |x| < 2**-57 */ if (ix < 0x3fc60000) /* |x| < 2**-57 */
{ {
if (fabsl (x) < LDBL_MIN) math_check_force_underflow (x);
{
long double force_underflow = x * x;
math_force_eval (force_underflow);
}
long double force_inexact = huge + x; long double force_inexact = huge + x;
math_force_eval (force_inexact); math_force_eval (force_inexact);
return x; /* return x with inexact if x!=0 */ return x; /* return x with inexact if x!=0 */

6
sysdeps/ieee754/ldbl-128/e_atanhl.c
View File

@ -60,11 +60,7 @@ __ieee754_atanhl(long double x)
} }
if(ix<0x3fc60000 && (huge+x)>zero) /* x < 2^-57 */ if(ix<0x3fc60000 && (huge+x)>zero) /* x < 2^-57 */
{ {
if (fabsl (x) < LDBL_MIN) math_check_force_underflow (x);
{
long double force_underflow = x * x;
math_force_eval (force_underflow);
}
return x; return x;
} }

6
sysdeps/ieee754/ldbl-128/e_expl.c
View File

@ -232,11 +232,7 @@ __ieee754_expl (long double x)
else else
{ {
result *= scale_u.d; result *= scale_u.d;
if (result < LDBL_MIN) math_check_force_underflow_nonneg (result);
{
long double force_underflow = result * result;
math_force_eval (force_underflow);
}
return result; return result;
} }
} }

6
sysdeps/ieee754/ldbl-128/e_gammal_r.c
View File

@ -194,11 +194,7 @@ __ieee754_gammal_r (long double x, int *signgamp)
ret = M_PIl / (-x * sinpix ret = M_PIl / (-x * sinpix
* gammal_positive (-x, &exp2_adj)); * gammal_positive (-x, &exp2_adj));
ret = __scalbnl (ret, -exp2_adj); ret = __scalbnl (ret, -exp2_adj);
if (ret < LDBL_MIN) math_check_force_underflow_nonneg (ret);
{
long double force_underflow = ret * ret;
math_force_eval (force_underflow);
}
} }
} }
} }

6
sysdeps/ieee754/ldbl-128/e_j1l.c
View File

@ -700,11 +700,7 @@ __ieee754_j1l (long double x)
if (xx <= 0x1p-58L) if (xx <= 0x1p-58L)
{ {
long double ret = x * 0.5L; long double ret = x * 0.5L;
if (fabsl (ret) < LDBL_MIN) math_check_force_underflow (ret);
{
long double force_underflow = ret * ret;
math_force_eval (force_underflow);
}
return ret; return ret;
} }
if (xx <= 2.0L) if (xx <= 2.0L)

7
sysdeps/ieee754/ldbl-128/e_jnl.c
View File

@ -297,11 +297,8 @@ __ieee754_jnl (int n, long double x)
} }
if (ret == 0) if (ret == 0)
ret = __copysignl (LDBL_MIN, ret) * LDBL_MIN; ret = __copysignl (LDBL_MIN, ret) * LDBL_MIN;
else if (fabsl (ret) < LDBL_MIN) else
{ math_check_force_underflow (ret);
long double force_underflow = ret * ret;
math_force_eval (force_underflow);
}
return ret; return ret;
} }
strong_alias (__ieee754_jnl, __jnl_finite) strong_alias (__ieee754_jnl, __jnl_finite)

6
sysdeps/ieee754/ldbl-128/e_sinhl.c
View File

@ -88,11 +88,7 @@ __ieee754_sinhl (long double x)
{ {
if (ix < 0x3fc60000) /* |x| < 2^-57 */ if (ix < 0x3fc60000) /* |x| < 2^-57 */
{ {
if (fabsl (x) < LDBL_MIN) math_check_force_underflow (x);
{
long double force_underflow = x * x;
math_force_eval (force_underflow);
}
if (shuge + x > one) if (shuge + x > one)
return x; /* sinh(tiny) = tiny with inexact */ return x; /* sinh(tiny) = tiny with inexact */
} }

6
sysdeps/ieee754/ldbl-128/k_sincosl.c
View File

@ -111,11 +111,7 @@ __kernel_sincosl(long double x, long double y, long double *sinx, long double *c
polynomial of degree 16(17). */ polynomial of degree 16(17). */
if (tix < 0x3fc60000) /* |x| < 2^-57 */ if (tix < 0x3fc60000) /* |x| < 2^-57 */
{ {
if (fabsl (x) < LDBL_MIN) math_check_force_underflow (x);
{
long double force_underflow = x * x;
math_force_eval (force_underflow);
}
if (!((int)x)) /* generate inexact */ if (!((int)x)) /* generate inexact */
{ {
*sinx = x; *sinx = x;

6
sysdeps/ieee754/ldbl-128/k_sinl.c
View File

@ -92,11 +92,7 @@ __kernel_sinl(long double x, long double y, int iy)
polynomial of degree 17. */ polynomial of degree 17. */
if (tix < 0x3fc60000) /* |x| < 2^-57 */ if (tix < 0x3fc60000) /* |x| < 2^-57 */
{ {
if (fabsl (x) < LDBL_MIN) math_check_force_underflow (x);
{
long double force_underflow = x * x;
math_force_eval (force_underflow);
}
if (!((int)x)) return x; /* generate inexact */ if (!((int)x)) return x; /* generate inexact */
} }
z = x * x; z = x * x;

6
sysdeps/ieee754/ldbl-128/k_tanl.c
View File

@ -101,11 +101,7 @@ __kernel_tanl (long double x, long double y, int iy)
return one / fabs (x); return one / fabs (x);
else if (iy == 1) else if (iy == 1)
{ {
if (fabsl (x) < LDBL_MIN) math_check_force_underflow (x);
{
long double force_underflow = x * x;
math_force_eval (force_underflow);
}
return x; return x;
} }
else else

6
sysdeps/ieee754/ldbl-128/s_asinhl.c
View File

@ -52,11 +52,7 @@ __asinhl (long double x)
return x + x; /* x is inf or NaN */ return x + x; /* x is inf or NaN */
if (ix < 0x3fc70000) if (ix < 0x3fc70000)
{ /* |x| < 2^ -56 */ { /* |x| < 2^ -56 */
if (fabsl (x) < LDBL_MIN) math_check_force_underflow (x);
{
long double force_underflow = x * x;
math_force_eval (force_underflow);
}
if (huge + x > one) if (huge + x > one)
return x; /* return x inexact except 0 */ return x; /* return x inexact except 0 */
} }

6
sysdeps/ieee754/ldbl-128/s_atanl.c
View File

@ -202,11 +202,7 @@ __atanl (long double x)
if (k <= 0x3fc50000) /* |x| < 2**-58 */ if (k <= 0x3fc50000) /* |x| < 2**-58 */
{ {
if (fabsl (x) < LDBL_MIN) math_check_force_underflow (x);
{
long double force_underflow = x * x;
math_force_eval (force_underflow);
}
/* Raise inexact. */ /* Raise inexact. */
if (huge + x > 0.0) if (huge + x > 0.0)
return x; return x;

6
sysdeps/ieee754/ldbl-128/s_erfl.c
View File

@ -785,11 +785,7 @@ __erfl (long double x)
{ {
/* Avoid spurious underflow. */ /* Avoid spurious underflow. */
long double ret = 0.0625 * (16.0 * x + (16.0 * efx) * x); long double ret = 0.0625 * (16.0 * x + (16.0 * efx) * x);
if (fabsl (ret) < LDBL_MIN) math_check_force_underflow (ret);
{
long double force_underflow = ret * ret;
math_force_eval (force_underflow);
}
return ret; return ret;
} }
return x + efx * x; return x + efx * x;

6
sysdeps/ieee754/ldbl-128/s_expm1l.c
View File

@ -128,11 +128,7 @@ __expm1l (long double x)
when the result does underflow. */ when the result does underflow. */
if (fabsl (x) < 0x1p-113L) if (fabsl (x) < 0x1p-113L)
{ {
if (fabsl (x) < LDBL_MIN) math_check_force_underflow (x);
{
long double force_underflow = x * x;
math_force_eval (force_underflow);
}
return x; return x;
} }

4
sysdeps/ieee754/ldbl-128/s_fmal.c
View File

@ -94,8 +94,8 @@ __fmal (long double x, long double y, long double z)
&& w.ieee.mantissa1 == 0 && w.ieee.mantissa1 == 0
&& w.ieee.mantissa0 == 0))) && w.ieee.mantissa0 == 0)))
{ {
volatile long double force_underflow = x * y; long double force_underflow = x * y;
(void) force_underflow; math_force_eval (force_underflow);
} }
return v.d * 0x1p-114L; return v.d * 0x1p-114L;
} }

6
sysdeps/ieee754/ldbl-128/s_log1pl.c
View File

@ -140,11 +140,7 @@ __log1pl (long double xm1)
if ((hx & 0x7fffffff) < 0x3f8e0000) if ((hx & 0x7fffffff) < 0x3f8e0000)
{ {
if (fabsl (xm1) < LDBL_MIN) math_check_force_underflow (xm1);
{
long double force_underflow = xm1 * xm1;
math_force_eval (force_underflow);
}
if ((int) xm1 == 0) if ((int) xm1 == 0)
return xm1; return xm1;
} }

6
sysdeps/ieee754/ldbl-128/s_tanhl.c
View File

@ -75,11 +75,7 @@ __tanhl (long double x)
return x; /* x == +- 0 */ return x; /* x == +- 0 */
if (ix < 0x3fc60000) /* |x| < 2^-57 */ if (ix < 0x3fc60000) /* |x| < 2^-57 */
{ {
if (fabsl (x) < LDBL_MIN) math_check_force_underflow (x);
{
long double force_underflow = x * x;
math_force_eval (force_underflow);
}
return x * (one + tiny); /* tanh(small) = small */ return x * (one + tiny); /* tanh(small) = small */
} }
u.parts32.w0 = ix; /* Absolute value of x. */ u.parts32.w0 = ix; /* Absolute value of x. */

6
sysdeps/ieee754/ldbl-128ibm/e_asinl.c
View File

@ -147,11 +147,7 @@ __ieee754_asinl (long double x)
{ {
if (a < 6.938893903907228e-18L) /* |x| < 2**-57 */ if (a < 6.938893903907228e-18L) /* |x| < 2**-57 */
{ {
if (fabsl (x) < LDBL_MIN) math_check_force_underflow (x);
{
long double force_underflow = x * x;
math_force_eval (force_underflow);
}
long double force_inexact = huge + x; long double force_inexact = huge + x;
math_force_eval (force_inexact); math_force_eval (force_inexact);
return x; /* return x with inexact if x!=0 */ return x; /* return x with inexact if x!=0 */

6
sysdeps/ieee754/ldbl-128ibm/e_atanhl.c
View File

@ -57,11 +57,7 @@ __ieee754_atanhl(long double x)
} }
if(ix<0x3c70000000000000LL&&(huge+x)>zero) /* x<2**-56 */ if(ix<0x3c70000000000000LL&&(huge+x)>zero) /* x<2**-56 */
{ {
if (fabsl (x) < LDBL_MIN) math_check_force_underflow (x);
{
long double force_underflow = x * x;
math_force_eval (force_underflow);
}
return x; return x;
} }
x = fabsl (x); x = fabsl (x);

6
sysdeps/ieee754/ldbl-128ibm/e_gammal_r.c
View File

@ -194,11 +194,7 @@ __ieee754_gammal_r (long double x, int *signgamp)
ret = M_PIl / (-x * sinpix ret = M_PIl / (-x * sinpix
* gammal_positive (-x, &exp2_adj)); * gammal_positive (-x, &exp2_adj));
ret = __scalbnl (ret, -exp2_adj); ret = __scalbnl (ret, -exp2_adj);
if (ret < LDBL_MIN) math_check_force_underflow_nonneg (ret);
{
long double force_underflow = ret * ret;
math_force_eval (force_underflow);
}
} }
} }
} }

7
sysdeps/ieee754/ldbl-128ibm/e_jnl.c
View File

@ -297,11 +297,8 @@ __ieee754_jnl (int n, long double x)
} }
if (ret == 0) if (ret == 0)
ret = __copysignl (LDBL_MIN, ret) * LDBL_MIN; ret = __copysignl (LDBL_MIN, ret) * LDBL_MIN;
else if (fabsl (ret) < LDBL_MIN) else
{ math_check_force_underflow (ret);
long double force_underflow = ret * ret;
math_force_eval (force_underflow);
}
return ret; return ret;
} }
strong_alias (__ieee754_jnl, __jnl_finite) strong_alias (__ieee754_jnl, __jnl_finite)

6
sysdeps/ieee754/ldbl-128ibm/e_sinhl.c
View File

@ -54,11 +54,7 @@ __ieee754_sinhl(long double x)
/* |x| in [0,40], return sign(x)*0.5*(E+E/(E+1))) */ /* |x| in [0,40], return sign(x)*0.5*(E+E/(E+1))) */
if (ix < 0x4044000000000000LL) { /* |x|<40 */ if (ix < 0x4044000000000000LL) { /* |x|<40 */
if (ix<0x3c90000000000000LL) { /* |x|<2**-54 */ if (ix<0x3c90000000000000LL) { /* |x|<2**-54 */
if (fabsl (x) < LDBL_MIN) math_check_force_underflow (x);
{
long double force_underflow = x * x;
math_force_eval (force_underflow);
}
if(shuge+x>one) return x;/* sinhl(tiny) = tiny with inexact */ if(shuge+x>one) return x;/* sinhl(tiny) = tiny with inexact */
} }
t = __expm1l(fabsl(x)); t = __expm1l(fabsl(x));

6
sysdeps/ieee754/ldbl-128ibm/k_sincosl.c
View File

@ -114,11 +114,7 @@ __kernel_sincosl(long double x, long double y, long double *sinx, long double *c
polynomial of degree 16(17). */ polynomial of degree 16(17). */
if (tix < 0x3c600000) /* |x| < 2^-57 */ if (tix < 0x3c600000) /* |x| < 2^-57 */
{ {
if (fabsl (x) < LDBL_MIN) math_check_force_underflow (x);
{
long double force_underflow = x * x;
math_force_eval (force_underflow);
}
if (!((int)x)) /* generate inexact */ if (!((int)x)) /* generate inexact */
{ {
*sinx = x; *sinx = x;

6
sysdeps/ieee754/ldbl-128ibm/k_sinl.c
View File

@ -95,11 +95,7 @@ __kernel_sinl(long double x, long double y, int iy)
polynomial of degree 17. */ polynomial of degree 17. */
if (tix < 0x3c600000) /* |x| < 2^-57 */ if (tix < 0x3c600000) /* |x| < 2^-57 */
{ {
if (fabsl (x) < LDBL_MIN) math_check_force_underflow (x);
{
long double force_underflow = x * x;
math_force_eval (force_underflow);
}
if (!((int)x)) return x; /* generate inexact */ if (!((int)x)) return x; /* generate inexact */
} }
z = x * x; z = x * x;

6
sysdeps/ieee754/ldbl-128ibm/k_tanl.c
View File

@ -101,11 +101,7 @@ __kernel_tanl (long double x, long double y, int iy)
return one / fabs (x); return one / fabs (x);
else if (iy == 1) else if (iy == 1)
{ {
if (fabsl (x) < LDBL_MIN) math_check_force_underflow (x);
{
long double force_underflow = x * x;
math_force_eval (force_underflow);
}
return x; return x;
} }
else else

6
sysdeps/ieee754/ldbl-128ibm/s_asinhl.c
View File

@ -46,11 +46,7 @@ long double __asinhl(long double x)
ix = hx&0x7fffffffffffffffLL; ix = hx&0x7fffffffffffffffLL;
if(ix>=0x7ff0000000000000LL) return x+x; /* x is inf or NaN */ if(ix>=0x7ff0000000000000LL) return x+x; /* x is inf or NaN */
if(ix< 0x3c70000000000000LL) { /* |x|<2**-56 */ if(ix< 0x3c70000000000000LL) { /* |x|<2**-56 */
if (fabsl (x) < LDBL_MIN) math_check_force_underflow (x);
{
long double force_underflow = x * x;
math_force_eval (force_underflow);
}
if(huge+x>one) return x; /* return x inexact except 0 */ if(huge+x>one) return x; /* return x inexact except 0 */
} }
if(ix>0x4370000000000000LL) { /* |x| > 2**56 */ if(ix>0x4370000000000000LL) { /* |x| > 2**56 */

6
sysdeps/ieee754/ldbl-128ibm/s_atanl.c
View File

@ -199,11 +199,7 @@ __atanl (long double x)
if (k <= 0x3c800000) /* |x| <= 2**-55. */ if (k <= 0x3c800000) /* |x| <= 2**-55. */
{ {
if (fabsl (x) < LDBL_MIN) math_check_force_underflow (x);
{
long double force_underflow = x * x;
math_force_eval (force_underflow);
}
/* Raise inexact. */ /* Raise inexact. */
if (1e300L + x > 0.0) if (1e300L + x > 0.0)
return x; return x;

6
sysdeps/ieee754/ldbl-128ibm/s_erfl.c
View File

@ -805,11 +805,7 @@ __erfl (long double x)
if (x == 0) if (x == 0)
return x; return x;
long double ret = 0.0625 * (16.0 * x + (16.0 * efx) * x); long double ret = 0.0625 * (16.0 * x + (16.0 * efx) * x);
if (fabsl (ret) < LDBL_MIN) math_check_force_underflow (ret);
{
long double force_underflow = ret * ret;
math_force_eval (force_underflow);
}
return ret; return ret;
} }
return x + efx * x; return x + efx * x;

Some files were not shown because too many files have changed in this diff Show More