Re: [eigen] Bug in the Cholesky/Core module? |
[ Thread Index |
Date Index
| More lists.tuxfamily.org/eigen Archives
]
- To: eigen@xxxxxxxxxxxxxxxxxxx
- Subject: Re: [eigen] Bug in the Cholesky/Core module?
- From: "Timothy Hunter" <tjhunter@xxxxxxxxxxxx>
- Date: Fri, 22 Aug 2008 10:58:20 -0700
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:received:received:message-id:date:from:sender :to:subject:in-reply-to:mime-version:content-type :content-transfer-encoding:content-disposition:references :x-google-sender-auth; bh=ij6WmcPp14m4Fjj+COdatGufaSmVXKxthbRIFCgJy3s=; b=IqHstgyKJ6RVmuDRHENHhI+1Pt7Zm2zlrTJl063IXvf7GeJE6u5bzuYTNBGEF3TVIJ MX+wsXGl7LwQ+WzW9mESjM2UF+D2DSvLLPUkxH/takHc01pxH8MoZplrLvJu9arer0Yx JzxlJgLbv1AmvOwcQXWslQnl/0Pgzw9vQNqN0=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=message-id:date:from:sender:to:subject:in-reply-to:mime-version :content-type:content-transfer-encoding:content-disposition :references:x-google-sender-auth; b=wz95Ub6o7xEXYG3lOh5oKKbyB1voukcNJiZq3uAAujzQ/ntqAH9Al4UTW2WTqDPgEC 8vBzUZXu1CRKzEZiO8I09cOJ7L5SGCbd1xH3kkSLdTcCKmoXJ1S25idWfugGAugysrli d42Sy+v6UrepuNqW4/fmslrrYJzrCk81FYhzs=
You are welcome. I am amazed by the rate at which bugs are fixed!
I am currently writing a convex optimization solver using Eigen. Do
you have any small tutorial/style recommendation on how to write code
for it? In my case, I want to write some matrix specializations
(tridiagonal, bock diagonal and block tridiagonal).
On Fri, Aug 22, 2008 at 10:49 AM, Gael Guennebaud
<gael.guennebaud@xxxxxxxxx> wrote:
> Hi,
>
> thanks for the report, bug fixed (rev 850982)
>
> gael.
>
>
> On Fri, Aug 22, 2008 at 7:34 PM, Timothy Hunter <tjhunter@xxxxxxxxxxxx> wrote:
>> Hi List,
>>
>> the following code:
>>
>> typedef Eigen::Matrix<float, 1, 1> MyMatrix;
>>
>> typedef Eigen::Matrix<float, 1, 1> MyVector;
>>
>> MyMatrix A1 = MyMatrix::Random();
>>
>> A1=(A1.transpose()*A1).eval();
>>
>> Eigen::Cholesky<MyMatrix> chol(A1);
>>
>> chol.solve(MyVector::Random());
>>
>> produces a 'you_mixed_matrices_of_different_sizes' error in Assign.h
>>
>> I know there is no need for the cholesky decomposition in this case, but I
>> had it as the result of deeply nested templates. Can anyone reproduce it?
>>
>> --
>>
>> Timothy Hunter
>>
>> Student (Stanford University)
>
>
>