Use Rcpp for covariance operations