Summary: | crash using pastix support with pastix 5.2.2.22 | ||||||
---|---|---|---|---|---|---|---|
Product: | Eigen | Reporter: | fabien chêne <fabien.chene> | ||||
Component: | Sparse | Assignee: | Nobody <eigen.nobody> | ||||
Status: | RESOLVED FIXED | ||||||
Severity: | Crash | CC: | chtz, gael.guennebaud | ||||
Priority: | Normal | ||||||
Version: | 3.3 (current stable) | ||||||
Hardware: | x86 - 64-bit | ||||||
OS: | Linux | ||||||
Whiteboard: | |||||||
Attachments: |
|
Could you try adding out.resize(matrix.rows(), matrix.cols()); 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 > src/PaStiXSupport/PaStiXSupport.h? 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), https://bitbucket.org/eigen/eigen/commits/b8fb5886e55f/ 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. |
Created attachment 645 [details] testcase 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.