| Re: [eigen] 2.0.16 and 3.0.1 for the end of the week. |
[ Thread Index |
Date Index
| More lists.tuxfamily.org/eigen Archives
]
- To: eigen@xxxxxxxxxxxxxxxxxxx
- Subject: Re: [eigen] 2.0.16 and 3.0.1 for the end of the week.
- From: Gael Guennebaud <gael.guennebaud@xxxxxxxxx>
- Date: Fri, 27 May 2011 09:26:14 +0200
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:mime-version:in-reply-to:references:from:date :message-id:subject:to:content-type:content-transfer-encoding; bh=ypsCoOqWvTwc0jpvVyU7fKZlGiFyCY6IY1e/8EPQgl0=; b=t6h7Oh0kR75xMw5kUv0h8L2uxY6bx1efK7DCG9r4WTBS3hVYIFOOeQ9ZVeUcLNDdid tseoGLwqum8HyoRuViNMQqCIAcmovW1CWOE27/qk9+/3lMAUR8CBGpRsjLaC1BTTdEhu xmbnWTROO9o0xrBTPLPR/WLDKy21TU/XL9D1I=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:in-reply-to:references:from:date:message-id:subject:to :content-type:content-transfer-encoding; b=eo83KKgC9Af6QcaXZszZzcCFUviqCsaDhQnaiozoz54PY8YRZS4Jb8RPb1fJGIMlhb VAIttt1Qo4I+o43dS/ZyLtSetROIp2MdNmRH1AtKn6wZV+02T3tyEJEjwDiM+zaajdyR vPl/bBnHLCABdeuhYokhwin44qpH2WZZN+2gk=
confirmed, and this is produced by the changeset a59984204f9d "fix
matrix product with EIGEN_DEFAULT_TO_ROW_MAJOR". The 3.0 branch does
not exhibit this regression.
gael
On Thu, May 26, 2011 at 9:28 PM, Robert Lupton the Good
<rhl@xxxxxxxxxxxxxxxxxxx> wrote:
> Compiled -O3 n a 64-bit Linux box the following code takes 9ms with 2.0.10 and 285ms with 2.0.15
>
> I don't know the hg root to clone the release candidate for 2.0.16 so I haven't checked.
>
> R
>
>
> #include <Eigen/Core>
> #include <Eigen/Array>
>
> int main(){
> int n = 50;
> Eigen::MatrixXd m1(n,n);
> Eigen::MatrixXd m2(n,n);
> Eigen::MatrixXd m3(n,n);
> m1.setRandom();
> m2.setRandom();
> m3.setRandom();
>
> Eigen::VectorXd v1(n);
> Eigen::VectorXd v2(n);
> Eigen::VectorXd v3(n);
> v1.setRandom();
> v2.setRandom();
> v3.setRandom();
>
> for (int i = 0; i != 100; ++i) {
> v3 = m1*(m2*v1+m3.transpose()*m2*v2);
> }
>
> return 0;
> }
>
>
>
>