|[eigen] part<SelfAdjoint> in Eigen3|
[ Thread Index | Date Index | More lists.tuxfamily.org/eigen Archives ]
Eigen2 documentation mentions the way to obtain an optimized computation for an _expression_ evaluating to a selfadjoint matrix:
n.part<SelfAdjoint>() = m+m.adjoint ; (1)
n.part<SelfAdjoint>() = (m*m.adjoint()).lazy(); (2)
I didn’t find the way to have that in Eigen3 (MatrixBase>Derived>::part is still there but flagged as deprecated). The Porting from Eigen2 to Eigen3 page does not mention this case (it gives only the translation for part<SelfAdjoint|Upper> and part<SelfAdjoint|Lower>). Is there a new direct way to perform such computations?
And what would be the best way to write (2) when m itself is selfadjoint?
|Mail converted by MHonArc 2.6.19+||http://listengine.tuxfamily.org/|