From: Ulrich Drepper Date: Wed, 26 Aug 2009 07:12:58 +0000 (-0700) Subject: Optimized isnan for 64-bit machines. X-Git-Tag: fedora/glibc-2.10.90-16~1^2~1 X-Git-Url: http://git.ipfire.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=15e7f0a44d6f17eb445f7d5f02ace7ec66306d7a;p=thirdparty%2Fglibc.git Optimized isnan for 64-bit machines. --- diff --git a/ChangeLog b/ChangeLog index e56e139dfd5..c28272c7357 100644 --- a/ChangeLog +++ b/ChangeLog @@ -1,3 +1,7 @@ +2009-08-26 Ulrich Drepper + + * sysdeps/ieee754/dbl-64/wordsize-64/s_isnan.c: New file. + 2009-08-25 Ulrich Drepper * sysdeps/ieee754/dbl-64/wordsize-64/s_lround.c: Add llround aliases. diff --git a/sysdeps/ieee754/dbl-64/wordsize-64/s_isnan.c b/sysdeps/ieee754/dbl-64/wordsize-64/s_isnan.c new file mode 100644 index 00000000000..3b08c54dd45 --- /dev/null +++ b/sysdeps/ieee754/dbl-64/wordsize-64/s_isnan.c @@ -0,0 +1,39 @@ +/* @(#)s_isnan.c 5.1 93/09/24 */ +/* + * ==================================================== + * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. + * + * Developed at SunPro, a Sun Microsystems, Inc. business. + * Permission to use, copy, modify, and distribute this + * software is freely granted, provided that this notice + * is preserved. + * ==================================================== + */ + +/* + * isnan(x) returns 1 is x is nan, else 0; + * no branching! + */ + +#include "math.h" +#include "math_private.h" + +#ifdef __STDC__ + int __isnan(double x) +#else + int __isnan(x) + double x; +#endif +{ + int64_t hx; + EXTRACT_WORDS64(hx,x); + hx &= UINT64_C(0x7fffffffffffffff); + hx = UINT64_C(0x7ff0000000000000) - hx; + return (int)(((uint64_t)hx)>>63); +} +hidden_def (__isnan) +weak_alias (__isnan, isnan) +#ifdef NO_LONG_DOUBLE +strong_alias (__isnan, __isnanl) +weak_alias (__isnan, isnanl) +#endif