sysdeps/ieee754/ldbl-128/ fixes

This commit is contained in:
Jakub Jelinek 2010-01-16 17:55:37 -08:00 committed by Ulrich Drepper
parent 893549c5a0
commit 64f388285e
6 changed files with 25 additions and 15 deletions

View File

@ -1,3 +1,15 @@
2010-01-16 Jakub Jelinek <jakub@redhat.com>
* sysdeps/ieee754/ldbl-128/s_ceill.c (huge): Add L suffix to the
constant.
* sysdeps/ieee754/ldbl-128/s_floorl.c (huge): Likewise.
* sysdeps/ieee754/ldbl-128/s_expm1l.c (big): Decrease to avoid
overflow.
* sysdeps/ieee754/ldbl-128/s_log1pl.c (big): Remove unused variable.
* sysdeps/ieee754/ldbl-128/s_nexttowardf.c (__nexttowardf): Use
math_opt_barrier and math_force_eval macros.
2010-01-15 Ulrich Drepper <drepper@redhat.com> 2010-01-15 Ulrich Drepper <drepper@redhat.com>
* posix/regcomp.c: Fix a few more cases of ignored return values. * posix/regcomp.c: Fix a few more cases of ignored return values.

View File

@ -30,9 +30,9 @@ static char rcsid[] = "$NetBSD: $";
#include "math_private.h" #include "math_private.h"
#ifdef __STDC__ #ifdef __STDC__
static const long double huge = 1.0e4930; static const long double huge = 1.0e4930L;
#else #else
static long double huge = 1.0e4930; static long double huge = 1.0e4930L;
#endif #endif
#ifdef __STDC__ #ifdef __STDC__

View File

@ -86,7 +86,7 @@ static const long double
/* ln (2^16384 * (1 - 2^-113)) */ /* ln (2^16384 * (1 - 2^-113)) */
maxlog = 1.1356523406294143949491931077970764891253E4L, maxlog = 1.1356523406294143949491931077970764891253E4L,
/* ln 2^-114 */ /* ln 2^-114 */
minarg = -7.9018778583833765273564461846232128760607E1L, big = 2e4932L; minarg = -7.9018778583833765273564461846232128760607E1L, big = 1e4932L;
long double long double

View File

@ -30,9 +30,9 @@ static char rcsid[] = "$NetBSD: $";
#include "math_private.h" #include "math_private.h"
#ifdef __STDC__ #ifdef __STDC__
static const long double huge = 1.0e4930; static const long double huge = 1.0e4930L;
#else #else
static long double huge = 1.0e4930; static long double huge = 1.0e4930L;
#endif #endif
#ifdef __STDC__ #ifdef __STDC__

View File

@ -117,7 +117,6 @@ static const long double C2 = 1.428606820309417232121458176568075500134E-6L;
static const long double sqrth = 0.7071067811865475244008443621048490392848L; static const long double sqrth = 0.7071067811865475244008443621048490392848L;
/* ln (2^16384 * (1 - 2^-113)) */ /* ln (2^16384 * (1 - 2^-113)) */
static const long double maxlog = 1.1356523406294143949491931077970764891253E4L; static const long double maxlog = 1.1356523406294143949491931077970764891253E4L;
static const long double big = 2e4932L;
static const long double zero = 0.0L; static const long double zero = 0.0L;
long double long double

View File

@ -44,10 +44,12 @@ static char rcsid[] = "$NetBSD: $";
return x+y; return x+y;
if((long double) x==y) return y; /* x=y, return y */ if((long double) x==y) return y; /* x=y, return y */
if(ix==0) { /* x == 0 */ if(ix==0) { /* x == 0 */
float x2; float u;
SET_FLOAT_WORD(x,(u_int32_t)((hy>>32)&0x80000000)|1);/* return +-minsub*/ SET_FLOAT_WORD(x,(u_int32_t)((hy>>32)&0x80000000)|1);/* return +-minsub*/
x2 = x*x; u = math_opt_barrier (x);
if(x2==x) return x2; else return x; /* raise underflow flag */ u = u * u;
math_force_eval (u); /* raise underflow flag */
return x;
} }
if(hx>=0) { /* x > 0 */ if(hx>=0) { /* x > 0 */
if(hy<0||(ix>>23)>(iy>>48)-0x3f80 if(hy<0||(ix>>23)>(iy>>48)-0x3f80
@ -68,12 +70,9 @@ static char rcsid[] = "$NetBSD: $";
} }
hy = hx&0x7f800000; hy = hx&0x7f800000;
if(hy>=0x7f800000) return x+x; /* overflow */ if(hy>=0x7f800000) return x+x; /* overflow */
if(hy<0x00800000) { /* underflow */ if(hy<0x00800000) {
float x2 = x*x; float u = x*x; /* underflow */
if(x2!=x) { /* raise underflow flag */ math_force_eval (u); /* raise underflow flag */
SET_FLOAT_WORD(x2,hx);
return x2;
}
} }
SET_FLOAT_WORD(x,hx); SET_FLOAT_WORD(x,hx);
return x; return x;