| 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 |
|
} |