Local quartic C^2 spline quasi-interpolation on 3D bounded domains