diff --git a/src/algebra/CRSMatrixDescriptor.hpp b/src/algebra/CRSMatrixDescriptor.hpp
index 8f123896d486b7d6b67be1350a96183423b4235e..e367db354e356035d18a8742b2410a8781dde09a 100644
--- a/src/algebra/CRSMatrixDescriptor.hpp
+++ b/src/algebra/CRSMatrixDescriptor.hpp
@@ -35,7 +35,7 @@ class CRSMatrixDescriptor
     Array<IndexType> row_map(m_nb_rows + 1);
     row_map[0] = 0;
     for (IndexType i = 0; i < m_nb_rows; ++i) {
-      Assert(non_zeros[i] >= 0);
+      Assert((non_zeros[i] >= 0) and (non_zeros[i] <= m_nb_columns));
       row_map[i + 1] = row_map[i] + non_zeros[i];
     }