--- trunk/OOPSE/libmdtools/Mat3x3d.cpp	2004/06/09 16:16:33	1254
+++ trunk/OOPSE/libmdtools/Mat3x3d.cpp	2004/06/11 17:16:21	1268
@@ -186,7 +186,6 @@ Euler3 Mat3x3d::toEuler(){
   //z-axis (again). 
 
   Euler3 e;
-  Mat3x3d m;
   double cosTheta;
   double sinTheta;
   const double eps = 1.0e-8;
@@ -226,3 +225,14 @@ Euler3 Mat3x3d::toEuler(){
 
   return e;
 }
+
+
+Vector3d operator*(const Mat3x3d& m, const Vector3d& v){
+  Vector3d result;
+  
+  result.x = m.element[0][0] * v.x + m.element[0][1] * v.y + m.element[0][2]*v.z;
+  result.x = m.element[1][0] * v.x + m.element[1][1] * v.y + m.element[1][2]*v.z;
+  result.x = m.element[2][0] * v.x + m.element[2][1] * v.y + m.element[2][2]*v.z;
+
+  return result;
+}