2  * Double-precision erf(x) function.
24 /* Top 32 bits of a double.  */
26 top32 (double x)  in top32()
34 double
35 erf (double x)  in erf()
50 	      double y =  fma (TwoOverSqrtPiMinusOne, x, x);  in erf()
56       double x2 = x * x;  in erf()
60 	  double r1 = fma (x2, PA[1], PA[0]);  in erf()
61 	  double r2 = fma (x2, PA[3], PA[2]);  in erf()
62 	  double r3 = fma (x2, PA[5], PA[4]);  in erf()
63 	  double r4 = fma (x2, PA[7], PA[6]);  in erf()
64 	  double r5 = fma (x2, PA[9], PA[8]);  in erf()
65 	  double x4 = x2 * x2;  in erf()
66 	  double r = r5;  in erf()
75 	  double x4, x8, r1n, r2n, r1d, r2d, r3d;  in erf()
84 	  double P = r1n + x4 * r2n + x8 * NA[4];  in erf()
85 	  double Q = r1d + x4 * r2d + x8 * r3d;  in erf()
91       double a2, a4, a6, r1n, r2n, r3n, r4n, r1d, r2d, r3d, r4d;  in erf()
92       double a = fabs (x) - 1.0;  in erf()
104       double P = r1n + a2 * r2n + a4 * r3n + a6 * r4n;  in erf()
105       double Q = r1d + a2 * r2d + a4 * r3d + a6 * r4d;  in erf()
113       double a = fabs (x);  in erf()
116       double r1 = fma (a, PC[1], PC[0]);  in erf()
117       double r2 = fma (a, PC[3], PC[2]);  in erf()
118       double r3 = fma (a, PC[5], PC[4]);  in erf()
119       double r4 = fma (a, PC[7], PC[6]);  in erf()
120       double r5 = fma (a, PC[9], PC[8]);  in erf()
121       double r6 = fma (a, PC[11], PC[10]);  in erf()
122       double r7 = fma (a, PC[13], PC[12]);  in erf()
123       double r8 = fma (a, PC[15], PC[14]);  in erf()
125       double a2 = a * a;  in erf()
127       double r = r8;  in erf()
143       double a = fabs (x);  in erf()
146       double r1 = fma (a, PD[1], PD[0]);  in erf()
147       double r2 = fma (a, PD[3], PD[2]);  in erf()
148       double r3 = fma (a, PD[5], PD[4]);  in erf()
149       double r4 = fma (a, PD[7], PD[6]);  in erf()
150       double r5 = fma (a, PD[9], PD[8]);  in erf()
151       double r6 = fma (a, PD[11], PD[10]);  in erf()
152       double r7 = fma (a, PD[13], PD[12]);  in erf()
153       double r8 = fma (a, PD[15], PD[14]);  in erf()
154       double r9 = fma (a, PD[17], PD[16]);  in erf()
156       double a2 = a * a;  in erf()
158       double r = r9;  in erf()
175       double a = fabs (x);  in erf()
178       double r1 = fma (a, PE[1], PE[0]);  in erf()
179       double r2 = fma (a, PE[3], PE[2]);  in erf()
180       double r3 = fma (a, PE[5], PE[4]);  in erf()
181       double r4 = fma (a, PE[7], PE[6]);  in erf()
182       double r5 = fma (a, PE[9], PE[8]);  in erf()
183       double r6 = fma (a, PE[11], PE[10]);  in erf()
184       double r7 = fma (a, PE[13], PE[12]);  in erf()
186       double a2 = a * a;  in erf()
188       double r = r7;  in erf()
203       double a = fabs (x);  in erf()
206       double r1 = fma (a, PF[1], PF[0]);  in erf()
207       double r2 = fma (a, PF[3], PF[2]);  in erf()
208       double r3 = fma (a, PF[5], PF[4]);  in erf()
209       double r4 = fma (a, PF[7], PF[6]);  in erf()
210       double r5 = fma (a, PF[9], PF[8]);  in erf()
211       double r6 = fma (a, PF[11], PF[10]);  in erf()
212       double r7 = fma (a, PF[13], PF[12]);  in erf()
213       double r8 = fma (a, PF[15], PF[14]);  in erf()
214       double r9 = PF[16];  in erf()
216       double a2 = a * a;  in erf()
218       double r = r9;  in erf()
237 	return (double) (1.0 - (sign << 1)) + 1.0 / x;  in erf()