From 2c75859cbcfcd7b9ee667a7fb5436a96e14cb192 Mon Sep 17 00:00:00 2001 From: Anton Khirnov Date: Fri, 8 Feb 2019 16:44:54 +0100 Subject: ell_grid_solve: do not read from possibly-uninitialized boundary layer --- ell_grid_solve.c | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/ell_grid_solve.c b/ell_grid_solve.c index f992186..9d50c82 100644 --- a/ell_grid_solve.c +++ b/ell_grid_solve.c @@ -419,7 +419,7 @@ static int solve_exact(EGSContext *ctx) if (!is_bnd[bnd_loc] || bnd->type != MG2D_BC_TYPE_FIXVAL) continue; - for (ptrdiff_t bnd_layer = 1; bnd_layer <= ctx->fd_stencil; bnd_layer++) + for (ptrdiff_t bnd_layer = 1; bnd_layer < ctx->fd_stencil; bnd_layer++) for (ptrdiff_t bnd_idx = -(ptrdiff_t)ctx->fd_stencil; bnd_idx < (ptrdiff_t)(bnd_len[0] + ctx->fd_stencil); bnd_idx++) { e->rhs[mat_row_idx] -= bnd->val[bnd_layer * bnd->val_stride + bnd_idx] * dst[dir * bnd_layer * bnd_stride[1] + bnd_idx * bnd_stride[0]]; dst[dir * bnd_layer * bnd_stride[1] + bnd_idx * bnd_stride[0]] = 0.0; -- cgit v1.2.3