Implemented specialised version of jacobian for vectorial base function set

[[Imported from SVN: r3510]]

fem/common/basefunctionsets.cc

@@ -66,6 +66,37 @@ namespace Dune {
     evaluate(baseFunct, diffVar, quad.point(quadPoint), phi);
   }
 
+  template <class FunctionSpaceImp>
+  void VectorialBaseFunctionSet<FunctionSpaceImp>::
+  jacobian(int baseFunct, const DomainType& xLocal,
+           JacobianRangeType& gradPhi) const
+  {
+    gradPhi *= 0.;
+
+    for (int i = 0; i < FunctionSpaceImp::DimDomain; ++i) {
+      diffVar1_[0] = i;
+      baseFunctions_[util_.containedDof(baseFunct)]->
+      evaluate(diffVar1_, xLocal, tmp_);
+      gradPhi[util_.component(baseFunct)][i] = tmp_[0];
+    }
+  }
+
+  template <class FunctionSpaceImp>
+  template <class QuadratureImp>
+  void VectorialBaseFunctionSet<FunctionSpaceImp>::
+  jacobian(int baseFunct, QuadratureImp& quad, int quadPoint,
+           JacobianRangeType& gradPhi) const
+  {
+    gradPhi *= 0.;
+
+    for (int i = 0; i < FunctionSpaceImp::DimDomain; ++i) {
+      diffVar1_[0] = i;
+      baseFunctions_[util_.containedDof(baseFunct)]->
+      evaluate(diffVar1_, quad.point(quadPoint), tmp_);
+      gradPhi[util_.component(baseFunct)][i] = tmp_[0];
+    }
+  }
+
   template <class FunctionSpaceImp>
   typename VectorialBaseFunctionSet<FunctionSpaceImp>::DofType
   VectorialBaseFunctionSet<FunctionSpaceImp>::

fem/common/basefunctionsets.hh

@@ -126,8 +126,6 @@ namespace Dune {
       }
     }
 
-    // * override other functions as well!!!!
-
     int numBaseFunctions() const;
 
     template <int diffOrd>
@@ -142,6 +140,13 @@ namespace Dune {
                   QuadratureType & quad,
                   int quadPoint, RangeType & phi ) const;
 
+    void jacobian(int baseFunct, const DomainType& xLocal,
+                  JacobianRangeType& gradPhi) const;
+
+    template <class QuadratureImp>
+    void jacobian(int baseFunct, QuadratureImp& quad, int quadPoint,
+                  JacobianRangeType& gradPhi) const;
+
     // * add those methods with quadratures as well
     DofType evaluateSingle(int baseFunct,
                            const DomainType& xLocal,