On trivariate blending sums of univariate and bivariate quadratic spline quasi-interpolants on bounded domains