https://github.com/geodynamics/citcoms
Revision 110bd1cb786a3753c46228e6da40f7beedcd71c3 authored by Thorsten Becker on 16 January 2010, 02:47:04 UTC, committed by Thorsten Becker on 16 January 2010, 02:47:04 UTC
search of why certain prescribed velocity models do not converge like
they should.

- allow for each inner loop solution to have a higher accuracy than
the outer, Uzawa loop. This is set by new parameter
"inner_accuracy_scale", which pre-multiplies the inner accuracy (imp)
setting (defaults to unity)

- remove any rigid body rotation at each Uzawa iteration
step. Parmaeter "inner_remove_rigid_rotation", by default off. For
this, I also added an experimental routine assign_v_to_vector() in
order to reasign the V velocities to the U solution vector after
removing the net rotations from V (how is supposed to work for the
regular operational mode?)




1 parent 33e0d0b
History
Tip revision: 110bd1cb786a3753c46228e6da40f7beedcd71c3 authored by Thorsten Becker on 16 January 2010, 02:47:04 UTC
Implementation of two, perhaps entirely unnecessary features, in
Tip revision: 110bd1c

back to top