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
Tip revision: 110bd1cb786a3753c46228e6da40f7beedcd71c3 authored by Thorsten Becker on 16 January 2010, 02:47:04 UTC
Implementation of two, perhaps entirely unnecessary features, in
Implementation of two, perhaps entirely unnecessary features, in
Tip revision: 110bd1c
File | Mode | Size |
---|---|---|
checkpoint | ||
stokes | ||
Makefile.am | -rw-r--r-- | 1.0 KB |
array2d.cc | -rw-r--r-- | 3.9 KB |
coupled.cfg | -rw-r--r-- | 1.8 KB |
exchange.py | -rwxr-xr-x | 5.0 KB |
gnomonic.c | -rw-r--r-- | 2.2 KB |
multicoupled.cfg | -rw-r--r-- | 2.4 KB |
signon.py | -rwxr-xr-x | 2.4 KB |
test1.sh | -rwxr-xr-x | 2.3 KB |
test2.sh | -rwxr-xr-x | 4.0 KB |
test5.sh | -rwxr-xr-x | 2.4 KB |
Computing file changes ...