diff --git a/src/fix_ave_correlate.cpp b/src/fix_ave_correlate.cpp
index 6002dc7f4a..f9c3cac235 100644
--- a/src/fix_ave_correlate.cpp
+++ b/src/fix_ave_correlate.cpp
@@ -541,7 +541,7 @@ void FixAveCorrelate::accumulate()
     for (k = 0; k < nsample; k++) {
       ipair = 0;
       for (i = 0; i < nvalues; i++)
-        for (j = 0; j < i-1; j++)
+        for (j = 0; j < i; j++)
           corr[k][ipair++] += values[m][i]*values[n][j];
       m--;
       if (m < 0) m = nrepeat-1;
@@ -561,7 +561,7 @@ void FixAveCorrelate::accumulate()
     for (k = 0; k < nsample; k++) {
       ipair = 0;
       for (i = 0; i < nvalues; i++)
-        for (j = 0; j < i; j++)
+        for (j = 0; j <= i; j++)
           corr[k][ipair++] += values[m][i]*values[n][j];
       m--;
       if (m < 0) m = nrepeat-1;