mirror/eigen
2
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-04-06 19:10:36 +08:00
Code Issues Packages Projects Releases Wiki Activity

+5% optimization in 4x4 inverse:

Browse Source 
-only evaluate block expressions for which that is beneficial
-don't check for invertibility unless requested
This commit is contained in:
Benoit Jacob 2008-04-15 20:39:27 +00:00
parent 6747b45ae7
commit 43e2bc14fe
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
Eigen/src/LU/Inverse.h
View File

@ -196,13 +196,13 @@ void Inverse<ExpressionType, CheckExistence>::_compute_in_size4_case(const Expre
{
const Block22 Q = xpr.template block<2,2>(0,2);
const Block22 P_inverse_times_Q = P_inverse * Q;
const Block22 R = xpr.template block<2,2>(2,0);
const XprBlock22 R = xpr.template block<2,2>(2,0);
const Block22 R_times_P_inverse = R * P_inverse;
const Block22 R_times_P_inverse_times_Q = R_times_P_inverse * Q;
const Block22 S = xpr.template block<2,2>(2,2);
const XprBlock22 S = xpr.template block<2,2>(2,2);
const Block22 X = S - R_times_P_inverse_times_Q;
Block22 Y;
if(ei_compute_size2_inverse<Block22, Block22, true>(X, &Y))
if(ei_compute_size2_inverse<Block22, Block22, CheckExistence>(X, &Y))
{
m_inverse.template block<2,2>(2,2) = Y;
m_inverse.template block<2,2>(2,0) = - Y * R_times_P_inverse;