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