Commit 0398c3f3 authored by Dennis Gläser's avatar Dennis Gläser
Browse files

[fluxvars] return outflow fluxes directly in upwind scheme for branching points

This should reduce the computational effort
parent a9ec57b7
......@@ -113,9 +113,9 @@ public:
Scalar branchingPointUpwindTerm = 0.0;
Scalar sumUpwindFluxes = 0.0;
// the inside flux
// if the inside flux is positive (outflow) do fully upwind and return flux
if (!std::signbit(flux))
branchingPointUpwindTerm += upwindTerm(insideVolVars)*flux;
return upwindTerm(insideVolVars)*flux;
else
sumUpwindFluxes += flux;
......
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