See also https://github.com/timholy/PositiveFactorizations.jl

--Tim

On Thursday, June 9, 2016 7:56:41 PM CDT Viral Shah wrote:
> 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


Reply via email to