Migrate from dune-fufem to dune-functions bases

Migrate from dune-fufem to dune-functions bases.

dune/gfe/surfacecosseratenergy.hh

@@ -203,12 +203,7 @@ public:
     // Length scale parameter
     L_c_ = parameters.template get<double>("L_c");
 
-    // Curvature exponent
-    q_ = parameters.template get<double>("q");
-
-    // Shear correction factor
-    kappa_ = parameters.template get<double>("kappa");
-
+    // Curvature parameters
     b1_ = parameters.template get<double>("b1");
     b2_ = parameters.template get<double>("b2");
     b3_ = parameters.template get<double>("b3");
@@ -278,7 +273,7 @@ RT energy(const typename Basis::LocalView& localView,
     auto id = idSet.subId(it.inside(), it.indexInInside(), 1);
     auto boundaryGeometry = geometriesOnShellBoundary_.at(id);
     auto quadOrder = (it.type().isSimplex()) ? deformationLocalFiniteElement.localBasis().order()
-                                                  : deformationLocalFiniteElement.localBasis().order() * gridDim;
+                                                  : deformationLocalFiniteElement.localBasis().order() * boundaryDim;
 
     const auto& quad = Dune::QuadratureRules<DT, boundaryDim>::rule(it.type(), quadOrder);
     for (size_t pt=0; pt<quad.size(); pt++) {
@@ -505,12 +500,6 @@ private:
   /** \brief Length scale parameter */
   double L_c_;
 
-  /** \brief Curvature exponent */
-  double q_;
-
-  /** \brief Shear correction factor */
-  double kappa_;
-
   /** \brief Curvature parameters */
   double b1_, b2_, b3_;