Re: [eigen] Instability in LLT and LDLT methods. |
[ Thread Index |
Date Index
| More lists.tuxfamily.org/eigen Archives
]
- To: eigen@xxxxxxxxxxxxxxxxxxx
- Subject: Re: [eigen] Instability in LLT and LDLT methods.
- From: Gael Guennebaud <gael.guennebaud@xxxxxxxxx>
- Date: Wed, 28 Jan 2009 00:05:56 +0100
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:mime-version:received:in-reply-to:references :date:message-id:subject:from:to:content-type :content-transfer-encoding; bh=GIpMheCtrWM4Upq9rIKSqC1Ox8ftFsAVm5w0Z6QY1Ps=; b=UGRPOQIFPC6XazbK/xzwZ7774jg2wTaDPsZz82CFMUpn4jSqiFbCVZPHhmPqpR2J13 5f0DaNIi3+d4pORjAmDvckHU4H1HZrOK6slpNQu5pvEnNwojhXYFItVIl6QSmQr/uZ26 fsCvHNPHQ6dMBsMLxS3L87RUgYzF0HfSlx/Wo=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:in-reply-to:references:date:message-id:subject:from:to :content-type:content-transfer-encoding; b=UMhIVwE3z4Aq1gAWeClrqOLOEMKnUQQTP+lAz4R8PdCOr4iwXN+USZ0wkGAYN+tmpL BgTABHFTcS2fCDXjU6jkxS3TaJOdTZRDlNfRZrL7qxVLDs7lU0RsM0hmZbl9fzoVQNRi wzmLd9n0/JokpDpDOi5GGxd+oB/LjTJGr/crc=
Hi,
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...
Gael.
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
>