1d_fd_simple_shear_transient

commit 58cbd5fd997a4f217da227a1336b991edebbf3df
parent a6079cd8f16165bafbdd72d9ed5f456b9069430e
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Mon, 29 Apr 2019 19:51:45 +0200

Increase solver tolerance criteria

Diffstat:
M
1d_fd_simple_shear.png
 | 
0

M
1d_fd_simple_shear_rheology.png
 | 
0

M
1d_fd_simple_shear_rheology_iverson.png
 | 
0

M
1d_fd_simple_shear_rheology_kamb.png
 | 
0

M
1d_fd_simple_shear_rheology_tulaczyk.png
 | 
0

M
diurnal.gif
 | 
0

M
main.c
 | 
4
++--

7 files changed, 2 insertions(+), 2 deletions(-)
diff --git a/1d_fd_simple_shear.png b/1d_fd_simple_shear.png
Binary files differ.
diff --git a/1d_fd_simple_shear_rheology.png b/1d_fd_simple_shear_rheology.png
Binary files differ.
diff --git a/1d_fd_simple_shear_rheology_iverson.png b/1d_fd_simple_shear_rheology_iverson.png
Binary files differ.
diff --git a/1d_fd_simple_shear_rheology_kamb.png b/1d_fd_simple_shear_rheology_kamb.png
Binary files differ.
diff --git a/1d_fd_simple_shear_rheology_tulaczyk.png b/1d_fd_simple_shear_rheology_tulaczyk.png
Binary files differ.
diff --git a/diurnal.gif b/diurnal.gif
Binary files differ.
diff --git a/main.c b/main.c
@@ -236,7 +236,7 @@ int main(int argc, char* argv[])
     while (sim.t <= sim.t_end) {
 
         if (sim.fluid) {
-            if (darcy_solver_1d(&sim, 10000, 1e-3))
+            if (darcy_solver_1d(&sim, 10000, 1e-5))
                 exit(1);
 #ifdef DEBUG
             puts(".. p_f_ghost:"); print_array(sim.p_f_ghost, sim.nz+2);
@@ -257,7 +257,7 @@ int main(int argc, char* argv[])
         puts(".. mu:"); print_array(sim.mu, sim.nz);
 #endif
 
-        if (implicit_1d_jacobian_poisson_solver(&sim, 10000, 1e-3))
+        if (implicit_1d_jacobian_poisson_solver(&sim, 10000, 1e-5))
             exit(1);
 
         compute_shear_strain_rate_plastic(&sim);