1d_fd_simple_shear

Continuum model for granular flows with pore-pressure dynamics
git clone git://src.adamsgaard.dk/1d_fd_simple_shear
Log | Files | Refs | README | LICENSE

commit 6b4638289f974b82cb8da373267c1ba86770bae2
parent 76960076e66459c6f146cc0a3b33b9bdccc2c6d6
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Sat,  6 Apr 2019 12:26:52 +0200

Clean-up script

Diffstat:
Mjulia/1d_fd_simple_shear.jl | 10+---------
1 file changed, 1 insertion(+), 9 deletions(-)

diff --git a/julia/1d_fd_simple_shear.jl b/julia/1d_fd_simple_shear.jl @@ -180,11 +180,6 @@ function implicit_1d_jacobian_poisson_solver(g, p, μ, Δz, set_bc_dirichlet(g, "-z") set_bc_neumann(g, "+z") - # Damsgaard et al 2013 include fixed grains in plot - # for i=1:9 - # set_bc_dirichlet(g, "-z", idx_offset=i) - # end - if verbose println("g after BC: ") println(g) @@ -192,7 +187,6 @@ function implicit_1d_jacobian_poisson_solver(g, p, μ, Δz, # perform a single jacobi iteration in each cell for iz=1:length(p) - #for iz=10:length(p) poisson_solver_1d_iteration(g, g_out, r_norm, μ, p, iz, Δz) end @@ -209,9 +203,7 @@ function implicit_1d_jacobian_poisson_solver(g, p, μ, Δz, # stop iterating if the relative tolerance is satisfied if r_norm_max <= rel_tol - #if verbose - @info ".. Solution converged after $iter iterations" - #end + @info ".. Solution converged after $iter iterations" return end end