Changeset ef8667 for src/Parser/PdbParser.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/PdbParser.cpp (modified) (2 diffs)

src/Parser/PdbParser.cpp

@@ -198 +198 @@
 
   // check for maximum number of time steps
-  size_t max_timesteps = 0;
-  for (vector<const atom *>::const_iterator atomIt = AtomList.begin();
-      atomIt != AtomList.end(); atomIt++) {
-    const atom * _atom = *atomIt;
-    LOG(4, "INFO: Atom " << _atom->getName() << " "
-        << *dynamic_cast<const AtomInfo *>(_atom) << ".");
-    if (_atom->getTrajectorySize() > max_timesteps)
-      max_timesteps = _atom->getTrajectorySize();
-  }
-  LOG(2,"INFO: Found a maximum of " << max_timesteps << " time steps to store.");
+  std::pair<size_t, size_t> minmax_trajectories =
+      getMinMaxTrajectories(AtomList);
+  LOG(2, "INFO: There are " << minmax_trajectories.second << " steps to save.");
 
   // re-distribute serials
@@ -222 +215 @@
 
   // store all time steps (always do first step)
-  for (size_t step = 0; (step == 0) || (step < max_timesteps); ++step) {
+  for (size_t step = 0; (step == 0) || (step < minmax_trajectories.second); ++step) {
     {
       // add initial remark