From e741b7becae92e87d9851e41858b892cfeb02755 Mon Sep 17 00:00:00 2001 From: Benoit Jacob Date: Mon, 4 Aug 2008 23:47:09 +0000 Subject: [PATCH] further big perf improvement in Inverse --- Eigen/src/LU/Inverse.h | 10 ++++------ 1 file changed, 4 insertions(+), 6 deletions(-) diff --git a/Eigen/src/LU/Inverse.h b/Eigen/src/LU/Inverse.h index 1a1fc2334..dd34d65e2 100644 --- a/Eigen/src/LU/Inverse.h +++ b/Eigen/src/LU/Inverse.h @@ -65,9 +65,8 @@ struct ei_compute_inverse_in_general_case inverse.row(size-1) /= matrix(size-1,size-1); for(int k = size-1; k >= 1; k--) - { - inverse.block(0,0,k,size) -= matrix.col(k).start(k) * inverse.row(k); - } + for(int row = 0; row < k; row++) + inverse.row(row) -= inverse.row(k) * matrix.coeff(row,k); } }; @@ -103,9 +102,8 @@ struct ei_compute_inverse_in_general_case inverse.col(size-1) /= matrix(size-1,size-1); for(int k = size-1; k >= 1; k--) - { - inverse.block(0,0,size,k) -= inverse.col(k) * matrix.row(k).start(k); - } + for(int col = 0; col < k; col++) + inverse.col(col) -= inverse.col(k) * matrix.coeff(k,col); } };