clmulr
tmp = clmulr(msgpoly, consts->barrett_reduction_const_1);
return clmulr(tmp, consts->barrett_reduction_const_2);