saberbud/ceres-solver
1
0
Fork 0
mirror of https://github.com/ceres-solver/ceres-solver.git synced 2026-08-30 00:50:37 +08:00
Code Issues Packages Projects Releases Wiki Activity

Use more performant, less conservative Eigen solvers.

Browse Source 
colPivHouseholderQR -> householderQR
ldlt -> llt.

The resulting performance differences are significant enough
to justify switching.

LAPACK's dgels routine used for solving linear least squares
problems does not use pivoting either.

Similarly, we are not actually using the fact that the matrix
being factorized can be indefinite when using LDLT factorization, so
its not clear that the performance hit is worth it.

These two changes result in Eigen being able to use blocking
algorithms, which for Cholesky factorization, brings the performance
closer to hardware optimized LAPACK. Similarly for dense QR
factorization, on intel there is a 2x speedup.

Change-Id: I4459ee0fc8eb87d58e2b299dfaa9e656d539dc5e
This commit is contained in:
Sameer Agarwal
2013-08-12 16:28:37 -07:00
committed by Keir Mierle
parent fb465a03b8
commit 080d1d04bd
9 changed files with 12 additions and 12 deletions
Download Patch File Download Diff File Expand all files Collapse all files
internal/ceres/implicit_schur_complement_test.cc
+2 -2
View File
@@ -109,7 +109,7 @@ class ImplicitSchurComplementTest : public ::testing::Test {
solution->setZero();
VectorRef schur_solution(solution->data() + num_cols_ - num_schur_rows,
num_schur_rows);
schur_solution = lhs->selfadjointView<Eigen::Upper>().ldlt().solve(*rhs);
schur_solution = lhs->selfadjointView<Eigen::Upper>().llt().solve(*rhs);
eliminator->BackSubstitute(A_.get(), b_.get(), D,
schur_solution.data(), solution->data());
}
@@ -156,7 +156,7 @@ class ImplicitSchurComplementTest : public ::testing::Test {
// Reference solution to the f_block.
const Vector reference_f_sol =
lhs.selfadjointView<Eigen::Upper>().ldlt().solve(rhs);
lhs.selfadjointView<Eigen::Upper>().llt().solve(rhs);
// Backsubstituted solution from the implicit schur solver using the
// reference solution to the f_block.