From 980db06e311fee62bc5b9c2253417f642627ed56 Mon Sep 17 00:00:00 2001
From: DennisGlaeser <dennis.glaeser@iws.uni-stuttgart.de>
Date: Fri, 12 Jan 2018 15:31:31 +0100
Subject: [PATCH] [cclocalassembler] use more generic computation of
maxNeighbors
The previous implementation failed for mpfa on surface/unstructured grids.
---
dumux/assembly/cclocalassembler.hh | 5 ++++-
1 file changed, 4 insertions(+), 1 deletion(-)
diff --git a/dumux/assembly/cclocalassembler.hh b/dumux/assembly/cclocalassembler.hh
index 4724492c4f..09ae189798 100644
--- a/dumux/assembly/cclocalassembler.hh
+++ b/dumux/assembly/cclocalassembler.hh
@@ -35,6 +35,7 @@
#include <dumux/assembly/diffmethod.hh>
#include <dumux/assembly/numericdifferentiation.hh>
#include <dumux/assembly/fvlocalassemblerbase.hh>
+#include <dumux/discretization/fluxstencil.hh>
namespace Dumux {
@@ -137,14 +138,16 @@ class CCLocalAssembler<TypeTag, Assembler, DiffMethod::numeric, /*implicit=*/tru
using LocalResidualValues = typename GET_PROP_TYPE(TypeTag, NumEqVector);
using Element = typename GET_PROP_TYPE(TypeTag, GridView)::template Codim<0>::Entity;
using ElementSolutionVector = typename GET_PROP_TYPE(TypeTag, ElementSolutionVector);
+ using FVElementGeometry = typename GET_PROP_TYPE(TypeTag, FVElementGeometry);
using GridVariables = typename GET_PROP_TYPE(TypeTag, GridVariables);
using JacobianMatrix = typename GET_PROP_TYPE(TypeTag, JacobianMatrix);
enum { numEq = GET_PROP_VALUE(TypeTag, NumEq) };
enum { dim = GET_PROP_TYPE(TypeTag, GridView)::dimension };
+ using FluxStencil = Dumux::FluxStencil<TypeTag>;
+ static constexpr int maxNeighbors = FluxStencil::maxFluxStencilSize*FVElementGeometry::maxNumElementScvfs;
static constexpr bool enableGridFluxVarsCache = GET_PROP_VALUE(TypeTag, EnableGridFluxVariablesCache);
- static constexpr int maxNeighbors = 4*(2*dim);
public:
--
GitLab