ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/group/trunk/OOPSE/libmdtools/Integrator.cpp
(Generate patch)

Comparing trunk/OOPSE/libmdtools/Integrator.cpp (file contents):
Revision 726 by tim, Tue Aug 26 20:37:30 2003 UTC vs.
Revision 763 by tim, Mon Sep 15 16:52:02 2003 UTC

# Line 153 | Line 153 | template<typename T> void Integrator<T>::integrate(voi
153    double sampleTime = info->sampleTime;
154    double statusTime = info->statusTime;
155    double thermalTime = info->thermalTime;
156 +  double resetTime = info->resetTime;
157  
158 +
159    double currSample;
160    double currThermal;
161    double currStatus;
162 <
162 >  double currReset;
163 >  
164    int calcPot, calcStress;
165    int isError;
166  
# Line 174 | Line 177 | template<typename T> void Integrator<T>::integrate(voi
177    // initialize the forces before the first step
178  
179    calcForce(1, 1);
180 <  // myFF->doForces(1,1);
178 <
180 >  
181    if (info->setTemp){
182      thermalize();
183    }
184  
183 –  calcPot = 0;
184 –  calcStress = 0;
185 –  currSample = sampleTime;
186 –  currThermal = thermalTime;
187 –  currStatus = statusTime;
188 –  
185    calcPot     = 0;
186    calcStress  = 0;
187    currSample  = sampleTime + info->getTime();
188    currThermal = thermalTime+ info->getTime();
189    currStatus  = statusTime + info->getTime();
190 +  currReset   = resetTime  + info->getTime();
191  
192    dumpOut->writeDump(info->getTime());
193    statOut->writeStat(info->getTime());
# Line 231 | Line 228 | template<typename T> void Integrator<T>::integrate(voi
228        currStatus += statusTime;
229      }
230  
231 +    if (info->resetIntegrator){
232 +      if (info->getTime() >= currReset){
233 +        this->resetIntegrator();
234 +        currReset += resetTime;
235 +      }
236 +    }
237 +
238   #ifdef IS_MPI
239      strcpy(checkPointMsg, "successfully took a time step.");
240      MPIcheckPoint();
# Line 353 | Line 357 | template<typename T> void Integrator<T>::moveA(void){
357        angle = dt2 * ji[0] / I[0][0];
358        this->rotate(1, 2, angle, ji, A);
359  
356 –
360        dAtom->setJ(ji);
361        dAtom->setA(A);
362      }
# Line 750 | Line 753 | template<typename T> void Integrator<T>::thermalize(){
753   template<typename T> void Integrator<T>::thermalize(){
754    tStats->velocitize();
755   }
756 +
757 + template<typename T> double Integrator<T>::getConservedQuantity(void){
758 +  return tStats->getTotalE();
759 + }

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines