Skip to content
Snippets Groups Projects
Commit 2a89b1e8 authored by Lambert Theisen's avatar Lambert Theisen
Browse files

Change minor things without effect in solver.py

parent ef6cc5b4
No related branches found
No related tags found
3 merge requests!5Update to version 1.1,!3Feat/symmetrization,!2Feat/symmetrization
# pylint: disable=invalid-name # pylint: disable=invalid-name,too-many-lines
""" """
Solver module, contains the Solver class. Solver module, contains the Solver class.
...@@ -321,6 +321,15 @@ class Solver: ...@@ -321,6 +321,15 @@ class Solver:
(p, u, sigma) = df.TrialFunctions(w_stress) (p, u, sigma) = df.TrialFunctions(w_stress)
(q, v, psi) = df.TestFunctions(w_stress) (q, v, psi) = df.TestFunctions(w_stress)
# Setup source functions
f_heat = self.heat_source
f_mass = self.mass_source
if self.mode == "r13":
cpl = 1.0
else:
cpl = 0.0
# Setup projections # Setup projections
s_n = df.dot(s, n) s_n = df.dot(s, n)
r_n = df.dot(r, n) r_n = df.dot(r, n)
...@@ -333,15 +342,6 @@ class Solver: ...@@ -333,15 +342,6 @@ class Solver:
sigma_nt = df.dot(sigma*n, t) sigma_nt = df.dot(sigma*n, t)
psi_nt = df.dot(psi*n, t) psi_nt = df.dot(psi*n, t)
# Setup source functions
f_heat = self.heat_source
f_mass = self.mass_source
if self.mode == "r13":
cpl = 1.0
else:
cpl = 0.0
# Sub functionals: # Sub functionals:
def b(scalar, vector): def b(scalar, vector):
return 1 * scalar * df.div(vector) * df.dx return 1 * scalar * df.div(vector) * df.dx
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment