int main()
{
slong iter;
flint_rand_t state;
flint_printf("set_d....");
fflush(stdout);
flint_randinit(state);
for (iter = 0; iter < 100000 * arb_test_multiplier(); iter++)
{
fmpr_t a, b, c;
mag_t m;
double x;
fmpr_init(a);
fmpr_init(b);
fmpr_init(c);
mag_init(m);
x = d_randtest2(state);
x = ldexp(x, 100 - n_randint(state, 200));
if (n_randint(state, 100) == 0)
x = 0.0;
fmpr_set_d(a, x);
mag_set_d(m, x);
mag_get_fmpr(b, m);
fmpr_set(c, a);
fmpr_mul_ui(c, c, 1025, MAG_BITS, FMPR_RND_UP);
fmpr_mul_2exp_si(c, c, -10);
MAG_CHECK_BITS(m)
if (!(fmpr_cmpabs(a, b) <= 0 && fmpr_cmpabs(b, c) <= 0))
{
flint_printf("FAIL\n\n");
flint_printf("a = "); fmpr_print(a); flint_printf("\n\n");
flint_printf("b = "); fmpr_print(b); flint_printf("\n\n");
flint_printf("c = "); fmpr_print(c); flint_printf("\n\n");
abort();
}
fmpr_clear(a);
fmpr_clear(b);
fmpr_clear(c);
mag_clear(m);
}
flint_randclear(state);
flint_cleanup();
flint_printf("PASS\n");
return EXIT_SUCCESS;
}
int
fmpcb_cmpabs_approx(const fmpcb_t x, const fmpcb_t y)
{
const fmpr_struct *xm, *ym;
if (fmpr_cmpabs(fmprb_midref(fmpcb_realref(x)), fmprb_midref(fmpcb_imagref(x))) >= 0)
xm = fmprb_midref(fmpcb_realref(x));
else
xm = fmprb_midref(fmpcb_imagref(x));
if (fmpr_cmpabs(fmprb_midref(fmpcb_realref(y)), fmprb_midref(fmpcb_imagref(y))) >= 0)
ym = fmprb_midref(fmpcb_realref(y));
else
ym = fmprb_midref(fmpcb_imagref(y));
return fmpr_cmpabs(xm, ym);
}
int main()
{
slong iter;
flint_rand_t state;
flint_printf("cmpabs....");
fflush(stdout);
flint_randinit(state);
for (iter = 0; iter < 100000; iter++)
{
slong bits;
fmpr_t x, y;
mpfr_t X, Y;
int cmp1, cmp2;
bits = 2 + n_randint(state, 200);
fmpr_init(x);
fmpr_init(y);
mpfr_init2(X, bits);
mpfr_init2(Y, bits);
fmpr_randtest_special(x, state, bits, 10);
fmpr_randtest_special(y, state, bits, 10);
fmpr_get_mpfr(X, x, MPFR_RNDN);
fmpr_get_mpfr(Y, y, MPFR_RNDN);
mpfr_abs(X, X, MPFR_RNDN);
mpfr_abs(Y, Y, MPFR_RNDN);
cmp1 = fmpr_cmpabs(x, y);
cmp2 = mpfr_cmp(X, Y);
if (cmp1 != cmp2)
{
flint_printf("FAIL\n\n");
flint_printf("x = "); fmpr_print(x); flint_printf("\n\n");
flint_printf("y = "); fmpr_print(y); flint_printf("\n\n");
flint_printf("cmp1 = %d, cmp2 = %d\n\n", cmp1, cmp2);
abort();
}
fmpr_clear(x);
fmpr_clear(y);
mpfr_clear(X);
mpfr_clear(Y);
}
flint_randclear(state);
flint_cleanup();
flint_printf("PASS\n");
return EXIT_SUCCESS;
}
int main()
{
slong iter;
flint_rand_t state;
flint_printf("rfac_ui....");
fflush(stdout);
flint_randinit(state);
for (iter = 0; iter < 2000 * arb_test_multiplier(); iter++)
{
fmpr_t x, y;
fmpz_t f;
mag_t xb;
ulong n;
fmpr_init(x);
fmpr_init(y);
fmpz_init(f);
mag_init(xb);
mag_randtest_special(xb, state, 80);
n = n_randtest(state) % 2000;
mag_rfac_ui(xb, n);
fmpz_fac_ui(f, n);
fmpr_one(x);
fmpr_div_fmpz(x, x, f, 2 * MAG_BITS, FMPR_RND_UP);
mag_get_fmpr(y, xb);
MAG_CHECK_BITS(xb)
if (!(fmpr_cmpabs(y, x) >= 0))
{
flint_printf("FAIL\n\n");
flint_printf("n = %wu\n\n", n);
flint_printf("x = "); fmpr_print(x); flint_printf("\n\n");
flint_printf("y = "); fmpr_print(y); flint_printf("\n\n");
abort();
}
fmpr_clear(x);
fmpr_clear(y);
fmpz_clear(f);
mag_clear(xb);
}
flint_randclear(state);
flint_cleanup();
flint_printf("PASS\n");
return EXIT_SUCCESS;
}
int main()
{
slong iter;
flint_rand_t state;
flint_printf("binpow_uiui....");
fflush(stdout);
flint_randinit(state);
for (iter = 0; iter < 10000 * arb_test_multiplier(); iter++)
{
fmpr_t x, y;
mag_t b;
ulong m, n;
fmpr_init(x);
fmpr_init(y);
mag_init(b);
m = n_randtest(state);
n = n_randtest(state);
fmpr_one(x);
fmpr_div_ui(x, x, m, 128, FMPR_RND_UP);
fmpr_add_ui(x, x, 1, 128, FMPR_RND_UP);
fmpr_pow_sloppy_ui(x, x, n, 128, FMPR_RND_UP);
mag_binpow_uiui(b, m, n);
mag_get_fmpr(y, b);
MAG_CHECK_BITS(b)
if (!(fmpr_cmpabs(x, y) <= 0))
{
flint_printf("FAIL\n\n");
flint_printf("m = %wu\n\n", m);
flint_printf("n = %wu\n\n", n);
flint_printf("x = "); fmpr_printd(x, 10); flint_printf("\n\n");
flint_printf("y = "); fmpr_printd(y, 10); flint_printf("\n\n");
flint_abort();
}
fmpr_clear(x);
fmpr_clear(y);
mag_clear(b);
}
flint_randclear(state);
flint_cleanup();
flint_printf("PASS\n");
return EXIT_SUCCESS;
}
