51 pp9  = -3.37925756196555959454018189718117864e-11L,	/* 0xbfdc293e, 0x1d9060cb, 0xd043204a, 0x314cd7…
60 qq9  =  3.37405581964478060434410167262684979e-11L;	/* 0x3fdc28c8, 0xfb8fa1be, 0x10e57eec, 0xaa19e4…
192 	long double ax,R,S,P,Q,s,y,z,r;  in erfl()  local
204 	ax = fabsl(x);  in erfl()
205 	if(ax < 0.84375) {  in erfl()
206 	    if(ax < 0x1p-40L) {  in erfl()
207 	        if(ax < 0x1p-16373L)	  in erfl()
219 	if(ax < 1.25) {  in erfl()
220 	    s = ax-one;  in erfl()
227 	if (ax >= 9) {			/* inf>|x|>= 9 */  in erfl()
230 	s = one/(ax*ax);  in erfl()
231 	if(ax < 2.85715) {	/* |x| < 2.85715 */  in erfl()
246 	z = (float)ax;  in erfl()
247 	r = expl(-z*z-0.5625)*expl((z-ax)*(z+ax)+R/S);  in erfl()
248 	if(x>=0) return (one-r/ax); else return (r/ax-one);  in erfl()
254 	long double ax,R,S,P,Q,s,y,z,r;  in erfcl()  local
265 	ax = fabsl(x);  in erfcl()
266 	if(ax < 0.84375L) {  in erfcl()
267 	    if(ax < 0x1p-34L)  in erfcl()
275 	    if(ax < 0.25L) {  	/* x<1/4 */  in erfcl()
283 	if(ax < 1.25L) {  in erfcl()
284 	    s = ax-one;  in erfcl()
296 	if(ax < 108) {			/* |x| < 108 */  in erfcl()
297  	    s = one/(ax*ax);  in erfcl()
298 	    if(ax < 2.85715) {		/* |x| < 2.85715 */  in erfcl()
305 	    } else if(ax < 9) {  in erfcl()
319 	    z = (float)ax;  in erfcl()
320 	    r = expl(-z*z-0.5625)*expl((z-ax)*(z+ax)+R/S);  in erfcl()
321 	    if(x>0) return r/ax; else return two-r/ax;  in erfcl()