Implement faster covariance matrix sampling