`11210 libm should be cstyle(1ONBLD) clean`

```@@ -20,19 +20,20 @@
*/

/*
*/
+
/*
* Use is subject to license terms.
*/

#pragma weak __erf = erf
#pragma weak __erfc = erfc

-/* INDENT OFF */
+
/*
* double erf(double x)
* double erfc(double x)
*                           x
*                    2      |\
```

```@@ -123,21 +124,22 @@
*      7. Special case:
*              erf(0)  = 0, erf(inf)  = 1, erf(-inf) = -1,
*              erfc(0) = 1, erfc(inf) = 0, erfc(-inf) = 2,
*      erfc/erf(NaN) is NaN
*/
-/* INDENT ON */

#include "libm_macros.h"
#include <math.h>

static const double xxx[] = {
-/* tiny */      1e-300,
+/* tiny */
+        1e-300,
/* half */      5.00000000000000000000e-01,     /* 3FE00000, 00000000 */
/* one */       1.00000000000000000000e+00,     /* 3FF00000, 00000000 */
/* two */       2.00000000000000000000e+00,     /* 40000000, 00000000 */
/* erx */       8.45062911510467529297e-01,     /* 3FEB0AC1, 60000000 */
+
/*
* Coefficients for approximation to  erf on [0,0.84375]
*/
/* efx */        1.28379167095512586316e-01,    /* 3FC06EBA, 8214DB69 */
/* efx8 */       1.02703333676410069053e+00,    /* 3FF06EBA, 8214DB69 */
```

```@@ -149,10 +151,11 @@
/* qq1 */        3.97917223959155352819e-01,    /* 3FD97779, CDDADC09 */
/* qq2 */        6.50222499887672944485e-02,    /* 3FB0A54C, 5536CEBA */
/* qq3 */        5.08130628187576562776e-03,    /* 3F74D022, C4D36B0F */
/* qq4 */        1.32494738004321644526e-04,    /* 3F215DC9, 221C1A10 */
/* qq5 */       -3.96022827877536812320e-06,    /* BED09C43, 42A26120 */
+
/*
* Coefficients for approximation to  erf  in [0.84375,1.25]
*/
/* pa0 */       -2.36211856075265944077e-03,    /* BF6359B8, BEF77538 */
/* pa1 */        4.14856118683748331666e-01,    /* 3FDA8D00, AD92B34D */
```

```@@ -165,10 +168,11 @@
/* qa2 */        5.40397917702171048937e-01,    /* 3FE14AF0, 92EB6F33 */
/* qa3 */        7.18286544141962662868e-02,    /* 3FB2635C, D99FE9A7 */
/* qa4 */        1.26171219808761642112e-01,    /* 3FC02660, E763351F */
/* qa5 */        1.36370839120290507362e-02,    /* 3F8BEDC2, 6B51DD1C */
/* qa6 */        1.19844998467991074170e-02,    /* 3F888B54, 5735151D */
+
/*
* Coefficients for approximation to  erfc in [1.25,1/0.35]
*/
/* ra0 */       -9.86494403484714822705e-03,    /* BF843412, 600D6435 */
/* ra1 */       -6.93858572707181764372e-01,    /* BFE63416, E4BA7360 */
```

```@@ -184,10 +188,11 @@
/* sa4 */        6.45387271733267880336e+02,    /* 40842B19, 21EC2868 */
/* sa5 */        4.29008140027567833386e+02,    /* 407AD021, 57700314 */
/* sa6 */        1.08635005541779435134e+02,    /* 405B28A3, EE48AE2C */
/* sa7 */        6.57024977031928170135e+00,    /* 401A47EF, 8E484A93 */
/* sa8 */       -6.04244152148580987438e-02,    /* BFAEEFF2, EE749A62 */
+
/*
* Coefficients for approximation to  erfc in [1/.35,28]
*/
/* rb0 */       -9.86494292470009928597e-03,    /* BF843412, 39E86F4A */
/* rb1 */       -7.99283237680523006574e-01,    /* BFE993BA, 70C285DE */
```

```@@ -208,10 +213,11 @@
#define tiny    xxx[0]
#define half    xxx[1]
#define one     xxx[2]
#define two     xxx[3]
#define erx     xxx[4]
+
/*
* Coefficients for approximation to  erf on [0,0.84375]
*/
#define efx     xxx[5]
#define efx8    xxx[6]
```

