tFix writes for Dirichlet BC and fix allocation of extra fluidity array - 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 fd7d2ebbea9d3079fd4cbd210b7bf081eb3f7764
DIR parent 60815244104ff136da823b895613c0979ecbffcf
HTML Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Wed, 10 Apr 2019 17:25:37 +0200
Fix writes for Dirichlet BC and fix allocation of extra fluidity array
Diffstat:
M simulation.c | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
---
DIR diff --git a/simulation.c b/simulation.c
t@@ -108,7 +108,7 @@ void set_bc_dirichlet(
if (boundary == -1)
g_ghost[0] = value;
else if (boundary == +1)
- g_ghost[idx1g(nz)+1] = value;
+ g_ghost[nz+1] = value;
else {
fprintf(stderr, "set_bc_dirichlet: Unknown boundary %d\n", boundary);
exit(1);
t@@ -153,7 +153,7 @@ int implicit_1d_jacobian_poisson_solver(
const int max_iter,
const double rel_tol)
{
- double* g_ghost_out = empty(sim->nz);
+ double* g_ghost_out = empty(sim->nz+2);
double* r_norm = empty(sim->nz);
int iter;