chol fails on Hermitian positive definite matrix?

1 view (last 30 days)
I get the following code output on the attached matrix:
K>> max(max(abs(A-A')))
ans =
0
K>> min(eig(A))
ans =
0.0900
K>> chol(A)
Error using chol
Matrix must be positive definite.

Accepted Answer

Christine Tobler
Christine Tobler on 20 May 2019
That is astonishing. The matrix is not badly conditioned, I would expect CHOL to work fine here.
I tried on my machine and didn't get an error message from CHOL, so this seems related to your version or machine. Could you tell me what MATLAB version and what OS you're using?
  2 Comments
Marco Duarte
Marco Duarte on 20 May 2019
K>> version
ans =
8.3.0.532 (R2014a)
OS is MacOS High Sierra 10.13.6
I tried the same operation with a newer version of Matlab (2018b) and chol(A) works. So it seems to be specific to R2014a.
Christine Tobler
Christine Tobler on 22 May 2019
I could reproduce this bug for MATLAB R2014a on a Mac. The bug did not appear on MATLAB R2014b on that Mac, and neither did it happen on MATLAB R2014a on Linux. This seems to be a version and OS-specific bug.

Sign in to comment.

More Answers (0)

Categories

Find more on Linear Algebra in Help Center and File Exchange

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!