254 |
|
|
255 |
|
moveA(); |
256 |
|
|
257 |
< |
if (nConstrained){ |
258 |
< |
constrainA(); |
259 |
< |
} |
257 |
> |
|
258 |
|
|
259 |
|
|
260 |
|
#ifdef IS_MPI |
277 |
|
|
278 |
|
moveB(); |
279 |
|
|
282 |
– |
if (nConstrained){ |
283 |
– |
constrainB(); |
284 |
– |
} |
280 |
|
|
281 |
+ |
|
282 |
|
#ifdef IS_MPI |
283 |
|
strcpy(checkPointMsg, "Succesful moveB\n"); |
284 |
|
MPIcheckPoint(); |
357 |
|
dAtom->setA(A); |
358 |
|
} |
359 |
|
} |
360 |
+ |
|
361 |
+ |
if (nConstrained){ |
362 |
+ |
constrainA(); |
363 |
+ |
} |
364 |
|
} |
365 |
|
|
366 |
|
|
402 |
|
dAtom->setJ(ji); |
403 |
|
} |
404 |
|
} |
405 |
+ |
|
406 |
+ |
if (nConstrained){ |
407 |
+ |
constrainB(); |
408 |
+ |
} |
409 |
|
} |
410 |
|
|
411 |
|
template<typename T> void Integrator<T>::preMove(void){ |
564 |
|
painCave.isFatal = 1; |
565 |
|
simError(); |
566 |
|
} |
567 |
+ |
|
568 |
|
} |
569 |
|
|
570 |
|
template<typename T> void Integrator<T>::constrainB(void){ |
761 |
|
|
762 |
|
template<typename T> double Integrator<T>::getConservedQuantity(void){ |
763 |
|
return tStats->getTotalE(); |
764 |
< |
} |
764 |
> |
} |