| 63 |
|
} |
| 64 |
|
|
| 65 |
|
FluctuatingChargePropagator::~FluctuatingChargePropagator() { |
| 66 |
– |
if (fqConstraints_ != NULL) delete fqConstraints_; |
| 66 |
|
} |
| 67 |
|
|
| 68 |
|
void FluctuatingChargePropagator::setForceManager(ForceManager* forceMan) { |
| 92 |
|
mol = info_->nextMolecule(i)) { |
| 93 |
|
for (atom = mol->beginFluctuatingCharge(j); atom != NULL; |
| 94 |
|
atom = mol->nextFluctuatingCharge(j)) { |
| 95 |
< |
atom->setFlucQPos(0.0); |
| 96 |
< |
atom->setFlucQVel(0.0); |
| 95 |
> |
cerr << " starts with = " << atom->getFlucQPos() << "\n"; |
| 96 |
> |
//atom->setFlucQPos(0.0); |
| 97 |
> |
//atom->setFlucQVel(0.0); |
| 98 |
|
} |
| 99 |
|
} |
| 100 |
|
|
| 112 |
|
DumpStatusFunction dsf(info_); // we want a dump file written |
| 113 |
|
// every iteration |
| 114 |
|
minim->minimize(problem, endCriteria); |
| 115 |
+ |
cerr << "back from minim\n"; |
| 116 |
|
initialized_ = true; |
| 117 |
|
} |
| 118 |
|
|