find and fix lower-bound bug in c/ch12/solns/elasto.c

bueler · May 7, 2023 · 38e49bf · 38e49bf
1 parent b5cb6b3
commit 38e49bf
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/c/ch12/solns/elasto.c b/c/ch12/solns/elasto.c
@@ -132,7 +132,7 @@ PetscErrorCode FormBounds(SNES snes, Vec Xl, Vec Xu) {
   DMDALocalInfo  info;
   PetscInt       i, j;
   PetscReal      **aXu, dx, dy, x, y;
-  PetscCall(VecSet(Xu,PETSC_NINFINITY));
+  PetscCall(VecSet(Xl,PETSC_NINFINITY));
   PetscCall(SNESGetDM(snes,&da));
   PetscCall(DMDAGetLocalInfo(da,&info));
   dx = 1.0 / (PetscReal)(info.mx-1);