Skip to content
Snippets Groups Projects
Commit 37dacbab authored by Timo Koch's avatar Timo Koch
Browse files

[upwind] Get reference to elemvolvars first

parent aa96a4a3
No related branches found
No related tags found
1 merge request!922Fix/box elemvolvars access
...@@ -55,11 +55,12 @@ public: ...@@ -55,11 +55,12 @@ public:
// TODO: pass this from outside? // TODO: pass this from outside?
static const Scalar upwindWeight = getParam<Scalar>("Implicit.UpwindWeight"); static const Scalar upwindWeight = getParam<Scalar>("Implicit.UpwindWeight");
const auto& elemVolVars = fluxVars.elemVolVars();
const auto& scvf = fluxVars.scvFace(); const auto& scvf = fluxVars.scvFace();
const auto& insideScv = fluxVars.fvGeometry().scv(scvf.insideScvIdx()); const auto& insideScv = fluxVars.fvGeometry().scv(scvf.insideScvIdx());
const auto& outsideScv = fluxVars.fvGeometry().scv(scvf.outsideScvIdx()); const auto& outsideScv = fluxVars.fvGeometry().scv(scvf.outsideScvIdx());
const auto& insideVolVars = fluxVars.elemVolVars()[insideScv]; const auto& insideVolVars = elemVolVars[insideScv];
const auto& outsideVolVars = fluxVars.elemVolVars()[outsideScv]; const auto& outsideVolVars = elemVolVars[outsideScv];
if (std::signbit(flux)) // if sign of flux is negative if (std::signbit(flux)) // if sign of flux is negative
return flux*(upwindWeight*upwindTerm(outsideVolVars) return flux*(upwindWeight*upwindTerm(outsideVolVars)
......
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