77 |
|
Vector3d trq; |
78 |
|
Vector3d EFfrc; |
79 |
|
Vector3d pos; |
80 |
– |
RealType chrg; |
81 |
– |
RealType pot, fieldPot; |
82 |
– |
RealType chrgToKcal = 23.0609; |
83 |
– |
RealType debyeToKcal = 4.8018969509; |
80 |
|
|
85 |
– |
|
81 |
|
if (doElectricField) { |
82 |
< |
fieldPot = 0.0; |
82 |
> |
const RealType chrgToKcal = 23.0609; |
83 |
> |
const RealType debyeToKcal = 4.8018969509; |
84 |
> |
RealType pot; |
85 |
> |
RealType fieldPot = 0.0; |
86 |
|
|
87 |
|
for (mol = info_->beginMolecule(i); mol != NULL; |
88 |
|
mol = info_->nextMolecule(i)) { |
91 |
|
atom = mol->nextAtom(j)) { |
92 |
|
|
93 |
|
bool isCharge = false; |
94 |
< |
chrg = 0.0; |
94 |
> |
RealType chrg = 0.0; |
95 |
|
|
96 |
|
AtomType* atype = atom->getAtomType(); |
97 |
|
|