Skip to content
Snippets Groups Projects
Commit 4ce0f391 authored by Robert K's avatar Robert K
Browse files

[bugfix][femDGOperator] Fix case where limiter is set but parameters are

empty for simplex grids.
parent 8f1d8306
No related branches found
No related tags found
No related merge requests found
......@@ -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
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment