This bugzilla service is closed. All entries have been migrated to https://gitlab.com/libeigen/eigen
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: 2019-12-04 15:19 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)
Comment 4 Nobody 2019-12-04 15:19:55 UTC
-- 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.

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