| 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."); |