Revert exp reimplementation (causes test failures).

	Revert:

	2017-12-19  Joseph Myers  <joseph@codesourcery.com>

	* sysdeps/x86_64/fpu/libm-test-ulps: Update.

	2017-12-19  Patrick McGehearty  <patrick.mcgehearty@oracle.com>

	* sysdeps/ieee754/dbl-64/e_exp.c: Include <math-svid-compat.h> and
	<errno.h>.  Include "eexp.tbl".
	(half): New constant.
	(one): Likewise.
	(__ieee754_exp): Rewrite.
	(__slowexp): Remove prototype.
	* sysdeps/ieee754/dbl-64/eexp.tbl: New file.
	* sysdeps/ieee754/dbl-64/slowexp.c: Remove file.
	* sysdeps/i386/fpu/slowexp.c: Likewise.
	* sysdeps/ia64/fpu/slowexp.c: Likewise.
	* sysdeps/m68k/m680x0/fpu/slowexp.c: Likewise.
	* sysdeps/x86_64/fpu/multiarch/slowexp-avx.c: Likewise.
	* sysdeps/x86_64/fpu/multiarch/slowexp-fma.c: Likewise.
	* sysdeps/x86_64/fpu/multiarch/slowexp-fma4.c: Likewise.
	* sysdeps/generic/math_private.h (__slowexp): Remove prototype.
	* sysdeps/ieee754/dbl-64/e_pow.c: Remove mention of slowexp.c in
	comment.
	* sysdeps/powerpc/power4/fpu/Makefile [$(subdir) = math]
	(CPPFLAGS-slowexp.c): Remove variable.
	* sysdeps/x86_64/fpu/multiarch/Makefile (libm-sysdep_routines):
	Remove slowexp-fma, slowexp-fma4 and slowexp-avx.
	(CFLAGS-slowexp-fma.c): Remove variable.
	(CFLAGS-slowexp-fma4.c): Likewise.
	(CFLAGS-slowexp-avx.c): Likewise.
	* sysdeps/x86_64/fpu/multiarch/e_exp-avx.c (__slowexp): Do not
	define as macro.
	* sysdeps/x86_64/fpu/multiarch/e_exp-fma.c (__slowexp): Likewise.
	* sysdeps/x86_64/fpu/multiarch/e_exp-fma4.c (__slowexp): Likewise.
	* math/Makefile (type-double-routines): Remove slowexp.
	* manual/probes.texi (slowexp_p6): Remove.
	(slowexp_p32): Likewise.

17 files changed, 319 insertions, 484 deletions

diff --git a/sysdeps/generic/math_private.h b/sysdeps/generic/math_private.h
index 689dc54767..f29898c19c 100644
--- a/sysdeps/generic/math_private.h
+++ b/sysdeps/generic/math_private.h
@@ -262,6 +262,7 @@ extern double __sin32 (double __x, double __res, double __res1);
 extern double __cos32 (double __x, double __res, double __res1);
 extern double __mpsin (double __x, double __dx, bool __range_reduce);
 extern double __mpcos (double __x, double __dx, bool __range_reduce);
+extern double __slowexp (double __x);
 extern double __slowpow (double __x, double __y, double __z);
 extern void __docos (double __x, double __dx, double __v[]);
 
diff --git a/sysdeps/i386/fpu/slowexp.c b/sysdeps/i386/fpu/slowexp.c
new file mode 100644
index 0000000000..1cc8931700
--- /dev/null
+++ b/sysdeps/i386/fpu/slowexp.c
@@ -0,0 +1 @@
+/* Not needed.  */
diff --git a/sysdeps/ia64/fpu/slowexp.c b/sysdeps/ia64/fpu/slowexp.c
new file mode 100644
index 0000000000..1cc8931700
--- /dev/null
+++ b/sysdeps/ia64/fpu/slowexp.c
@@ -0,0 +1 @@
+/* Not needed.  */
diff --git a/sysdeps/ieee754/dbl-64/e_exp.c b/sysdeps/ieee754/dbl-64/e_exp.c
index 6a7122f585..6757a14ce1 100644
--- a/sysdeps/ieee754/dbl-64/e_exp.c
+++ b/sysdeps/ieee754/dbl-64/e_exp.c
@@ -1,4 +1,3 @@
-/* EXP function - Compute double precision exponential */
 /*
  * IBM Accurate Mathematical Library
  * written by International Business Machines Corp.
@@ -24,7 +23,7 @@
 /*           exp1                                                          */
 /*                                                                         */
 /* FILES NEEDED:dla.h endian.h mpa.h mydefs.h uexp.h                       */
