Re: [eigen] SelfAdjointView<> as LHS expression for optimizations? |

[ Thread Index | Date Index | More lists.tuxfamily.org/eigen Archives ]

*To*: eigen@xxxxxxxxxxxxxxxxxxx*Subject*: Re: [eigen] SelfAdjointView<> as LHS expression for optimizations?*From*: Jose Luis Blanco <joseluisblancoc@xxxxxxxxx>*Date*: Thu, 11 Nov 2010 01:00:35 +0100*Dkim-signature*: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:mime-version:received:received:in-reply-to :references:date:message-id:subject:from:to:content-type; bh=kc59ca6q1xLuX1H8zhQidSIjTYNbACjusM4m+OMGnOo=; b=PJ0tM4B+HQ4xVBsPmV7f/aoH/eLGis0by6hNmw1H7t3do42Mjgr9Bcv9dpMmu3pBJE pfEdZntz+BXYZCfXcxkPo1bdLdRYcHTzzuv/pf3WS1pFp4N9IwiF6cCfe8t003cCgomR DOWPwK8aqhoNTe2eEgZFfhDHXWhYdf2fWB9fs=*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; b=DA7pRDm8EhZwTJW5CfaI13cz3Mo4ASuO/PcP3/IZOKzxQQ5W9PEh8hZLgQsQ1zzdOc ALIgKfzRLRlCicMsChMaqMO1PVXMDlWwHyiDnm9c8/JwcRfD7/rTMimDxNjbYDc+biNJ V0KCYam3cZINqa8C3D7PLl/iRCvNuwZGhaxm4=

Thanks a lot for the help and the work, I tested it and it's OK! Just to complete this thread, here's the code for the case one also wants the resulting matrix with both triangular parts filled instead of just one lo/up triangular part: template <typename MAT_H, typename MAT_C, typename MAT_R> inline void multiply_HCHt( const MAT_H &H, const MAT_C &C, MAT_R &R) { R = ( (H * C.template selfadjointView<Eigen::Lower>()).eval() * H.adjoint()).eval().template selfadjointView<Eigen::Lower>(); } Best, JL

**References**:**[eigen] SelfAdjointView<> as LHS expression for optimizations?***From:*Jose Luis Blanco

**Re: [eigen] SelfAdjointView<> as LHS expression for optimizations?***From:*Gael Guennebaud

**Re: [eigen] SelfAdjointView<> as LHS expression for optimizations?***From:*Benoit Jacob

**Re: [eigen] SelfAdjointView<> as LHS expression for optimizations?***From:*Jose Luis Blanco

**Re: [eigen] SelfAdjointView<> as LHS expression for optimizations?***From:*Gael Guennebaud

**Re: [eigen] SelfAdjointView<> as LHS expression for optimizations?***From:*Gael Guennebaud

**Messages sorted by:**[ date | thread ]- Prev by Date:
**Re: [eigen] SelfAdjointView<> as LHS expression for optimizations?** - Next by Date:
**[eigen] How to obtain the matrices Q and R of a QR decomposition** - Previous by thread:
**Re: [eigen] SelfAdjointView<> as LHS expression for optimizations?** - Next by thread:
**Re: [eigen] SelfAdjointView<> as LHS expression for optimizations?**

Mail converted by MHonArc 2.6.19+ | http://listengine.tuxfamily.org/ |