127 |
|
} |
128 |
|
|
129 |
|
void FluctuatingChargeLangevin::applyConstraints() { |
130 |
+ |
|
131 |
+ |
if (!hasFlucQ_) return; |
132 |
+ |
|
133 |
|
SimInfo::MoleculeIterator i; |
134 |
|
Molecule::FluctuatingChargeIterator j; |
135 |
|
Molecule* mol; |
136 |
|
Atom* atom; |
137 |
< |
RealType cvel, cpos, cfrc, cmass, randomForce, frictionForce, velstep; |
137 |
> |
RealType cvel, cfrc, cmass, randomForce, frictionForce; |
138 |
|
RealType velStep, oldFF; // used to test for convergence |
139 |
|
|
140 |
|
for (mol = info_->beginMolecule(i); mol != NULL; |
143 |
|
atom = mol->nextFluctuatingCharge(j)) { |
144 |
|
|
145 |
|
cvel = atom->getFlucQVel(); |
143 |
– |
cpos = atom->getFlucQPos(); |
146 |
|
cfrc = atom->getFlucQFrc(); |
147 |
|
cmass = atom->getChargeMass(); |
148 |
|
|