-/*              mpa.c mpexp.x                                              */
+/*              mpa.c mpexp.x slowexp.c                                    */
 /*                                                                         */
 /* An ultimate exp routine. Given an IEEE double machine number x          */
 /* it computes the correctly rounded (to nearest) value of e^x             */
@@ -33,238 +32,207 @@
 /*                                                                         */
 /***************************************************************************/
 
-/*  IBM exp(x) replaced by following exp(x) in 2017. IBM exp1(x,xx) remains.  */
-/* exp(x)
-   Hybrid algorithm of Peter Tang's Table driven method (for large
-   arguments) and an accurate table (for small arguments).
-   Written by K.C. Ng, November 1988.
-   Revised by Patrick McGehearty, Nov 2017 to use j/64 instead of j/32
-   Method (large arguments):
-	1. Argument Reduction: given the input x, find r and integer k
-	   and j such that
-	             x = (k+j/64)*(ln2) + r,  |r| <= (1/128)*ln2
-
-	2. exp(x) = 2^k * (2^(j/64) + 2^(j/64)*expm1(r))
-	   a. expm1(r) is approximated by a polynomial:
-	      expm1(r) ~ r + t1*r^2 + t2*r^3 + ... + t5*r^6
-	      Here t1 = 1/2 exactly.
-	   b. 2^(j/64) is represented to twice double precision
-	      as TBL[2j]+TBL[2j+1].
-
-   Note: If divide were fast enough, we could use another approximation
-	 in 2.a:
-	      expm1(r) ~ (2r)/(2-R), R = r - r^2*(t1 + t2*r^2)
-	      (for the same t1 and t2 as above)
-
-   Special cases:
-	exp(INF) is INF, exp(NaN) is NaN;
-	exp(-INF)=  0;
-	for finite argument, only exp(0)=1 is exact.
-
-   Accuracy:
-	According to an error analysis, the error is always less than
-	an ulp (unit in the last place).  The largest errors observed
-	are less than 0.55 ulp for normal results and less than 0.75 ulp
-	for subnormal results.
-
-   Misc. info.
-	For IEEE double
-		if x >  7.09782712893383973096e+02 then exp(x) overflow
-		if x < -7.45133219101941108420e+02 then exp(x) underflow.  */
-
 #include <math.h>
-#include <math-svid-compat.h>
-#include <math_private.h>
-#include <errno.h>
 #include "endian.h"
 #include "uexp.h"
-#include "uexp.tbl"
 #include "mydefs.h"
 #include "MathLib.h"
+#include "uexp.tbl"
+#include <math_private.h>
 #include <fenv.h>
 #include <float.h>
 
-extern double __ieee754_exp (double);
-
-#include "eexp.tbl"
-
-static const double
-  half = 0.5,
-  one = 1.0;
+#ifndef SECTION
+# define SECTION
+#endif
 
+double __slowexp (double);
 
+/* An ultimate exp routine. Given an IEEE double machine number x it computes
+   the correctly rounded (to nearest) value of e^x.  */
 double
