isposdef() or catch the exception from cholfact. -viral
On Thursday, June 9, 2016 at 10:16:29 PM UTC-4, [email protected] wrote: > > 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 > >
