and use functions from Matrix to reset covariances.