hessian warning

No replies
carey's picture
Offline
Joined: 10/19/2009

congrats on getting the correct gradient and calculated hessian implemented!

suggest a warning to the user when the calculated hessian is not positive definite, something that should be easy to implement because the you invert the hessian to get the standard errors. see attached code which shows that openmx can appear to have converged when the model is unidentified.

also a warning when a gradient element is exactly 0 can help a user diagnose a problem.

greg