do not refactor mass matrices in SIQN transport calculation

firedrakeproject · Nov 7, 2024 · dfae07d · dfae07d
1 parent bccf8b9
commit dfae07d
Showing 1 changed file with 4 additions and 2 deletions.
diff --git a/gusto/timestepping/semi_implicit_quasi_newton.py b/gusto/timestepping/semi_implicit_quasi_newton.py
@@ -531,11 +531,13 @@ def __init__(self, equation, alpha):
 
         # now we can set up the explicit and implicit problems
         explicit_forcing_problem = LinearVariationalProblem(
-            a.form, L_explicit.form, self.xF, bcs=bcs
+            a.form, L_explicit.form, self.xF, bcs=bcs,
+            constant_jacobian=True
         )
 
         implicit_forcing_problem = LinearVariationalProblem(
-            a.form, L_implicit.form, self.xF, bcs=bcs
+            a.form, L_implicit.form, self.xF, bcs=bcs,
+            constant_jacobian=True
         )
 
         self.solvers = {}