Cubature formulas based on bivariate quadratic spline quasi-interpolants