-__ieee754_exp (double x_arg)
+SECTION
+__ieee754_exp (double x)
 {
-  double z, t;
+  double bexp, t, eps, del, base, y, al, bet, res, rem, cor;
+  mynumber junk1, junk2, binexp = {{0, 0}};
+  int4 i, j, m, n, ex;
   double retval;
-  int hx, ix, k, j, m;
-  int fe_val;
-  union
-  {
-    int i_part[2];
-    double x;
-  } xx;
-  union
-  {
-    int y_part[2];
-    double y;
-  } yy;
-  xx.x = x_arg;
-
-  ix = xx.i_part[HIGH_HALF];
-  hx = ix & ~0x80000000;
-
-  if (hx < 0x3ff0a2b2)
-    {				/* |x| < 3/2 ln 2 */
-      if (hx < 0x3f862e42)
-	{			/* |x| < 1/64 ln 2 */
-	  if (hx < 0x3ed00000)
-	    {			/* |x| < 2^-18 */
-	      if (hx < 0x3e300000)
-		{
-		  retval = one + xx.x;
-		  return retval;
-		}
-	      retval = one + xx.x * (one + half * xx.x);
-	      return retval;
-	    }
-	  /* Use FE_TONEAREST rounding mode for computing yy.y.
-	     Avoid set/reset of rounding mode if in FE_TONEAREST mode.  */
-	  fe_val = get_rounding_mode ();
-	  if (fe_val == FE_TONEAREST)
-	    {
-	      t = xx.x * xx.x;
-	      yy.y = xx.x + (t * (half + xx.x * t2)
-			     + (t * t) * (t3 + xx.x * t4 + t * t5));
-	      retval = one + yy.y;
-	    }
-	  else
-	    {
-	      libc_fesetround (FE_TONEAREST);
-	      t = xx.x * xx.x;
-	      yy.y = xx.x + (t * (half + xx.x * t2)
-			     + (t * t) * (t3 + xx.x * t4 + t * t5));
-	      retval = one + yy.y;
-	      libc_fesetround (fe_val);
-	    }
-	  return retval;
-	}
-
-      /* Find the multiple of 2^-6 nearest x.  */
-      k = hx >> 20;
-      j = (0x00100000 | (hx & 0x000fffff)) >> (0x40c - k);
-      j = (j - 1) & ~1;
-      if (ix < 0)
-	j += 134;
-      /* Use FE_TONEAREST rounding mode for computing yy.y.
-	 Avoid set/reset of rounding mode if in FE_TONEAREST mode.  */
-      fe_val = get_rounding_mode ();
-      if (fe_val == FE_TONEAREST)
-	{
-	  z = xx.x - TBL2[j];
-	  t = z * z;
-	  yy.y = z + (t * (half + (z * t2))
-		      + (t * t) * (t3 + z * t4 + t * t5));
-	  retval = TBL2[j + 1] + TBL2[j + 1] * yy.y;
-	}
-      else
-	{
-	  libc_fesetround (FE_TONEAREST);
-	  z = xx.x - TBL2[j];
-	  t = z * z;
-	  yy.y = z + (t * (half + (z * t2))
-		      + (t * t) * (t3 + z * t4 + t * t5));
-	  retval = TBL2[j + 1] + TBL2[j + 1] * yy.y;
-	  libc_fesetround (fe_val);
-	}
-      return retval;
-    }
-
-  if (hx >= 0x40862e42)
-    {				/* x is large, infinite, or nan.  */
-      if (hx >= 0x7ff00000)
-	{
-	  if (ix == 0xfff00000 && xx.i_part[LOW_HALF] == 0)
-	    return zero;	/* exp(-inf) = 0.  */
-	  return (xx.x * xx.x);	/* exp(nan/inf) is nan or inf.  */
-	}
-      if (xx.x > threshold1)
-	{			/* Set overflow error condition.  */
-	  retval = hhuge * hhuge;
-	  return retval;
-	}
-      if (-xx.x > threshold2)
-	{			/* Set underflow error condition.  */
-	  double force_underflow = tiny * tiny;
-	  math_force_eval (force_underflow);
-	  retval = force_underflow;
-	  return retval;
-	}
-    }
-
-  /* Use FE_TONEAREST rounding mode for computing yy.y.
-     Avoid set/reset of rounding mode if already in FE_TONEAREST mode.  */
-  fe_val = get_rounding_mode ();
-  if (fe_val == FE_TONEAREST)
-    {
-      t = invln2_64 * xx.x;
-      if (ix < 0)
-	t -= half;
-      else
-	t += half;
-      k = (int) t;
-      j = (k & 0x3f) << 1;
-      m = k >> 6;
-      z = (xx.x - k * ln2_64hi) - k * ln2_64lo;
-
-      /* z is now in primary range.  */
-      t = z * z;
-      yy.y = z + (t * (half + z * t2) + (t * t) * (t3 + z * t4 + t * t5));
-      yy.y = TBL[j] + (TBL[j + 1] + TBL[j] * yy.y);
-    }
-  else
-    {
-      libc_fesetround (FE_TONEAREST);
-      t = invln2_64 * xx.x;
-      if (ix < 0)
-	t -= half;
-      else
-	t += half;
-      k = (int) t;
-      j = (k & 0x3f) << 1;
-      m = k >> 6;
-      z = (xx.x - k * ln2_64hi) - k * ln2_64lo;
-
-      /* z is now in primary range.  */
-      t = z * z;
-      yy.y = z + (t * (half + z * t2) + (t * t) * (t3 + z * t4 + t * t5));
-      yy.y = TBL[j] + (TBL[j + 1] + TBL[j] * yy.y);
-      libc_fesetround (fe_val);
-    }
 
-  if (m < -1021)
-    {
-      yy.y_part[HIGH_HALF] += (m + 54) << 20;
-      retval = twom54 * yy.y;
-      if (retval < DBL_MIN)
-	{
-	  double force_underflow = tiny * tiny;
-	  math_force_eval (force_underflow);
-	}
-      return retval;
-    }
-  yy.y_part[HIGH_HALF] += m << 20;
-  return yy.y;
+  {
+    SET_RESTORE_ROUND (FE_TONEAREST);
+
+    junk1.x = x;
+    m = junk1.i[HIGH_HALF];
+    n = m & hugeint;
+
+    if (n > smallint && n < bigint)
+      {
+	y = x * log2e.x + three51.x;
+	bexp = y - three51.x;	/*  multiply the result by 2**bexp        */
+
+	junk1.x = y;
+
+	eps = bexp * ln_two2.x;	/* x = bexp*ln(2) + t - eps               */
+	t = x - bexp * ln_two1.x;
+
+	y = t + three33.x;
+	base = y - three33.x;	/* t rounded to a multiple of 2**-18      */
+	junk2.x = y;
+	del = (t - base) - eps;	/*  x = bexp*ln(2) + base + del           */
+	eps = del + del * del * (p3.x * del + p2.x);
+
+	binexp.i[HIGH_HALF] = (junk1.i[LOW_HALF] + 1023) << 20;
+
+	i = ((junk2.i[LOW_HALF] >> 8) & 0xfffffffe) + 356;
+	j = (junk2.i[LOW_HALF] & 511) << 1;
+
+	al = coar.x[i] * fine.x[j];
+	bet = ((coar.x[i] * fine.x[j + 1] + coar.x[i + 1] * fine.x[j])
+	       + coar.x[i + 1] * fine.x[j + 1]);
+
+	rem = (bet + bet * eps) + al * eps;
+	res = al + rem;
+	cor = (al - res) + rem;
+	if (res == (res + cor * err_0))
+	  {
+	    retval = res * binexp.x;
+	    goto ret;
+	  }
+	else
+	  {
+	    retval = __slowexp (x);
+	    goto ret;
+	  }			/*if error is over bound */
+      }
+
+    if (n <= smallint)
+      {
+	retval = 1.0;
+	goto ret;
+      }
+
+    if (n >= badint)
+      {
+	if (n > infint)
+	  {
+	    retval = x + x;
+	    goto ret;
+	  }			/* x is NaN */
+	if (n < infint)
+	  {
+	    if (x > 0)
+	      goto ret_huge;
+	    else
+	      goto ret_tiny;
+	  }
+	/* x is finite,  cause either overflow or underflow  */
+	if (junk1.i[LOW_HALF] != 0)
+	  {
+	    retval = x + x;
+	    goto ret;
+	  }			/*  x is NaN  */
+	retval = (x > 0) ? inf.x : zero;	/* |x| = inf;  return either inf or 0 */
+	goto ret;
+      }
+
+    y = x * log2e.x + three51.x;
+    bexp = y - three51.x;
+    junk1.x = y;
+    eps = bexp * ln_two2.x;
+    t = x - bexp * ln_two1.x;
+    y = t + three33.x;
+    base = y - three33.x;
+    junk2.x = y;
+    del = (t - base) - eps;
+    eps = del + del * del * (p3.x * del + p2.x);
+    i = ((junk2.i[LOW_HALF] >> 8) & 0xfffffffe) + 356;
+    j = (junk2.i[LOW_HALF] & 511) << 1;
+    al = coar.x[i] * fine.x[j];
+    bet = ((coar.x[i] * fine.x[j + 1] + coar.x[i + 1] * fine.x[j])
+	   + coar.x[i + 1] * fine.x[j + 1]);
+    rem = (bet + bet * eps) + al * eps;
+    res = al + rem;
+    cor = (al - res) + rem;
+    if (m >> 31)
+      {
+	ex = junk1.i[LOW_HALF];
+	if (res < 1.0)
+	  {
+	    res += res;
+	    cor += cor;
+	    ex -= 1;
+	  }
+	if (ex >= -1022)
+	  {
+	    binexp.i[HIGH_HALF] = (1023 + ex) << 20;
+	    if (res == (res + cor * err_0))
+	      {
+		retval = res * binexp.x;
+		goto ret;
+	      }
+	    else
+	      {
+		retval = __slowexp (x);
+		goto check_uflow_ret;
+	      }			/*if error is over bound */
+	  }
+	ex = -(1022 + ex);
+	binexp.i[HIGH_HALF] = (1023 - ex) << 20;
+	res *= binexp.x;
+	cor *= binexp.x;
+	eps = 1.0000000001 + err_0 * binexp.x;
+	t = 1.0 + res;
+	y = ((1.0 - t) + res) + cor;
+	res = t + y;
+	cor = (t - res) + y;
+	if (res == (res + eps * cor))
+	  {
+	    binexp.i[HIGH_HALF] = 0x00100000;
+	    retval = (res - 1.0) * binexp.x;
+	    goto check_uflow_ret;
+	  }
+	else
+	  {
+	    retval = __slowexp (x);
+	    goto check_uflow_ret;
+	  }			/*   if error is over bound    */
+      check_uflow_ret:
+	if (retval < DBL_MIN)
+	  {
+	    double force_underflow = tiny * tiny;
+	    math_force_eval (force_underflow);
+	  }
+	if (retval == 0)
+	  goto ret_tiny;
+	goto ret;
+      }
+    else
+      {
+	binexp.i[HIGH_HALF] = (junk1.i[LOW_HALF] + 767) << 20;
+	if (res == (res + cor * err_0))
+	  retval = res * binexp.x * t256.x;
+	else
+	  retval = __slowexp (x);
+	if (isinf (retval))
+	  goto ret_huge;
+	else
+	  goto ret;
+      }
+  }
+ret:
+  return retval;
+
+ ret_huge:
+  return hhuge * hhuge;
+
+ ret_tiny:
+  return tiny * tiny;
 }
 #ifndef __ieee754_exp
 strong_alias (__ieee754_exp, __exp_finite)
 #endif
 
