160 |
|
double thermalTime = info->thermalTime; |
161 |
|
double resetTime = info->resetTime; |
162 |
|
|
163 |
< |
|
163 |
> |
double difference; |
164 |
|
double currSample; |
165 |
|
double currThermal; |
166 |
|
double currStatus; |
179 |
|
|
180 |
|
readyCheck(); |
181 |
|
|
182 |
+ |
// remove center of mass drift velocity (in case we passed in a configuration |
183 |
+ |
// that was drifting |
184 |
+ |
tStats->removeCOMdrift(); |
185 |
+ |
|
186 |
|
// initialize the forces before the first step |
187 |
|
|
188 |
|
calcForce(1, 1); |
215 |
|
#endif // is_mpi |
216 |
|
|
217 |
|
while (info->getTime() < runTime && !stopIntegrator()){ |
218 |
< |
if ((info->getTime() + dt) >= currStatus){ |
218 |
> |
difference = info->getTime() + dt - currStatus; |
219 |
> |
if (difference > 0 || fabs(difference) < 1e-4 ){ |
220 |
|
calcPot = 1; |
221 |
|
calcStress = 1; |
222 |
|
} |