New user self-registration is disabled due to spam. Please email eigen-core-team @ lists.tuxfamily.org if you need an account.
Bug 1145 - crash using pastix support with pastix 5.2.2.22
Summary: crash using pastix support with pastix 5.2.2.22
Status: RESOLVED FIXED
Alias: None
Product: Eigen
Classification: Unclassified
Component: Sparse (show other bugs)
Version: 3.3 (current stable)
Hardware: x86 - 64-bit Linux
: Normal Crash
Assignee: Nobody
URL:
Whiteboard:
Keywords:
Depends on:
Blocks:
 
Reported: 2016-01-13 22:00 UTC by fabien chêne
Modified: 2016-01-20 17:52 UTC (History)
2 users (show)



Attachments
testcase (450 bytes, text/x-c++src)
2016-01-13 22:00 UTC, fabien chêne
no flags Details

Description fabien chêne 2016-01-13 22:00:04 UTC
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.
Comment 1 Christoph Hertzberg 2016-01-20 10:08:35 UTC
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 ...
Comment 2 fabien chêne 2016-01-20 12:20:48 UTC
(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...
Comment 3 Gael Guennebaud 2016-01-20 17:52:53 UTC
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)

Note You need to log in before you can comment on or make changes to this bug.