--- trunk/OOPSE/libmdtools/Vector3d.cpp 2004/08/09 14:50:35 1451 +++ trunk/OOPSE/libmdtools/Vector3d.cpp 2004/08/23 15:11:36 1452 @@ -12,6 +12,24 @@ Vector3d& Vector3d::operator*= ( const Mat3x3d & m){ return *this; } +Mat3x3d Vector3d::makeSkewMat(){ + Mat3x3d m; + + m.element[0][0] = 0; + m.element[0][1] = z; + m.element[0][2] = -y; + + m.element[1][0] = -z; + m.element[1][1] = 0; + m.element[1][2] = x; + + m.element[2][0] = y; + m.element[2][1] = -x; + m.element[2][2] = 0; + + return m; +} + Mat3x3d outProduct(const Vector3d& v1, const Vector3d& v2){ Mat3x3d m; @@ -29,3 +47,4 @@ Mat3x3d outProduct(const Vector3d& v1, const Vector3d& return m; } +