New function covariance_calculate_unnormalized