Skip to content
Snippets Groups Projects
Commit d08424eb authored by Stephan Hilb's avatar Stephan Hilb
Browse files

fix missing term in residual error estimator

parent bf78de49
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
scripts/run_experiments.jl
+ 2
− 1
View file @ d08424eb
...@@ -619,7 +619,8 @@ function estimate_res!(st::L1L2TVState) ...@@ -619,7 +619,8 @@ function estimate_res!(st::L1L2TVState)
cellf(hcell; g, u, nablau, p1, p2, tdata) = cellf(hcell; g, u, nablau, p1, p2, tdata) =
st.alpha2 * st.alpha2 *
st.T(adjoint, tdata, st.alpha2 * st.T(tdata, u) - g) + 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) = facetf(hfacet, n; g, u, nablau, p1, p2, tdata) =
#display(n' * p2) #display(n' * p2)
# norm2 is calculated later after both facet contributions are # norm2 is calculated later after both facet contributions are
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment