$L^2 $-Approximation Theory of Even Order Multivariate Splines