Smoothing splines via the penalized least squares method provide versatile and effective nonparametric models for regression with Gaussian responses. The computation of smoothing splines is generally of the order "O"("n"-super-3), "n" being the sample size, which severely limits its practical...