```@@ -223,10 +229,11 @@
#define qq1     xxx[12]
#define qq2     xxx[13]
#define qq3     xxx[14]
#define qq4     xxx[15]
#define qq5     xxx[16]
+
/*
* Coefficients for approximation to  erf  in [0.84375,1.25]
*/
#define pa0     xxx[17]
#define pa1     xxx[18]
```

```@@ -239,10 +246,11 @@
#define qa2     xxx[25]
#define qa3     xxx[26]
#define qa4     xxx[27]
#define qa5     xxx[28]
#define qa6     xxx[29]
+
/*
* Coefficients for approximation to  erfc in [1.25,1/0.35]
*/
#define ra0     xxx[30]
#define ra1     xxx[31]
```

```@@ -258,10 +266,11 @@
#define sa4     xxx[41]
#define sa5     xxx[42]
#define sa6     xxx[43]
#define sa7     xxx[44]
#define sa8     xxx[45]
+
/*
* Coefficients for approximation to  erfc in [1/.35,28]
*/
#define rb0     xxx[46]
#define rb1     xxx[47]
```

```@@ -277,126 +286,145 @@
#define sb5     xxx[57]
#define sb6     xxx[58]
#define sb7     xxx[59]

double
-erf(double x) {
+erf(double x)
+{
int hx, ix, i;
double R, S, P, Q, s, y, z, r;

-        hx = ((int *) &x)[HIWORD];
+        hx = ((int *)&x)[HIWORD];
ix = hx & 0x7fffffff;
+
if (ix >= 0x7ff00000) { /* erf(nan)=nan */
if (ix >= 0x7ff80000)           /* assumes sparc-like QNaN */
return (x);
#endif
-                i = ((unsigned) hx >> 31) << 1;
-                return ((double) (1 - i) + one / x);    /* erf(+-inf)=+-1 */
+                i = ((unsigned)hx >> 31) << 1;
+                return ((double)(1 - i) + one / x);     /* erf(+-inf)=+-1 */
}

if (ix < 0x3feb0000) {  /* |x|<0.84375 */
if (ix < 0x3e300000) {  /* |x|<2**-28 */
if (ix < 0x00800000)    /* avoid underflow */
return (0.125 * (8.0 * x + efx8 * x));
+
return (x + efx * x);
}
+
z = x * x;
r = pp0 + z * (pp1 + z * (pp2 + z * (pp3 + z * pp4)));
-                s = one +
-                        z *(qq1 + z * (qq2 + z * (qq3 + z * (qq4 + z * qq5))));
+                s = one + z * (qq1 + z * (qq2 + z * (qq3 + z * (qq4 + z *
+                    qq5))));
y = r / s;
return (x + x * y);
}
+
if (ix < 0x3ff40000) {  /* 0.84375 <= |x| < 1.25 */
s = fabs(x) - one;
-                P = pa0 + s * (pa1 + s * (pa2 + s * (pa3 + s * (pa4 +
-                        s * (pa5 + s * pa6)))));
-                Q = one + s * (qa1 + s * (qa2 + s * (qa3 + s * (qa4 +
-                        s * (qa5 + s * qa6)))));
+                P = pa0 + s * (pa1 + s * (pa2 + s * (pa3 + s * (pa4 + s * (pa5 +
+                    s * pa6)))));
+                Q = one + s * (qa1 + s * (qa2 + s * (qa3 + s * (qa4 + s * (qa5 +
+                    s * qa6)))));
+
if (hx >= 0)
return (erx + P / Q);
else
return (-erx - P / Q);
}
+
if (ix >= 0x40180000) { /* inf > |x| >= 6 */
if (hx >= 0)
return (one - tiny);
else
return (tiny - one);
}
+
x = fabs(x);
s = one / (x * x);
+
if (ix < 0x4006DB6E) {  /* |x| < 1/0.35 */
-                R = ra0 + s * (ra1 + s * (ra2 + s * (ra3 + s * (ra4 +
-                        s * (ra5 + s * (ra6 + s * ra7))))));
-                S = one + s * (sa1 + s * (sa2 + s * (sa3 + s * (sa4 +
-                        s * (sa5 + s * (sa6 + s * (sa7 + s * sa8)))))));
+                R = ra0 + s * (ra1 + s * (ra2 + s * (ra3 + s * (ra4 + s * (ra5 +
+                    s * (ra6 + s * ra7))))));
+                S = one + s * (sa1 + s * (sa2 + s * (sa3 + s * (sa4 + s * (sa5 +
+                    s * (sa6 + s * (sa7 + s * sa8)))))));
} else {                        /* |x| >= 1/0.35 */
-                R = rb0 + s * (rb1 + s * (rb2 + s * (rb3 + s * (rb4 +
-                        s * (rb5 + s * rb6)))));
-                S = one + s * (sb1 + s * (sb2 + s * (sb3 + s * (sb4 +
-                        s * (sb5 + s * (sb6 + s * sb7))))));
+                R = rb0 + s * (rb1 + s * (rb2 + s * (rb3 + s * (rb4 + s * (rb5 +
+                    s * rb6)))));
+                S = one + s * (sb1 + s * (sb2 + s * (sb3 + s * (sb4 + s * (sb5 +
+                    s * (sb6 + s * sb7))))));
}
+
z = x;
-        ((int *) &z)[LOWORD] = 0;
+        ((int *)&z)[LOWORD] = 0;
r = exp(-z * z - 0.5625) * exp((z - x) * (z + x) + R / S);
+
if (hx >= 0)
return (one - r / x);
else
return (r / x - one);
}

double
-erfc(double x) {
+erfc(double x)
+{
int hx, ix;
double R, S, P, Q, s, y, z, r;

-        hx = ((int *) &x)[HIWORD];
+        hx = ((int *)&x)[HIWORD];
ix = hx & 0x7fffffff;
+
if (ix >= 0x7ff00000) { /* erfc(nan)=nan */
if (ix >= 0x7ff80000)           /* assumes sparc-like QNaN */
return (x);
#endif
/* erfc(+-inf)=0,2 */
-                return ((double) (((unsigned) hx >> 31) << 1) + one / x);
+                return ((double)(((unsigned)hx >> 31) << 1) + one / x);
}

if (ix < 0x3feb0000) {  /* |x| < 0.84375 */
if (ix < 0x3c700000)    /* |x| < 2**-56 */
return (one - x);
+
z = x * x;
r = pp0 + z * (pp1 + z * (pp2 + z * (pp3 + z * pp4)));
-                s = one +
-                        z * (qq1 + z * (qq2 + z * (qq3 + z * (qq4 + z * qq5))));
+                s = one + z * (qq1 + z * (qq2 + z * (qq3 + z * (qq4 + z *
+                    qq5))));
y = r / s;
+
if (hx < 0x3fd00000) {  /* x < 1/4 */
return (one - (x + x * y));
} else {
r = x * y;
r += (x - half);
return (half - r);
}
}
+
if (ix < 0x3ff40000) {  /* 0.84375 <= |x| < 1.25 */
s = fabs(x) - one;
-                P = pa0 + s * (pa1 + s * (pa2 + s * (pa3 + s * (pa4 +
-                        s * (pa5 + s * pa6)))));
-                Q = one + s * (qa1 + s * (qa2 + s * (qa3 + s * (qa4 +
-                        s * (qa5 + s * qa6)))));
+                P = pa0 + s * (pa1 + s * (pa2 + s * (pa3 + s * (pa4 + s * (pa5 +
+                    s * pa6)))));
+                Q = one + s * (qa1 + s * (qa2 + s * (qa3 + s * (qa4 + s * (qa5 +
+                    s * qa6)))));
+
if (hx >= 0) {
z = one - erx;
return (z - P / Q);
} else {
z = erx + P / Q;
return (one + z);
}
}
+
if (ix < 0x403c0000) {  /* |x|<28 */
x = fabs(x);
s = one / (x * x);
+
if (ix < 0x4006DB6D) {  /* |x| < 1/.35 ~ 2.857143 */
R = ra0 + s * (ra1 + s * (ra2 + s * (ra3 + s * (ra4 +
s * (ra5 + s * (ra6 + s * ra7))))));
S = one + s * (sa1 + s * (sa2 + s * (sa3 + s * (sa4 +
s * (sa5 + s * (sa6 + s * (sa7 + s * sa8)))))));
```

```@@ -408,13 +436,15 @@
R = rb0 + s * (rb1 + s * (rb2 + s * (rb3 + s * (rb4 +
s * (rb5 + s * rb6)))));
S = one + s * (sb1 + s * (sb2 + s * (sb3 + s * (sb4 +
s * (sb5 + s * (sb6 + s * sb7))))));
}
+
z = x;
-                ((int *) &z)[LOWORD] = 0;
+                ((int *)&z)[LOWORD] = 0;
r = exp(-z * z - 0.5625) * exp((z - x) * (z + x) + R / S);
+
if (hx > 0)
return (r / x);
else
return (two - r / x);
} else {
```