Let's assume that I have a symmetric matrix $A$. What's the fastest way to find its nearest positive definite matrix in Mathematica?

Any help would be appreciated.

  • 2
    $\begingroup$ Diagonalize, zero out negative values on the diagonal, reverse, and you have the closest positive semidefinite matrix. $\endgroup$ – Daniel Lichtblau Aug 3 '17 at 18:57
  • 2
    $\begingroup$ @Anoldmaninthesea. Asking for a a positive definite matrix is like asking which number in the open interval (0, 1) is nearest to 2 $\endgroup$ – Coolwater Aug 3 '17 at 19:29
  • 3
    $\begingroup$ What people are trying to say is that there is no "nearest" PD matrix, only PSD. $\endgroup$ – Daniel Lichtblau Aug 3 '17 at 21:01
  • 2
    $\begingroup$ Positively definitely... $\endgroup$ – Daniel Lichtblau Nov 17 '18 at 22:53
  • 2
    $\begingroup$ By the way, it is Eigensystem that I had in mind, not Eigenvalues alone. $\endgroup$ – Daniel Lichtblau Nov 18 '18 at 14:16

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.