solverMLE(G,U,ZeroTolerance=>n)
After computing the variety of the zero-dimensional ideal of the score equations, solverMLE needs to determine which points lie in the PD cone. A matrix is assumed to be positive definite if for all eigenvalues e:
- realPart e > ZeroTolerance
- abs imaginaryPart e <= ZeroTolerance
If an MLE matrix contains an imaginary part below the tolerance level, then only the real part is reported in the output. (See checkPD(...,ZeroTolerance=>...))