From 4ce0f3919755d191acf54d786a65e69f9e35cd32 Mon Sep 17 00:00:00 2001
From: Robert K <robertk@posteo.org>
Date: Sun, 14 Feb 2021 13:57:59 +0100
Subject: [PATCH] [bugfix][femDGOperator] Fix case where limiter is set but
parameters are empty for simplex grids.
---
python/dune/femdg/_operators.py | 12 ++++++++++--
1 file changed, 10 insertions(+), 2 deletions(-)
diff --git a/python/dune/femdg/_operators.py b/python/dune/femdg/_operators.py
index b182b6c5..1fff8b12 100644
--- a/python/dune/femdg/_operators.py
+++ b/python/dune/femdg/_operators.py
@@ -217,12 +217,20 @@ def femDGOperator(Model, space,
limiter = "minmod"
limiterstr = limiter if space.grid.type.isSimplex else "lp"
# force default values for how reconstruction is done
- parameters["femdg.limiter.admissiblefunctions"] = "default"
+ if parameters is None:
+ from dune.fem import parameter
+ parameter.append({"femdg.limiter.admissiblefunctions":"default"})
+ else:
+ parameters["femdg.limiter.admissiblefunctions"] = "default"
if limiter.lower() == "lp":
limiter = "minmod"
# force default values for how reconstruction is done
- parameters["femdg.limiter.admissiblefunctions"] = "lp"
+ if parameters is None:
+ from dune.fem import parameter
+ parameter.append({"femdg.limiter.admissiblefunctions":"lp"})
+ else:
+ parameters["femdg.limiter.admissiblefunctions"] = "lp"
if limiter.lower() == "scaling":
# check for scaling limiter interface
--
GitLab