Changeset ef8667 for src/Parser/XyzParser.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/XyzParser.cpp (modified) (1 diff)

src/Parser/XyzParser.cpp

@@ -169 +169 @@
 
   // 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< xyz >::save() - not all atoms have same number of trajectories: "
-      +toString(min_trajectories)+" != "+toString(max_trajectories)+".");
-  LOG(2, "INFO: There are " << max_trajectories << " steps to save.");
+  std::pair<size_t, size_t> minmax_trajectories =
+      getMinMaxTrajectories(atoms);
+  LOG(2, "INFO: There are " << minmax_trajectories.second << " steps to save.");
 
   // always store at least one step
-  for (size_t step = 0; (step < max_trajectories) || (step == 0); ++step) {
+  for (size_t step = 0; (step < minmax_trajectories.second) || (step == 0); ++step) {
     if (step != 0)
       *file << "\n";