MUL_PQR_XR
MUL(mul[MUL_PQR_XR], REC_PQR_Z); /* Xr = Rxyz * xr */
coeff[MUL_PQR_XR] = gf_div(ny ^ nz, xd);