-#ifndef SECTION
-# define SECTION
-#endif
-
 /* Compute e^(x+xx).  The routine also receives bound of error of previous
    calculation.  If after computing exp the error exceeds the allowed bounds,
    the routine returns a non-positive number.  Otherwise it returns the
diff --git a/sysdeps/ieee754/dbl-64/e_pow.c b/sysdeps/ieee754/dbl-64/e_pow.c
index 2eb8dbfd5f..9f6439ee42 100644
--- a/sysdeps/ieee754/dbl-64/e_pow.c
+++ b/sysdeps/ieee754/dbl-64/e_pow.c
@@ -25,7 +25,7 @@
 /*             log1                                                        */
 /*             checkint                                                    */
 /* FILES NEEDED: dla.h endian.h mpa.h mydefs.h                             */
-/*               halfulp.c mpexp.c mplog.c slowpow.c mpa.c                 */
+/*               halfulp.c mpexp.c mplog.c slowexp.c slowpow.c mpa.c       */
 /*                          uexp.c  upow.c				   */
 /*               root.tbl uexp.tbl upow.tbl                                */
 /* An ultimate power routine. Given two IEEE double machine numbers y,x    */
diff --git a/sysdeps/ieee754/dbl-64/eexp.tbl b/sysdeps/ieee754/dbl-64/eexp.tbl
deleted file mode 100644
index 5941b9522b..0000000000
--- a/sysdeps/ieee754/dbl-64/eexp.tbl
+++ /dev/null
@@ -1,255 +0,0 @@
-/* EXP function tables - for use in computing double precision exponential
-   Copyright (C) 2017 Free Software Foundation, Inc.
-   This file is part of the GNU C Library.
-
-   The GNU C Library is free software; you can redistribute it and/or
-   modify it under the terms of the GNU Lesser General Public
-   License as published by the Free Software Foundation; either
-   version 2.1 of the License, or (at your option) any later version.
-
-   The GNU C Library is distributed in the hope that it will be useful,
-   but WITHOUT ANY WARRANTY; without even the implied warranty of
-   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
-   Lesser General Public License for more details.
-
-   You should have received a copy of the GNU Lesser General Public
-   License along with the GNU C Library; if not, see
-   <http://www.gnu.org/licenses/>.  */
-
-
-/*
-   TBL[2*j] is 2**(j/64), rounded to nearest.
-   TBL[2*j+1] is 2**(j/64) - TBL[2*j], rounded to nearest.
-   These values are used to approximate exp(x) using the formula
-   given in the comments for e_exp.c.  */
-
-static const double TBL[128] = {
-    0x1.0000000000000p+0,  0x0.0000000000000p+0,
-    0x1.02c9a3e778061p+0, -0x1.19083535b085dp-56,
-    0x1.059b0d3158574p+0,  0x1.d73e2a475b465p-55,
-    0x1.0874518759bc8p+0,  0x1.186be4bb284ffp-57,
-    0x1.0b5586cf9890fp+0,  0x1.8a62e4adc610bp-54,
-    0x1.0e3ec32d3d1a2p+0,  0x1.03a1727c57b52p-59,
-    0x1.11301d0125b51p+0, -0x1.6c51039449b3ap-54,
-    0x1.1429aaea92de0p+0, -0x1.32fbf9af1369ep-54,
-    0x1.172b83c7d517bp+0, -0x1.19041b9d78a76p-55,
-    0x1.1a35beb6fcb75p+0,  0x1.e5b4c7b4968e4p-55,
-    0x1.1d4873168b9aap+0,  0x1.e016e00a2643cp-54,
-    0x1.2063b88628cd6p+0,  0x1.dc775814a8495p-55,
-    0x1.2387a6e756238p+0,  0x1.9b07eb6c70573p-54,
-    0x1.26b4565e27cddp+0,  0x1.2bd339940e9d9p-55,
-    0x1.29e9df51fdee1p+0,  0x1.612e8afad1255p-55,
-    0x1.2d285a6e4030bp+0,  0x1.0024754db41d5p-54,
-    0x1.306fe0a31b715p+0,  0x1.6f46ad23182e4p-55,
-    0x1.33c08b26416ffp+0,  0x1.32721843659a6p-54,
-    0x1.371a7373aa9cbp+0, -0x1.63aeabf42eae2p-54,
-    0x1.3a7db34e59ff7p+0, -0x1.5e436d661f5e3p-56,
-    0x1.3dea64c123422p+0,  0x1.ada0911f09ebcp-55,
-    0x1.4160a21f72e2ap+0, -0x1.ef3691c309278p-58,
-    0x1.44e086061892dp+0,  0x1.89b7a04ef80d0p-59,
-    0x1.486a2b5c13cd0p+0,  0x1.3c1a3b69062f0p-56,
-    0x1.4bfdad5362a27p+0,  0x1.d4397afec42e2p-56,
-    0x1.4f9b2769d2ca7p+0, -0x1.4b309d25957e3p-54,
-    0x1.5342b569d4f82p+0, -0x1.07abe1db13cadp-55,
-    0x1.56f4736b527dap+0,  0x1.9bb2c011d93adp-54,
-    0x1.5ab07dd485429p+0,  0x1.6324c054647adp-54,
-    0x1.5e76f15ad2148p+0,  0x1.ba6f93080e65ep-54,
-    0x1.6247eb03a5585p+0, -0x1.383c17e40b497p-54,
-    0x1.6623882552225p+0, -0x1.bb60987591c34p-54,
-    0x1.6a09e667f3bcdp+0, -0x1.bdd3413b26456p-54,
-    0x1.6dfb23c651a2fp+0, -0x1.bbe3a683c88abp-57,
-    0x1.71f75e8ec5f74p+0, -0x1.16e4786887a99p-55,
-    0x1.75feb564267c9p+0, -0x1.0245957316dd3p-54,
-    0x1.7a11473eb0187p+0, -0x1.41577ee04992fp-55,
-    0x1.7e2f336cf4e62p+0,  0x1.05d02ba15797ep-56,
-    0x1.82589994cce13p+0, -0x1.d4c1dd41532d8p-54,
-    0x1.868d99b4492edp+0, -0x1.fc6f89bd4f6bap-54,
-    0x1.8ace5422aa0dbp+0,  0x1.6e9f156864b27p-54,
-    0x1.8f1ae99157736p+0,  0x1.5cc13a2e3976cp-55,
-    0x1.93737b0cdc5e5p+0, -0x1.75fc781b57ebcp-57,
-    0x1.97d829fde4e50p+0, -0x1.d185b7c1b85d1p-54,
-    0x1.9c49182a3f090p+0,  0x1.c7c46b071f2bep-56,
-    0x1.a0c667b5de565p+0, -0x1.359495d1cd533p-54,
-    0x1.a5503b23e255dp+0, -0x1.d2f6edb8d41e1p-54,
-    0x1.a9e6b5579fdbfp+0,  0x1.0fac90ef7fd31p-54,
-    0x1.ae89f995ad3adp+0,  0x1.7a1cd345dcc81p-54,
-