| 88 |
|
Molecule* mol; |
| 89 |
|
Atom* atom; |
| 90 |
|
|
| 91 |
< |
for (mol = info_->beginMolecule(i); mol != NULL; |
| 92 |
< |
mol = info_->nextMolecule(i)) { |
| 93 |
< |
for (atom = mol->beginFluctuatingCharge(j); atom != NULL; |
| 94 |
< |
atom = mol->nextFluctuatingCharge(j)) { |
| 95 |
< |
cerr << " starts with = " << atom->getFlucQPos() << "\n"; |
| 96 |
< |
//atom->setFlucQPos(0.0); |
| 97 |
< |
//atom->setFlucQVel(0.0); |
| 98 |
< |
} |
| 99 |
< |
} |
| 91 |
> |
// For single-minima flucq, this ensures a net neutral system, but |
| 92 |
> |
// for multiple minima, this is no longer the right thing to do: |
| 93 |
> |
// |
| 94 |
> |
// for (mol = info_->beginMolecule(i); mol != NULL; |
| 95 |
> |
// mol = info_->nextMolecule(i)) { |
| 96 |
> |
// for (atom = mol->beginFluctuatingCharge(j); atom != NULL; |
| 97 |
> |
// atom = mol->nextFluctuatingCharge(j)) { |
| 98 |
> |
// atom->setFlucQPos(0.0); |
| 99 |
> |
// atom->setFlucQVel(0.0); |
| 100 |
> |
// } |
| 101 |
> |
// } |
| 102 |
|
|
| 103 |
|
FluctuatingChargeObjectiveFunction flucQobjf(info_, forceMan_, |
| 104 |
|
fqConstraints_); |
