| 25 |
|
if (info->the_integrator != NULL){ |
| 26 |
|
delete info->the_integrator; |
| 27 |
|
} |
| 28 |
< |
info->the_integrator = this; |
| 29 |
< |
|
| 28 |
> |
|
| 29 |
|
nAtoms = info->n_atoms; |
| 30 |
|
|
| 31 |
|
// check for constraints |
| 173 |
|
dt = info->dt; |
| 174 |
|
dt2 = 0.5 * dt; |
| 175 |
|
|
| 176 |
+ |
readyCheck(); |
| 177 |
+ |
|
| 178 |
|
// initialize the forces before the first step |
| 179 |
|
|
| 180 |
|
calcForce(1, 1); |
| 181 |
+ |
|
| 182 |
+ |
if (nConstrained){ |
| 183 |
+ |
preMove(); |
| 184 |
+ |
constrainA(); |
| 185 |
+ |
calcForce(1, 1); |
| 186 |
+ |
constrainB(); |
| 187 |
+ |
} |
| 188 |
|
|
| 189 |
|
if (info->setTemp){ |
| 190 |
|
thermalize(); |
| 200 |
|
dumpOut->writeDump(info->getTime()); |
| 201 |
|
statOut->writeStat(info->getTime()); |
| 202 |
|
|
| 203 |
< |
readyCheck(); |
| 203 |
> |
|
| 204 |
|
|
| 205 |
|
#ifdef IS_MPI |
| 206 |
|
strcpy(checkPointMsg, "The integrator is ready to go."); |