about summary refs log tree commit diff
path: root/sysdeps/ieee754/ldbl-128ibm/s_fmal.c
diff options
context:
space:
mode:
authorPaul E. Murphy <murphyp@linux.vnet.ibm.com>2016-08-08 15:58:28 -0500
committerPaul E. Murphy <murphyp@linux.vnet.ibm.com>2016-08-19 11:29:43 -0500
commit4482ff226e4b286ab171f3c5841ae1f7e61780cd (patch)
tree69bf06744f290b3656ad238c0b5dd700c850f9e9 /sysdeps/ieee754/ldbl-128ibm/s_fmal.c
parent01ee387015a2075c45a4e1ad45d39e50b5a6d40b (diff)
downloadglibc-4482ff226e4b286ab171f3c5841ae1f7e61780cd.tar.gz
glibc-4482ff226e4b286ab171f3c5841ae1f7e61780cd.tar.xz
glibc-4482ff226e4b286ab171f3c5841ae1f7e61780cd.zip
Merge common usage of mul_split function
A number of files share identical code for the
mul_split function.

This moves the duplicated function mul_split into its
own header, and refactors the fma usage into a single
selection macro.  Likewise, mul_split when used by a
long double implementation is renamed mul_splitl for
clarity.
Diffstat (limited to 'sysdeps/ieee754/ldbl-128ibm/s_fmal.c')
-rw-r--r--sysdeps/ieee754/ldbl-128ibm/s_fmal.c27
1 files changed, 1 insertions, 26 deletions
diff --git a/sysdeps/ieee754/ldbl-128ibm/s_fmal.c b/sysdeps/ieee754/ldbl-128ibm/s_fmal.c
index 177a04817b..140576315a 100644
--- a/sysdeps/ieee754/ldbl-128ibm/s_fmal.c
+++ b/sysdeps/ieee754/ldbl-128ibm/s_fmal.c
@@ -22,6 +22,7 @@
 #include <math.h>
 #include <math_private.h>
 #include <math_ldbl_opt.h>
+#include <mul_split.h>
 #include <stdlib.h>
 
 /* Calculate X + Y exactly and store the result in *HI + *LO.  It is
@@ -36,32 +37,6 @@ add_split (double *hi, double *lo, double x, double y)
   *lo = (x - *hi) + y;
 }
 
-/* Calculate X * Y exactly and store the result in *HI + *LO.  It is
-   given that the values are small enough that no overflow occurs and
-   large enough (or zero) that no underflow occurs.  */
-
-static void
-mul_split (double *hi, double *lo, double x, double y)
-{
-#ifdef __FP_FAST_FMA
-  /* Fast built-in fused multiply-add.  */
-  *hi = x * y;
-  *lo = __builtin_fma (x, y, -*hi);
-#else
-  /* Apply Dekker's algorithm.  */
-  *hi = x * y;
-# define C ((1 << (DBL_MANT_DIG + 1) / 2) + 1)
-  double x1 = x * C;
-  double y1 = y * C;
-# undef C
-  x1 = (x - x1) + x1;
-  y1 = (y - y1) + y1;
-  double x2 = x - x1;
-  double y2 = y - y1;
-  *lo = (((x1 * y1 - *hi) + x1 * y2) + x2 * y1) + x2 * y2;
-#endif
-}
-
 /* Value with extended range, used in intermediate computations.  */
 typedef struct
 {