| 35 | 
  | 
 *                                                                       | 
| 36 | 
  | 
 * [1]  Meineke, et al., J. Comp. Chem. 26, 252-271 (2005).              | 
| 37 | 
  | 
 * [2]  Fennell & Gezelter, J. Chem. Phys. 124, 234104 (2006).           | 
| 38 | 
< | 
 * [3]  Sun, Lin & Gezelter, J. Chem. Phys. 128, 24107 (2008).           | 
| 38 | 
> | 
 * [3]  Sun, Lin & Gezelter, J. Chem. Phys. 128, 234107 (2008).           | 
| 39 | 
  | 
 * [4]  Kuang & Gezelter,  J. Chem. Phys. 133, 164101 (2010). | 
| 40 | 
  | 
 * [5]  Vardeman, Stocker & Gezelter, J. Chem. Theory Comput. 7, 834 (2011). | 
| 41 | 
  | 
 */ | 
| 290 | 
  | 
 | 
| 291 | 
  | 
    for (unsigned int i = 0; i < atoms_.size(); i++) { | 
| 292 | 
  | 
       | 
| 293 | 
+ | 
      atype = atoms_[i]->getAtomType(); | 
| 294 | 
+ | 
 | 
| 295 | 
  | 
      afrc = atoms_[i]->getFrc(); | 
| 296 | 
  | 
      apos = atoms_[i]->getPos(); | 
| 297 | 
  | 
      rpos = apos - pos; | 
| 309 | 
  | 
        atrq = atoms_[i]->getTrq(); | 
| 310 | 
  | 
        trq += atrq; | 
| 311 | 
  | 
      } | 
| 312 | 
+ | 
 | 
| 313 | 
  | 
      if ((sl & DataStorage::dslElectricField) && (atype->isElectrostatic())) { | 
| 314 | 
  | 
        ef += atoms_[i]->getElectricField(); | 
| 315 | 
  | 
        eCount++; |