46 |
|
Integrator::Integrator(SimInfo* info) |
47 |
|
: info_(info), forceMan_(NULL) , needPotential(false), needStress(false), |
48 |
|
needReset(false), velocitizer_(NULL), needVelocityScaling(false), |
49 |
+ |
rnemd_(NULL), useRNEMD(false), |
50 |
|
dumpWriter(NULL), statWriter(NULL), thermo(info), |
51 |
|
currentSnapshot_(info->getSnapshotManager()->getCurrentSnapshot()) { |
52 |
|
|
96 |
|
needReset = true; |
97 |
|
resetTime = simParams->getResetTime(); |
98 |
|
} |
99 |
+ |
|
100 |
+ |
|
101 |
+ |
if (simParams->haveUseRNEMD()) { |
102 |
+ |
if (simParams->getUseRNEMD()) { |
103 |
+ |
useRNEMD = simParams->getUseRNEMD(); |
104 |
+ |
if (simParams->haveRNEMD_swapTime()) { |
105 |
+ |
RNEMD_swapTime = simParams->getRNEMD_swapTime(); |
106 |
+ |
} |
107 |
+ |
} |
108 |
+ |
} |
109 |
|
|
110 |
|
// Create a default ForceManager: If the subclass wants to use |
111 |
|
// a different ForceManager, use setForceManager |
130 |
|
// Create a default a velocitizer: If the subclass wants to use |
131 |
|
// a different velocitizer, use setVelocitizer |
132 |
|
velocitizer_ = new Velocitizer(info); |
133 |
+ |
// Create a default a RNEMD. |
134 |
+ |
rnemd_ = new RNEMD(info); |
135 |
|
|
136 |
|
} |
137 |
|
|
138 |
|
Integrator::~Integrator(){ |
139 |
|
delete forceMan_; |
140 |
|
delete velocitizer_; |
141 |
+ |
delete rnemd_; |
142 |
|
|
143 |
|
delete dumpWriter; |
144 |
|
delete statWriter; |