diff --git a/Eigen/src/IterativeLinearSolvers/IncompleteLUT.h b/Eigen/src/IterativeLinearSolvers/IncompleteLUT.h
index fd2a188e5..5ff97d6d3 100644
--- a/Eigen/src/IterativeLinearSolvers/IncompleteLUT.h
+++ b/Eigen/src/IterativeLinearSolvers/IncompleteLUT.h
@@ -381,7 +381,7 @@ void IncompleteLUT<Scalar>::factorize(const _MatrixType& amat)
           {
             newpos = sizel;
             sizel++;
-            eigen_internal_assert(sizel<ii);
+            eigen_internal_assert(sizel<=ii);
           }
           ju(newpos) = j;
           u(newpos) = -prod;