From d2b4aa318a764d8a6c46ee5f3034190e2e5371f3 Mon Sep 17 00:00:00 2001 From: Daniel Shapero Date: Mon, 17 May 2021 08:30:46 -0700 Subject: [PATCH] Increase the max solver iterations This was necessary to make the solver converge only on the first call after changes to icepack to regularize Glen's flow law around zero strain rate. --- PigInversion.py | 6 +++++- 1 file changed, 5 insertions(+), 1 deletion(-) diff --git a/PigInversion.py b/PigInversion.py index b8c91e5..6b26527 100755 --- a/PigInversion.py +++ b/PigInversion.py @@ -670,7 +670,11 @@ def main(): if inversionParams['initWithDeg1']: Q1 = firedrake.FunctionSpace(mesh, family='CG', degree=1) area = firedrake.assemble(firedrake.Constant(1) * firedrake.dx(mesh)) - opts = {'dirichlet_ids': meshOpts['dirichlet_ids']} # Opts from mesh + opts = { + 'dirichlet_ids': meshOpts['dirichlet_ids'], ## Opts from mesh + 'diagnostic_solver_parameters': {'max_iterations': 100} + } + # # Input model geometry and velocity zb, s, h, floating, grounded = \