diff --git a/dune/gfe/riemanniantrsolver.cc b/dune/gfe/riemanniantrsolver.cc index 22b8353202a10953251a02f55663d3f106f49178..1403c64af1c98c20547f90fef58d493269efff9b 100644 --- a/dune/gfe/riemanniantrsolver.cc +++ b/dune/gfe/riemanniantrsolver.cc @@ -473,7 +473,7 @@ void RiemannianTrustRegionSolver<GridType,TargetSpace>::solve() } if (energy >= oldEnergy && - (std::abs(oldEnergy-energy)/energy < 1e-9 || relativeModelDecrease < 1e-9)) { + (std::abs((oldEnergy-energy)/energy) < 1e-9 || relativeModelDecrease < 1e-9)) { if (this->verbosity_ == NumProc::FULL and rank==0) std::cout << "Suspecting rounding problems" << std::endl;