| Re: [eigen] request for comments / help |
[ Thread Index |
Date Index
| More lists.tuxfamily.org/eigen Archives
]
- To: eigen@xxxxxxxxxxxxxxxxxxx
- Subject: Re: [eigen] request for comments / help
- From: Gael Guennebaud <gael.guennebaud@xxxxxxxxx>
- Date: Wed, 13 Jun 2012 19:06:24 +0200
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=20120113; h=mime-version:in-reply-to:references:from:date:message-id:subject:to :content-type:content-transfer-encoding; bh=x9g1vDP4if94x8FUFpBRg3vqxGYCQ7lq5QJmaaxgkr0=; b=IyHoACySNT90QcRL63fczbIV2yyX8YtaUjojbWDTW3vhKXwB9bZzdpRwRG+n4BiF1f QFy19+YvzL9gl/4M/uatUJiXBD+UNnbgTCEJInUyw2MhfB2kfJUQNc5L7zyJT+WrJP5i m27qiEpoZCyD9VtRE1xlbMN4y2Ua3oqJg/E40OUReZ20D+/YOeHvljJmbfZqmChjq8CQ uwyoVXhP5Wib6selGhVZmrwCSczVCk+cH8LXuBg6pu/tbgRqV8pM84/j3MD4WYbcLd2O RPRMGlByx0GjTiwpRTVpbJc0blqiuaUMShCVC9SH0ZSo5YzGLlSa4/AhRB0u9GCvh3N8 b0dA==
Well, I have no doubt that in most cases QR-Cholesky works fine and is
pretty efficient, however Cholesky involves squaring thus reducing
numerical accuracy. So I'm tempted to favor a full orthogonalization
method. This is also the approach followed by Lapack (dgelsy).
Nevertheless your approach has the advantage of simplicity and being
already implemented! So if you're ok I propose in a first step to add
it in unsupported/, I'm sure it might be useful to many users.
Gael.
On Tue, Jun 12, 2012 at 8:49 AM, Helmut Jarausch
<jarausch@xxxxxxxxxxxxxxxxxxx> wrote:
> On 06/11/2012 06:31:22 PM, Gael Guennebaud wrote:
>>
>> On Mon, Jun 11, 2012 at 5:56 PM, Helmut Jarausch
>> <jarausch@xxxxxxxxxxxxxxxxxxx> wrote:
>> > I've attached the corrected version of QR_Cholesky.
>> > Is it impudent to suggest adding this to Eigen?
>> >
>> > I think it's the standard way of solving the minimum length least
>> > squares
>> > problem.
>>
>> I think it would be better to perform a full orthogonalization of the
>> form:
>>
>> A * P = Q * U * Z
>>
>> where U of the form:
>>
>> | U11 0 |
>> | 0 0 |
>>
>> This can be done by applying householder reflexions on the non zero
>> upper trapezoidal part of the matrix R.
>>
>
> I have been using that myself - more than 20 years ago.
> Recent references, like Peter Deuflhard's book on
> Numerical Methods of Nonlinear Problems, Springer 2004,
> recommendsthe use of QR-Cholesky (section 4.4.1).
> Since Peter Deuflhard and his group have a lot of experience in
> real life problems I am sure this suggestion has proven itself in practise.
>
> But if someone has an implementation of full orthogonalization we could do
> some benchmarking.
>
> Helmut.
>
>
>