From cdd2982589044c83f035b69142a3bad88ae40cbd Mon Sep 17 00:00:00 2001 From: Josh Hope-Collins Date: Thu, 7 Nov 2024 13:28:43 +0000 Subject: [PATCH] Do not refactor mass matrices in SIQN transport calculation (#572) --- gusto/timestepping/semi_implicit_quasi_newton.py | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/gusto/timestepping/semi_implicit_quasi_newton.py b/gusto/timestepping/semi_implicit_quasi_newton.py index 2c49304a..76c517a7 100644 --- 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 = {}