full matrixes

diff --git a/src/spacetimecorr.cpp b/src/spacetimecorr.cpp
index 68fa667235a879c47418e4655b3ec483730118db..507001d5773137f6e5c0e52e8ab89f9838975e9f 100644 (file)
--- a/src/spacetimecorr.cpp
+++ b/src/spacetimecorr.cpp
@@ -506,7 +506,7 @@ Domains::finishAnalysis(int nframes)
             }
             for (int j = 0; j < frankenstein_trajectory.size() - i; j++) {
                 for (int f1 = 0; f1 < index.size(); f1++) {
-                    for (int f2 = f1; f2 < index.size(); f2++) {
+                    for (int f2 = 0; f2 < index.size(); f2++) {
                         rvec temp1, temp2;
                         rvec_sub(frankenstein_trajectory[basic_frame][f1], frankenstein_trajectory[j][f1], temp1);
                         rvec_dec(temp1, medx[f1]);
@@ -519,7 +519,7 @@ Domains::finishAnalysis(int nframes)
                 }
             }
             for (int j = 0; j < index.size(); j++) {
-                for (int f = j; f < index.size(); f++) {
+                for (int f = 0; f < index.size(); f++) {
                     crltns[i][j][f] = a[j][f] / (std::sqrt(b[j][f] * c[j][f]));
                 }
             }