diff options
| author |  Desire NUENTSA <desire.nuentsa_wakam@inria.fr> | 2013-06-11 14:48:04 +0200 |
|---|---|---|
| committer | Desire NUENTSA <desire.nuentsa_wakam@inria.fr> | 2013-06-11 14:48:04 +0200 |
| commit | 1bf18bd57feed5d23c752c97457c2a9884e291bb (patch) | |
| tree | 77d60d1725c463f997c2f17eb3c51a4c54a86f2c /Eigen/src/SparseLU/SparseLU_column_dfs.h | |
| parent | 9266f65318381bd34ddc719f21f8fa9b4e1521e4 (diff) | |
Fix bug in SparseLU dfs for dense matrices
Diffstat (limited to 'Eigen/src/SparseLU/SparseLU_column_dfs.h')
| -rw-r--r-- | Eigen/src/SparseLU/SparseLU_column_dfs.h | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/Eigen/src/SparseLU/SparseLU_column_dfs.h b/Eigen/src/SparseLU/SparseLU_column_dfs.h index bd450ddc7..bc4cfbf37 100644 --- a/Eigen/src/SparseLU/SparseLU_column_dfs.h +++ b/Eigen/src/SparseLU/SparseLU_column_dfs.h @@ -101,7 +101,7 @@ Index SparseLUImpl<Scalar,Index>::column_dfs(const Index m, const Index jcol, In column_dfs_traits<IndexVector, ScalarVector> traits(jcol, jsuper, glu, *this); // For each nonzero in A(*,jcol) do dfs - for (Index k = 0; lsub_col[k] != emptyIdxLU; k++) + for (Index k = 0; ((k < m) ? lsub_col[k] != emptyIdxLU : false) ; k++) { Index krow = lsub_col(k); lsub_col(k) = emptyIdxLU; |