Commit 233bd1b9 authored by Lukas Riedel's avatar Lukas Riedel

Merge branch '206-flux-reconstruction-on-fv-does-not-work-with-alpha-volumes' into 'master'

Resolve "Flux reconstruction on FV does not work with alpha volumes"

Closes #206

See merge request !219
parents 9592528a 55c83f72
......@@ -27,6 +27,7 @@
* Ensure UTF-8 locale definition in Docker images !210
* `grid.extensions` values were erroneously read when loading GMSH grids #193
* Compiling failed on macOS because C++11 was selected as standard #198
* Flux reconstruction now considers volume elements without DOFs !219
## 2.0.0 (2020-05-14)
......
......@@ -385,8 +385,10 @@ public:
rl_view.setWeight(local_assembler.weight());
Dune::PDELab::LocalAssemblerCallSwitch<LOP,LOP::doAlphaVolume>::
alpha_volume(lop,eg,lfsw,pl,lfsv_volume,rl_view);
// Under certain circumstances volume element has no DOFs (e.g. finite volumes)
if (lfsv_volume.size() > 0)
Dune::PDELab::LocalAssemblerCallSwitch<LOP,LOP::doAlphaVolume>::
alpha_volume(lop,eg,lfsw,pl,lfsv_volume,rl_view);
// Load residual view into the local residual vector
for (unsigned int i = 0; i < lfsv_volume.size(); ++i)
......@@ -737,4 +739,4 @@ private:
} // namespace Dorie
} // namespace Dune
#endif // DUNE_DORIE_RAVIART_THOMAS_RESIDUALENGINE_HH
\ No newline at end of file
#endif // DUNE_DORIE_RAVIART_THOMAS_RESIDUALENGINE_HH
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