From 34164e8f6e63a2b32b22dc4bc23614ada2d20ac1 Mon Sep 17 00:00:00 2001
From: Mathis Kelm <mathis.kelm@iws.uni-stuttgart.de>
Date: Tue, 27 Sep 2022 17:57:59 +0200
Subject: [PATCH] [assembly] Use fvGeometry.numScvf() to size residual/jacobian
---
dumux/assembly/fclocalassembler.hh | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/dumux/assembly/fclocalassembler.hh b/dumux/assembly/fclocalassembler.hh
index f09a175c57..33af2d0bcb 100644
--- a/dumux/assembly/fclocalassembler.hh
+++ b/dumux/assembly/fclocalassembler.hh
@@ -358,7 +358,7 @@ public:
//////////////////////////////////////////////////////////////////////////////////////////////////
// one residual per element facet
- const auto numElementResiduals = element.subEntities(1);
+ const auto numElementResiduals = fvGeometry.numScv();
// create the vector storing the partial derivatives
ElementResidualVector partialDerivs(numElementResiduals);
@@ -573,7 +573,7 @@ public:
auto elemSol = elementSolution(element, curSol, fvGeometry.gridGeometry());
// create the vector storing the partial derivatives
- ElementResidualVector partialDerivs(element.subEntities(1));
+ ElementResidualVector partialDerivs(fvGeometry.numScv());
// calculation of the derivatives
for (const auto& scv : scvs(fvGeometry))
--
GitLab