diff --git a/dune/elasticity/assemblers/neohookeoperatorassembler.hh b/dune/elasticity/assemblers/neohookeoperatorassembler.hh
index 7b4062d8838dad91c643d3dd0d6561fd3269202b..192ece1f8f55cd01985dd56ec5537903da74cb79 100644
--- a/dune/elasticity/assemblers/neohookeoperatorassembler.hh
+++ b/dune/elasticity/assemblers/neohookeoperatorassembler.hh
@@ -192,7 +192,7 @@ private:
      *  Note that it is assumed the defGrad is the deformation gradient and not the displacement gradient   
      */
     Dune::FieldMatrix<ctype,dim,dim> hessDefDet(const Dune::FieldMatrix<ctype,dim,dim>& defGrad,
-                            const Dune::FieldVector<ctype,dim>& testGrad, const Dune::FieldVector<ctype,dim>& ansatzGrad)
+                            const Dune::FieldVector<ctype,dim>& testGrad, const Dune::FieldVector<ctype,dim>& ansatzGrad) const
     {
         Dune::FieldMatrix<ctype, dim, dim> res(0);