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