Lines Matching refs:y
281 long double s,y,t;
288 y = fabsl(x);
289 if (y <= 0.84375L) {
290 if (y<=tiny) return x+P[0]*x;
291 s = y*y;
295 if (y<=1.25L) {
296 s = y-one;
299 } else if (y<=1.75L) {
300 s = y-onehalf;
304 if (y<=9.0L) t = erfcl(y); else t = tiny;
311 long double s,y,t;
337 y = __poly_libmq(x,15,R2);
338 t = y/__poly_libmq(x,16,S2);
340 y = __poly_libmq(x,14,R1);
341 t = y/__poly_libmq(x,15,S1);
346 * -x*x with a small absolute error. To this end, we set y
348 * zeros that y*y can be computed exactly and we rewrite x*x
349 * as y*y + (x-y)*(x+y), distributing the latter expression
352 * We could construct y in a portable way by setting
355 * y = (long double)i * 1/ptwo;
357 * where ptwo is some power of two large enough to make x-y
362 y = x;
363 *(2+(int*)&y) = *(3+(int*)&y) = 0;
364 t *= expl(-y*y)*expl(-(x-y)*(x+y));