| 5 |
|
#include "Thermo.hpp" |
| 6 |
|
#include "ExtendedSystem.hpp" |
| 7 |
|
|
| 8 |
< |
ExtendedSystem::ExtendedSystem( SimInfo &info ) { |
| 8 |
> |
ExtendedSystem::ExtendedSystem( SimInfo* the_entry_plug ) { |
| 9 |
|
|
| 10 |
|
// get what information we need from the SimInfo object |
| 11 |
|
|
| 12 |
< |
entry_plug = &info; |
| 12 |
> |
entry_plug = the_entry_plug; |
| 13 |
|
nAtoms = entry_plug->n_atoms; |
| 14 |
|
atoms = entry_plug->atoms; |
| 15 |
|
nMols = entry_plug->n_mol; |
| 33 |
|
// amu*Ang^2*fs^-2/K |
| 34 |
|
DirectionalAtom* dAtom; |
| 35 |
|
|
| 36 |
– |
|
| 36 |
|
ke_temp = ke * e_convert; |
| 37 |
|
NkBT = (double)ndf * kB * targetTemp; |
| 38 |
|
|
| 40 |
|
// qmass is set in the parameter file |
| 41 |
|
|
| 42 |
|
zeta += dt * ( (ke_temp*2.0 - NkBT) / qmass ); |
| 43 |
+ |
std::cerr << "ke_temp = " << ke_temp << "\n"; |
| 44 |
+ |
|
| 45 |
|
zetaScale = zeta * dt; |
| 46 |
+ |
|
| 47 |
|
|
| 48 |
+ |
|
| 49 |
|
// perform thermostat scaling on linear velocities and angular momentum |
| 50 |
|
for(i = 0; i < nAtoms; i++){ |
| 51 |
|
|