| 174 |
|
|
| 175 |
|
// initialize the forces before the first step |
| 176 |
|
|
| 177 |
+ |
std::cerr << "Before initial Force calc\n"; |
| 178 |
+ |
|
| 179 |
|
calcForce(1, 1); |
| 180 |
|
|
| 181 |
|
if (nConstrained){ |
| 183 |
|
constrainA(); |
| 184 |
|
calcForce(1, 1); |
| 185 |
|
constrainB(); |
| 186 |
+ |
std::cerr << "premove done\n"; |
| 187 |
|
} |
| 188 |
|
|
| 189 |
+ |
|
| 190 |
+ |
|
| 191 |
|
if (info->setTemp){ |
| 192 |
|
thermalize(); |
| 193 |
|
} |
| 245 |
|
} |
| 246 |
|
} |
| 247 |
|
|
| 248 |
+ |
std::cerr << "done with time = " << info->getTime() << "\n"; |
| 249 |
+ |
|
| 250 |
|
#ifdef IS_MPI |
| 251 |
|
strcpy(checkPointMsg, "successfully took a time step."); |
| 252 |
|
MPIcheckPoint(); |