- fixed a bug in matrix filling

diff --git a/src/domaintype.cpp b/src/domaintype.cpp
index 86c0f6f21fa3db97b3af76eecc6c83db6db56d19..9659d70a306b72b4492f9e158eb96393c5a44d32 100644 (file)
--- a/src/domaintype.cpp
+++ b/src/domaintype.cpp
@@ -53,7 +53,7 @@ void domainType::update(const std::vector< std::vector< RVec > > frame, const in
         #pragma omp parallel for
         for (unsigned long j = 0; j < graph.front().size(); j++) {
             for (unsigned int k1 = 0; k1 < graph[i][j].size(); k1++) {
-                for (unsigned int k2 = 0; k2 < graph[i][j].size(); k2++) {
+                for (unsigned int k2 = k1; k2 < graph[i][j].size(); k2++) {
                     if ((frame[j][k1] - frame[j][k2] - graph[i][j][k1][k2].radius).norm() <= static_cast< float >(eps)) {
                         if (k1 == k2) {
                             graph[i][j][k1][k2].num++;