From d08424eb6e797e5ab646f6614cfcb963fc59f854 Mon Sep 17 00:00:00 2001
From: Stephan Hilb <stephan@ecshi.net>
Date: Fri, 25 Feb 2022 23:35:52 +0100
Subject: [PATCH] fix missing term in residual error estimator
---
scripts/run_experiments.jl | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/scripts/run_experiments.jl b/scripts/run_experiments.jl
index ce1cc2c..62ecfe4 100644
--- a/scripts/run_experiments.jl
+++ b/scripts/run_experiments.jl
@@ -619,7 +619,8 @@ function estimate_res!(st::L1L2TVState)
cellf(hcell; g, u, nablau, p1, p2, tdata) =
st.alpha2 *
st.T(adjoint, tdata, st.alpha2 * st.T(tdata, u) - g) +
- st.T(adjoint, tdata, p1)
+ st.T(adjoint, tdata, p1) +
+ (isSnabla ? zero(p1) : st.beta * u)
facetf(hfacet, n; g, u, nablau, p1, p2, tdata) =
#display(n' * p2)
# norm2 is calculated later after both facet contributions are
--
GitLab