int main()
{
slong iter;
flint_rand_t state;
flint_printf("bell_sum_taylor....");
fflush(stdout);
flint_randinit(state);
for (iter = 0; iter < 1000 * arb_test_multiplier(); iter++)
{
arb_t s1, s2;
fmpz_t a, b, n;
slong prec;
arb_init(s1);
arb_init(s2);
fmpz_init(a);
fmpz_init(b);
fmpz_init(n);
prec = 2 + n_randint(state, 300);
fmpz_randtest_unsigned(n, state, 1 + n_randint(state, 100));
fmpz_randtest_unsigned(a, state, 1 + n_randint(state, 100));
fmpz_add_ui(b, a, n_randint(state, 100));
arb_bell_sum_bsplit(s1, n, a, b, NULL, prec);
arb_bell_sum_taylor(s2, n, a, b, NULL, prec);
if (!arb_overlaps(s1, s2) || (arb_rel_accuracy_bits(s1) < prec - 4)
|| (arb_rel_accuracy_bits(s2) < prec - 4))
{
flint_printf("FAIL: overlap or accuracy\n\n");
flint_printf("prec = %wd\n\n", prec);
flint_printf("n = "); fmpz_print(n); flint_printf("\n\n");
flint_printf("a = "); fmpz_print(a); flint_printf("\n\n");
flint_printf("b = "); fmpz_print(b); flint_printf("\n\n");
flint_printf("s1 = "); arb_printn(s1, 100, 0); flint_printf("\n\n");
flint_printf("s2 = "); arb_printn(s2, 100, 0); flint_printf("\n\n");
flint_abort();
}
arb_clear(s1);
arb_clear(s2);
fmpz_clear(a);
fmpz_clear(b);
fmpz_clear(n);
}
flint_randclear(state);
flint_cleanup();
flint_printf("PASS\n");
return EXIT_SUCCESS;
}
void
fmpz_randtest(fmpz_t f, flint_rand_t state, mp_bitcnt_t bits)
{
ulong m;
fmpz_randtest_unsigned(f, state, bits);
m = n_randlimb(state);
if (m & UWORD(1))
fmpz_neg(f, f);
}
void
fmpz_mat_randtest_unsigned(fmpz_mat_t mat, flint_rand_t state, mp_bitcnt_t bits)
{
slong r, c, i, j;
r = mat->r;
c = mat->c;
for (i = 0; i < r; i++)
for (j = 0; j < c; j++)
fmpz_randtest_unsigned(mat->rows[i] + j, state, bits);
}
void
fmpz_poly_randtest_unsigned(fmpz_poly_t f, flint_rand_t state,
long len, mp_bitcnt_t bits)
{
long i;
fmpz_poly_fit_length(f, len);
for (i = 0; i < len; i++)
fmpz_randtest_unsigned(f->coeffs + i, state, bits);
_fmpz_poly_set_length(f, len);
_fmpz_poly_normalise(f);
}
int
main(void)
{
int i, result;
FLINT_TEST_INIT(state);
flint_printf("swap....");
fflush(stdout);
for (i = 0; i < 10000; i++)
{
fmpz_t p;
fmpz_mod_poly_t a, b, c;
fmpz_init(p);
fmpz_randtest_unsigned(p, state, 2 * FLINT_BITS);
fmpz_add_ui(p, p, 2);
fmpz_mod_poly_init(a, p);
fmpz_mod_poly_init(b, p);
fmpz_mod_poly_init(c, p);
fmpz_mod_poly_randtest(a, state, n_randint(state, 100));
fmpz_mod_poly_randtest(b, state, n_randint(state, 100));
fmpz_mod_poly_set(c, b);
fmpz_mod_poly_swap(a, b);
result = (fmpz_mod_poly_equal(a, c));
if (!result)
{
flint_printf("FAIL:\n");
fmpz_mod_poly_print(a), flint_printf("\n\n");
fmpz_mod_poly_print(b), flint_printf("\n\n");
fmpz_mod_poly_print(c), flint_printf("\n\n");
abort();
}
fmpz_mod_poly_clear(a);
fmpz_mod_poly_clear(b);
fmpz_mod_poly_clear(c);
fmpz_clear(p);
}
FLINT_TEST_CLEANUP(state);
flint_printf("PASS\n");
return 0;
}
int
main(void)
{
int i, result;
flint_rand_t state;
printf("neg....");
fflush(stdout);
flint_randinit(state);
for (i = 0; i < 10000; i++)
{
fmpz_t p;
fmpz_mod_poly_t a, b, c;
fmpz_init(p);
fmpz_randtest_unsigned(p, state, 2 * FLINT_BITS);
fmpz_add_ui(p, p, 2);
fmpz_mod_poly_init(a, p);
fmpz_mod_poly_init(b, p);
fmpz_mod_poly_init(c, p);
fmpz_mod_poly_randtest(a, state, n_randint(state, 100));
fmpz_mod_poly_neg(b, a);
fmpz_mod_poly_neg(c, b);
result = (fmpz_mod_poly_equal(a, c));
if (!result)
{
printf("FAIL:\n");
fmpz_mod_poly_print(a), printf("\n\n");
fmpz_mod_poly_print(b), printf("\n\n");
fmpz_mod_poly_print(c), printf("\n\n");
abort();
}
fmpz_mod_poly_clear(a);
fmpz_mod_poly_clear(b);
fmpz_mod_poly_clear(c);
fmpz_clear(p);
}
flint_randclear(state);
_fmpz_cleanup();
printf("PASS\n");
return 0;
}
void
fmpz_randtest_mod(fmpz_t f, flint_rand_t state, const fmpz_t m)
{
fmpz_t t;
fmpz_init(t);
fmpz_randtest_unsigned(t, state, fmpz_bits(m) + 2);
fmpz_mod(t, t, m);
if (n_randlimb(state) & UWORD(1))
{
fmpz_sub(t, m, t);
fmpz_sub_ui(t, t, UWORD(1));
}
fmpz_set(f, t);
fmpz_clear(t);
}
int
main(void)
{
int i, result;
FLINT_TEST_INIT(state);
flint_printf("zero....");
fflush(stdout);
for (i = 0; i < 10000; i++)
{
fmpz_t p;
fmpz_mod_poly_t a;
fmpz_init(p);
fmpz_randtest_unsigned(p, state, 2 * FLINT_BITS);
fmpz_add_ui(p, p, 2);
fmpz_mod_poly_init(a, p);
fmpz_mod_poly_randtest(a, state, n_randint(state, 100));
fmpz_mod_poly_zero(a);
result = (fmpz_mod_poly_is_zero(a));
if (!result)
{
flint_printf("FAIL:\n");
flint_printf("a = "), fmpz_mod_poly_print(a), flint_printf("\n\n");
abort();
}
fmpz_mod_poly_clear(a);
fmpz_clear(p);
}
FLINT_TEST_CLEANUP(state);
flint_printf("PASS\n");
return 0;
}
void fmpq_poly_randtest_unsigned(fmpq_poly_t poly, flint_rand_t state,
long len, mp_bitcnt_t bits)
{
ulong m;
m = n_randlimb(NULL);
fmpq_poly_fit_length(poly, len);
_fmpq_poly_set_length(poly, len);
if (m & 1UL)
{
_fmpz_vec_randtest_unsigned(poly->coeffs, state, len, bits);
}
else
{
fmpz_t x;
fmpz_init(x);
fmpz_randtest_unsigned(x, state, bits / 2);
_fmpz_vec_randtest_unsigned(poly->coeffs, state, len, (bits + 1) / 2);
_fmpz_vec_scalar_mul_fmpz(poly->coeffs, poly->coeffs, len, x);
fmpz_clear(x);
}
if (m & 2UL)
{
fmpz_randtest_not_zero(poly->den, state, FLINT_MAX(bits, 1));
fmpz_abs(poly->den, poly->den);
fmpq_poly_canonicalise(poly);
}
else
{
fmpz_set_ui(poly->den, 1);
_fmpq_poly_normalise(poly);
}
}
int
main(void)
{
int i, result;
FLINT_TEST_INIT(state);
flint_printf("divrem_f....");
fflush(stdout);
/* Check q*b + r = a when gcd(lead(B),p) = 1, no aliasing */
for (i = 0; i < 5000; i++)
{
fmpz_t f, p;
fmpz_mod_poly_t a, b, q, r, t;
fmpz_init(f);
fmpz_init(p);
fmpz_randtest_unsigned(p, state, 2 * FLINT_BITS);
fmpz_add_ui(p, p, 2);
fmpz_mod_poly_init(a, p);
fmpz_mod_poly_init(b, p);
fmpz_mod_poly_init(q, p);
fmpz_mod_poly_init(r, p);
fmpz_mod_poly_init(t, p);
fmpz_mod_poly_randtest(a, state, n_randint(state, 100));
fmpz_mod_poly_randtest_not_zero(b, state, n_randint(state, 100) + 1);
{
fmpz_t d;
fmpz *leadB = fmpz_mod_poly_lead(b);
fmpz_init(d);
fmpz_gcd(d, p, leadB);
while (!fmpz_is_one(d))
{
fmpz_divexact(leadB, leadB, d);
fmpz_gcd(d, p, leadB);
}
fmpz_clear(d);
}
fmpz_mod_poly_divrem_f(f, q, r, a, b);
fmpz_mod_poly_mul(t, q, b);
fmpz_mod_poly_add(t, t, r);
result = (fmpz_is_one(f) && fmpz_mod_poly_equal(a, t));
if (!result)
{
flint_printf("FAIL (divrem):\n");
flint_printf("p = "), fmpz_print(p), flint_printf("\n\n");
flint_printf("f = "), fmpz_print(f), flint_printf("\n\n");
flint_printf("a = "), fmpz_mod_poly_print(a), flint_printf("\n\n");
flint_printf("b = "), fmpz_mod_poly_print(b), flint_printf("\n\n");
flint_printf("q = "), fmpz_mod_poly_print(q), flint_printf("\n\n");
flint_printf("r = "), fmpz_mod_poly_print(r), flint_printf("\n\n");
flint_printf("t = "), fmpz_mod_poly_print(t), flint_printf("\n\n");
abort();
}
fmpz_mod_poly_clear(a);
fmpz_mod_poly_clear(b);
fmpz_mod_poly_clear(q);
fmpz_mod_poly_clear(r);
fmpz_mod_poly_clear(t);
fmpz_clear(f);
fmpz_clear(p);
}
/* Check f | p when gcd(lead(B),p) > 1 */
for (i = 0; i < 5000; i++)
{
fmpz_t f, p, q1, q2;
fmpz_mod_poly_t a, b, q, r, t;
fmpz_init(f);
fmpz_init(p);
fmpz_init(q1);
fmpz_init(q2);
fmpz_randtest_unsigned(q1, state, 2 * FLINT_BITS);
fmpz_randtest_unsigned(q2, state, 2 * FLINT_BITS);
fmpz_add_ui(q1, q1, 2);
fmpz_add_ui(q2, q2, 2);
fmpz_mul(p, q1, q2);
fmpz_mod_poly_init(a, p);
fmpz_mod_poly_init(b, p);
fmpz_mod_poly_init(q, p);
fmpz_mod_poly_init(r, p);
fmpz_mod_poly_init(t, p);
fmpz_mod_poly_randtest(a, state, n_randint(state, 100));
fmpz_mod_poly_randtest_not_zero(b, state, n_randint(state, 100) + 1);
{
fmpz_t d;
fmpz *leadB = fmpz_mod_poly_lead(b);
fmpz_init(d);
fmpz_gcd(d, p, leadB);
if (fmpz_is_one(d))
fmpz_set(leadB, q1);
fmpz_clear(d);
}
fmpz_mod_poly_divrem_f(f, q, r, a, b);
fmpz_mod_poly_mul(t, q, b);
fmpz_mod_poly_add(t, t, r);
result = (fmpz_cmp_ui(f, 1) > 0 && fmpz_cmp(f, p) < 0 && fmpz_divisible(p, f));
if (!result)
{
flint_printf("FAIL (factor):\n");
flint_printf("p = "), fmpz_print(p), flint_printf("\n\n");
flint_printf("f = "), fmpz_print(f), flint_printf("\n\n");
flint_printf("q1 = "), fmpz_print(q1), flint_printf("\n\n");
flint_printf("q2 = "), fmpz_print(q2), flint_printf("\n\n");
flint_printf("a = "), fmpz_mod_poly_print(a), flint_printf("\n\n");
flint_printf("b = "), fmpz_mod_poly_print(b), flint_printf("\n\n");
flint_printf("q = "), fmpz_mod_poly_print(q), flint_printf("\n\n");
flint_printf("r = "), fmpz_mod_poly_print(r), flint_printf("\n\n");
flint_printf("t = "), fmpz_mod_poly_print(t), flint_printf("\n\n");
abort();
}
fmpz_mod_poly_clear(a);
fmpz_mod_poly_clear(b);
fmpz_mod_poly_clear(q);
fmpz_mod_poly_clear(r);
fmpz_mod_poly_clear(t);
fmpz_clear(f);
fmpz_clear(p);
fmpz_clear(q1);
fmpz_clear(q2);
}
FLINT_TEST_CLEANUP(state);
flint_printf("PASS\n");
return 0;
}
int main()
{
slong iter;
flint_rand_t state;
flint_printf("zeta_nzeros_gram....");
fflush(stdout);
flint_randinit(state);
for (iter = 0; iter < 130 + 20 * arb_test_multiplier(); iter++)
{
arb_t t, x;
fmpz_t N, n;
slong prec1, prec2;
arb_init(t);
arb_init(x);
fmpz_init(n);
fmpz_init(N);
if (iter < 130)
{
fmpz_set_si(n, iter - 1);
}
else
{
fmpz_randtest_unsigned(n, state, 20);
fmpz_add_ui(n, n, 129);
}
prec1 = 2 + n_randtest(state) % 100;
prec2 = 2 + n_randtest(state) % 100;
acb_dirichlet_zeta_nzeros_gram(N, n);
acb_dirichlet_gram_point(t, n, NULL, NULL, prec1);
acb_dirichlet_zeta_nzeros(x, t, prec2);
if (!arb_contains_fmpz(x, N))
{
flint_printf("FAIL: containment\n\n");
flint_printf("n = "); fmpz_print(n);
flint_printf(" prec1 = %wd prec2 = %wd\n\n", prec1, prec2);
flint_printf("N = "); fmpz_print(N); flint_printf("\n\n");
flint_printf("t = "); arb_printn(t, 100, 0); flint_printf("\n\n");
flint_printf("x = "); arb_printn(x, 100, 0); flint_printf("\n\n");
flint_abort();
}
arb_clear(t);
arb_clear(x);
fmpz_clear(n);
fmpz_clear(N);
}
flint_randclear(state);
flint_cleanup();
flint_printf("PASS\n");
return EXIT_SUCCESS;
}
int
main(void)
{
int i, result;
FLINT_TEST_INIT(state);
flint_printf("pow_trunc....");
fflush(stdout);
/* Check aliasing */
for (i = 0; i < 10000; i++)
{
fmpz_mod_poly_t a, b, c;
fmpz_t p;
slong e, trunc;
fmpz_init(p);
fmpz_randtest_unsigned(p, state, 2 * FLINT_BITS);
fmpz_add_ui(p, p, 2);
fmpz_mod_poly_init(a, p);
fmpz_mod_poly_init(b, p);
fmpz_mod_poly_init(c, p);
fmpz_mod_poly_randtest(a, state, n_randint(state, 30));
e = n_randint(state, 20);
trunc = n_randint(state, 30);
fmpz_mod_poly_set(c, a);
fmpz_mod_poly_pow_trunc(b, a, e, trunc);
fmpz_mod_poly_pow_trunc(c, c, e, trunc);
result = (fmpz_mod_poly_equal(b, c));
if (!result)
{
flint_printf("FAIL aliasing:\n");
flint_printf("a->length = %wd, p = %wu, exp = %wd, trunc = %wd\n",
a->length, a->p, e, trunc);
flint_printf("a:\n"); fmpz_mod_poly_print(a), flint_printf("\n\n");
flint_printf("b:\n"); fmpz_mod_poly_print(b), flint_printf("\n\n");
flint_printf("c:\n"); fmpz_mod_poly_print(c), flint_printf("\n\n");
abort();
}
fmpz_clear(p);
fmpz_mod_poly_clear(a);
fmpz_mod_poly_clear(b);
fmpz_mod_poly_clear(c);
}
/* Check powering against naive method */
for (i = 0; i < 10000; i++)
{
fmpz_mod_poly_t a, b, c;
fmpz_t p;
slong e, trunc;
fmpz_init(p);
fmpz_randtest_unsigned(p, state, 2 * FLINT_BITS);
fmpz_add_ui(p, p, 2);
fmpz_mod_poly_init(a, p);
fmpz_mod_poly_init(b, p);
fmpz_mod_poly_init(c, p);
fmpz_mod_poly_randtest(a, state, n_randint(state, 30));
e = n_randint(state, 20);
trunc = n_randint(state, 30);
fmpz_mod_poly_pow_trunc(b, a, e, trunc);
fmpz_mod_poly_pow(c, a, e);
fmpz_mod_poly_truncate(c, trunc);
result = (fmpz_mod_poly_equal(b, c)
|| (a->length == 0 && e == 0 && c->length == 1 && c->coeffs[0] == 1));
if (!result)
{
flint_printf("FAIL:\n");
flint_printf("a->length = %wd, p = %wu, exp = %wd, trunc = %wd\n",
a->length, a->p, e, trunc);
flint_printf("a:\n"); fmpz_mod_poly_print(a), flint_printf("\n\n");
flint_printf("b:\n"); fmpz_mod_poly_print(b), flint_printf("\n\n");
flint_printf("c:\n"); fmpz_mod_poly_print(c), flint_printf("\n\n");
abort();
}
fmpz_clear(p);
fmpz_mod_poly_clear(a);
fmpz_mod_poly_clear(b);
fmpz_mod_poly_clear(c);
}
FLINT_TEST_CLEANUP(state);
flint_printf("PASS\n");
return 0;
}
int
main(void)
{
int i, result;
flint_rand_t state;
printf("derivative....");
fflush(stdout);
flint_randinit(state);
/* Check aliasing */
for (i = 0; i < 10000; i++)
{
fmpz_t p;
fmpz_mod_poly_t a, b, c;
fmpz_init(p);
fmpz_randtest_unsigned(p, state, 2 * FLINT_BITS);
fmpz_add_ui(p, p, 2);
fmpz_mod_poly_init(a, p);
fmpz_mod_poly_init(b, p);
fmpz_mod_poly_init(c, p);
fmpz_mod_poly_randtest(a, state, n_randint(state, 100));
fmpz_mod_poly_set(b, a);
fmpz_mod_poly_derivative(c, b);
fmpz_mod_poly_derivative(b, b);
result = (fmpz_mod_poly_equal(b, c));
if (!result)
{
printf("FAIL (alias):\n");
fmpz_mod_poly_print(a), printf("\n\n");
fmpz_mod_poly_print(b), printf("\n\n");
fmpz_mod_poly_print(c), printf("\n\n");
abort();
}
fmpz_mod_poly_clear(a);
fmpz_mod_poly_clear(b);
fmpz_mod_poly_clear(c);
fmpz_clear(p);
}
/* Check constants have derivative zero */
for (i = 0; i < 10000; i++)
{
fmpz_t p;
fmpz_mod_poly_t a, b;
fmpz_init(p);
fmpz_randtest_unsigned(p, state, 2 * FLINT_BITS);
fmpz_add_ui(p, p, 2);
fmpz_mod_poly_init(a, p);
fmpz_mod_poly_init(b, p);
fmpz_mod_poly_randtest(a, state, n_randint(state, 2));
fmpz_mod_poly_derivative(b, a);
result = (b->length == 0);
if (!result)
{
printf("FAIL (da == 0):\n");
fmpz_mod_poly_print(a), printf("\n\n");
fmpz_mod_poly_print(b), printf("\n\n");
abort();
}
fmpz_mod_poly_clear(a);
fmpz_mod_poly_clear(b);
fmpz_clear(p);
}
/* Check (f g)' = f' g + f g' */
for (i = 0; i < 10000; i++)
{
fmpz_t p;
fmpz_mod_poly_t a, b, c, d, lhs, rhs;
fmpz_init(p);
fmpz_randtest_unsigned(p, state, 2 * FLINT_BITS);
fmpz_add_ui(p, p, 2);
fmpz_mod_poly_init(a, p);
fmpz_mod_poly_init(b, p);
fmpz_mod_poly_init(c, p);
fmpz_mod_poly_init(d, p);
fmpz_mod_poly_init(lhs, p);
fmpz_mod_poly_init(rhs, p);
fmpz_mod_poly_randtest(a, state, n_randint(state, 100));
fmpz_mod_poly_randtest(b, state, n_randint(state, 100));
fmpz_mod_poly_mul(lhs, a, b);
fmpz_mod_poly_derivative(lhs, lhs);
fmpz_mod_poly_derivative(c, a);
fmpz_mod_poly_derivative(d, b);
fmpz_mod_poly_mul(c, c, b);
fmpz_mod_poly_mul(d, a, d);
fmpz_mod_poly_add(rhs, c, d);
result = fmpz_mod_poly_equal(lhs, rhs);
if (!result)
{
printf("FAIL (Leibniz):\n");
printf("a = "), fmpz_mod_poly_print(a), printf("\n\n");
printf("b = "), fmpz_mod_poly_print(b), printf("\n\n");
printf("(ab)' = "), fmpz_mod_poly_print(lhs), printf("\n\n");
printf("a'b + ab' = "), fmpz_mod_poly_print(rhs), printf("\n\n");
abort();
}
fmpz_mod_poly_clear(a);
fmpz_mod_poly_clear(b);
fmpz_mod_poly_clear(c);
fmpz_mod_poly_clear(d);
fmpz_mod_poly_clear(lhs);
fmpz_mod_poly_clear(rhs);
fmpz_clear(p);
}
flint_randclear(state);
_fmpz_cleanup();
printf("PASS\n");
return 0;
}
int
main(void)
{
int i, result;
flint_rand_t state;
printf("div_basecase....");
fflush(stdout);
flint_randinit(state);
/* Compare to divrem_basecase */
for (i = 0; i < 5000; i++)
{
fmpz_t p;
fmpz_mod_poly_t a, b, q, q2, r2;
fmpz_init(p);
fmpz_randtest_unsigned(p, state, 2 * FLINT_BITS);
fmpz_add_ui(p, p, 2);
fmpz_mod_poly_init(a, p);
fmpz_mod_poly_init(b, p);
fmpz_mod_poly_init(q, p);
fmpz_mod_poly_init(q2, p);
fmpz_mod_poly_init(r2, p);
fmpz_mod_poly_randtest(a, state, n_randint(state, 100));
fmpz_mod_poly_randtest_not_zero(b, state, n_randint(state, 100) + 1);
{
fmpz_t d;
fmpz *leadB = fmpz_mod_poly_lead(b);
fmpz_init(d);
fmpz_gcd(d, p, leadB);
while (!fmpz_is_one(d))
{
fmpz_divexact(leadB, leadB, d);
fmpz_gcd(d, p, leadB);
}
fmpz_clear(d);
}
fmpz_mod_poly_div_basecase(q, a, b);
fmpz_mod_poly_divrem_basecase(q2, r2, a, b);
result = (fmpz_mod_poly_equal(q, q2));
if (!result)
{
printf("FAIL:\n");
printf("p = "), fmpz_print(p), printf("\n\n");
printf("a = "), fmpz_mod_poly_print(a), printf("\n\n");
printf("b = "), fmpz_mod_poly_print(b), printf("\n\n");
printf("q = "), fmpz_mod_poly_print(q), printf("\n\n");
printf("q2 = "), fmpz_mod_poly_print(q2), printf("\n\n");
printf("r2 = "), fmpz_mod_poly_print(r2), printf("\n\n");
abort();
}
fmpz_mod_poly_clear(a);
fmpz_mod_poly_clear(b);
fmpz_mod_poly_clear(q);
fmpz_mod_poly_clear(q2);
fmpz_mod_poly_clear(r2);
fmpz_clear(p);
}
/* Alias a and q */
for (i = 0; i < 500; i++)
{
fmpz_t p;
fmpz_mod_poly_t a, b, q;
fmpz_init(p);
fmpz_randtest_unsigned(p, state, 2 * FLINT_BITS);
fmpz_add_ui(p, p, 2);
fmpz_mod_poly_init(a, p);
fmpz_mod_poly_init(b, p);
fmpz_mod_poly_init(q, p);
fmpz_mod_poly_randtest(a, state, n_randint(state, 100));
fmpz_mod_poly_randtest_not_zero(b, state, n_randint(state, 100) + 1);
{
fmpz_t d;
fmpz *leadB = fmpz_mod_poly_lead(b);
fmpz_init(d);
fmpz_gcd(d, p, leadB);
while (!fmpz_is_one(d))
{
fmpz_divexact(leadB, leadB, d);
fmpz_gcd(d, p, leadB);
}
fmpz_clear(d);
}
fmpz_mod_poly_div_basecase(q, a, b);
fmpz_mod_poly_div_basecase(a, a, b);
result = (fmpz_mod_poly_equal(q, a));
if (!result)
{
printf("FAIL:\n");
printf("p = "), fmpz_print(p), printf("\n\n");
printf("a = "), fmpz_mod_poly_print(a), printf("\n\n");
printf("b = "), fmpz_mod_poly_print(b), printf("\n\n");
printf("q = "), fmpz_mod_poly_print(q), printf("\n\n");
abort();
}
fmpz_mod_poly_clear(a);
fmpz_mod_poly_clear(b);
fmpz_mod_poly_clear(q);
fmpz_clear(p);
}
/* Alias b and q */
for (i = 0; i < 500; i++)
{
fmpz_t p;
fmpz_mod_poly_t a, b, q;
fmpz_init(p);
fmpz_randtest_unsigned(p, state, 2 * FLINT_BITS);
fmpz_add_ui(p, p, 2);
fmpz_mod_poly_init(a, p);
fmpz_mod_poly_init(b, p);
fmpz_mod_poly_init(q, p);
fmpz_mod_poly_randtest(a, state, n_randint(state, 100));
fmpz_mod_poly_randtest_not_zero(b, state, n_randint(state, 100) + 1);
{
fmpz_t d;
fmpz *leadB = fmpz_mod_poly_lead(b);
fmpz_init(d);
fmpz_gcd(d, p, leadB);
while (!fmpz_is_one(d))
{
fmpz_divexact(leadB, leadB, d);
fmpz_gcd(d, p, leadB);
}
fmpz_clear(d);
}
fmpz_mod_poly_div_basecase(q, a, b);
fmpz_mod_poly_div_basecase(b, a, b);
result = (fmpz_mod_poly_equal(q, b));
if (!result)
{
printf("FAIL:\n");
printf("p = "), fmpz_print(p), printf("\n\n");
printf("a = "), fmpz_mod_poly_print(a), printf("\n\n");
printf("b = "), fmpz_mod_poly_print(b), printf("\n\n");
printf("q = "), fmpz_mod_poly_print(q), printf("\n\n");
abort();
}
fmpz_mod_poly_clear(a);
fmpz_mod_poly_clear(b);
fmpz_mod_poly_clear(q);
fmpz_clear(p);
}
flint_randclear(state);
_fmpz_cleanup();
printf("PASS\n");
return 0;
}
int
main(void)
{
int i, result;
flint_rand_t state;
printf("evaluate_fmpz....");
fflush(stdout);
flint_randinit(state);
/* Check aliasing */
for (i = 0; i < 10000; i++)
{
fmpz_t a, b, p;
fmpz_mod_poly_t f;
fmpz_init(p);
fmpz_randtest_unsigned(p, state, 2 * FLINT_BITS);
fmpz_add_ui(p, p, 2);
fmpz_init(a);
fmpz_init(b);
fmpz_randm(a, state, p);
fmpz_mod_poly_init(f, p);
fmpz_mod_poly_randtest(f, state, n_randint(state, 100));
fmpz_mod_poly_evaluate_fmpz(b, f, a);
fmpz_mod_poly_evaluate_fmpz(a, f, a);
result = (fmpz_equal(a, b));
if (!result)
{
printf("FAIL:\n");
fmpz_mod_poly_print(f), printf("\n\n");
fmpz_print(a), printf("\n\n");
fmpz_print(b), printf("\n\n");
abort();
}
fmpz_mod_poly_clear(f);
fmpz_clear(a);
fmpz_clear(b);
fmpz_clear(p);
}
/* Check that the result agrees with Z[X] */
for (i = 0; i < 10000; i++)
{
fmpz_t a, b, c, p;
fmpz_mod_poly_t f;
fmpz_poly_t g;
fmpz_init(p);
fmpz_randtest_unsigned(p, state, 2 * FLINT_BITS);
fmpz_add_ui(p, p, 2);
fmpz_init(a);
fmpz_init(b);
fmpz_init(c);
fmpz_mod_poly_init(f, p);
fmpz_poly_init(g);
fmpz_randm(a, state, p);
fmpz_mod_poly_randtest(f, state, n_randint(state, 100));
fmpz_mod_poly_get_fmpz_poly(g, f);
fmpz_mod_poly_evaluate_fmpz(b, f, a);
fmpz_poly_evaluate_fmpz(c, g, a);
fmpz_mod(c, c, p);
result = (fmpz_equal(b, c));
if (!result)
{
printf("FAIL (cmp with fmpz_poly):\n");
fmpz_mod_poly_print(f), printf("\n\n");
fmpz_poly_print(g), printf("\n\n");
fmpz_print(a), printf("\n\n");
fmpz_print(b), printf("\n\n");
fmpz_print(c), printf("\n\n");
abort();
}
fmpz_clear(a);
fmpz_clear(b);
fmpz_clear(c);
fmpz_clear(p);
fmpz_mod_poly_clear(f);
fmpz_poly_clear(g);
}
flint_randclear(state);
_fmpz_cleanup();
printf("PASS\n");
return 0;
}
int
main(void)
{
int i;
flint_rand_t state;
printf("init/init2/realloc/clear....");
fflush(stdout);
flint_randinit(state);
for (i = 0; i < 10000; i++)
{
fmpz_t p;
fmpz_mod_poly_t a;
fmpz_init(p);
fmpz_randtest_unsigned(p, state, 2 * FLINT_BITS);
fmpz_add_ui(p, p, 2);
fmpz_mod_poly_init2(a, p, n_randint(state, 100));
fmpz_mod_poly_clear(a);
fmpz_clear(p);
}
for (i = 0; i < 10000; i++)
{
fmpz_t p;
fmpz_mod_poly_t a;
fmpz_init(p);
fmpz_randtest_unsigned(p, state, 2 * FLINT_BITS);
fmpz_add_ui(p, p, 2);
fmpz_mod_poly_init2(a, p, n_randint(state, 100));
fmpz_mod_poly_realloc(a, n_randint(state, 100));
fmpz_mod_poly_clear(a);
fmpz_clear(p);
}
for (i = 0; i < 10000; i++)
{
fmpz_t p;
fmpz_mod_poly_t a;
fmpz_init(p);
fmpz_randtest_unsigned(p, state, 2 * FLINT_BITS);
fmpz_add_ui(p, p, 2);
fmpz_mod_poly_init(a, p);
fmpz_mod_poly_randtest(a, state, n_randint(state, 100));
fmpz_mod_poly_clear(a);
fmpz_clear(p);
}
flint_randclear(state);
_fmpz_cleanup();
printf("PASS\n");
return 0;
}
