MUL_PQR_XP
MUL(mul[MUL_PQR_XP], REC_PQR_X); /* Xp = Pxyz * xp */
coeff[MUL_PQR_XP] = gf_div(nyyz ^ nyzz, xd);