SparseLU factorization error: THE MATRIX IS STRUCTURALLY SINGULAR ... ZERO COLUMN

Submitted by cDc

Assigned to Nobody

Link to original bugzilla bug (#1053)
Version: 3.2

Description

Created attachment 597
test

SparseLU factorization step fails always (except for very small and particular patrices) with an error message similar to:

THE MATRIX IS STRUCTURALLY SINGULAR ... ZERO COLUMN AT 428

The problem appears only when compiling with RowMajor as define:

#define EIGEN_DEFAULT_TO_ROW_MAJOR

Please find attached a small code that reproduces the problem with a test file (which is not important, all problems give a similar error). Please run as:

SparseLU.exe coeff.bin

If the EIGEN_DEFAULT_TO_ROW_MAJOR define is commented the factorization completes as expected.

Attachment 597, "test":
SparseLU.zip