Skip to content

GitLab

Commit 7a28cfae authored by David Vincent's avatar David Vincent
Browse files

modifying the tolerances of newton solver for implicit time scheme: relative tolerance is now 1e-8

parent 06ed68ad
Pipeline #1496 passed with stage
in 29 minutes and 5 seconds
Hide whitespace changes
Inline Side-by-side
Showing with 2 additions and 2 deletions
dgpy/scripts/slim.py
View file @ 7a28cfae
......@@ -815,7 +815,7 @@ class Loop:
equation._temporal_solver = dgpy.dgDIRK(equation._equation, equation._sys.dof_manager, 2)
equation._temporal_solver.getNewton().setVerb(10)
equation._temporal_solver.getNewton().setAtol(1e-8)
equation._temporal_solver.getNewton().setRtol(1e-4)
equation._temporal_solver.getNewton().setRtol(1e-8)
elif equation._temporal_scheme == "multirate":
equation._temporal_solver = dgpy.dgMultirateERK(equation._domain._groups, equation._equation, dgpy.ERK_22_C)
equation._dt_mr = equation._temporal_solver.splitGroupsForMultirate(15, equation._solution, [equation._solution, equation._domain._bath])
......@@ -853,7 +853,7 @@ class Loop:
equation._temporal_solver = dgpy.dgMDIRK(equation._equation, equation._sys_tracer1.dof_manager, equation._sys_tracer2.dof_manager, 2)
equation._temporal_solver.getNewton().setVerb(10)
equation._temporal_solver.getNewton().setAtol(1.e-8)
equation._temporal_solver.getNewton().setRtol(1e-4)
equation._temporal_solver.getNewton().setRtol(1e-8)
else:
dgpy.Msg.Fatal("Temporal scheme '"+equation._temporal_scheme+"' not implemented!")
self._number_of_tracer += 1
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment