54 qq4u  = LD80C(0x87332f82cff4ff96, -11,  5.15746855583604912827e-4L),
126 sa4u  = LD80C(0xcaa83f403713e33e, 11,  3.24251544209971162003e+3L),
128 sa6u  = LD80C(0xb6ef97f9c753157b, 11,  2.92697460344182158454e+3L),
166 sb3u = LD80C(0x809c4ade46b927c7, 11,  2.05776827838541292848e+3L),
169 sb6u = LD80C(0x9d09a35988934631, 11,  2.51260238030767176221e+3L),
218 	long double ax,R,S,P,Q,s,y,z,r;  in erfl()  local
232 	ax = fabsl(x);  in erfl()
233 	if(ax < 0.84375) {  in erfl()
234 	    if(ax < 0x1p-34L) {  in erfl()
235 	        if(ax < 0x1p-16373L)	  in erfl()
245 	if(ax < 1.25) {  in erfl()
246 	    s = ax-one;  in erfl()
251 	if(ax >= 7) {			/* inf>|x|>= 7 */  in erfl()
254 	s = one/(ax*ax);  in erfl()
255 	if(ax < 2.85715) {	/* |x| < 2.85715 */  in erfl()
264 	z=(float)ax;  in erfl()
265 	r=expl(-z*z-0.5625)*expl((z-ax)*(z+ax)+R/S);  in erfl()
266 	if(x>=0) RETURNI(one-r/ax); else RETURNI(r/ax-one);  in erfl()
272 	long double ax,R,S,P,Q,s,y,z,r;  in erfcl()  local
285 	ax = fabsl(x);  in erfcl()
286 	if(ax < 0.84375L) {  in erfcl()
287 	    if(ax < 0x1p-34L)  in erfcl()
293 	    if(ax < 0.25L) {  	/* x<1/4 */  in erfcl()
301 	if(ax < 1.25L) {  in erfcl()
302 	    s = ax-one;  in erfcl()
312 	if(ax < 108) {			/* |x| < 108 */  in erfcl()
313  	    s = one/(ax*ax);  in erfcl()
314 	    if(ax < 2.85715) {		/* |x| < 2.85715 */  in erfcl()
319 	    } else if(ax < 7) {		/* | |x| < 7 */  in erfcl()
327 	    z = (float)ax;  in erfcl()
328 	    r = expl(-z*z-0.5625)*expl((z-ax)*(z+ax)+R/S);  in erfcl()
329 	    if(x>0) RETURNI(r/ax); else RETURNI(two-r/ax);  in erfcl()