From 7d37258b757b8a5fc820804bc11d361665bdab3e Mon Sep 17 00:00:00 2001
From: Timo Koch <timo.koch@iws.uni-stuttgart.de>
Date: Thu, 2 Sep 2021 12:11:46 +0200
Subject: [PATCH] [staggered][bugfix] Reinitialize maps in update
If update is called twice the connectivity map is wrong sinnce
we push additional entries in there, i.e. each node in the stencil
will appear twice.
This makes sure all vectors are correctly reinitalized in every update
---
.../staggered/freeflow/connectivitymap.hh | 10 ++++++----
1 file changed, 6 insertions(+), 4 deletions(-)
diff --git a/dumux/discretization/staggered/freeflow/connectivitymap.hh b/dumux/discretization/staggered/freeflow/connectivitymap.hh
index 75efceeeb6..26215028bb 100644
--- a/dumux/discretization/staggered/freeflow/connectivitymap.hh
+++ b/dumux/discretization/staggered/freeflow/connectivitymap.hh
@@ -63,10 +63,12 @@ public:
const auto numDofsCC = gridGeometry.gridView().size(0);
const auto numDofsFace = gridGeometry.gridView().size(1);
const auto numBoundaryFacets = gridGeometry.numBoundaryScvf();
- cellCenterToCellCenterMap_.resize(numDofsCC);
- cellCenterToFaceMap_.resize(numDofsCC);
- faceToCellCenterMap_.resize(2*numDofsFace - numBoundaryFacets);
- faceToFaceMap_.resize(2*numDofsFace - numBoundaryFacets);
+
+ // reinitialize maps
+ cellCenterToCellCenterMap_ = Map(numDofsCC);
+ cellCenterToFaceMap_ = Map(numDofsCC);
+ faceToCellCenterMap_ = Map(2*numDofsFace - numBoundaryFacets);
+ faceToFaceMap_ = Map(2*numDofsFace - numBoundaryFacets);
// restrict the FvGeometry locally
auto fvGeometry = localView(gridGeometry);
--
GitLab