160 |
|
double thermalTime = info->thermalTime; |
161 |
|
double resetTime = info->resetTime; |
162 |
|
|
163 |
< |
|
163 |
> |
double difference; |
164 |
|
double currSample; |
165 |
|
double currThermal; |
166 |
|
double currStatus; |
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 |
|
} |
347 |
|
integrableObjects[i]->getVel(vel); |
348 |
|
integrableObjects[i]->getPos(pos); |
349 |
|
integrableObjects[i]->getFrc(frc); |
349 |
– |
|
350 |
– |
std::cerr << "i =\t" << i << "\t" << frc[0] << "\t" << frc[1]<< "\t" << frc[2] << "\n"; |
350 |
|
|
351 |
|
mass = integrableObjects[i]->getMass(); |
352 |
|
|