[fusl] clang-format fusl

fusl/src/math/exp.c

 /* origin: FreeBSD /usr/src/lib/msun/src/e_exp.c */
 /*
  * ====================================================
  * Copyright (C) 2004 by Sun Microsystems, Inc. All rights reserved.
  *
  * Permission to use, copy, modify, and distribute this
  * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
  */
@@ skipping 49 matching lines @@
  *      1 ulp (unit in the last place).
  *
  * Misc. info.
  *      For IEEE double
  *          if x >  709.782712893383973096 then exp(x) overflows
  *          if x < -745.133219101941108420 then exp(x) underflows
  */
 
 #include "libm.h"
 
-static const double
-half[2] = {0.5,-0.5},
-ln2hi = 6.93147180369123816490e-01, /* 0x3fe62e42, 0xfee00000 */
-ln2lo = 1.90821492927058770002e-10, /* 0x3dea39ef, 0x35793c76 */
-invln2 = 1.44269504088896338700e+00, /* 0x3ff71547, 0x652b82fe */
-P1   =  1.66666666666666019037e-01, /* 0x3FC55555, 0x5555553E */
-P2   = -2.77777777770155933842e-03, /* 0xBF66C16C, 0x16BEBD93 */
-P3   =  6.61375632143793436117e-05, /* 0x3F11566A, 0xAF25DE2C */
-P4   = -1.65339022054652515390e-06, /* 0xBEBBBD41, 0xC5D26BF1 */
-P5   =  4.13813679705723846039e-08; /* 0x3E663769, 0x72BEA4D0 */
+static const double half[2] = {0.5, -0.5},
+                    ln2hi =
+                        6.93147180369123816490e-01, /* 0x3fe62e42, 0xfee00000 */
+    ln2lo = 1.90821492927058770002e-10,             /* 0x3dea39ef, 0x35793c76 */
+    invln2 = 1.44269504088896338700e+00,            /* 0x3ff71547, 0x652b82fe */
+    P1 = 1.66666666666666019037e-01,                /* 0x3FC55555, 0x5555553E */
+    P2 = -2.77777777770155933842e-03,               /* 0xBF66C16C, 0x16BEBD93 */
+    P3 = 6.61375632143793436117e-05,                /* 0x3F11566A, 0xAF25DE2C */
+    P4 = -1.65339022054652515390e-06,               /* 0xBEBBBD41, 0xC5D26BF1 */
+    P5 = 4.13813679705723846039e-08;                /* 0x3E663769, 0x72BEA4D0 */
 
-double exp(double x)
-{
-        double_t hi, lo, c, xx, y;
-        int k, sign;
-        uint32_t hx;
+double exp(double x) {
+  double_t hi, lo, c, xx, y;
+  int k, sign;
+  uint32_t hx;
 
-        GET_HIGH_WORD(hx, x);
-        sign = hx>>31;
-        hx &= 0x7fffffff;  /* high word of |x| */
+  GET_HIGH_WORD(hx, x);
+  sign = hx >> 31;
+  hx &= 0x7fffffff; /* high word of |x| */
 
-        /* special cases */
-        if (hx >= 0x4086232b) {  /* if |x| >= 708.39... */
-                if (isnan(x))
-                        return x;
-                if (x > 709.782712893383973096) {
-                        /* overflow if x!=inf */
-                        x *= 0x1p1023;
-                        return x;
-                }
-                if (x < -708.39641853226410622) {
-                        /* underflow if x!=-inf */
-                        FORCE_EVAL((float)(-0x1p-149/x));
-                        if (x < -745.13321910194110842)
-                                return 0;
-                }
-        }
+  /* special cases */
+  if (hx >= 0x4086232b) { /* if |x| >= 708.39... */
+    if (isnan(x))
+      return x;
+    if (x > 709.782712893383973096) {
+      /* overflow if x!=inf */
+      x *= 0x1p1023;
+      return x;
+    }
+    if (x < -708.39641853226410622) {
+      /* underflow if x!=-inf */
+      FORCE_EVAL((float)(-0x1p-149 / x));
+      if (x < -745.13321910194110842)
+        return 0;
+    }
+  }
 
-        /* argument reduction */
-        if (hx > 0x3fd62e42) {  /* if |x| > 0.5 ln2 */
-                if (hx >= 0x3ff0a2b2)  /* if |x| >= 1.5 ln2 */
-                        k = (int)(invln2*x + half[sign]);
-                else
-                        k = 1 - sign - sign;
-                hi = x - k*ln2hi;  /* k*ln2hi is exact here */
-                lo = k*ln2lo;
-                x = hi - lo;
-        } else if (hx > 0x3e300000)  {  /* if |x| > 2**-28 */
-                k = 0;
-                hi = x;
-                lo = 0;
-        } else {
-                /* inexact if x!=0 */
-                FORCE_EVAL(0x1p1023 + x);
-                return 1 + x;
-        }
+  /* argument reduction */
+  if (hx > 0x3fd62e42) {  /* if |x| > 0.5 ln2 */
+    if (hx >= 0x3ff0a2b2) /* if |x| >= 1.5 ln2 */
+      k = (int)(invln2 * x + half[sign]);
+    else
+      k = 1 - sign - sign;
+    hi = x - k * ln2hi; /* k*ln2hi is exact here */
+    lo = k * ln2lo;
+    x = hi - lo;
+  } else if (hx > 0x3e300000) { /* if |x| > 2**-28 */
+    k = 0;
+    hi = x;
+    lo = 0;
+  } else {
+    /* inexact if x!=0 */
+    FORCE_EVAL(0x1p1023 + x);
+    return 1 + x;
+  }
 
-        /* x is now in primary range */
-        xx = x*x;
-        c = x - xx*(P1+xx*(P2+xx*(P3+xx*(P4+xx*P5))));
-        y = 1 + (x*c/(2-c) - lo + hi);
-        if (k == 0)
-                return y;
-        return scalbn(y, k);
+  /* x is now in primary range */
+  xx = x * x;
+  c = x - xx * (P1 + xx * (P2 + xx * (P3 + xx * (P4 + xx * P5))));
+  y = 1 + (x * c / (2 - c) - lo + hi);
+  if (k == 0)
+    return y;
+  return scalbn(y, k);
 }