To maximize the probability of drawing the correct error given the estimated Gaussian parameters is equivalent to minimize the negative full log-likelihood function [1]}.
To ensure a positive definite covariance matrix we use the LDL decomposition as in [2]}.
