From 746b409fd79d9ddd39902d487207f41bb925d8f6 Mon Sep 17 00:00:00 2001
From: Kilian Weishaupt <kilian.weishaupt@iws.uni-stuttgart.de>
Date: Wed, 27 Oct 2021 13:13:34 +0200
Subject: [PATCH] [staggered][fvElementGeometry] Use fixed index map for
lateralOrthogonal scvf
---
.../facecentered/staggered/fvelementgeometry.hh | 4 +++-
1 file changed, 3 insertions(+), 1 deletion(-)
diff --git a/dumux/discretization/facecentered/staggered/fvelementgeometry.hh b/dumux/discretization/facecentered/staggered/fvelementgeometry.hh
index 1d7987f2c2..598e234290 100644
--- a/dumux/discretization/facecentered/staggered/fvelementgeometry.hh
+++ b/dumux/discretization/facecentered/staggered/fvelementgeometry.hh
@@ -79,7 +79,9 @@ public:
const SubControlVolumeFace& lateralOrthogonalScvf(const SubControlVolumeFace& scvf) const
{
assert(scvf.isLateral());
- return gridGeometry().scvf(gridGeometry().lateralOrthogonalScvf(scvf));
+ const auto otherGlobalIdx = scvfIndices_()[GridGeometry::GeometryHelper::lateralOrthogonalScvfLocalIndex(scvf.localIndex())];
+ return gridGeometry().scvf(otherGlobalIdx);
+
}
//! Return the frontal sub control volume face on a the boundary for a given sub control volume
--
GitLab