diff --git a/dune/gfe/bendingisometryhelper.hh b/dune/gfe/bendingisometryhelper.hh
index c3075959de53fb93281317b0e5702c6a055506c7..8db7bf173becd76960361798c0c037329eb91de5 100644
--- a/dune/gfe/bendingisometryhelper.hh
+++ b/dune/gfe/bendingisometryhelper.hh
@@ -332,7 +332,6 @@ namespace Dune::GFE::Impl
          */
         size_t edgeIndex =  representationLFE.localCoefficients().localKey(i).subEntity();
         auto edgeGeometry = element.template subEntity<1>(edgeIndex).geometry();
-        auto edgeEntity = element.template subEntity<1>(edgeIndex);
 
         /**
          * @brief Get the right orientation for tangent and normal vectors on current edge.
diff --git a/dune/gfe/functions/discretekirchhoffbendingisometry.hh b/dune/gfe/functions/discretekirchhoffbendingisometry.hh
index a7296ae5869fe5b81df2b661b0a2d5b4b8b314e6..98e022be979cb93b344bd35d94810cf605846fdf 100644
--- a/dune/gfe/functions/discretekirchhoffbendingisometry.hh
+++ b/dune/gfe/functions/discretekirchhoffbendingisometry.hh
@@ -61,8 +61,8 @@ namespace Dune::GFE
     DiscreteKirchhoffBendingIsometry(const DiscreteKirchhoffBasis& discreteKirchhoffBasis,
                                      const CoefficientBasis& coefficientBasis,
                                      Coefficients& globalIsometryCoefficients)
-      : basis_(discreteKirchhoffBasis),
-      coefficientBasis_(coefficientBasis),
+      : coefficientBasis_(coefficientBasis),
+      basis_(discreteKirchhoffBasis),
       localView_(basis_),
       localViewCoefficient_(coefficientBasis_),
       globalIsometryCoefficients_(globalIsometryCoefficients)
diff --git a/src/bending-plate.cc b/src/bending-plate.cc
index 6edfea3ada4f781068ac87e8da0248b4907b83da..e17d8c6ae16fd785f40b47c2fe5cd20afd1cad3d 100644
--- a/src/bending-plate.cc
+++ b/src/bending-plate.cc
@@ -261,7 +261,6 @@ int main(int argc, char *argv[])
   Dune::GFE::RiemannianTrustRegionSolver<CoefficientBasis, Coefficient> RTRsolver;
 
   std::string Solver_name = parameterSet.get<std::string>("Solver", "RNHM");
-  double numerical_energy; //final discrete energy
 
   if(Solver_name == "RNHM")
   {
@@ -273,7 +272,6 @@ int main(int argc, char *argv[])
 
     RNHMsolver.solve();
     isometryCoefficients = RNHMsolver.getSol();
-    numerical_energy = RNHMsolver.getStatistics().finalEnergy;
   } else if (Solver_name =="RiemannianTR")
   {
     std::cout << "Using Riemannian Trust-region method for energy minimization." << std::endl;
@@ -285,7 +283,6 @@ int main(int argc, char *argv[])
 
     RTRsolver.solve();
     isometryCoefficients = RTRsolver.getSol();
-    numerical_energy = RTRsolver.getStatistics().finalEnergy;
   } else
     DUNE_THROW(Dune::Exception, "Unknown Solver type for bending isometries.");