243 |
|
for( tl=0; tl < n_loops; tl++ ){ |
244 |
|
|
245 |
|
if (!strcasecmp( entry_plug->ensemble, "NVT")) |
246 |
< |
myES->NoseHooverNVT( dt , tStats->getKinetic() ); |
246 |
> |
myES->NoseHooverNVT( dt / 2.0 , tStats->getKinetic() ); |
247 |
|
|
248 |
|
for( j=0; j<nAtoms; j++ ){ |
249 |
|
|
411 |
|
|
412 |
|
|
413 |
|
if (!strcasecmp( entry_plug->ensemble, "NVT")) |
414 |
< |
myES->NoseHooverNVT( dt , tStats->getKinetic() ); |
414 |
> |
myES->NoseHooverNVT( dt / 2.0, tStats->getKinetic() ); |
415 |
|
|
416 |
|
if (!strcasecmp( entry_plug->ensemble, "NPT") ) |
417 |
|
myES->NoseHooverAndersonNPT( dt, |
446 |
|
trans_kE = 0.0; |
447 |
|
|
448 |
|
if (!strcasecmp( entry_plug->ensemble, "NVT")) |
449 |
< |
myES->NoseHooverNVT( dt , tStats->getKinetic() ); |
449 |
> |
myES->NoseHooverNVT( dt / 2.0, tStats->getKinetic() ); |
450 |
|
|
451 |
|
for( i=0; i<nAtoms; i++ ){ |
452 |
|
|
589 |
|
} |
590 |
|
|
591 |
|
if (!strcasecmp( entry_plug->ensemble, "NVT")) |
592 |
< |
myES->NoseHooverNVT( dt , tStats->getKinetic() ); |
592 |
> |
myES->NoseHooverNVT( dt / 2.0, tStats->getKinetic() ); |
593 |
|
|
594 |
|
if (!strcasecmp( entry_plug->ensemble, "NPT") ) |
595 |
|
myES->NoseHooverAndersonNPT( dt, |