this is important enough to need a regression test

Author: bueler

c/ch7/solns/bratu2D.c

@@ -264,7 +264,6 @@ PetscErrorCode NonlinearGS(SNES snes, Vec u, Vec b, void *ctx) {
                     //            + hxhy * (2 u - au[j-1][i] - au[j+1][i])
                     //            - darea * lambda * e^u - bij
                     uu = au[j][i];
-                    phi0 = 0.0;
                     for (k = 0; k < maxits; k++) {
                         phi =   hyhx * (2.0 * uu - au[j][i-1] - au[j][i+1])
                               + hxhy * (2.0 * uu - au[j-1][i] - au[j+1][i])

c/ch7/solns/makefile

@@ -5,10 +5,18 @@ bratu2D: bratu2D.o ../../ch6/poissonfunctions.o
 	-${CLINKER} -o bratu2D bratu2D.o ../../ch6/poissonfunctions.o ${PETSC_LIB}
 	${RM} bratu2D.o ../../ch6/poissonfunctions.o
 
+# testing
+
+runbratu2D_1:
+	-@../../testit.sh bratu2D "-lb_exact -snes_converged_reason -lb_showcounts -snes_type fas -snes_fas_type full -fas_levels_snes_type ngs -fas_levels_snes_ngs_sweeps 2 -fas_levels_snes_max_it 1 -fas_coarse_snes_type ngs -fas_coarse_snes_ngs_sweeps 2 -fas_coarse_snes_max_it 4 -da_refine 2 -snes_monitor_short" 1 1
+
+test_bratu2D: runbratu2D_1
+
+test: test_bratu2D
+
 # etc
 
-.PHONY: distclean
+.PHONY: distclean test test_bratu2D runbratu2D_1
 
 distclean:
-	@rm -f bratu2D
-
+	@rm -f bratu2D *tmp

c/ch7/solns/output/bratu2D.test1

@@ -0,0 +1,7 @@
+  0 SNES Function norm 9.04754 
+  1 SNES Function norm 0.000449564 
+  2 SNES Function norm 1.87245e-06 
+  3 SNES Function norm 8.93257e-09 
+Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 3
+flops = 8.763e+04,  residual calls = 69,  NGS calls = 58
+done on 9 x 9 grid:   error |u-uexact|_inf = 3.169e-04