diff --git a/dune/gfe/riemanniantrsolver.cc b/dune/gfe/riemanniantrsolver.cc
index 42dc81f283d96219d60d24d3c955f1947c85a26b..1f53e565ad5abfb5560b098ae973c9003d37804b 100644
--- a/dune/gfe/riemanniantrsolver.cc
+++ b/dune/gfe/riemanniantrsolver.cc
@@ -331,7 +331,12 @@ void RiemannianTrustRegionSolver<Basis,TargetSpace>::solve()
     //   Trust-Region Solver
     // /////////////////////////////////////////////////////
 
+    Dune::Timer energyTimer;
     double oldEnergy = assembler_->computeEnergy(x_);
+    if (this->verbosity_ == Solver::FULL)
+        std::cout << "Energy computation took " << energyTimer.elapsed() << " sec." << std::endl;
+
+
     oldEnergy = grid_->comm().sum(oldEnergy);
 
     bool recomputeGradientHessian = true;