From d96b819c587000b0036945dc2580c5efd87e18ad Mon Sep 17 00:00:00 2001
From: David Seus <david.seus@ians.uni-stuttgart.de>
Date: Fri, 30 Aug 2019 11:25:21 +0200
Subject: [PATCH] add FEM polynomial degree parameter
---
LDDsimulation/LDDsimulation.py | 7 +++++--
1 file changed, 5 insertions(+), 2 deletions(-)
diff --git a/LDDsimulation/LDDsimulation.py b/LDDsimulation/LDDsimulation.py
index baba899..9959bbf 100644
--- a/LDDsimulation/LDDsimulation.py
+++ b/LDDsimulation/LDDsimulation.py
@@ -33,7 +33,8 @@ class LDDsimulation(object):
tol: float = None,#
LDDsolver_tol: float = 1E-8,
# maximal number of L-iterations that the LDD solver uses.
- max_iter_num: int = 1000):
+ max_iter_num: int = 1000,
+ FEM_Lagrange_degree: int = 1):
hlp.print_once("\n ###############################################\n",
"############# Begin Simulation ################\n",
"###############################################\n")
@@ -48,6 +49,7 @@ class LDDsimulation(object):
# dimension of the simulation. This is by default 2 as I don't intend
# to implement 3D.
self.dim = dimension
+ self.FEM_Lagrange_degree = FEM_Lagrange_degree
#Parameters
# df.parameters["allow_extrapolation"] = True
# df.parameters["refinement_algorithm"] = "plaza_with_parent_facets"
@@ -1107,7 +1109,8 @@ class LDDsimulation(object):
saturation = self.saturation[subdom_num],#
timestep_size = self.timestep_size,#
interpolation_degree = self.interpolation_degree,#
- tol = self.tol)
+ tol = self.tol,
+ degree=self.FEM_Lagrange_degree)
}
)
# setup the linear solvers and set the solver parameters.
--
GitLab