Convergence of an efficient local least-squares fitting method for bases with compact support

An efficient local least-squares algorithm for non-orthogonal projection onto smooth piecewise-polynomial basis functions is analyzed. The algorithm runs in optimal time and delivers the same order of accuracy as the standard orthogonal projection. Numerical results indicate that in many computational situations, the new algorithm offers an effective alternative to global least-squares approximation.