covariance_accumulate: New function for one-pass computation of the