diff --git a/dune/gfe/riemanniantrsolver.cc b/dune/gfe/riemanniantrsolver.cc
index 532b0361b46e87267ec6f36490f9293f18285f69..7f57ff7afa7c60b7d3164bf43631310c4fb942cc 100644
--- a/dune/gfe/riemanniantrsolver.cc
+++ b/dune/gfe/riemanniantrsolver.cc
@@ -340,7 +340,9 @@ void RiemannianTrustRegionSolver<GridType,TargetSpace>::solve()
 
             std::cout << "Solve quadratic problem..." << std::endl;
 
+            Dune::Timer solutionTimer;
             innerSolver_->solve();
+            std::cout << "Solving the quadratic problem took " << solutionTimer.elapsed() << " seconds." << std::endl;
 
             if (mgStep)
                 corr_global = mgStep->getSol();