diff --git a/dune/gfe/rodassembler.hh b/dune/gfe/rodassembler.hh
index ce794c7138d6a8d5bc52be72496a2dc8167d35e3..32c5caf31d9e4e67992c411c8bad83b39c24426a 100644
--- a/dune/gfe/rodassembler.hh
+++ b/dune/gfe/rodassembler.hh
@@ -68,8 +68,8 @@ public:
         {   return  dynamic_cast<RodLocalStiffness<GridView, double>* >(this->localStiffness_)->referenceConfiguration_;
         }
 
-        void assembleGradient(const std::vector<RigidBodyMotion<double,3> >& sol,
-                              Dune::BlockVector<Dune::FieldVector<double, blocksize> >& grad) const;
+  virtual void assembleGradient(const std::vector<RigidBodyMotion<double,3> >& sol,
+                                Dune::BlockVector<Dune::FieldVector<double, blocksize> >& grad) const override;
 
         void getStrain(const std::vector<RigidBodyMotion<double,3> >& sol,
                        Dune::BlockVector<Dune::FieldVector<double, blocksize> >& strain) const;
@@ -137,11 +137,11 @@ public:
     void assembleMatrix(const std::vector<RigidBodyMotion<double,2> >& sol,
                         Dune::BCRSMatrix<MatrixBlock>& matrix);
 
-    void assembleGradient(const std::vector<RigidBodyMotion<double,2> >& sol,
-                          Dune::BlockVector<Dune::FieldVector<double, blocksize> >& grad) const;
+    virtual void assembleGradient(const std::vector<RigidBodyMotion<double,2> >& sol,
+                          Dune::BlockVector<Dune::FieldVector<double, blocksize> >& grad) const override;
 
     /** \brief Compute the energy of a deformation state */
-    double computeEnergy(const std::vector<RigidBodyMotion<double,2> >& sol) const;
+    virtual double computeEnergy(const std::vector<RigidBodyMotion<double,2> >& sol) const override;
 
 protected:
 
diff --git a/dune/gfe/rodlocalstiffness.hh b/dune/gfe/rodlocalstiffness.hh
index f922e16dc090e75000853b4fd9e3820c5ef97170..b37bdadcfd2759a8292fe97552c771a0e25b092b 100644
--- a/dune/gfe/rodlocalstiffness.hh
+++ b/dune/gfe/rodlocalstiffness.hh
@@ -101,7 +101,7 @@ public:
                        const std::array<RigidBodyMotion<RT,3>, dim+1>& localSolution) const;
 
     virtual RT energy (const typename Basis::LocalView& localView,
-                       const std::vector<RigidBodyMotion<RT,3> >& localSolution) const
+                       const std::vector<RigidBodyMotion<RT,3> >& localSolution) const override
     {
         assert(localSolution.size()==2);
         std::array<RigidBodyMotion<RT,3>, 2> localSolutionArray = {localSolution[0], localSolution[1]};