| 1340 |
|
void SimSetup::makeIntegrator(void){ |
| 1341 |
|
int k; |
| 1342 |
|
|
| 1343 |
+ |
NVE<RealIntegrator>* myNVE = NULL; |
| 1344 |
|
NVT<RealIntegrator>* myNVT = NULL; |
| 1345 |
|
NPTi<NPT<RealIntegrator> >* myNPTi = NULL; |
| 1346 |
|
NPTf<NPT<RealIntegrator> >* myNPTf = NULL; |
| 1350 |
|
case NVE_ENS: |
| 1351 |
|
if (globals->haveZconstraints()){ |
| 1352 |
|
setupZConstraint(info[k]); |
| 1353 |
< |
new ZConstraint<NVE<RealIntegrator> >(&(info[k]), the_ff); |
| 1353 |
> |
myNVE = new ZConstraint<NVE<RealIntegrator> >(&(info[k]), the_ff); |
| 1354 |
|
} |
| 1355 |
< |
else |
| 1356 |
< |
new NVE<RealIntegrator>(&(info[k]), the_ff); |
| 1355 |
> |
else{ |
| 1356 |
> |
myNVE = new NVE<RealIntegrator>(&(info[k]), the_ff); |
| 1357 |
> |
} |
| 1358 |
> |
|
| 1359 |
> |
info->the_integrator = myNVE; |
| 1360 |
|
break; |
| 1361 |
|
|
| 1362 |
|
case NVT_ENS: |
| 1378 |
|
painCave.isFatal = 1; |
| 1379 |
|
simError(); |
| 1380 |
|
} |
| 1381 |
+ |
|
| 1382 |
+ |
info->the_integrator = myNVT; |
| 1383 |
|
break; |
| 1384 |
|
|
| 1385 |
|
case NPTi_ENS: |
| 1421 |
|
painCave.isFatal = 1; |
| 1422 |
|
simError(); |
| 1423 |
|
} |
| 1424 |
+ |
|
| 1425 |
+ |
info->the_integrator = myNPTi; |
| 1426 |
|
break; |
| 1427 |
|
|
| 1428 |
|
case NPTf_ENS: |
| 1464 |
|
painCave.isFatal = 1; |
| 1465 |
|
simError(); |
| 1466 |
|
} |
| 1467 |
+ |
|
| 1468 |
+ |
info->the_integrator = myNPTf; |
| 1469 |
|
break; |
| 1470 |
|
|
| 1471 |
|
default: |