| 1 |
< |
/* |
| 1 |
> |
/* |
| 2 |
|
* Copyright (c) 2005 The University of Notre Dame. All Rights Reserved. |
| 3 |
|
* |
| 4 |
|
* The University of Notre Dame grants you ("Licensee") a |
| 43 |
|
|
| 44 |
|
namespace oopse { |
| 45 |
|
|
| 46 |
< |
GhostTorsion::GhostTorsion(Atom *atom1, Atom *atom2, DirectionalAtom* ghostAtom, |
| 47 |
< |
TorsionType *tt) : Torsion(atom1, atom2, ghostAtom, ghostAtom, tt) {} |
| 46 |
> |
GhostTorsion::GhostTorsion(Atom *atom1, Atom *atom2, DirectionalAtom* ghostAtom, |
| 47 |
> |
TorsionType *tt) : Torsion(atom1, atom2, ghostAtom, ghostAtom, tt) {} |
| 48 |
|
|
| 49 |
< |
void GhostTorsion::calcForce() { |
| 50 |
< |
DirectionalAtom* ghostAtom = static_cast<DirectionalAtom*>(atom2_); |
| 49 |
> |
void GhostTorsion::calcForce() { |
| 50 |
> |
DirectionalAtom* ghostAtom = static_cast<DirectionalAtom*>(atom3_); |
| 51 |
|
|
| 52 |
|
Vector3d pos1 = atom1_->getPos(); |
| 53 |
|
Vector3d pos2 = atom2_->getPos(); |
| 92 |
|
|
| 93 |
|
f3.negate(); |
| 94 |
|
ghostAtom->addTrq(cross(r43, f3)); |
| 95 |
< |
} |
| 95 |
> |
} |
| 96 |
|
|
| 97 |
|
} |
| 98 |
|
|