int main() { slong iter; flint_rand_t state; flint_printf("fast_mul_2exp_si...."); fflush(stdout); flint_randinit(state); for (iter = 0; iter < 100000 * arb_test_multiplier(); iter++) { fmpr_t x, y, z; mag_t xb, yb; slong e; fmpr_init(x); fmpr_init(y); fmpr_init(z); mag_init(xb); mag_init(yb); mag_randtest(xb, state, 15); e = n_randint(state, 10000) - n_randint(state, 10000); mag_get_fmpr(x, xb); mag_fast_mul_2exp_si(yb, xb, e); fmpr_mul_2exp_si(y, x, e); mag_get_fmpr(z, yb); MAG_CHECK_BITS(yb) if (!fmpr_equal(z, y)) { flint_printf("FAIL\n\n"); flint_printf("x = "); fmpr_printd(x, 15); flint_printf("\n\n"); flint_printf("y = "); fmpr_printd(y, 15); flint_printf("\n\n"); flint_printf("z = "); fmpr_printd(z, 15); flint_printf("\n\n"); flint_abort(); } fmpr_clear(x); fmpr_clear(y); fmpr_clear(z); mag_clear(xb); mag_clear(yb); } flint_randclear(state); flint_cleanup(); flint_printf("PASS\n"); return EXIT_SUCCESS; }
static void _acb_sqr_fast(acb_t z, const acb_t x, slong prec) { int inexact; mag_t am, bm, er, fr; mag_fast_init_set_arf(am, arb_midref(a)); mag_fast_init_set_arf(bm, arb_midref(b)); mag_init(er); mag_init(fr); mag_fast_addmul(er, am, ar); mag_fast_addmul(er, bm, br); mag_fast_mul_2exp_si(er, er, 1); mag_fast_addmul(er, ar, ar); mag_fast_addmul(er, br, br); mag_fast_addmul(fr, bm, ar); mag_fast_addmul(fr, am, br); mag_fast_addmul(fr, ar, br); mag_fast_mul_2exp_si(fr, fr, 1); inexact = arf_complex_sqr(arb_midref(e), arb_midref(f), arb_midref(a), arb_midref(b), prec, ARB_RND); if (inexact & 1) arf_mag_add_ulp(arb_radref(e), er, arb_midref(e), prec); else mag_set(arb_radref(e), er); if (inexact & 2) arf_mag_add_ulp(arb_radref(f), fr, arb_midref(f), prec); else mag_set(arb_radref(f), fr); }