Main bulk of stuff needed to update strain rate in momentum diagnostics. I've tested that the code compiles and runs, but need to tidy the code and verify that it also work in sync with the viscosity, as planned.
Minor comment clarification.
Some tidying up following my chat to Dr. Wilson. Need to test this next in a non-linear run, to verify that things are being updates as I think they are.
Merge trunk into branch.
Merge in trunk.
Merge in trunk.
Merge in trunk.
Merge remote-tracking branch 'fluidity/master' into fix_non_linear_nnewt_viscosity
Ensure that SecondInvariant routine is called for all materials during non lienar iterations. Also add some write statements.
Add a 1 time-step, multi-material test (case 1a of the slab detachment benchmark), which checks viscosity values at end of time-step. These values are updated every non-linear iteration, dependent upon the strain-rate second invariant.
Change from GAMG to MG due to petsc limitations of buildbot.
Merge in master.
Add viscous dissipation check to this test.
Fix python script and adjust test values.
Minor tweek. Silly me.
Exploit symmetry in system to reduce problem size by a factor of 2.
Make max viscosity check relative to a given value, as suggested by Cian.
Make max viscosity check relative to a given value, for all tests, as suggested by Cian.
Further simply test. Use checkpointed mesh and remove mesh adaptivity. Also - less non linear iterations so that test runs quicker.
Forgot makefile.