Lines Matching defs:ex
95 mpfr_exp_t ex, ey, e;
99 ex = mpfr_get_exp (mpc_realref (z));
100 MPC_ASSERT(ex <= -1);
103 e = (ex >= ey) ? ex : ey;
153 ex = mpfr_get_exp (mpc_realref (s));
154 /* ulp(Re(s)) = 2^(ex+1-p) */
157 if (ex + 1 > e) /* divide kx by 2^(ex+1-e) */
158 while (ex + 1 > e)
161 ex--;
163 else /* multiply the error by 2^(e-(ex+1)), thus add e-(ex+1) to err */
164 err += e - (ex + 1);
218 mpfr_exp_t ex, ey, es, err;
238 ex = mpfr_get_exp (x);
242 if (ey > - prec_im - 4 || ex < 2 * ey + 4 + prec_re)
425 mpfr_exp_t ex, ey;
448 ex = mpfr_get_exp (mpc_realref(z1));
454 ex = ex - mpfr_get_exp (mpc_realref(z1));
455 ex = (ex <= 0) ? 0 : ex;
456 /* err(x) <= 2^ex * ulp(x) */
457 ex = ex + mpfr_get_exp (mpc_realref(z1)) - p;
458 /* err(x) <= 2^ex */
461 ex = (ex >= ey) ? ex : ey; /* err(x), err(y) <= 2^ex, i.e., the norm
462 of the error is bounded by |h|<=2^(ex+1/2) */
468 ex = (2 * ex + 1) - 2 - (ey - 1); /* |h^2/4/|t| <= 2^ex */
469 ex = (ex + 1) / 2; /* ceil(ex/2) */
470 /* express ex in terms of ulp(z1) */
473 ex = ex - ey + p;
475 err = (ex <= 0) ? 1 : ex + 1;
478 ex = mpfr_get_exp (mpc_realref(z1));
484 ex -= mpfr_get_exp (mpc_realref(z1)); /* cancellation in x */
486 ex = (ex >= ey) ? ex : ey; /* maximum cancellation */
487 err += ex;
491 ex = mpfr_get_exp (mpc_realref(z1));
493 ex = (ex >= ey) ? ex : ey;
494 err += ex - p; /* revert to absolute error <= 2^err */
496 err -= ex - 1; /* 1/|t| <= 1/|z| <= 2^(1-ex) */