| 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: "Gael Guennebaud" <gael.guennebaud@xxxxxxxxx>
- Date: Fri, 22 Aug 2008 19:49:11 +0200
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:received:received:message-id:date:from:to :subject:in-reply-to:mime-version:content-type :content-transfer-encoding:content-disposition:references; bh=wp9T6otHci2ZjD0/bCH9DnmjgwotryUbTt/UEOk47O0=; b=DLUqBlDeix8tKMtrn5N0egl4dEsi9eqJ0vGYrbyAnbj7BnNcPxZPRJo69HsNSa7Khl ntDeaGWaD7nZvQZXDhMYnSqs9Wg0ekeXKZSZEhCW1Pppk02T+rLyd0kZwTMi/EqCi6pP FWfewnAvxztKrD85VsN1HwGhnA2K8MLXI6+kc=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=message-id:date:from:to:subject:in-reply-to:mime-version :content-type:content-transfer-encoding:content-disposition :references; b=QKLPDRnEkBKP6oc47R2be07KiisPweTQq918pj8qg2OBcUv4qmhv8E0fQ/95Ce67yz 3zY65wBD0jXElLmp04xd2ceDb3kKVNTKrrngoftoYToerUR5fUnEU0dpyb8lgb9Vzhne t3kxF04w+f8d5Dg1mveQWgqpJ16FOi/VPojp4=
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)