From 2dc4055c22ae1c0c2dad3b3888aadf01acbf6e46 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Thu, 28 Aug 2014 08:08:49 +0000
Subject: [PATCH] Only warn when energy increases

It appears that there may be very small energy increases due to
rounding errors.

[[Imported from SVN: r9867]]
---
 dune/gfe/mixedriemanniantrsolver.cc | 5 +++--
 1 file changed, 3 insertions(+), 2 deletions(-)

diff --git a/dune/gfe/mixedriemanniantrsolver.cc b/dune/gfe/mixedriemanniantrsolver.cc
index 05422631..32cf12df 100644
--- a/dune/gfe/mixedriemanniantrsolver.cc
+++ b/dune/gfe/mixedriemanniantrsolver.cc
@@ -474,8 +474,9 @@ void MixedRiemannianTrustRegionSolver<GridType,Basis0,TargetSpace0,Basis1,Target
 
               std::cout << "Energy: " << energy << std::endl;
 
-              if (energy >= oldEnergy)
-                DUNE_THROW(Dune::Exception, "energy increase!");
+              if (energy > oldEnergy)
+                //DUNE_THROW(Dune::Exception, "energy increase!");
+                std::cout << "Warning: energy increase!" << std::endl;
 
               oldEnergy = energy;
             }
-- 
GitLab