| 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, |