dune/gfe/surfacecosseratenergy.hh

@@ -228,12 +228,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");
@@ -282,7 +277,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()) ? localFiniteElement.localBasis().order()
-                                                  : localFiniteElement.localBasis().order() * gridDim;
+                                                  : localFiniteElement.localBasis().order() * boundaryDim;
 
     const auto& quad = Dune::QuadratureRules<DT, boundaryDim>::rule(it.type(), quadOrder);
     for (size_t pt=0; pt<quad.size(); pt++) {
@@ -485,12 +480,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_;