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
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);
}
