It would be good to specify in the code the type of covariance estimator used by the Normality class.
Further in the future it would be good to add the capability to switch types of covariance estimator to the Normality reconciliator with a covariance_type input.