about summary refs log tree commit diff
path: root/sysdeps/ieee754/dbl-64/mpa.c
diff options
context:
space:
mode:
authorSiddhesh Poyarekar <siddhesh@redhat.com>2012-12-31 16:05:13 +0530
committerSiddhesh Poyarekar <siddhesh@redhat.com>2012-12-31 16:05:13 +0530
commit18ea052c3e3fa91cf8a80e20c48023fd0ba5fab8 (patch)
tree3ec6ec05dba547cfeeb36b60a6f37b2e93ce90e2 /sysdeps/ieee754/dbl-64/mpa.c
parent9c89fca6e0bd8ecf72c793965b8bdf51b883a64d (diff)
downloadglibc-18ea052c3e3fa91cf8a80e20c48023fd0ba5fab8.tar.gz
glibc-18ea052c3e3fa91cf8a80e20c48023fd0ba5fab8.tar.xz
glibc-18ea052c3e3fa91cf8a80e20c48023fd0ba5fab8.zip
Favour normal numbers
Diffstat (limited to 'sysdeps/ieee754/dbl-64/mpa.c')
-rw-r--r--sysdeps/ieee754/dbl-64/mpa.c9
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