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} }