diff --git a/dune/gfe/densities/bulkcosseratdensity.hh b/dune/gfe/densities/bulkcosseratdensity.hh
index d2a6d9757b20e79abb00a894a56e7b9bb8079211..4067354914e20d7bba238799c930f1bfaa711d27 100644
--- a/dune/gfe/densities/bulkcosseratdensity.hh
+++ b/dune/gfe/densities/bulkcosseratdensity.hh
@@ -231,7 +231,7 @@ namespace Dune::GFE
     }
 
     // Construct a copy of this density but using 'adouble' as the number type
-    virtual std::unique_ptr<LocalDensity<Element,ATargetSpace> > makeActiveDensity() const
+    virtual std::unique_ptr<LocalDensity<Element,ATargetSpace> > makeActiveDensity() const override
     {
       // curvatureType_ is a local enum type, and therefore its type changes
       // together with the type of the surrounding class.
diff --git a/dune/gfe/densities/cosseratvolumeloaddensity.hh b/dune/gfe/densities/cosseratvolumeloaddensity.hh
index 99c6a62928072e4cfc8498be4df3e883eeafacee..8c6a7ca071b3c60354b4635eb6cc1ebb8ce9c107 100644
--- a/dune/gfe/densities/cosseratvolumeloaddensity.hh
+++ b/dune/gfe/densities/cosseratvolumeloaddensity.hh
@@ -52,7 +52,7 @@ namespace Dune::GFE
     }
 
     // Construct a copy of this density but using 'adouble' as the number type
-    virtual std::unique_ptr<LocalDensity<ElementOrIntersection,ATargetSpace> > makeActiveDensity() const
+    virtual std::unique_ptr<LocalDensity<ElementOrIntersection,ATargetSpace> > makeActiveDensity() const override
     {
       auto result = std::make_unique<CosseratVolumeLoadDensity<ElementOrIntersection,adouble> >(volumeLoad_);
 
diff --git a/dune/gfe/densities/harmonicdensity.hh b/dune/gfe/densities/harmonicdensity.hh
index 455360b8582af90454c5c264a56f9a5a63bd9864..945849a020d93a97485ef56616fd64232424289d 100644
--- a/dune/gfe/densities/harmonicdensity.hh
+++ b/dune/gfe/densities/harmonicdensity.hh
@@ -63,7 +63,7 @@ namespace Dune::GFE
     }
 
     // Construct a copy of this density but using 'adouble' as the number type
-    virtual std::unique_ptr<LocalDensity<ElementOrIntersection,ATargetSpace> > makeActiveDensity() const
+    virtual std::unique_ptr<LocalDensity<ElementOrIntersection,ATargetSpace> > makeActiveDensity() const override
     {
       auto result = std::make_unique<HarmonicDensity<ElementOrIntersection,ATargetSpace> >();
       if (this->elementOrIntersection_)