[Richards] Use same local for flux reconstruction

parent 90acc43d
......@@ -530,16 +530,15 @@ public:
std::shared_ptr<GFFluxReconstruction> gf_ptr;
if constexpr (enable_rt_engine)
{
// copy local operator
SLOP _slop = *slop;
// create a flux reconstruction grid function
gf_ptr = std::make_shared<GFFluxReconstruction>(this->_log,
gv,inifile.sub("fluxReconstruction"));
_slop.setTime(state.time);
slop->setTime(state.time);
// update it with the state
gf_ptr->update(*(state.coefficients),*(state.grid_function_space),_slop);
gf_ptr->update(*(state.coefficients),*(state.grid_function_space),*slop);
slop->setTime(current_time());
} else {
DUNE_THROW(NotImplemented,
"Flux reconstruction is not implemented for the selected "
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment