55 |
|
currSnapshot_ = info_->getSnapshotManager()->getCurrentSnapshot(); |
56 |
|
simParam = info_->getSimParams(); |
57 |
|
|
58 |
< |
if (simParam->haveThermIntLambda()){ |
59 |
< |
tIntLambda_ = simParam->getThermIntLambda(); |
58 |
> |
if (simParam->haveThermodynamicIntegrationLambda()){ |
59 |
> |
tIntLambda_ = simParam->getThermodynamicIntegrationLambda(); |
60 |
|
} |
61 |
|
else{ |
62 |
|
tIntLambda_ = 1.0; |
69 |
|
simError(); |
70 |
|
} |
71 |
|
|
72 |
< |
if (simParam->haveThermIntK()){ |
73 |
< |
tIntK_ = simParam->getThermIntK(); |
72 |
> |
if (simParam->haveThermodynamicIntegrationK()){ |
73 |
> |
tIntK_ = simParam->getThermodynamicIntegrationK(); |
74 |
|
} |
75 |
|
else{ |
76 |
|
tIntK_ = 1.0; |
144 |
|
tempTau = curSnapshot->statData.getTau(); |
145 |
|
tempTau *= factor_; |
146 |
|
curSnapshot->statData.setTau(tempTau); |
147 |
< |
|
147 |
> |
|
148 |
|
// do crystal restraint forces for thermodynamic integration |
149 |
|
if (simParam->getUseSolidThermInt()) { |
150 |
|
|