--- trunk/src/primitives/GhostTorsion.cpp 2010/06/28 20:07:12 1458 +++ trunk/src/primitives/GhostTorsion.cpp 2010/07/06 21:23:46 1459 @@ -56,7 +56,7 @@ namespace OpenMD { Vector3d r21 = pos1 - pos2; Vector3d r32 = pos2 - pos3; - Vector3d r43 = ghostAtom->getA().getColumn(2); + Vector3d r43 = ghostAtom->getA().transpose().getColumn(2); // Calculate the cross products and distances Vector3d A = cross(r21, r32);