void nmod_poly_mat_right_tmul(nmod_poly_mat_t C, const nmod_poly_mat_t L, const nmod_poly_mat_t Bin, long m, long k){
long qq = nmod_poly_mat_nrows(Bin);
long rr = nmod_poly_mat_ncols(Bin);
long pp = nmod_poly_mat_nrows(L);
nmod_poly_mat_t B;
nmod_poly_mat_init(B, rr, qq, nmod_poly_mat_modulus(Bin));
long i, j;
for (i = 0; i < rr; i++)
for (j = 0; j < qq; j++)
nmod_poly_reverse(nmod_poly_mat_entry(B, i, j), nmod_poly_mat_entry(Bin, j, i), m+1);
nmod_poly_mat_mul(C, L, B);
for (i = 0; i < pp; i++)
for (j = 0; j < qq; j++){
nmod_poly_truncate(nmod_poly_mat_entry(C, i, j), k+m);
nmod_poly_shift_right(nmod_poly_mat_entry(C, i, j), nmod_poly_mat_entry(C, i, j), m);
}
nmod_poly_mat_clear(B);
}
int main()
{
setbuf(stdout, NULL);
for (long l = 5; l <= 60; l += 5) {
// for (long n = 1024; n <= 65536; n *= 2) {
for (long idx = 0; idx < 13; idx ++) {
long n = 1024*(1L << idx/2);
if (idx & 1) n += n/2;
SetSeed((ZZ(l) << 64) + ZZ(n));
long p;
RandomLen(p, l);
if (p % 2 == 0) p++;
zz_p::init(p);
zz_pX a, b, c, f;
random(a, n);
random(b, n);
random(f, n);
SetCoeff(f, n);
zz_pXModulus F(f);
zz_pXMultiplier B(b, F);
double t;
MulMod(c, a, B, F);
long iter = 1;
do {
t = GetTime();
for (long i = 0; i < iter; i++) MulMod(c, a, B, F);
t = GetTime() - t;
iter *= 2;
} while (t < 3);
iter /= 2;
t = GetTime();
for (long i = 0; i < iter; i++) MulMod(c, a, B, F);
t = GetTime()-t;
double NTLTime = t;
Flintzz_pX f_a(a), f_b(b), f_c, f_f(f), f_finv;
nmod_poly_reverse(f_finv.value, f_f.value, f_f.value->length);
nmod_poly_inv_series_newton(f_finv.value, f_finv.value, f_f.value->length);
nmod_poly_mulmod_preinv(f_c.value, f_a.value, f_b.value, f_f.value, f_finv.value);
t = GetTime();
for (long i = 0; i < iter; i++)
nmod_poly_mulmod_preinv(f_c.value, f_a.value, f_b.value, f_f.value, f_finv.value);
t = GetTime()-t;
double FlintTime = t;
printf("%8.2f", FlintTime/NTLTime);
}
printf("\n");
}
}