Changeset ef8667 for src/Parser/FormatParser_common.cpp

Timestamp:

Jul 3, 2017, 3:07:28 PM (8 years ago)

Author:

Frederik Heber <frederik.heber@…>

Branches:

ForceAnnealing_oldresults, IndependentFragmentGrids_IntegrationTest

Children:

67e885

Parents:

efc80e

git-author:

Frederik Heber <frederik.heber@…> (06/17/17 23:09:19)

git-committer:

Frederik Heber <frederik.heber@…> (07/03/17 15:07:28)

Message:

Extracted getMinMaxTrajectories() into FormatParser_common.

File:

• src/Parser/FormatParser_common.cpp (modified) (2 diffs)

src/Parser/FormatParser_common.cpp

@@ -36 +36 @@
 
 #include <iostream>
-
+#include <limits>
+
+#include "CodePatterns/Assert.hpp"
 #include "CodePatterns/Observer/Notification.hpp"
 #include "World.hpp"
@@ -179 +181 @@
   return iter->second;
 }
+
+std::pair<size_t, size_t> FormatParser_common::getMinMaxTrajectories(
+    const std::vector<const atom *> &_atoms) const
+{
+  // get max and min trajectories
+  size_t min_trajectories = std::numeric_limits<size_t>::max();
+  size_t max_trajectories = std::numeric_limits<size_t>::min();
+  for (std::vector<const atom *>::const_iterator iter = _atoms.begin();
+      iter != _atoms.end();
+      ++iter) {
+    if (max_trajectories < (*iter)->getTrajectorySize())
+      max_trajectories = (*iter)->getTrajectorySize();
+    if (min_trajectories > (*iter)->getTrajectorySize())
+      min_trajectories = (*iter)->getTrajectorySize();
+  }
+  // no atoms? Then, they all have same amount
+  if (_atoms.size() == 0)
+    min_trajectories = max_trajectories = 1;
+  ASSERT(min_trajectories == max_trajectories,
+      "FormatParser_common::getMinMaxTrajectories() - not all atoms have same number of trajectories: "
+      +toString(min_trajectories)+" != "+toString(max_trajectories)+".");
+
+  return std::make_pair(min_trajectories, max_trajectories);
+}