tOnly restore mu_wall if rate-limited - cngf-pf - continuum model for granular flows with pore-pressure dynamics (renamed from 1d_fd_simple_shear)
HTML git clone git://src.adamsgaard.dk/cngf-pf
DIR Log
DIR Files
DIR Refs
DIR README
DIR LICENSE
---
DIR commit 274f48009ff4183a39733aa2fec335097f3da033
DIR parent 62bcffb0514d7c3517e4d43ea22648e78d2f4763
HTML Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Mon, 26 Aug 2019 14:19:26 +0200
Only restore mu_wall if rate-limited
Diffstat:
M main.c | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
---
DIR diff --git a/main.c b/main.c
t@@ -338,7 +338,8 @@ main(int argc, char* argv[])
} while ((!isnan(sim.v_x_fix) || !isnan(sim.v_x_limit))
&& fabs(res_norm) > RTOL);
- sim.mu_wall = mu_wall_orig;
+ if (!isnan(sim.v_x_limit))
+ sim.mu_wall = mu_wall_orig;
sim.t += sim.dt;
filetimeclock += sim.dt;
iter++;