Skip to content
Snippets Groups Projects
Commit 1a895655 authored by Tobias Leibner's avatar Tobias Leibner
Browse files

[gamm-2019-talk-on-conservative-rb.cc] fix compilation

parent 42a4d40a
No related branches found
No related tags found
No related merge requests found
......@@ -282,7 +282,9 @@ compute_estimate(const GP& grid,
const auto pressure_grad = p->jacobian(xx)[0];
const auto t_val = t->evaluate(xx);
auto difference = diff * pressure_grad + t_val;
return (diff_inv * difference) * difference;
auto diff_inv_times_diff = difference;
diff_inv.mv(difference, diff_inv_times_diff);
return diff_inv_times_diff * difference;
},
std::max(df->order() + std::max(p->order() - 1, 0), t->order()) + over_integrate);
// compute indicators and estimator
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment