Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • src/Atom/atom.cpp

    r7e51e1 r8cc22f  
    6969    mol(0)
    7070{
    71   AtomicPosition = pointer->AtomicPosition; // copy trajectory of coordination
    72   AtomicVelocity = pointer->AtomicVelocity; // copy trajectory of velocity
    73   AtomicForce = pointer->AtomicForce;
    7471  // sign on to global atom change tracker
    7572  AtomObserver::getInstance().AtomInserted(this);
     
    9390
    9491
    95 void atom::UpdateSteps()
    96 {
    97   LOG(4,"atom::UpdateSteps() called.");
     92void atom::UpdateStep(const unsigned int _step)
     93{
     94  LOG(4,"atom::UpdateStep() called.");
    9895  // append to position, velocity and force vector
    99   AtomInfo::AppendTrajectoryStep();
     96  AtomInfo::AppendTrajectoryStep(WorldTime::getTime()+1);
    10097  // append to ListOfBonds vector
    101   BondedParticleInfo::AppendTrajectoryStep();
    102 }
    103 
    104 void atom::removeSteps()
    105 {
    106   LOG(4,"atom::removeSteps() called.");
     98  BondedParticleInfo::AppendTrajectoryStep(WorldTime::getTime()+1);
     99}
     100
     101void atom::removeStep(const unsigned int _step)
     102{
     103  LOG(4,"atom::removeStep() called.");
    107104  // append to position, velocity and force vector
    108   AtomInfo::removeTrajectoryStep();
     105  AtomInfo::removeTrajectoryStep(_step);
    109106  // append to ListOfBonds vector
    110   BondedParticleInfo::removeTrajectoryStep();
     107  BondedParticleInfo::removeTrajectoryStep(_step);
    111108}
    112109
     
    313310atom* NewAtom(atomId_t _id){
    314311  atom * res = new atom();
    315   // extent trajectory to current time step
    316   const size_t CurrentTime = WorldTime::getTime();
    317   for (size_t step = res->getTrajectorySize(); step <= CurrentTime; ++step)
    318     res->UpdateSteps();
    319312  res->setId(_id);
    320313  return res;
Note: See TracChangeset for help on using the changeset viewer.