added multTransposed to helper functions.

[[Imported from SVN: r2809]]

common/fmatrix.hh

@@ -1424,6 +1424,23 @@ namespace Dune {
       return ret;
     }
 
+    //! calculates ret = matrix^T * x
+    template <typename K, int dim>
+    static FieldVector<K,dim> multTransposed(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];
+        }
+      }
+      return ret;
+    }
+
   } // end namespace FMatrixHelp
 
 #ifdef DUNE_EXPRESSIONTEMPLATES