diff options
| author | Siddhesh Poyarekar <siddhesh@redhat.com> | 2012-12-31 16:05:13 +0530 |
|---|---|---|
| committer | Siddhesh Poyarekar <siddhesh@redhat.com> | 2012-12-31 16:05:13 +0530 |
| commit | 18ea052c3e3fa91cf8a80e20c48023fd0ba5fab8 (patch) | |
| tree | 3ec6ec05dba547cfeeb36b60a6f37b2e93ce90e2 /sysdeps | |
| parent | 9c89fca6e0bd8ecf72c793965b8bdf51b883a64d (diff) | |
Favour normal numbers
Diffstat (limited to 'sysdeps')
| -rw-r--r-- | sysdeps/ieee754/dbl-64/mpa.c | 9 |
1 files changed, 5 insertions, 4 deletions
diff --git a/sysdeps/ieee754/dbl-64/mpa.c b/sysdeps/ieee754/dbl-64/mpa.c index cf4f6f7c53..9ab0a48c11 100644 --- a/sysdeps/ieee754/dbl-64/mpa.c +++ b/sysdeps/ieee754/dbl-64/mpa.c @@ -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 |