Re: [eigen] symmetrize the eigne way |
[ Thread Index |
Date Index
| More lists.tuxfamily.org/eigen Archives
]
- To: eigen@xxxxxxxxxxxxxxxxxxx
- Subject: Re: [eigen] symmetrize the eigne way
- From: Gael Guennebaud <gael.guennebaud@xxxxxxxxx>
- Date: Sat, 30 Jan 2010 09:03:21 +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 :from:date:message-id:subject:to:content-type :content-transfer-encoding; bh=BzfPo+/4sOpR6w7dV/5QCxTelii3nIbOODCEiSBvyQQ=; b=n5gsUy15u6aEu5zCmLK8rmfhEgIq8FMkKspQ/lWS2R4YS+e+JuzAfElnSUIhzYs2iE 8Odq+H6OmnI3zGZd0FNPljB5Jjcvc1grpKLWc35tas200W13n7rBZPOU06xGL/O+gfjI DediNzQuqxhoZ+CDcDbAsD0OqJrt2rU/zbSog=
- 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=U8L/E/xsxqX3FxhHJrNiA6HYoD/sc+VuYHiK4Q+JLsjEAF7Kg3OnortOQWW+2upvbd ogfjJB6CMFMdVloLQdquJuscxEN9g+UUDqKNpwHKQpJPc6QW5fF43wXoq6qxDYkllUXc azfOomzU0zmTzTcdtKDNIMzTYIoJB867m1AoE=
that was a bug. both versions were supposed to work. problem solved.
gael.
On Mon, Jan 25, 2010 at 4:06 PM, Benoit Jacob <jacob.benoit.1@xxxxxxxxx> wrote:
> 2010/1/25 Thomas Capricelli <orzel@xxxxxxxxxxxxxxx>:
>>
>> Suppose i want to create a selfadjoint matrix this way:
>> for (j = 0; j < n; ++j)
>> for (i = j+1; i < n; ++i)
>> s(i,j) = s(j,i);
>>
>> How could i do that using Eigen? I've thought about :
>> s.template triangularView<StrictlyLower>() = s.template triangularView<StrictlyUpper>().transpose();
>>
>> But this wont compile:
>> Eigen/src/Core/TriangularMatrix.h:425: error: no matching function for call to ‘Eigen::Matrix<double, 33331, 33331, 0, 33331, 33331>::copyCoeff(int&, int&, const Eigen::TriangularView<Eigen::Transpose<Eigen::Matrix<double, 33331, 33331, 0, 33331, 33331> >, 9u>&)’
>
> OK how about:
>
> s.template triangularView<StrictlyLower>() = s.transpose();
>
> I think the more general problem with what you tried is that we don't
> support syntax using triangular views on both sides. IIRC?
>
> Benoit>
>>
>> ++
>> --
>> Thomas Capricelli <orzel@xxxxxxxxxxxxxxx>
>> http://www.freehackers.org/thomas
>>
>>
>>
>>
>
>
>