New user self-registration is currently disabled. Please email eigen-core-team @ lists.tuxfamily.org if you need an account.
Bug 1026 - factorize() von SparseLU solver never returns
factorize() von SparseLU solver never returns
Status: RESOLVED FIXED
Product: Eigen
Classification: Unclassified
Component: LU
3.2
x86 - 64-bit Linux
: Normal Wrong Result
Assigned To: Nobody
:
Depends on:
Blocks:
  Show dependency treegraph
 
Reported: 2015-06-23 08:20 UTC by Fabian
Modified: 2015-06-26 12:04 UTC (History)
2 users (show)



Attachments
compile and run (890 bytes, text/plain)
2015-06-23 08:20 UTC, Fabian
no flags Details

Description Fabian 2015-06-23 08:20:03 UTC
Created attachment 585 [details]
compile and run

We run a large simulation program, which makes use of the sparse capabilities of Eigen. One unexpected simulation scenario lead to the case that we  tried to invert a Matrix part, which did not include any components different to zero. It took us a while to figure out, that out program got stuck at the attempt to factorize a zero matrix. Of course there is no sense in factorizing a Zero-Matrix.

If there is no factorization possible, is there any hint given from eigen to detect this? In the actual implementation, the factorize method just never returns, although it should be easy to detect the zero matrix.

The attachment gives a minimal example of this behavior.

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