Replace naïve covariance algorithm to avoid numerical instability and arithmetic overflow by the stable one from ...

[mbok]

@TechReport{pebay2008formulas,
title={Formulas for robust, one-pass parallel computation of covariances and arbitrary-order statistical moments.},
author={Pebay, Philippe Pierre},
year={2008},
institution={Sandia National Laboratories}
}