mirror of
https://sourceware.org/git/glibc.git
synced 2024-11-24 02:03:35 +08:00
9946e7a949
The wrapper implementations of ldexp / scalbn / scalbln (architecture-independent), and their float / long double variants, return sNaN for sNaN input. This patch fixes them to add relevant arguments to themselves so that qNaN is returned in this case. Tested for x86_64 and x86. [BZ #20225] * math/s_ldexp.c (__ldexp): Add non-finite or zero argument to itself. * math/s_ldexpf.c (__ldexpf): Likewise. * math/s_ldexpl.c (__ldexpl): Likewise. * math/w_scalbln.c (__w_scalbln): Likewise. * math/w_scalblnf.c (__w_scalblnf): Likewise. * math/w_scalblnl.c (__w_scalblnl): Likewise. * math/libm-test.inc (scalbn_test_data): Add sNaN tests. (scalbln_test_data): Likewise.
35 lines
984 B
C
35 lines
984 B
C
/* @(#)s_ldexp.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.
|
|
* ====================================================
|
|
*/
|
|
|
|
#if defined(LIBM_SCCS) && !defined(lint)
|
|
static char rcsid[] = "$NetBSD: s_ldexp.c,v 1.6 1995/05/10 20:47:40 jtc Exp $";
|
|
#endif
|
|
|
|
#include <math.h>
|
|
#include <math_private.h>
|
|
#include <errno.h>
|
|
|
|
double __ldexp(double value, int exp)
|
|
{
|
|
if(!isfinite(value)||value==0.0) return value + value;
|
|
value = __scalbn(value,exp);
|
|
if(!isfinite(value)||value==0.0) __set_errno (ERANGE);
|
|
return value;
|
|
}
|
|
weak_alias (__ldexp, ldexp)
|
|
weak_alias (__ldexp, scalbn)
|
|
#ifdef NO_LONG_DOUBLE
|
|
strong_alias (__ldexp, __ldexpl)
|
|
weak_alias (__ldexp, ldexpl)
|
|
weak_alias (__ldexp, scalbnl)
|
|
#endif
|