| Re: Crash on 1x1 matrix (was Re: [eigen] 2.0.1 soon?) |
[ Thread Index |
Date Index
| More lists.tuxfamily.org/eigen Archives
]
- To: eigen@xxxxxxxxxxxxxxxxxxx
- Subject: Re: Crash on 1x1 matrix (was Re: [eigen] 2.0.1 soon?)
- From: Gael Guennebaud <gael.guennebaud@xxxxxxxxx>
- Date: Mon, 23 Mar 2009 16:02:07 +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=riNIr2wcL67BWFTGeHz47/oA7CZWhIvoyv29fp2w/rg=; b=wZsh7pP05XRGHRNw2LvYjMAgLYlBOt+qlYDC6Ax39YZh2i9+e25ByS4yFufVHE1stk EtjwFghHlxb6k7V29dtr0BSwsZE1LBsPAwHtBI1ktPwdFiopG4Coib5FAgzWlghA4Xjn Hgx95ZHv/KsFC1f13b0KSkGjXJvn9xJZg6RqQ=
- 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=tmc/3AHtC2DT+HwGKnscRsqSGi/iA2wxi9L0rwQK3XgY2axJ+uUAkb9byptjko/dN9 v1wrkmWGjDh6aB40dzmvEtf0mJV4M5c8zRk6jl4M5x2mTDVHwMzby6JNthLWdAxQKe6Y ccPavrC4F9IKEuUgDIdDeX4Il54Iskq913Otk=
should be fixed in trunk. backporting it to the 2.0.x branch is another story...
On Mon, Mar 23, 2009 at 1:43 PM, Robert Lupton the Good
<rhl@xxxxxxxxxxxxxxxxxxx> wrote:
>> <rhl@xxxxxxxxxxxxxxxxxxx> wrote:
>>>
>>> There are some problems with 0- and 1- sized matrices that I'd
>>> like to get reported before 2.0.1! eigen shouldn't crash on the
>>> trivial cases.
>>
>> can you be more specific ? I don't remember we still have such problems.
>
> I was going to! Here you go; this crashes with
>
> Assertion failed: (dim > 0), function Matrix, file
> /u/lsst/products/external/eigen/2.0.0/include/Eigen/src/Core/Matrix.h, line
> 312.
>
> And no, nImage isn't really a const int in the application, so I couldn't
> use
> Matrix1d even if it existed!
>
> R
>
> #include "Eigen/Core"
> #include "Eigen/QR"
>
> int main() {
> const int nImage = 1;
>
> Eigen::MatrixXd R(nImage, nImage);
>
> R(0, 0) = 1.0;
>
> Eigen::SelfAdjointEigenSolver<Eigen::MatrixXd> eVecValues(R);
> }
>
>
>
>