1
0
mirror of https://sourceware.org/git/glibc.git synced 2025-07-29 11:41:21 +03:00

Optimize use of isnan, isinf, finite

This commit is contained in:
Ulrich Drepper
2011-10-08 10:18:26 -04:00
parent 187da0aedc
commit 7edb55ce06
35 changed files with 236 additions and 63 deletions

View File

@ -1,5 +1,41 @@
2011-10-08 Ulrich Drepper <drepper@gmail.com> 2011-10-08 Ulrich Drepper <drepper@gmail.com>
* include/math.h: Declare __isinf_ns, __isinf_nsf, __isinf_nsl.
* sysdeps/ieee754/dbl-64/s_isinf_ns.c: New file.
* sysdeps/ieee754/dbl-64/wordsize-64/s_isinf_ns.c: New file.
* sysdeps/ieee754/flt-32/s_isinf_nsf.c: New file.
* sysdeps/ieee754/ldbl-128/s_isinf_nsl.c: New file.
* sysdeps/ieee754/ldbl-96/s_isinf_nsl.c: New file.
* math/Makefile (libm-calls): Add s_isinf_ns.
* math/divtc3.c: Use __isinf_nsl instead of isinf.
* math/multc3.c: Likewise.
* math/s_casin.c: Likewise.
* math/s_casinf.c: Likewise.
* math/s_casinl.c: Likewise.
* math/s_ccos.c: Likewise.
* math/s_ccosf.c: Likewise.
* math/s_ccosl.c: Likewise.
* math/s_ctan.c: Likewise.
* math/s_ctanf.c: Likewise.
* math/s_ctanh.c: Likewise.
* math/s_ctanhf.c: Likewise.
* math/s_ctanhl.c: Likewise.
* math/s_ctanl.c: Likewise.
* math/w_fmod.c: Likewise.
* math/w_fmodf.c: Likewise.
* math/w_fmodl.c: Likewise.
* math/w_remainder.c: Likewise.
* math/w_remainderf.c: Likewise.
* math/w_remainderl.c: Likewise.
* sysdeps/ieee754/dbl-64/s_finite.c: Undefine __finite.
* sysdeps/ieee754/dbl-64/s_isnan.c: Undefine __isnan.
* sysdeps/ieee754/dbl-64/wordsize-64/s_finite.c: Undefine __finite.
* sysdeps/ieee754/dbl-64/wordsize-64/s_isnan.c: Undefine __nan.
* sysdeps/ieee754/flt-32/s_finitef.c: Undefine __finitef.
* sysdeps/ieee754/flt-32/s_isnanf.c: Undefine __nan.
* sysdeps/x86_64/fpu/math_private.h: Add optimized versions of __isnsn,
__isnanf, __isinf_ns, __isinf_nsf, __finite, and __finitef.
* stdio-common/printf_fp.c: Use the fact that isinf returns the sign * stdio-common/printf_fp.c: Use the fact that isinf returns the sign
of the number. of the number.
* stdio-common/printf_fphex.c: Likewise. * stdio-common/printf_fphex.c: Likewise.

View File

@ -31,4 +31,8 @@ libm_hidden_proto (__expl)
libm_hidden_proto (__expm1l) libm_hidden_proto (__expm1l)
# endif # endif
extern int __isinf_ns (double);
extern int __isinf_nsf (float);
extern int __isinf_nsl (long double);
#endif #endif

View File

@ -1,4 +1,4 @@
# Copyright (C) 1996-2001,2002,2003,2004,2005,2006 # Copyright (C) 1996-2001,2002,2003,2004,2005,2006,2011
# Free Software Foundation, Inc. # Free Software Foundation, Inc.
# This file is part of the GNU C Library. # This file is part of the GNU C Library.
@ -60,7 +60,7 @@ libm-calls = e_acos e_acosh e_asin e_atan2 e_atanh e_cosh e_exp e_fmod \
s_catan s_casin s_ccos s_csin s_ctan s_ctanh s_cacos \ s_catan s_casin s_ccos s_csin s_ctan s_ctanh s_cacos \
s_casinh s_cacosh s_catanh s_csqrt s_cpow s_cproj s_clog10 \ s_casinh s_cacosh s_catanh s_csqrt s_cpow s_cproj s_clog10 \
s_fma s_lrint s_llrint s_lround s_llround e_exp10 w_log2 \ s_fma s_lrint s_llrint s_lround s_llround e_exp10 w_log2 \
$(calls:s_%=m_%) s_isinf_ns $(calls:s_%=m_%)
include ../Makeconfig include ../Makeconfig

