Fq_freeze
for (i = 0; i < p; ++i) h[i] = Fq_freeze(fg[i]);
for (i = 0; i < p; ++i) h[i] = Fq_freeze(3 * f[i]);
ai = Fq_freeze(a1 * (int32_t)ai);
for (i = 0; i < p + 1; ++i) g[i] = Fq_freeze(f0 * g[i] - g0 * f[i]);
for (i = 0; i < p + 1; ++i) r[i] = Fq_freeze(f0 * r[i] - g0 * v[i]);
for (i = 0; i < p; ++i) out[i] = Fq_freeze(scale * (int32_t)v[p - 1 - i]);