Created attachment 645 [details]
Using the attached testcase, a crash occurs:
/home/fabien/EigenDev/Eigen3.3/eigen-eigen-09a8e2186610/Eigen/src/SparseCore/SparseSelfAdjointView.h:156: void Eigen::SparseSelfAdjointView<MatrixType, UpLo>::resize(Eigen::Index, Eigen::Index) [with MatrixType = Eigen::SparseMatrix<double, 0>; unsigned int _Mode = 1u; Eigen::Index = long int]: Assertion `rows == this->rows() && cols == this->cols() && "SparseSelfadjointView::resize() does not actually allow to resize."' failed.
I am using the 3.3_beta version.
Could you try adding
at the start of both grabMatrix(...) methods in src/PaStiXSupport/PaStiXSupport.h?
I'm afraid, PaStiX is barely tested ...
(In reply to Christoph Hertzberg from comment #1)
> Could you try adding
> out.resize(matrix.rows(), matrix.cols());
> at the start of both grabMatrix(...) methods in
Thanks, it works for the three decompositions, not sure however this is needed for LU.
> I'm afraid, PaStiX is barely tested ...
I've got some pending bug reports on PaStiX support...
Thanks (I've a working pastix installation, all test passes fine now),
Summary: Bug 1145: fix PastixSupport LLT/LDLT wrappers (missing resize prior to calls to selfAdjointView)
-- GitLab Migration Automatic Message --
This bug has been migrated to gitlab.com's GitLab instance and has been closed from further activity.
You can subscribe and participate further through the new bug through this link to our GitLab instance: https://gitlab.com/libeigen/eigen/issues/1145.