View File

@ -1,4 +1,4 @@
/* Copyright (C) 2005, 2006 Free Software Foundation, Inc. /* Copyright (C) 2005, 2006, 2011 Free Software Foundation, Inc.
This file is part of the GNU C Library. This file is part of the GNU C Library.
Contributed by Richard Henderson <rth@redhat.com>, 2005. Contributed by Richard Henderson <rth@redhat.com>, 2005.
@ -55,17 +55,19 @@ __divtc3 (long double a, long double b, long double c, long double d)
x = __copysignl (INFINITY, c) * a; x = __copysignl (INFINITY, c) * a;
y = __copysignl (INFINITY, c) * b; y = __copysignl (INFINITY, c) * b;
} }
else if ((isinf (a) || isinf (b)) && isfinite (c) && isfinite (d)) else if ((__isinf_nsl (a) || __isinf_nsl (b))
&& isfinite (c) && isfinite (d))
{ {
a = __copysignl (isinf (a) ? 1 : 0, a); a = __copysignl (__isinf_nsl (a) ? 1 : 0, a);
b = __copysignl (isinf (b) ? 1 : 0, b); b = __copysignl (__isinf_nsl (b) ? 1 : 0, b);
x = INFINITY * (a * c + b * d); x = INFINITY * (a * c + b * d);
y = INFINITY * (b * c - a * d); y = INFINITY * (b * c - a * d);
} }
else if ((isinf (c) || isinf (d)) && isfinite (a) && isfinite (b)) else if ((__isinf_nsl (c) || __isinf_nsl (d))
&& isfinite (a) && isfinite (b))
{ {
c = __copysignl (isinf (c) ? 1 : 0, c); c = __copysignl (__isinf_nsl (c) ? 1 : 0, c);
d = __copysignl (isinf (d) ? 1 : 0, d); d = __copysignl (__isinf_nsl (d) ? 1 : 0, d);
x = 0.0 * (a * c + b * d); x = 0.0 * (a * c + b * d);
y = 0.0 * (b * c - a * d); y = 0.0 * (b * c - a * d);
} }

View File

@ -1,4 +1,4 @@
/* Copyright (C) 2005, 2006 Free Software Foundation, Inc. /* Copyright (C) 2005, 2006, 2011 Free Software Foundation, Inc.
This file is part of the GNU C Library. This file is part of the GNU C Library.
Contributed by Richard Henderson <rth@redhat.com>, 2005. Contributed by Richard Henderson <rth@redhat.com>, 2005.
@ -39,28 +39,29 @@ __multc3 (long double a, long double b, long double c, long double d)
{ {
/* Recover infinities that computed as NaN + iNaN. */ /* Recover infinities that computed as NaN + iNaN. */
bool recalc = 0; bool recalc = 0;
if (isinf (a) || isinf (b)) if (__isinf_nsl (a) || __isinf_nsl (b))
{ {
/* z is infinite. "Box" the infinity and change NaNs in /* z is infinite. "Box" the infinity and change NaNs in
the other factor to 0. */ the other factor to 0. */
a = __copysignl (isinf (a) ? 1 : 0, a); a = __copysignl (__isinf_nsl (a) ? 1 : 0, a);
b = __copysignl (isinf (b) ? 1 : 0, b); b = __copysignl (__isinf_nsl (b) ? 1 : 0, b);
if (isnan (c)) c = __copysignl (0, c); if (isnan (c)) c = __copysignl (0, c);
if (isnan (d)) d = __copysignl (0, d); if (isnan (d)) d = __copysignl (0, d);
recalc = 1; recalc = 1;
} }
if (isinf (c) || isinf (d)) if (__isinf_nsl (c) || __isinf_nsl (d))
{ {
/* w is infinite. "Box" the infinity and change NaNs in /* w is infinite. "Box" the infinity and change NaNs in
the other factor to 0. */ the other factor to 0. */
c = __copysignl (isinf (c) ? 1 : 0, c); c = __copysignl (__isinf_nsl (c) ? 1 : 0, c);
d = __copysignl (isinf (d) ? 1 : 0, d); d = __copysignl (__isinf_nsl (d) ? 1 : 0, d);
if (isnan (a)) a = __copysignl (0, a); if (isnan (a)) a = __copysignl (0, a);
if (isnan (b)) b = __copysignl (0, b); if (isnan (b)) b = __copysignl (0, b);
recalc = 1; recalc = 1;
} }
if (!recalc if (!recalc
&& (isinf (ac) || isinf (bd) || isinf (ad) || isinf (bc))) && (__isinf_nsl (ac) || __isinf_nsl (bd)
|| __isinf_nsl (ad) || __isinf_nsl (bc)))
{ {
/* Recover infinities from overflow by changing NaNs to 0. */ /* Recover infinities from overflow by changing NaNs to 0. */
if (isnan (a)) a = __copysignl (0, a); if (isnan (a)) a = __copysignl (0, a);

View File

@ -1,5 +1,5 @@
/* Return arc sine of complex double value. /* Return arc sine of complex double value.
Copyright (C) 1997 Free Software Foundation, Inc. Copyright (C) 1997, 2011 Free Software Foundation, Inc.
This file is part of the GNU C Library. This file is part of the GNU C Library.
Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997. Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
@ -20,6 +20,7 @@
#include <complex.h> #include <complex.h>
#include <math.h> #include <math.h>
#include <math_private.h>
__complex__ double __complex__ double
@ -33,7 +34,7 @@ __casin (__complex__ double x)
{ {
res = x; res = x;
} }
else if (__isinf (__real__ x) || __isinf (__imag__ x)) else if (__isinf_ns (__real__ x) || __isinf_ns (__imag__ x))
{ {
__real__ res = __nan (""); __real__ res = __nan ("");
__imag__ res = __copysign (HUGE_VAL, __imag__ x); __imag__ res = __copysign (HUGE_VAL, __imag__ x);

View File

@ -1,5 +1,5 @@
/* Return arc sine of complex float value. /* Return arc sine of complex float value.
Copyright (C) 1997 Free Software Foundation, Inc. Copyright (C) 1997, 2011 Free Software Foundation, Inc.
This file is part of the GNU C Library. This file is part of the GNU C Library.
Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997. Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
@ -20,6 +20,7 @@
#include <complex.h> #include <complex.h>
#include <math.h> #include <math.h>
#include <math_private.h>
__complex__ float __complex__ float
@ -33,7 +34,7 @@ __casinf (__complex__ float x)
{ {
res = x; res = x;
} }
else if (__isinff (__real__ x) || __isinff (__imag__ x)) else if (__isinf_nsf (__real__ x) || __isinf_nsf (__imag__ x))
{ {
__real__ res = __nanf (""); __real__ res = __nanf ("");
__imag__ res = __copysignf (HUGE_VALF, __imag__ x); __imag__ res = __copysignf (HUGE_VALF, __imag__ x);

View File

@ -1,5 +1,5 @@
/* Return arc sine of complex long double value. /* Return arc sine of complex long double value.
Copyright (C) 1997 Free Software Foundation, Inc. Copyright (C) 1997, 2011 Free Software Foundation, Inc.
This file is part of the GNU C Library. This file is part of the GNU C Library.
Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997. Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
@ -20,6 +20,7 @@
#include <complex.h> #include <complex.h>
#include <math.h> #include <math.h>
#include <math_private.h>
__complex__ long double __complex__ long double
@ -33,7 +34,7 @@ __casinl (__complex__ long double x)
{ {
res = x; res = x;
} }
else if (__isinfl (__real__ x) || __isinfl (__imag__ x)) else if (__isinf_nsl (__real__ x) || __isinf_nsl (__imag__ x))
{ {
__real__ res = __nanl (""); __real__ res = __nanl ("");
__imag__ res = __copysignl (HUGE_VALL, __imag__ x); __imag__ res = __copysignl (HUGE_VALL, __imag__ x);

View File

@ -1,5 +1,5 @@
/* Return cosine of complex double value. /* Return cosine of complex double value.
Copyright (C) 1997 Free Software Foundation, Inc. Copyright (C) 1997, 2011 Free Software Foundation, Inc.
This file is part of the GNU C Library. This file is part of the GNU C Library.
Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997. Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
@ -21,6 +21,7 @@
#include <complex.h> #include <complex.h>
#include <fenv.h> #include <fenv.h>
#include <math.h> #include <math.h>
#include <math_private.h>
__complex__ double __complex__ double
@ -36,17 +37,17 @@ __ccos (__complex__ double x)
__imag__ res = 0.0; __imag__ res = 0.0;
#ifdef FE_INVALID #ifdef FE_INVALID
if (__isinf (__real__ x)) if (__isinf_ns (__real__ x))
feraiseexcept (FE_INVALID); feraiseexcept (FE_INVALID);
#endif #endif
} }
else if (__isinf (__imag__ x)) else if (__isinf_ns (__imag__ x))
{ {
__real__ res = HUGE_VAL; __real__ res = HUGE_VAL;
__imag__ res = __nan (""); __imag__ res = __nan ("");
#ifdef FE_INVALID #ifdef FE_INVALID
if (__isinf (__real__ x)) if (__isinf_ns (__real__ x))
feraiseexcept (FE_INVALID); feraiseexcept (FE_INVALID);
#endif #endif
} }

View File

@ -1,5 +1,5 @@
/* Return cosine of complex float value. /* Return cosine of complex float value.
Copyright (C) 1997 Free Software Foundation, Inc. Copyright (C) 1997, 2011 Free Software Foundation, Inc.
This file is part of the GNU C Library. This file is part of the GNU C Library.
Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997. Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
@ -21,6 +21,7 @@
#include <complex.h> #include <complex.h>
#include <fenv.h> #include <fenv.h>
#include <math.h> #include <math.h>
#include <math_private.h>
__complex__ float __complex__ float
@ -36,17 +37,17 @@ __ccosf (__complex__ float x)
__imag__ res = 0.0; __imag__ res = 0.0;
#ifdef FE_INVALID #ifdef FE_INVALID
if (__isinff (__real__ x)) if (__isinf_nsf (__real__ x))
feraiseexcept (FE_INVALID); feraiseexcept (FE_INVALID);
#endif #endif
} }
else if (__isinff (__imag__ x)) else if (__isinf_nsf (__imag__ x))
{ {
__real__ res = HUGE_VALF; __real__ res = HUGE_VALF;
__imag__ res = __nanf (""); __imag__ res = __nanf ("");
#ifdef FE_INVALID #ifdef FE_INVALID
if (__isinff (__real__ x)) if (__isinf_nsf (__real__ x))
feraiseexcept (FE_INVALID); feraiseexcept (FE_INVALID);
#endif #endif
} }

View File

@ -1,5 +1,5 @@
/* Return cosine of complex long double value. /* Return cosine of complex long double value.
Copyright (C) 1997 Free Software Foundation, Inc. Copyright (C) 1997, 2011 Free Software Foundation, Inc.
This file is part of the GNU C Library. This file is part of the GNU C Library.
Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997. Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
@ -21,6 +21,7 @@
#include <complex.h> #include <complex.h>
#include <fenv.h> #include <fenv.h>
#include <math.h> #include <math.h>
#include <math_private.h>
__complex__ long double __complex__ long double
@ -36,17 +37,17 @@ __ccosl (__complex__ long double x)
__imag__ res = 0.0; __imag__ res = 0.0;
#ifdef FE_INVALID #ifdef FE_INVALID
if (__isinfl (__real__ x)) if (__isinf_nsl (__real__ x))
feraiseexcept (FE_INVALID); feraiseexcept (FE_INVALID);
#endif #endif
} }
else if (__isinfl (__imag__ x)) else if (__isinf_nsl (__imag__ x))
{ {
__real__ res = HUGE_VALL; __real__ res = HUGE_VALL;
__imag__ res = __nanl (""); __imag__ res = __nanl ("");
#ifdef FE_INVALID #ifdef FE_INVALID
if (__isinfl (__real__ x)) if (__isinf_nsl (__real__ x))
feraiseexcept (FE_INVALID); feraiseexcept (FE_INVALID);
#endif #endif
} }

View File

@ -1,5 +1,5 @@
/* Complex tangent function for double. /* Complex tangent function for double.
Copyright (C) 1997, 2005 Free Software Foundation, Inc. Copyright (C) 1997, 2005, 2011 Free Software Foundation, Inc.
This file is part of the GNU C Library. This file is part of the GNU C Library.
Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997. Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
@ -32,7 +32,7 @@ __ctan (__complex__ double x)
if (!isfinite (__real__ x) || !isfinite (__imag__ x)) if (!isfinite (__real__ x) || !isfinite (__imag__ x))
{ {
if (__isinf (__imag__ x)) if (__isinf_ns (__imag__ x))
{ {
__real__ res = __copysign (0.0, __real__ x); __real__ res = __copysign (0.0, __real__ x);
__imag__ res = __copysign (1.0, __imag__ x); __imag__ res = __copysign (1.0, __imag__ x);
@ -47,7 +47,7 @@ __ctan (__complex__ double x)
__imag__ res = __nan (""); __imag__ res = __nan ("");
#ifdef FE_INVALID #ifdef FE_INVALID
if (__isinf (__real__ x)) if (__isinf_ns (__real__ x))
feraiseexcept (FE_INVALID); feraiseexcept (FE_INVALID);
#endif #endif
} }

View File

@ -1,5 +1,5 @@
/* Complex tangent function for float. /* Complex tangent function for float.
Copyright (C) 1997, 2005 Free Software Foundation, Inc. Copyright (C) 1997, 2005, 2011 Free Software Foundation, Inc.
This file is part of the GNU C Library. This file is part of the GNU C Library.
Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997. Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
@ -32,7 +32,7 @@ __ctanf (__complex__ float x)
if (!isfinite (__real__ x) || !isfinite (__imag__ x)) if (!isfinite (__real__ x) || !isfinite (__imag__ x))
{ {
if (__isinff (__imag__ x)) if (__isinf_nsf (__imag__ x))
{ {
__real__ res = __copysignf (0.0, __real__ x); __real__ res = __copysignf (0.0, __real__ x);
__imag__ res = __copysignf (1.0, __imag__ x); __imag__ res = __copysignf (1.0, __imag__ x);
@ -47,7 +47,7 @@ __ctanf (__complex__ float x)
__imag__ res = __nanf (""); __imag__ res = __nanf ("");
#ifdef FE_INVALID #ifdef FE_INVALID
if (__isinff (__real__ x)) if (__isinf_nsf (__real__ x))
feraiseexcept (FE_INVALID); feraiseexcept (FE_INVALID);
#endif #endif
} }

View File

@ -1,5 +1,5 @@
/* Complex hyperbole tangent for double. /* Complex hyperbole tangent for double.
Copyright (C) 1997, 2005 Free Software Foundation, Inc. Copyright (C) 1997, 2005, 2011 Free Software Foundation, Inc.
This file is part of the GNU C Library. This file is part of the GNU C Library.
Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997. Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
@ -32,7 +32,7 @@ __ctanh (__complex__ double x)
if (!isfinite (__real__ x) || !isfinite (__imag__ x)) if (!isfinite (__real__ x) || !isfinite (__imag__ x))
{ {
if (__isinf (__real__ x)) if (__isinf_ns (__real__ x))
{ {
__real__ res = __copysign (1.0, __real__ x); __real__ res = __copysign (1.0, __real__ x);
__imag__ res = __copysign (0.0, __imag__ x); __imag__ res = __copysign (0.0, __imag__ x);
@ -47,7 +47,7 @@ __ctanh (__complex__ double x)
__imag__ res = __nan (""); __imag__ res = __nan ("");
#ifdef FE_INVALID #ifdef FE_INVALID
if (__isinf (__imag__ x)) if (__isinf_ns (__imag__ x))
feraiseexcept (FE_INVALID); feraiseexcept (FE_INVALID);
#endif #endif
} }

View File

@ -1,5 +1,5 @@
/* Complex hyperbole tangent for float. /* Complex hyperbole tangent for float.
Copyright (C) 1997, 2005 Free Software Foundation, Inc. Copyright (C) 1997, 2005, 2011 Free Software Foundation, Inc.
This file is part of the GNU C Library. This file is part of the GNU C Library.
Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997. Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
@ -32,7 +32,7 @@ __ctanhf (__complex__ float x)
if (!isfinite (__real__ x) || !isfinite (__imag__ x)) if (!isfinite (__real__ x) || !isfinite (__imag__ x))
{ {
if (__isinff (__real__ x)) if (__isinf_nsf (__real__ x))
{ {
__real__ res = __copysignf (1.0, __real__ x); __real__ res = __copysignf (1.0, __real__ x);
__imag__ res = __copysignf (0.0, __imag__ x); __imag__ res = __copysignf (0.0, __imag__ x);
@ -47,7 +47,7 @@ __ctanhf (__complex__ float x)
__imag__ res = __nanf (""); __imag__ res = __nanf ("");
#ifdef FE_INVALID #ifdef FE_INVALID
if (__isinff (__imag__ x)) if (__isinf_nsf (__imag__ x))
feraiseexcept (FE_INVALID); feraiseexcept (FE_INVALID);
#endif #endif
} }

View File

@ -1,5 +1,5 @@
/* Complex hyperbole tangent for long double. /* Complex hyperbole tangent for long double.
Copyright (C) 1997, 2005 Free Software Foundation, Inc. Copyright (C) 1997, 2005, 2011 Free Software Foundation, Inc.
This file is part of the GNU C Library. This file is part of the GNU C Library.
Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997. Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
@ -32,7 +32,7 @@ __ctanhl (__complex__ long double x)
if (!isfinite (__real__ x) || !isfinite (__imag__ x)) if (!isfinite (__real__ x) || !isfinite (__imag__ x))
{ {
if (__isinfl (__real__ x)) if (__isinf_nsl (__real__ x))
{ {
__real__ res = __copysignl (1.0, __real__ x); __real__ res = __copysignl (1.0, __real__ x);
__imag__ res = __copysignl (0.0, __imag__ x); __imag__ res = __copysignl (0.0, __imag__ x);
@ -47,7 +47,7 @@ __ctanhl (__complex__ long double x)
__imag__ res = __nanl (""); __imag__ res = __nanl ("");
#ifdef FE_INVALID #ifdef FE_INVALID
if (__isinfl (__imag__ x)) if (__isinf_nsl (__imag__ x))
feraiseexcept (FE_INVALID); feraiseexcept (FE_INVALID);
#endif #endif
} }

View File

@ -1,5 +1,5 @@
/* Complex tangent function for long double. /* Complex tangent function for long double.
Copyright (C) 1997, 2005 Free Software Foundation, Inc. Copyright (C) 1997, 2005, 2011 Free Software Foundation, Inc.
This file is part of the GNU C Library. This file is part of the GNU C Library.
Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997. Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
@ -32,7 +32,7 @@ __ctanl (__complex__ long double x)
if (!isfinite (__real__ x) || !isfinite (__imag__ x)) if (!isfinite (__real__ x) || !isfinite (__imag__ x))
{ {
if (__isinfl (__imag__ x)) if (__isinf_nsl (__imag__ x))
{ {
__real__ res = __copysignl (0.0, __real__ x); __real__ res = __copysignl (0.0, __real__ x);
__imag__ res = __copysignl (1.0, __imag__ x); __imag__ res = __copysignl (1.0, __imag__ x);
@ -47,7 +47,7 @@ __ctanl (__complex__ long double x)
__imag__ res = __nanl (""); __imag__ res = __nanl ("");
#ifdef FE_INVALID #ifdef FE_INVALID
if (__isinfl (__real__ x)) if (__isinf_nsl (__real__ x))
feraiseexcept (FE_INVALID); feraiseexcept (FE_INVALID);
#endif #endif
} }

View File

@ -35,9 +35,9 @@ static char rcsid[] = "$NetBSD: w_fmod.c,v 1.6 1995/05/10 20:48:55 jtc Exp $";
double z; double z;
z = __ieee754_fmod(x,y); z = __ieee754_fmod(x,y);
if(_LIB_VERSION == _IEEE_ ||__isnan(y)||__isnan(x)) return z; if(_LIB_VERSION == _IEEE_ ||__isnan(y)||__isnan(x)) return z;
if(__isinf(x)||y==0.0) { if(__isinf_ns(x)||y==0.0) {
/* fmod(+-Inf,y) or fmod(x,0) */ /* fmod(+-Inf,y) or fmod(x,0) */
return __kernel_standard(x,y,27); return __kernel_standard(x,y,27);
} else } else
return z; return z;
#endif #endif

