| 119 | 
  | 
    std::vector<SDTuple3>::iterator i; | 
| 120 | 
  | 
    Vector3d normal(0.0, 0.0, 1.0); | 
| 121 | 
  | 
    double scd = 0.0; | 
| 122 | 
< | 
    for (i = tuples_.begin(); i != tuples_.end(); ++i) { | 
| 123 | 
< | 
        Vector3d zAxis = i->third->getPos() - i->first->getPos(); | 
| 122 | 
> | 
    for (i = tuples_.begin(); i != tuples_.end(); ++i) {         | 
| 123 | 
> | 
        //Egberts B. and Berendsen H.J.C, J.Chem.Phys. 89(6), 3718-3732, 1988 | 
| 124 | 
> | 
         | 
| 125 | 
> | 
        Vector3d zAxis =   i->third->getPos() - i->first->getPos(); | 
| 126 | 
  | 
        snapshot->wrapVector(zAxis); | 
| 127 | 
  | 
        Vector3d v12 =  i->second->getPos() - i->first->getPos(); | 
| 128 | 
  | 
        snapshot->wrapVector(v12); |