Skip to content

Commit

Permalink
set constant_jacobian for mass matrices in SIQN forcing solvers
Browse files Browse the repository at this point in the history
  • Loading branch information
JHopeCollins committed Nov 4, 2024
1 parent 7c6a607 commit 9b2a545
Showing 1 changed file with 4 additions and 2 deletions.
6 changes: 4 additions & 2 deletions gusto/timestepping/semi_implicit_quasi_newton.py
Original file line number Diff line number Diff line change
Expand Up @@ -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 = {}
Expand Down

0 comments on commit 9b2a545

Please sign in to comment.