diff --git a/dune/gfe/riemanniantrsolver.cc b/dune/gfe/riemanniantrsolver.cc index 236b732c7cca7b891b1254ca168e35a9e41b2bb1..b86fd19bdedc1749466dcb0b55b8a16c41bd2b31 100644 --- a/dune/gfe/riemanniantrsolver.cc +++ b/dune/gfe/riemanniantrsolver.cc @@ -240,21 +240,16 @@ void RiemannianTrustRegionSolver<GridType,TargetSpace>::solve() if (recomputeGradientHessian) { - double oldClock = omp_get_wtime(); - assembler_->assembleGradient(x_, rhs); rhs *= -1; // The right hand side is the _negative_ gradient std::cout << "gradient assembly took " << gradientTimer.elapsed() << " sec." << std::endl; - std::cout << "gradient assembly took " << (omp_get_wtime() - oldClock) << " sec." << std::endl; gradientTimer.reset(); - oldClock = omp_get_wtime(); assembler_->assembleMatrix(x_, *hessianMatrix_, i==0 // assemble occupation pattern only for the first call ); std::cout << "hessian assembly took " << gradientTimer.elapsed() << " sec." << std::endl; - std::cout << "hessian assembly took " << (omp_get_wtime() - oldClock) << " sec." << std::endl; recomputeGradientHessian = false; }