int
main(void)
{
int i, result;
flint_rand_t state;
printf("is_zero....");
fflush(stdout);
flint_randinit(state);
/* Check zero vector */
for (i = 0; i < 10000; i++)
{
fmpz *a;
long len = n_randint(state, 100);
a = _fmpz_vec_init(len);
_fmpz_vec_randtest(a, state, len, 200);
_fmpz_vec_zero(a, len);
result = (_fmpz_vec_is_zero(a, len));
if (!result)
{
printf("FAIL1:\n");
_fmpz_vec_print(a, len), printf("\n\n");
abort();
}
_fmpz_vec_clear(a, len);
}
/* Check non-zero vector */
for (i = 0; i < 10000; i++)
{
fmpz *a;
long len = n_randint(state, 100) + 1;
a = _fmpz_vec_init(len);
_fmpz_vec_randtest(a, state, len, 200);
fmpz_set_ui(a + (len - 1), 1UL);
result = (!_fmpz_vec_is_zero(a, len));
if (!result)
{
printf("FAIL2:\n");
_fmpz_vec_print(a, len), printf("\n\n");
abort();
}
_fmpz_vec_clear(a, len);
}
flint_randclear(state);
_fmpz_cleanup();
printf("PASS\n");
return 0;
}
void _qadic_exp_balanced(fmpz *rop, const fmpz *x, slong v, slong len,
const fmpz *a, const slong *j, slong lena,
const fmpz_t p, slong N, const fmpz_t pN)
{
const slong d = j[lena - 1];
fmpz_t pw;
fmpz *r, *s, *t;
slong i, w;
r = _fmpz_vec_init(d);
s = _fmpz_vec_init(2*d - 1);
t = _fmpz_vec_init(d);
fmpz_init(pw);
fmpz_pow_ui(pw, p, v);
_fmpz_vec_scalar_mul_fmpz(t, x, len, pw);
_fmpz_vec_scalar_mod_fmpz(t, t, len, pN);
_fmpz_vec_zero(t + len, d - len);
fmpz_set(pw, p);
fmpz_one(rop + 0);
_fmpz_vec_zero(rop + 1, d - 1);
w = 1;
while (!_fmpz_vec_is_zero(t, d))
{
fmpz_mul(pw, pw, pw);
for (i = 0; i < d; i++)
{
fmpz_fdiv_r(r + i, t + i, pw);
fmpz_sub(t + i, t + i, r + i);
}
if (!_fmpz_vec_is_zero(r, d))
{
_qadic_exp_bsplit(r, r, w, d, a, j, lena, p, N);
_fmpz_poly_mul(s, rop, d, r, d);
_fmpz_poly_reduce(s, 2*d - 1, a, j, lena);
_fmpz_vec_scalar_mod_fmpz(rop, s, d, pN);
}
w *= 2;
}
_fmpz_vec_clear(r, d);
_fmpz_vec_clear(s, 2*d - 1);
_fmpz_vec_clear(t, d);
fmpz_clear(pw);
}
void _fmpq_poly_scalar_div_mpq(fmpz * rpoly, fmpz_t rden,
const fmpz * poly, const fmpz_t den, long len,
const fmpz_t r, const fmpz_t s)
{
fmpz_t gcd1; /* GCD( poly, r ) */
fmpz_t gcd2; /* GCD( s, den ) */
fmpz_init(gcd1);
fmpz_init(gcd2);
fmpz_set_ui(gcd1, 1);
fmpz_set_ui(gcd2, 1);
if (*r != 1L)
{
_fmpz_vec_content(gcd1, poly, len);
if (*gcd1 != 1L)
fmpz_gcd(gcd1, gcd1, r);
}
if (*den != 1L && *s != 1L)
fmpz_gcd(gcd2, s, den);
if (*gcd1 == 1L)
{
if (*gcd2 == 1L)
{
_fmpz_vec_scalar_mul_fmpz(rpoly, poly, len, s);
fmpz_mul(rden, den, r);
}
else
{
fmpz_t s2;
fmpz_init(s2);
fmpz_divexact(s2, s, gcd2);
_fmpz_vec_scalar_mul_fmpz(rpoly, poly, len, s2);
fmpz_divexact(rden, den, gcd2);
fmpz_mul(rden, rden, r);
fmpz_clear(s2);
}
}
else
{
fmpz_t r2;
fmpz_init(r2);
fmpz_divexact(r2, r, gcd1);
if (*gcd2 == 1L)
{
_fmpz_vec_scalar_divexact_fmpz(rpoly, poly, len, gcd1);
_fmpz_vec_scalar_mul_fmpz(rpoly, rpoly, len, s);
fmpz_mul(rden, den, r2);
}
else
{
fmpz_t s2;
fmpz_init(s2);
fmpz_divexact(s2, s, gcd2);
_fmpz_vec_scalar_divexact_fmpz(rpoly, poly, len, gcd1);
_fmpz_vec_scalar_mul_fmpz(rpoly, rpoly, len, s2);
fmpz_divexact(rden, den, gcd2);
fmpz_mul(rden, rden, r2);
fmpz_clear(s2);
}
fmpz_clear(r2);
}
if (_fmpz_vec_is_zero(rpoly, len))
fmpz_set_ui(rden, 1);
if (fmpz_sgn(rden) < 0)
{
_fmpz_vec_neg(rpoly, rpoly, len);
fmpz_neg(rden, rden);
}
fmpz_clear(gcd1);
fmpz_clear(gcd2);
}