48 |
|
needVelocityScaling(false), dumpWriter(NULL), statWriter(NULL), thermo(info), |
49 |
|
currentSnapshot_(info->getSnapshotManager()->getCurrentSnapshot()) { |
50 |
|
|
51 |
< |
Globals* simParams = info->getSimParams(); |
51 |
> |
simParams = info->getSimParams(); |
52 |
|
|
53 |
|
if (simParams->haveDt()) { |
54 |
|
dt = simParams->getDt(); |
90 |
|
currentSnapshot_->setTime(0.0); |
91 |
|
} |
92 |
|
|
93 |
< |
//create a default a ForceManager |
94 |
< |
//if the subclass want to using different ForceManager, use setForceManager |
93 |
> |
//create a default ForceManager |
94 |
> |
//if the subclass wants to use a different ForceManager, use setForceManager |
95 |
|
forceMan_ = new ForceManager(info); |
96 |
|
|
97 |
+ |
//set the force manager for thermodynamic integration if specified |
98 |
+ |
if (simParams->getUseSolidThermInt() || simParams->getUseLiquidThermInt()){ |
99 |
+ |
ThermoIntegrationForceManager* thermoForce_ |
100 |
+ |
= new ThermoIntegrationForceManager(info); |
101 |
+ |
setForceManager(thermoForce_); |
102 |
+ |
} |
103 |
+ |
|
104 |
|
// check for the temperature set flag (velocity scaling) |
105 |
|
if (simParams->haveTempSet()) { |
106 |
|
needVelocityScaling = simParams->getTempSet(); |