Re: [eigen] Instability in LLT and LDLT methods.

[ Thread Index | Date Index | More Archives ]


yes, it seems the test to check whether the matrix is positive
definite was too strict. I changed the absolute tolerance a bit, but
we still need something better. Basically, in Cholesky we compute at
each iteration 1/sqrt(x), and so x must be >0 with some epsilon...


On Tue, Jan 27, 2009 at 7:35 PM, Keir Mierle <mierle@xxxxxxxxx> wrote:
> Here is a testcase that fails with LLT and LDLT but works fine with
> all of LU, SVD, and QR solving. Depends on my previous patch for QR
> solver (or comment out the qr().solve line).
> Keir

Mail converted by MHonArc 2.6.19+