added multAssignTransposed.

[[Imported from SVN: r3046]]

common/fmatrix.hh

@@ -1466,6 +1466,22 @@ namespace Dune {
       }
     }
 
+    //! calculates ret = matrix * x
+    template <typename K, int dim>
+    static void multAssignTransposed( const FieldMatrix<K,dim,dim> &matrix, const FieldVector<K,dim> & x, FieldVector<K,dim> & ret)
+    {
+      typedef typename FieldMatrix<K,dim,dim>::size_type size_type;
+
+      for(size_type i=0; i<dim; i++)
+      {
+        ret[i] = 0.0;
+        for(size_type j=0; j<dim; j++)
+        {
+          ret[i] += matrix[j][i]*x[j];
+        }
+      }
+    }
+
     //! calculates ret = matrix * x
     template <typename K, int dim>
     static FieldVector<K,dim> mult(const FieldMatrix<K,dim,dim> &matrix, const FieldVector<K,dim> & x)