GitLab

There were a couple of issues observed with the Newton solver:

1. AbsoluteLimit lower than 1e-8 is not observable in the (single precision) VTK output.
   • ...but in mass conservation, see #39 (closed).
2. If AbsoluteLimit is too high, the solver will not perform an iteration for small time steps and weak dynamics.
3. We don't use Reduction.

Propositions:

• AbsoluteLimit = 1e-10
• Reduction = 1e-4: To save time, allow a lower defect if the reduction is sufficient
• ForceIteration = true: make sure that we don't skip dynamics