In Matlab to check if a symmetric sparse matrix is positive definite, I can say [R,p]=chol(A) and then if p>0 etc. Is this functionality available in Julia? The cholfact standard routine throws an exception if its argument is not positive definite rather than returning any helpful information.
I looked at the code for cholfact in cholmod.jl in Base; it appears that I can write a modified version of cholfact that exposes this functionality. But it would be better if the functionality were available in the library so that my code is not sensitive to changes in undocumented low-level routines. Thanks, Steve Vavasis
