Abstract: To solve geodynamic problem, the computational geodynamic community has to solve large systems of nonlinear equations. This is mostly done by the Picard fixed-point iteration, but the Newton iteration is now starting to get used more. We found some fundamental problems with the Newton iteration for the Stokes equation, and have found ways to deal with them efficiently. We will furthermore show how this plays out for real use case examples.