View File

@ -38,9 +38,9 @@ static char rcsid[] = "$NetBSD: w_fmodf.c,v 1.3 1995/05/10 20:48:57 jtc Exp $";
float z; float z;
z = __ieee754_fmodf(x,y); z = __ieee754_fmodf(x,y);
if(_LIB_VERSION == _IEEE_ ||__isnanf(y)||__isnanf(x)) return z; if(_LIB_VERSION == _IEEE_ ||__isnanf(y)||__isnanf(x)) return z;
if(__isinff(x)||y==(float)0.0) { if(__isinf_nsf(x)||y==(float)0.0) {
/* fmodf(+-Inf,y) or fmodf(x,0) */ /* fmodf(+-Inf,y) or fmodf(x,0) */
return (float)__kernel_standard((double)x,(double)y,127); return (float)__kernel_standard((double)x,(double)y,127);
} else } else
return z; return z;
#endif #endif

View File

@ -39,9 +39,9 @@ static char rcsid[] = "$NetBSD: $";
long double z; long double z;
z = __ieee754_fmodl(x,y); z = __ieee754_fmodl(x,y);
if(_LIB_VERSION == _IEEE_ ||__isnanl(y)||__isnanl(x)) return z; if(_LIB_VERSION == _IEEE_ ||__isnanl(y)||__isnanl(x)) return z;
if(__isinfl(x)||y==0.0) { if(__isinf_nsl(x)||y==0.0) {
/* fmodl(+-Inf,y) or fmodl(x,0) */ /* fmodl(+-Inf,y) or fmodl(x,0) */
return __kernel_standard(x,y,227); return __kernel_standard(x,y,227);
} else } else
return z; return z;
#endif #endif

