pl/math/erfc_1u8.c

101       double p2 = fma (P21, r2, -P20);
102       double p3 = -r * fma (P20, r2, -0.5);
103       double p4 = fma (fma (P42, r2, -P41), r2, P40);
104       double p5 = -r * fma (fma (P52, r2, -P51), r2, P50);
107       double p6 = fma (Q5 * r, p5, p4) * R5;
108       double p7 = fma (Q6 * r, p6, p5) * R6;
109       double p8 = fma (Q7 * r, p7, p6) * R7;
110       double p9 = fma (Q8 * r, p8, p7) * R8;
111       double p10 = fma (Q9 * r, p9, p8) * R9;
113       double p90 = fma (p10, d, p9);
114       double p78 = fma (p8, d, p7);
115       double p56 = fma (p6, d, p5);
116       double p34 = fma (p4, d, p3);
117       double p12 = fma (p2, d, p1);
118       double y = fma (p90, d2, p78);
119       y = fma (y, d2, p56);
120       y = fma (y, d2, p34);
121       y = fma (y, d2, p12);
123       y = fma (-fma (y, d2, d), scale, erfcr);
125       /* Handle sign and scale back in a single fma.  */
128       y = fma (y, fac, off);