diff --git a/dune/solvers/solvers/quadraticipopt.hh b/dune/solvers/solvers/quadraticipopt.hh
index ff32c6955381affb74c21410b20c6cb37d66f966..41837e4ab2453824a482016629fb80c8d0a69c4c 100644
--- a/dune/solvers/solvers/quadraticipopt.hh
+++ b/dune/solvers/solvers/quadraticipopt.hh
@@ -574,6 +574,7 @@ void QuadraticIPOptSolver<MatrixType,VectorType>::solve()
   app->Options()->SetIntegerValue("max_iter", this->maxIterations_);
   app->Options()->SetStringValue("mu_strategy", "adaptive");
   app->Options()->SetStringValue("output_file", "ipopt.out");
+  app->Options()->SetStringValue("hessian_constant", "yes");
   switch (this->verbosity_) {
   case NumProc::QUIET:
       app->Options()->SetIntegerValue("print_level", 0);