Favour normal numbers

This commit is contained in:
Siddhesh Poyarekar 2012-12-31 16:05:13 +05:30
parent 9c89fca6e0
commit 18ea052c3e
2 changed files with 10 additions and 4 deletions

View File

@ -1,3 +1,8 @@
2012-12-31 Siddhesh Poyarekar <siddhesh@redhat.com>
* sysdeps/ieee754/dbl-64/mpa.c (__mp_dbl): Favour normal
numbers.
2012-12-30 Mike Frysinger <vapier@gentoo.org>
* math/bits/mathcalls.h (modf): Use __nonnull.

View File

@ -1,7 +1,7 @@
/*
* IBM Accurate Mathematical Library
* written by International Business Machines Corp.
* Copyright (C) 2001, 2011 Free Software Foundation
* Copyright (C) 2001-2012 Free Software Foundation
*
* This program is free software; you can redistribute it and/or modify
* it under the terms of the GNU Lesser General Public License as published by
@ -247,9 +247,10 @@ void __mp_dbl(const mp_no *x, double *y, int p) {
if (X[0] == ZERO) {*y = ZERO; return; }
if (EX> -42) norm(x,y,p);
else if (EX==-42 && X[1]>=TWO10) norm(x,y,p);
else denorm(x,y,p);
if (__glibc_likely (EX > -42 || (EX == -42 && X[1] >= TWO10)))
norm(x,y,p);
else
denorm(x,y,p);
}
#endif