mult_q
term = mult_q(term, div_q(x, 1 << PSC_Q_FRACTION));
term = mult_q(term, div_q(x, 2 << PSC_Q_FRACTION));
term = mult_q(term, div_q(x, 3 << PSC_Q_FRACTION));
term = mult_q(term, div_q(x, 4 << PSC_Q_FRACTION));
mult_q(sigma_q, sigma_q)));
g2_q = mult_q(g1_q, g1_q);
g0_q = mult_q(g0_q, g1_q);
g1_q = mult_q(g1_q, g2_q);