add Covariance1DFunction
add CovarianceMatrix
split operator() argument in two