| 1281 | 
  | 
 | 
| 1282 | 
  | 
void SimSetup::makeIntegrator( void ){ | 
| 1283 | 
  | 
 | 
| 1284 | 
< | 
  NVT*  myNVT = NULL; | 
| 1285 | 
< | 
  NPTi* myNPTi = NULL; | 
| 1286 | 
< | 
  NPTf* myNPTf = NULL; | 
| 1287 | 
< | 
  NPTim* myNPTim = NULL; | 
| 1288 | 
< | 
  NPTfm* myNPTfm = NULL; | 
| 1284 | 
> | 
  NVT<RealIntegrator>*  myNVT = NULL; | 
| 1285 | 
> | 
  NPTi<RealIntegrator>* myNPTi = NULL; | 
| 1286 | 
> | 
  NPTf<RealIntegrator>* myNPTf = NULL; | 
| 1287 | 
> | 
  NPTim<RealIntegrator>* myNPTim = NULL; | 
| 1288 | 
> | 
  NPTfm<RealIntegrator>* myNPTfm = NULL; | 
| 1289 | 
  | 
 | 
| 1290 | 
  | 
  switch( ensembleCase ){ | 
| 1291 | 
  | 
 | 
| 1292 | 
  | 
  case NVE_ENS: | 
| 1293 | 
< | 
    new NVE( info, the_ff ); | 
| 1293 | 
> | 
    new NVE<RealIntegrator>( info, the_ff ); | 
| 1294 | 
  | 
    break; | 
| 1295 | 
  | 
 | 
| 1296 | 
  | 
  case NVT_ENS: | 
| 1297 | 
< | 
    myNVT = new NVT( info, the_ff ); | 
| 1297 | 
> | 
    myNVT = new NVT<RealIntegrator>( info, the_ff ); | 
| 1298 | 
  | 
    myNVT->setTargetTemp(globals->getTargetTemp()); | 
| 1299 | 
  | 
 | 
| 1300 | 
  | 
    if (globals->haveTauThermostat())  | 
| 1310 | 
  | 
    break; | 
| 1311 | 
  | 
 | 
| 1312 | 
  | 
  case NPTi_ENS: | 
| 1313 | 
< | 
    myNPTi = new NPTi( info, the_ff ); | 
| 1313 | 
> | 
    myNPTi = new NPTi<RealIntegrator>( info, the_ff ); | 
| 1314 | 
  | 
    myNPTi->setTargetTemp( globals->getTargetTemp() ); | 
| 1315 | 
  | 
 | 
| 1316 | 
  | 
    if (globals->haveTargetPressure()) | 
| 1345 | 
  | 
    break; | 
| 1346 | 
  | 
 | 
| 1347 | 
  | 
  case NPTf_ENS: | 
| 1348 | 
< | 
    myNPTf = new NPTf( info, the_ff ); | 
| 1348 | 
> | 
    myNPTf = new NPTf<RealIntegrator>( info, the_ff ); | 
| 1349 | 
  | 
    myNPTf->setTargetTemp( globals->getTargetTemp()); | 
| 1350 | 
  | 
 | 
| 1351 | 
  | 
    if (globals->haveTargetPressure()) | 
| 1380 | 
  | 
    break; | 
| 1381 | 
  | 
     | 
| 1382 | 
  | 
  case NPTim_ENS: | 
| 1383 | 
< | 
    myNPTim = new NPTim( info, the_ff ); | 
| 1383 | 
> | 
    myNPTim = new NPTim<RealIntegrator>( info, the_ff ); | 
| 1384 | 
  | 
    myNPTim->setTargetTemp( globals->getTargetTemp()); | 
| 1385 | 
  | 
 | 
| 1386 | 
  | 
    if (globals->haveTargetPressure()) | 
| 1415 | 
  | 
    break; | 
| 1416 | 
  | 
 | 
| 1417 | 
  | 
  case NPTfm_ENS: | 
| 1418 | 
< | 
    myNPTfm = new NPTfm( info, the_ff ); | 
| 1418 | 
> | 
    myNPTfm = new NPTfm<RealIntegrator>( info, the_ff ); | 
| 1419 | 
  | 
    myNPTfm->setTargetTemp( globals->getTargetTemp()); | 
| 1420 | 
  | 
 | 
| 1421 | 
  | 
    if (globals->haveTargetPressure()) |