void
fmpz_poly_mat_mullow(fmpz_poly_mat_t C, const fmpz_poly_mat_t A,
const fmpz_poly_mat_t B, long len)
{
long ar, bc, br;
long i, j, k;
fmpz_poly_t t;
ar = A->r;
br = B->r;
bc = B->c;
if (br == 0 || len < 1)
{
fmpz_poly_mat_zero(C);
return;
}
if (C == A || C == B)
{
fmpz_poly_mat_t T;
fmpz_poly_mat_init(T, ar, bc);
fmpz_poly_mat_mullow(T, A, B, len);
fmpz_poly_mat_swap(C, T);
fmpz_poly_mat_clear(T);
return;
}
fmpz_poly_init(t);
for (i = 0; i < ar; i++)
{
for (j = 0; j < bc; j++)
{
fmpz_poly_mullow(fmpz_poly_mat_entry(C, i, j),
fmpz_poly_mat_entry(A, i, 0),
fmpz_poly_mat_entry(B, 0, j), len);
for (k = 1; k < br; k++)
{
fmpz_poly_mullow(t, fmpz_poly_mat_entry(A, i, k),
fmpz_poly_mat_entry(B, k, j), len);
fmpz_poly_add(fmpz_poly_mat_entry(C, i, j),
fmpz_poly_mat_entry(C, i, j), t);
}
}
}
fmpz_poly_clear(t);
}
void gmde_convert_soln(fmpz_poly_mat_t A, long *vA,
const padic_mat_struct *C, long N, const fmpz_t p)
{
long i, j, k;
fmpz_t s, t;
assert(N > 0);
assert(A->r == padic_mat(C)->r && A->c == padic_mat(C)->c);
fmpz_init(s);
fmpz_init(t);
/* Find valuation */
*vA = LONG_MAX;
for (k = 0; k < N; k++)
*vA = FLINT_MIN(*vA, padic_mat_val(C + k));
fmpz_poly_mat_zero(A);
for (k = N - 1; k >= 0; k--)
{
if (padic_mat_val(C + k) == *vA)
{
for (i = 0; i < A->r; i++)
for (j = 0; j < A->c; j++)
if (!fmpz_is_zero(padic_mat_entry(C + k, i, j)))
fmpz_poly_set_coeff_fmpz(fmpz_poly_mat_entry(A, i, j), k,
padic_mat_entry(C + k, i, j));
}
else
{
fmpz_pow_ui(s, p, padic_mat_val(C + k) - *vA);
for (i = 0; i < A->r; i++)
for (j = 0; j < A->c; j++)
if (!fmpz_is_zero(padic_mat_entry(C + k, i, j)))
{
fmpz_mul(t, s, padic_mat_entry(C + k, i, j));
fmpz_poly_set_coeff_fmpz(fmpz_poly_mat_entry(A, i, j),
k, t);
}
}
}
fmpz_clear(s);
fmpz_clear(t);
}
void
fmpz_poly_mat_sqrlow(fmpz_poly_mat_t B, const fmpz_poly_mat_t A, slong len)
{
slong n = A->r;
if (n == 0)
return;
if (len < 1)
{
fmpz_poly_mat_zero(B);
return;
}
if (n == 1)
{
fmpz_poly_sqrlow(E(B, 0, 0), E(A, 0, 0), len);
return;
}
if (n == 2)
{
fmpz_poly_t t, u;
fmpz_poly_init(t);
fmpz_poly_init(u);
fmpz_poly_addlow(t, E(A, 0, 0), E(A, 1, 1), len);
fmpz_poly_mullow(u, E(A, 0, 1), E(A, 1, 0), len);
fmpz_poly_sqrlow(E(B, 0, 0), E(A, 0, 0), len);
fmpz_poly_addlow(E(B, 0, 0), E(B, 0, 0), u, len);
fmpz_poly_sqrlow(E(B, 1, 1), E(A, 1, 1), len);
fmpz_poly_addlow(E(B, 1, 1), E(B, 1, 1), u, len);
fmpz_poly_mullow(E(B, 0, 1), E(A, 0, 1), t, len);
fmpz_poly_mullow(E(B, 1, 0), E(A, 1, 0), t, len);
fmpz_poly_clear(t);
fmpz_poly_clear(u);
return;
}
fmpz_poly_mat_mullow(B, A, A, len);
}
