Commit 7d40e0d6 authored by Dennis Gläser's avatar Dennis Gläser Committed by Timo Koch
Browse files

[mpfao][iavol] save neumann flux calculation in interaction volumes that do not touch boundary

parent b28f810f
......@@ -297,7 +297,7 @@ public:
Scalar getNeumannFlux(const LocalIndexPair& localIndexPair) const
{
if (fluxScvfIndexSet_().size() == 0 || GET_PROP_VALUE(TypeTag, UseTpfaBoundary))
if (!onBoundary() || fluxScvfIndexSet_().size() == 0 || GET_PROP_VALUE(TypeTag, UseTpfaBoundary))
return 0.0;
auto flux = CAinv_[localIndexPair.first] * neumannFluxes_;
......
Supports Markdown
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