crash when using IncompleteCholesky preconditioner
Submitted by fabien chêne
Assigned to Nobody
Link to original bugzilla bug (#1150)
Version: 3.3 (current stable)
Description
Created attachment 646
testcase
The attached testcase crashes when compiled without NDEBUG on this assertion:
IncompleteCholesky.h:155: void Eigen::IncompleteCholesky<Scalar, _UpLo, _OrderingType>::_solve_impl(const Rhs&, Dest&) const [with Rhs = Eigen::Matrix<double, -1, 1>; Dest = Eigen::Matrix<double, -1, 1>; Scalar = double; int _UpLo = 1; _OrderingType = Eigen::AMDOrdering<int>]: Assertion `m_factorizationIsOk && "factorize() should be called first"' failed.
With N < 10, the assertion does not trigger.
Attachment 646, "testcase":
TestIncompleteCholesky.cc
Edited by Eigen Bugzilla