anne-glerum

Hi @laila-16 , thank you for your PR. Could you explain why you think the default value should be changed? Is it not enough to be able to change the...

> according to the master thesis that was mentioned here #4496 , the growth of the time step should be limited by a factor of 1.86, changing the default value...

Thanks for the comments @naliboff , I've addressed them where possible. Maybe @bobmyhill , @jdannberg or @gassmoeller can also give this a look?

> 1. Can the stress_component_statistics still be removed, per a comment in the .cc file? If it does need to be included, I will review that as well. I think...

I'm not sure it's the same problem, but I just ran the gmg_mesh_def_ghost_nodes.prm test for #4248 and this also shows significant differences in RMS and max vel and iteration counts...

> Yes, that looks like the same, although the differences are not as large, but still look significant. What was your residual? The test just does a single Stokes solve,...

Thank you Valentina!

@naliboff For completeness, did you also check the Newton solver?

During t0 the reaction terms are not computed, so no changes in the strain would be accumulated anyway. If the Picard Stokes solver does several nonlinear iterations, is there a...

Where needed rheologies already set masks for the fields they know to _not_ be used as material fractions, but those rheologies do control the evolution of those fields. If there...