Lines Matching full:double

2  * Double-precision erf(x) function.
26 /* Top 32 bits of a double.  */
28 top32 (double x)
36 double
37 erf (double x)
52 	      double y =  fma (TwoOverSqrtPiMinusOne, x, x);
58       double x2 = x * x;
62 	  double r1 = fma (x2, PA[1], PA[0]);
63 	  double r2 = fma (x2, PA[3], PA[2]);
64 	  double r3 = fma (x2, PA[5], PA[4]);
65 	  double r4 = fma (x2, PA[7], PA[6]);
66 	  double r5 = fma (x2, PA[9], PA[8]);
67 	  double x4 = x2 * x2;
68 	  double r = r5;
77 	  double x4, x8, r1n, r2n, r1d, r2d, r3d;
86 	  double P = r1n + x4 * r2n + x8 * NA[4];
87 	  double Q = r1d + x4 * r2d + x8 * r3d;
93       double a2, a4, a6, r1n, r2n, r3n, r4n, r1d, r2d, r3d, r4d;
94       double a = fabs (x) - 1.0;
106       double P = r1n + a2 * r2n + a4 * r3n + a6 * r4n;
107       double Q = r1d + a2 * r2d + a4 * r3d + a6 * r4d;
115       double a = fabs (x);
118       double r1 = fma (a, PC[1], PC[0]);
119       double r2 = fma (a, PC[3], PC[2]);
120       double r3 = fma (a, PC[5], PC[4]);
121       double r4 = fma (a, PC[7], PC[6]);
122       double r5 = fma (a, PC[9], PC[8]);
123       double r6 = fma (a, PC[11], PC[10]);
124       double r7 = fma (a, PC[13], PC[12]);
125       double r8 = fma (a, PC[15], PC[14]);
127       double a2 = a * a;
129       double r = r8;
145       double a = fabs (x);
148       double r1 = fma (a, PD[1], PD[0]);
149       double r2 = fma (a, PD[3], PD[2]);
150       double r3 = fma (a, PD[5], PD[4]);
151       double r4 = fma (a, PD[7], PD[6]);
152       double r5 = fma (a, PD[9], PD[8]);
153       double r6 = fma (a, PD[11], PD[10]);
154       double r7 = fma (a, PD[13], PD[12]);
155       double r8 = fma (a, PD[15], PD[14]);
156       double r9 = fma (a, PD[17], PD[16]);
158       double a2 = a * a;
160       double r = r9;
177       double a = fabs (x);
180       double r1 = fma (a, PE[1], PE[0]);
181       double r2 = fma (a, PE[3], PE[2]);
182       double r3 = fma (a, PE[5], PE[4]);
183       double r4 = fma (a, PE[7], PE[6]);
184       double r5 = fma (a, PE[9], PE[8]);
185       double r6 = fma (a, PE[11], PE[10]);
186       double r7 = fma (a, PE[13], PE[12]);
188       double a2 = a * a;
190       double r = r7;
205       double a = fabs (x);
208       double r1 = fma (a, PF[1], PF[0]);
209       double r2 = fma (a, PF[3], PF[2]);
210       double r3 = fma (a, PF[5], PF[4]);
211       double r4 = fma (a, PF[7], PF[6]);
212       double r5 = fma (a, PF[9], PF[8]);
213       double r6 = fma (a, PF[11], PF[10]);
214       double r7 = fma (a, PF[13], PF[12]);
215       double r8 = fma (a, PF[15], PF[14]);
216       double r9 = PF[16];
218       double a2 = a * a;
220       double r = r9;
239 	return (double) (1.0 - (sign << 1)) + 1.0 / x;