diff --git a/Eigen/src/SVD/BDCSVD.h b/Eigen/src/SVD/BDCSVD.h
index 2c1cbf922..17bd28869 100644
--- a/Eigen/src/SVD/BDCSVD.h
+++ b/Eigen/src/SVD/BDCSVD.h
@@ -1099,7 +1099,7 @@ void BDCSVD<MatrixType, Options>::perturbCol0(const ArrayRef& col0, const ArrayR
 #endif
           // Avoid index out of bounds.
           // Will end up setting zhat(k) = 0.
-          if (l == 0) {
+          if (i >= k && l == 0) {
             m_info = NumericalIssue;
             prod = 0;
             break;