Constructing Good Coefficient Functionals for Bivariate $C^1$ Quadratic Spline Quasi-Interpolants