MUL_PQR_XQ
MUL(mul[MUL_PQR_XQ], REC_PQR_Y); /* Xq = Qxyz * xq */
coeff[MUL_PQR_XQ] = gf_div(nyy ^ nzz, xd);