From e436bb5f4177f0ae9d4c35cbe5af56b66d8eeeb5 Mon Sep 17 00:00:00 2001
From: Ned Coltman <edward.coltman@iws.uni-stuttgart.de>
Date: Wed, 16 Jan 2019 16:39:38 +0100
Subject: [PATCH] [scvf] fix the coupled ghost face index issue.
---
dumux/discretization/staggered/freeflow/subcontrolvolumeface.hh | 2 +-
dumux/freeflow/navierstokes/staggered/fluxvariables.hh | 2 +-
2 files changed, 2 insertions(+), 2 deletions(-)
diff --git a/dumux/discretization/staggered/freeflow/subcontrolvolumeface.hh b/dumux/discretization/staggered/freeflow/subcontrolvolumeface.hh
index d9f06de4a8..0fc4fc8e9e 100644
--- a/dumux/discretization/staggered/freeflow/subcontrolvolumeface.hh
+++ b/dumux/discretization/staggered/freeflow/subcontrolvolumeface.hh
@@ -108,7 +108,7 @@ public:
selfToOppositeDistance_(0.0),
dirIdx_(dirIdx),
isGhostFace_(true)
- {}
+ { axisData_.selfDof = dofIdx; }
//! The center of the sub control volume face
const GlobalPosition& center() const
diff --git a/dumux/freeflow/navierstokes/staggered/fluxvariables.hh b/dumux/freeflow/navierstokes/staggered/fluxvariables.hh
index 46ca4213af..8bc1dc638e 100644
--- a/dumux/freeflow/navierstokes/staggered/fluxvariables.hh
+++ b/dumux/freeflow/navierstokes/staggered/fluxvariables.hh
@@ -654,7 +654,7 @@ private:
//! helper function to conveniently create a ghost face used to retrieve boundary values from the problem
SubControlVolumeFace makeGhostFace_(const SubControlVolumeFace& ownScvf, const GlobalPosition& pos) const
{
- return SubControlVolumeFace(pos, std::vector<unsigned int>{ownScvf.insideScvIdx(), ownScvf.outsideScvIdx()}, ownScvf.directionIndex(), ownScvf.dofIndex(), ownScvf.index());
+ return SubControlVolumeFace(pos, std::vector<unsigned int>{ownScvf.insideScvIdx(), ownScvf.outsideScvIdx()}, ownScvf.directionIndex(), ownScvf.axisData().selfDof, ownScvf.index());
};
//! helper function to conveniently create a ghost face which is outside the domain, parallel to the scvf of interest
--
GitLab