Fix bug in SparseLU dfs for dense matrices

This commit is contained in:
Desire NUENTSA 2013-06-11 14:48:04 +02:00
parent 9266f65318
commit 1bf18bd57f

View File

@ -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;