From 48f8ec3aaeed5a239b699516e7305db73b70d502 Mon Sep 17 00:00:00 2001
From: Lisa Julia Nebel <lisa_julia.nebel@tu-dresden.de>
Date: Fri, 4 Sep 2020 18:52:31 +0200
Subject: [PATCH] Close the ADOL-C tape even if the energy calculation throws
an error in mixedlocalgfeadolcstiffness
---
dune/gfe/mixedlocalgfeadolcstiffness.hh | 8 ++++++--
1 file changed, 6 insertions(+), 2 deletions(-)
diff --git a/dune/gfe/mixedlocalgfeadolcstiffness.hh b/dune/gfe/mixedlocalgfeadolcstiffness.hh
index ddf89f96..c668f3ac 100644
--- a/dune/gfe/mixedlocalgfeadolcstiffness.hh
+++ b/dune/gfe/mixedlocalgfeadolcstiffness.hh
@@ -122,10 +122,14 @@ energy(const typename Basis::LocalView& localView,
}
using namespace Dune::TypeTree::Indices;
- energy = localEnergy_->energy(localView,
+ try {
+ energy = localEnergy_->energy(localView,
localAConfiguration0,
localAConfiguration1);
-
+ } catch (Dune::Exception &e) {
+ trace_off();
+ throw e;
+ }
energy >>= pureEnergy;
trace_off();
--
GitLab