| 146 |
|
|
| 147 |
|
|
| 148 |
|
template<typename T> void Integrator<T>::integrate(void){ |
| 149 |
– |
int i, j; // loop counters |
| 149 |
|
|
| 150 |
|
double runTime = info->run_time; |
| 151 |
|
double sampleTime = info->sampleTime; |
| 160 |
|
double currReset; |
| 161 |
|
|
| 162 |
|
int calcPot, calcStress; |
| 164 |
– |
int isError; |
| 163 |
|
|
| 164 |
|
tStats = new Thermo(info); |
| 165 |
|
statOut = new StatWriter(info); |
| 166 |
|
dumpOut = new DumpWriter(info); |
| 167 |
|
|
| 168 |
|
atoms = info->atoms; |
| 171 |
– |
DirectionalAtom* dAtom; |
| 169 |
|
|
| 170 |
|
dt = info->dt; |
| 171 |
|
dt2 = 0.5 * dt; |
| 402 |
|
} |
| 403 |
|
|
| 404 |
|
template<typename T> void Integrator<T>::constrainA(){ |
| 405 |
< |
int i, j, k; |
| 405 |
> |
int i, j; |
| 406 |
|
int done; |
| 407 |
|
double posA[3], posB[3]; |
| 408 |
|
double velA[3], velB[3]; |
| 546 |
|
} |
| 547 |
|
|
| 548 |
|
template<typename T> void Integrator<T>::constrainB(void){ |
| 549 |
< |
int i, j, k; |
| 549 |
> |
int i, j; |
| 550 |
|
int done; |
| 551 |
|
double posA[3], posB[3]; |
| 552 |
|
double velA[3], velB[3]; |
| 555 |
|
int a, b, ax, ay, az, bx, by, bz; |
| 556 |
|
double rma, rmb; |
| 557 |
|
double dx, dy, dz; |
| 558 |
< |
double rabsq, pabsq, rvab; |
| 562 |
< |
double diffsq; |
| 558 |
> |
double rvab; |
| 559 |
|
double gab; |
| 560 |
|
int iteration; |
| 561 |
|
|