View File

@ -34,7 +34,7 @@ static char rcsid[] = "$NetBSD: w_remainder.c,v 1.6 1995/05/10 20:49:44 jtc Exp
double z; double z;
z = __ieee754_remainder(x,y); z = __ieee754_remainder(x,y);
if(_LIB_VERSION == _IEEE_ || __isnan(y) || __isnan(x)) return z; if(_LIB_VERSION == _IEEE_ || __isnan(y) || __isnan(x)) return z;
if(y==0.0 || __isinf(x)) if(y==0.0 || __isinf_ns(x))
return __kernel_standard(x,y,28); /* remainder(x,0) */ return __kernel_standard(x,y,28); /* remainder(x,0) */
else else
return z; return z;

View File

@ -37,7 +37,7 @@ static char rcsid[] = "$NetBSD: w_remainderf.c,v 1.3 1995/05/10 20:49:46 jtc Exp
float z; float z;
z = __ieee754_remainderf(x,y); z = __ieee754_remainderf(x,y);
if(_LIB_VERSION == _IEEE_ || __isnanf(y) || __isnanf(x)) return z; if(_LIB_VERSION == _IEEE_ || __isnanf(y) || __isnanf(x)) return z;
if(y==(float)0.0 || __isinff(x)) if(y==(float)0.0 || __isinf_nsf(x))
/* remainder(x,0) */ /* remainder(x,0) */
return (float)__kernel_standard((double)x,(double)y,128); return (float)__kernel_standard((double)x,(double)y,128);
else else

