diff --git a/sPuReMD/src/lin_alg.c b/sPuReMD/src/lin_alg.c
index 110c43a3f747ac15e20bb0dc7958a9354963dbfd..93e833a73a7d2170c75bd4f379d60fc829e9c114 100644
--- a/sPuReMD/src/lin_alg.c
+++ b/sPuReMD/src/lin_alg.c
@@ -1699,7 +1699,7 @@ real sparse_approx_inverse( const sparse_matrix * const A,
                 }
                 // change the value if any of the column indices is seen
                 for ( d_j = A->start[pos_x[d_i]];
-                        d_j < A->start[pos_x[d_i + 1]]; ++d_j )
+                        d_j < A->start[pos_x[d_i] + 1]; ++d_j )
                 {
                     if ( Y[A->j[d_j]] == 1 )
                     {