I have contributed several functions to the R package
refund: Regression with Functional
Data. In particular, the R function fbps provides the sandwich smoother for the fast bivariate penalized smoothing of matrix data; the R function fpca.face provides fast covariance estimation of
functional data. An important feature of these functions is that they are fast and scalable for high dimensional data.