View File

@ -39,7 +39,7 @@ static char rcsid[] = "$NetBSD: $";
long double z; long double z;
z = __ieee754_remainderl(x,y); z = __ieee754_remainderl(x,y);
if(_LIB_VERSION == _IEEE_ || __isnanl(y) || __isnanl(x)) return z; if(_LIB_VERSION == _IEEE_ || __isnanl(y) || __isnanl(x)) return z;
if(y==0.0 || __isinfl(x)) if(y==0.0 || __isinf_nsl(x))
return __kernel_standard(x,y,228); /* remainder(x,0) */ return __kernel_standard(x,y,228); /* remainder(x,0) */
else else
return z; return z;

View File

@ -22,6 +22,7 @@ static char rcsid[] = "$NetBSD: s_finite.c,v 1.8 1995/05/10 20:47:17 jtc Exp $";
#include "math.h" #include "math.h"
#include "math_private.h" #include "math_private.h"
#undef __finite
#ifdef __STDC__ #ifdef __STDC__
int __finite(double x) int __finite(double x)
#else #else

View File

@ -0,0 +1,20 @@
/*
* Written by Ulrich Drepper <drepper@gmail.com>.
*/
/*
* __isinf_ns(x) returns != 0 if x is ±inf, else 0;
* no branching!
*/
#include "math.h"
#include "math_private.h"
#undef __isinf_ns
int
__isinf_ns (double x)
{
int32_t hx,lx;
EXTRACT_WORDS(hx,lx,x);
return !(lx | ((hx & 0x7fffffff) ^ 0x7ff00000));
}

