| Re: [eigen] Bug in cholesky module |
[ Thread Index |
Date Index
| More lists.tuxfamily.org/eigen Archives
]
- To: eigen@xxxxxxxxxxxxxxxxxxx
- Subject: Re: [eigen] Bug in cholesky module
- From: "Gael Guennebaud" <gael.guennebaud@xxxxxxxxx>
- Date: Sat, 16 Aug 2008 09:10: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=ruu9AUr9FOwrqn7nYqW9V4Z4L81hnCII/neNm9MR0oQ=; b=hJu8ajBboQ0HF5x5GUywc3tv0BxRxnxcv3M2ipYXXfh8Rf+vr+y0DZx62TPmNrkCYc JW1PiDpugxJigzUYzC0IdiZqWQRklj8NBj/Zy0EBAVPQehkq9LaqmeZUZg5jz++HJB3v D1hv5U+nmd+gqbtiH37FFBv1SkvWeoUo/sSnQ=
- 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=cASbI7+8nsbOR3ZaL53oo/jrcSqq7Zd538OzrmGeQz5UZg6D53CaPNhibir4a0mv06 CaZFgSgOiMYxbunVlWyRmu07M/QG9Ugo496aFLmGPBQqh0u1cGiDHDpACFzmqInZcdmU z3dsyh3KRFm6+C5E9/gi3Mm5hf0u4GFEI94E4=
Hi Timothy,
thanks for the report, bug fixed (in kdesupport)
gael.
On Sat, Aug 16, 2008 at 3:06 AM, Timothy Hunter <tjhunter@xxxxxxxxxxxx> wrote:
> Hi all,
>
> the CholeskyWithoutSquareRoot module does not work with the following code
> snippet:
>
> Eigen::Matrix4f m = Eigen::Matrix4f::Random();
>
> Eigen::CholeskyWithoutSquareRoot<Eigen::Matrix4f> chol =
> (m*m.transpose()).eval().choleskyNoSqrt();
>
> Eigen::Matrix<float, 4, 3> mat = Eigen::Matrix<float, 4, 3>::Random();
>
> chol.solve(mat);
>
> Error message:
>
> eigen2/Eigen/src/Core/CwiseBinaryOp.h:92: Eigen::CwiseBinaryOp<BinaryOp,
> Lhs, Rhs>::CwiseBinaryOp(const Lhs&, const Rhs&, const BinaryOp&) [with
> BinaryOp = Eigen::ei_scalar_quotient_op<float>, Lhs = Eigen::Matrix<float,
> 4, 3, 4, 3, 120u>, Rhs = Eigen::DiagonalCoeffs<Eigen::Matrix<float, 4, 4, 4,
> 4, 120u> >]: Assertion `lhs.rows() == rhs.rows() && lhs.cols() ==
> rhs.cols()' failed.
>
> It seems do happen when doing the component-wise division by the diagonal
> matrix. It works fine with a vector only though.
>
> Best regards
>
> Tim Hunter
>
> --
>
> Timothy Hunter
>
> Student (Stanford University)
>
> T. 404 421 3075