View File

@ -22,6 +22,7 @@ static char rcsid[] = "$NetBSD: s_isnan.c,v 1.8 1995/05/10 20:47:36 jtc Exp $";
#include "math.h" #include "math.h"
#include "math_private.h" #include "math_private.h"
#undef __isnan
#ifdef __STDC__ #ifdef __STDC__
int __isnan(double x) int __isnan(double x)
#else #else

View File

@ -17,6 +17,7 @@
#include "math.h" #include "math.h"
#include "math_private.h" #include "math_private.h"
#undef __finite
int int
__finite(double x) __finite(double x)
{ {

View File

@ -0,0 +1,20 @@
/*
* Written by Ulrich Drepper <drepper@gmail.com>.
*/
/*
* __isinf_ns(x) returns != 0 if x is ±inf, else 0;
* no branching!
*/
#include "math.h"
#include "math_private.h"
#undef __isinf_ns
int
__isinf_ns (double x)
{
int64_t ix;
EXTRACT_WORDS64(ix,x);
return (ix & UINT64_C(0x7fffffffffffffff)) == UINT64_C(0x7ff0000000000000);
}

View File

@ -18,6 +18,7 @@
#include "math.h" #include "math.h"
#include "math_private.h" #include "math_private.h"
#undef __isnan
#ifdef __STDC__ #ifdef __STDC__
int __isnan(double x) int __isnan(double x)
#else #else

View File

@ -25,6 +25,7 @@ static char rcsid[] = "$NetBSD: s_finitef.c,v 1.4 1995/05/10 20:47:18 jtc Exp $"
#include "math.h" #include "math.h"
#include "math_private.h" #include "math_private.h"
#undef __finitef
#ifdef __STDC__ #ifdef __STDC__
int __finitef(float x) int __finitef(float x)
#else #else

View File

@ -0,0 +1,20 @@
/*
* Written by Ulrich Drepper <drepper@gmail.com>.
*/
/*
* __isinf_nsf(x) returns != 0 if x is ±inf, else 0;
* no branching!
*/
#include "math.h"
#include "math_private.h"
#undef __isinf_nsf
int
__isinf_nsf (float x)
{
int32_t ix,t;
GET_FLOAT_WORD(ix,x);
return (ix & 0x7fffffff) == 0x7f800000;
}

View File

@ -4,7 +4,7 @@
/* /*
* ==================================================== * ====================================================
* Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. * Copyright (C) 1993, 2011 by Sun Microsystems, Inc. All rights reserved.
* *
* Developed at SunPro, a Sun Microsystems, Inc. business. * Developed at SunPro, a Sun Microsystems, Inc. business.
* Permission to use, copy, modify, and distribute this * Permission to use, copy, modify, and distribute this
@ -25,6 +25,7 @@ static char rcsid[] = "$NetBSD: s_isnanf.c,v 1.4 1995/05/10 20:47:38 jtc Exp $";
#include "math.h" #include "math.h"
#include "math_private.h" #include "math_private.h"
#undef __isnanf
#ifdef __STDC__ #ifdef __STDC__
int __isnanf(float x) int __isnanf(float x)
#else #else

View File

@ -0,0 +1,19 @@
/*
* Written by Ulrich Drepper <drepper@gmail.com>
*/
/*
* __isinf_nsl(x) returns != 0 if x is ±inf, else 0;
* no branching!
*/
#include "math.h"
#include "math_private.h"
int
__isinf_nsl (long double x)
{
int64_t hx,lx;
GET_LDOUBLE_WORDS64(hx,lx,x);
return !(lx | ((hx & 0x7fffffffffffffffLL) ^ 0x7fff000000000000LL));
}

View File

@ -0,0 +1,18 @@
/*
* Written by Ulrich Drepper <drepper@gmail.com>.
*/
/*
* __isinf_nsl(x) returns != 0 if x is ±inf, else 0;
*/
#include "math.h"
#include "math_private.h"
int
__isinf_nsl (long double x)
{
int32_t se,hx,lx;
GET_LDOUBLE_WORDS(se,hx,lx,x);
return !(((se & 0x7fff) ^ 0x7fff) | lx | (hx & 0x7fffffff));
}

View File

@ -56,3 +56,24 @@ do { \
} while (0) } while (0)
#endif #endif
#define __isnan(d) \
({ long int __di; EXTRACT_WORDS64 (__di, d); \
(__di & 0x7fffffffffffffffl) > 0x7ff0000000000000l; })
#define __isnanf(d) \
({ int __di; GET_FLOAT_WORD (__di, d); \
(__di & 0x7fffffff) > 0x7f800000; })
#define __isinf_ns(d) \
({ long int __di; EXTRACT_WORDS64 (__di, d); \
(__di & 0x7fffffffffffffffl) == 0x7ff0000000000000l; })
#define __isinf_nsf(d) \
({ int __di; GET_FLOAT_WORD (__di, d); \
(__di & 0x7fffffff) == 0x7f800000; })
#define __finite(d) \
({ long int __di; EXTRACT_WORDS64 (__di, d); \
(__di & 0x7fffffffffffffffl) < 0x7ff0000000000000l; })
#define __finitef(d) \
({ int __di; GET_FLOAT_WORD (__di, d); \
(__di & 0x7fffffff) < 0x7f800000; })