These are timing and accuracy results to solve the system
$$
\gathered
 -\nabla\cdot(a\nabla u)+b\cdot\nabla u+cu=f\quad\text{on }\Omega,
\\
a\frac{\partial u}{\partial n}+\gamma u=g\quad\text{on }\partial Omega.
\endgathered
$$
with $b=0$ with $\Omega=[0,1]^2$ and a regular
grid of size $129 \times 129$ with one set of diagonals and
linear finite elements.  Multigrid was used, with m1 presmoothing
steps, m2 postsmoothing steps, and p coarse-correction steps in each
MG-step, and r iterations of MG-step at each refinement level.
Richardson, Jacobi, or Gauss-Seidel smoothing was used.
At the coarsest level, only the smoothers were employed, i.e., no
special solver was used at the coarsest level.

The four columns of numbers are L2 error, H1 error, User ms, and GC ms.

The L2 and H1 errors were estimated by injecting the approximate
solution into a $257 \times 257$ grid and comparing it with the
interpolant of the true solution on the fine grid.  The problems
are roughly in increasing order of difficulty; in the last
problem, for example, the stiffness coefficient $a$ varies from
1.0 to 21.0 over the domain.

User ms is the amount of user+system time to solve the linear
system only; every once in a while there was a garbage collection
during the solution process, whose time is measured separately.

The tests were run on a box:

% cat /proc/cpu
model name	: Intel(R) Core(TM)2 Quad  CPU   Q8200  @ 2.33GHz

% uname -a
Linux heine 3.0.0-17-generic #30-Ubuntu SMP Thu Mar 8 20:45:39 UTC 2012 x86_64 x86_64 x86_64 GNU/Linux

% gsi -v
v4.6.4 20120222164157 x86_64-unknown-linux-gnu "./configure 'CC=gcc -march=native -fschedule-insns -frename-registers' '--enable-single-host' '--enable-shared' '--enable-multiple-versions'"

% gcc -v
Using built-in specs.
COLLECT_GCC=gcc
COLLECT_LTO_WRAPPER=/usr/lib/gcc/x86_64-linux-gnu/4.6.1/lto-wrapper
Target: x86_64-linux-gnu
Configured with: ../src/configure -v --with-pkgversion='Ubuntu/Linaro 4.6.1-9ubuntu3' --with-bugurl=file:///usr/share/doc/gcc-4.6/README.Bugs --enable-languages=c,c++,fortran,objc,obj-c++,go --prefix=/usr --program-suffix=-4.6 --enable-shared --enable-linker-build-id --with-system-zlib --libexecdir=/usr/lib --without-included-gettext --enable-threads=posix --with-gxx-include-dir=/usr/include/c++/4.6 --libdir=/usr/lib --enable-nls --with-sysroot=/ --enable-clocale=gnu --enable-libstdcxx-debug --enable-libstdcxx-time=yes --enable-plugin --enable-objc-gc --disable-werror --with-arch-32=i686 --with-tune=generic --enable-checking=release --build=x86_64-linux-gnu --host=x86_64-linux-gnu --target=x86_64-linux-gnu
Thread model: posix
gcc version 4.6.1 (Ubuntu/Linaro 4.6.1-9ubuntu3) 

I have not yet analyzed these results to see what the "best"
combination of m1, m2, r, and p is.

heine:~/programs/615project/2012/code> gsc++ -:m1000000
[ Meroon V3 Paques2001+1 $Revision: 1.2 $ ]
Gambit v4.6.4

> (load "all")
"/home/lucier/programs/615project/2012/code/all.o1"
> (load "final-tests")
"/home/lucier/programs/615project/2012/code/final-tests.o1"
> (run-multigrid-tests symmetric-tests)

a=1, c=1, u=1, Robin, gamma=x (error = 0 with exact linear algebra)
MG:   Richardson, r = 1, m1 = 1, m2 = 1, p = 1.         .0055458676     .0088262396     8.      0.
MG:   Richardson, r = 2, m1 = 1, m2 = 1, p = 1.         .0000494526     .0001162031     12.     0.
MG:   Richardson, r = 3, m1 = 1, m2 = 1, p = 1.         .0000005617     .0000018489     20.     0.
MG:   Richardson, r = 4, m1 = 1, m2 = 1, p = 1.         .0000000075     .0000000314     24.     0.
MG:   Richardson, r = 1, m1 = 2, m2 = 2, p = 1.         .0000576949     .0001674695     12.     0.
MG:   Richardson, r = 2, m1 = 2, m2 = 2, p = 1.         .0000000163     .0000000821     20.     0.
MG:   Richardson, r = 3, m1 = 2, m2 = 2, p = 1.         .0000000000     .0000000000     28.     0.
MG:   Richardson, r = 4, m1 = 2, m2 = 2, p = 1.         .0000000000     .0000000000     32.     0.
MG:   Richardson, r = 1, m1 = 3, m2 = 3, p = 1.         .0000016058     .0000077971     12.     0.
MG:   Richardson, r = 2, m1 = 3, m2 = 3, p = 1.         .0000000000     .0000000002     24.     0.
MG:   Richardson, r = 3, m1 = 3, m2 = 3, p = 1.         .0000000000     .0000000000     32.     0.
MG:   Richardson, r = 4, m1 = 3, m2 = 3, p = 1.         .0000000000     .0000000000     48.     0.
MG:   Richardson, r = 1, m1 = 4, m2 = 4, p = 1.         .0000001116     .0000006230     16.     0.
MG:   Richardson, r = 2, m1 = 4, m2 = 4, p = 1.         .0000000000     .0000000000     32.     0.
MG:   Richardson, r = 3, m1 = 4, m2 = 4, p = 1.         .0000000000     .0000000000     48.     0.
MG:   Richardson, r = 4, m1 = 4, m2 = 4, p = 1.         .0000000000     .0000000000     60.     0.
MG:   Richardson, r = 1, m1 = 1, m2 = 1, p = 2.         .0000000279     .0000012650     12.     0.
MG:   Richardson, r = 2, m1 = 1, m2 = 1, p = 2.         .0000000000     .0000000025     32.     0.
MG:   Richardson, r = 3, m1 = 1, m2 = 1, p = 2.         .0000000000     .0000000000     40.     0.
MG:   Richardson, r = 4, m1 = 1, m2 = 1, p = 2.         .0000000000     .0000000000     60.     0.
MG:   Richardson, r = 1, m1 = 2, m2 = 2, p = 2.         .0000000008     .0000000465     24.     0.
MG:   Richardson, r = 2, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     44.     0.
MG:   Richardson, r = 3, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     64.     0.
MG:   Richardson, r = 4, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     84.     0.
MG:   Richardson, r = 1, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000048     32.     0.
MG:   Richardson, r = 2, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     56.     0.
MG:   Richardson, r = 3, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     88.     0.
MG:   Richardson, r = 4, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     112.    0.
MG:   Richardson, r = 1, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000007     32.     0.
MG:   Richardson, r = 2, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     68.     0.
MG:   Richardson, r = 3, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     104.    0.
MG:   Richardson, r = 4, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     144.    0.
MG:       Jacobi, r = 1, m1 = 1, m2 = 1, p = 1.         .0000062710     .0000112146     4.      0.
MG:       Jacobi, r = 2, m1 = 1, m2 = 1, p = 1.         .0000000000     .0000000030     12.     0.
MG:       Jacobi, r = 3, m1 = 1, m2 = 1, p = 1.         .0000000000     .0000000000     16.     0.
MG:       Jacobi, r = 4, m1 = 1, m2 = 1, p = 1.         .0000000000     .0000000000     24.     0.
MG:       Jacobi, r = 1, m1 = 2, m2 = 2, p = 1.         .0000000019     .0000000103     8.      0.
MG:       Jacobi, r = 2, m1 = 2, m2 = 2, p = 1.         .0000000000     .0000000000     16.     0.
MG:       Jacobi, r = 3, m1 = 2, m2 = 2, p = 1.         .0000000000     .0000000000     28.     0.
MG:       Jacobi, r = 4, m1 = 2, m2 = 2, p = 1.         .0000000000     .0000000000     36.     0.
MG:       Jacobi, r = 1, m1 = 3, m2 = 3, p = 1.         .0000000000     .0000000001     12.     0.
MG:       Jacobi, r = 2, m1 = 3, m2 = 3, p = 1.         .0000000000     .0000000000     24.     0.
MG:       Jacobi, r = 3, m1 = 3, m2 = 3, p = 1.         .0000000000     .0000000000     36.     0.
MG:       Jacobi, r = 4, m1 = 3, m2 = 3, p = 1.         .0000000000     .0000000000     44.     0.
MG:       Jacobi, r = 1, m1 = 4, m2 = 4, p = 1.         .0000000000     .0000000000     16.     0.
MG:       Jacobi, r = 2, m1 = 4, m2 = 4, p = 1.         .0000000000     .0000000000     24.     0.
MG:       Jacobi, r = 3, m1 = 4, m2 = 4, p = 1.         .0000000000     .0000000000     44.     0.
MG:       Jacobi, r = 4, m1 = 4, m2 = 4, p = 1.         .0000000000     .0000000000     52.     0.
MG:       Jacobi, r = 1, m1 = 1, m2 = 1, p = 2.         .0000000023     .0000005152     16.     0.
MG:       Jacobi, r = 2, m1 = 1, m2 = 1, p = 2.         .0000000000     .0000000022     24.     0.
MG:       Jacobi, r = 3, m1 = 1, m2 = 1, p = 2.         .0000000000     .0000000000     40.     0.
MG:       Jacobi, r = 4, m1 = 1, m2 = 1, p = 2.         .0000000000     .0000000000     52.     0.
MG:       Jacobi, r = 1, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000005     20.     0.
MG:       Jacobi, r = 2, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     40.     0.
MG:       Jacobi, r = 3, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     60.     0.
MG:       Jacobi, r = 4, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     76.     0.
MG:       Jacobi, r = 1, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     28.     0.
MG:       Jacobi, r = 2, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     52.     0.
MG:       Jacobi, r = 3, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     76.     0.
MG:       Jacobi, r = 4, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     104.    0.
MG:       Jacobi, r = 1, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     32.     0.
MG:       Jacobi, r = 2, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     68.     0.
MG:       Jacobi, r = 3, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     96.     0.
MG:       Jacobi, r = 4, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     128.    0.
MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 1.         .0000000001     .0000000006     8.      0.
MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 1.         .0000000000     .0000000000     12.     0.
MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 1.         .0000000000     .0000000000     24.     0.
MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 1.         .0000000000     .0000000000     28.     0.
MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 1.         .0000000000     .0000000000     12.     0.
MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 1.         .0000000000     .0000000000     24.     0.
MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 1.         .0000000000     .0000000000     36.     0.
MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 1.         .0000000000     .0000000000     44.     0.
MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 1.         .0000000000     .0000000000     16.     0.
MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 1.         .0000000000     .0000000000     36.     0.
MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 1.         .0000000000     .0000000000     48.     0.
MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 1.         .0000000000     .0000000000     64.     0.
MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 1.         .0000000000     .0000000000     20.     0.
MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 1.         .0000000000     .0000000000     40.     0.
MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 1.         .0000000000     .0000000000     64.     0.
MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 1.         .0000000000     .0000000000     84.     0.
MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 2.         .0000000000     .0000000000     16.     0.
MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 2.         .0000000000     .0000000000     28.     0.
MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 2.         .0000000000     .0000000000     44.     0.
MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 2.         .0000000000     .0000000000     56.     0.
MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     24.     0.
MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     44.     0.
MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     64.     0.
MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     88.     0.
MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     28.     0.
MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     60.     0.
MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     84.     0.
MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     116.    0.
MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     36.     0.
MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     76.     0.
MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     108.    0.
MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     144.    0.

a=1, c=1, u=x, Neumann (error = 0 with exact linear algebra)
MG:   Richardson, r = 1, m1 = 1, m2 = 1, p = 1.         .0534908688     .0542735698     8.      0.
MG:   Richardson, r = 2, m1 = 1, m2 = 1, p = 1.         .0059216097     .0060273622     12.     0.
MG:   Richardson, r = 3, m1 = 1, m2 = 1, p = 1.         .0006581236     .0006706749     16.     0.
MG:   Richardson, r = 4, m1 = 1, m2 = 1, p = 1.         .0000732589     .0000746927     20.     0.
MG:   Richardson, r = 1, m1 = 2, m2 = 2, p = 1.         .0059941401     .0060858562     12.     0.
MG:   Richardson, r = 2, m1 = 2, m2 = 2, p = 1.         .0000748119     .0000761172     16.     0.
MG:   Richardson, r = 3, m1 = 2, m2 = 2, p = 1.         .0000009374     .0000009543     24.     0.
MG:   Richardson, r = 4, m1 = 2, m2 = 2, p = 1.         .0000000117     .0000000119     36.     0.
MG:   Richardson, r = 1, m1 = 3, m2 = 3, p = 1.         .0006804990     .0006910016     12.     0.
MG:   Richardson, r = 2, m1 = 3, m2 = 3, p = 1.         .0000009681     .0000009849     24.     0.
MG:   Richardson, r = 3, m1 = 3, m2 = 3, p = 1.         .0000000013     .0000000013     36.     0.
MG:   Richardson, r = 4, m1 = 3, m2 = 3, p = 1.         .0000000000     .0000000000     48.     0.
MG:   Richardson, r = 1, m1 = 4, m2 = 4, p = 1.         .0000783080     .0000795766     16.     0.
MG:   Richardson, r = 2, m1 = 4, m2 = 4, p = 1.         .0000000128     .0000000131     28.     0.
MG:   Richardson, r = 3, m1 = 4, m2 = 4, p = 1.         .0000000000     .0000000000     44.     0.
MG:   Richardson, r = 4, m1 = 4, m2 = 4, p = 1.         .0000000000     .0000000000     60.     0.
MG:   Richardson, r = 1, m1 = 1, m2 = 1, p = 2.         .0000000378     .0000016948     16.     0.
MG:   Richardson, r = 2, m1 = 1, m2 = 1, p = 2.         .0000000001     .0000000084     32.     0.
MG:   Richardson, r = 3, m1 = 1, m2 = 1, p = 2.         .0000000000     .0000000001     44.     0.
MG:   Richardson, r = 4, m1 = 1, m2 = 1, p = 2.         .0000000000     .0000000000     60.     0.
MG:   Richardson, r = 1, m1 = 2, m2 = 2, p = 2.         .0000000012     .0000000715     20.     0.
MG:   Richardson, r = 2, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     44.     0.
MG:   Richardson, r = 3, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     64.     0.
MG:   Richardson, r = 4, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     88.     0.
MG:   Richardson, r = 1, m1 = 3, m2 = 3, p = 2.         .0000000001     .0000000084     28.     0.
MG:   Richardson, r = 2, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     60.     0.
MG:   Richardson, r = 3, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     84.     0.
MG:   Richardson, r = 4, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     112.    0.
MG:   Richardson, r = 1, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000013     36.     0.
MG:   Richardson, r = 2, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     68.     0.
MG:   Richardson, r = 3, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     100.    0.
MG:   Richardson, r = 4, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     136.    0.
MG:       Jacobi, r = 1, m1 = 1, m2 = 1, p = 1.         .0031581475     .0031611772     4.      0.
MG:       Jacobi, r = 2, m1 = 1, m2 = 1, p = 1.         .0000199977     .0000200173     12.     0.
MG:       Jacobi, r = 3, m1 = 1, m2 = 1, p = 1.         .0000001266     .0000001267     16.     0.
MG:       Jacobi, r = 4, m1 = 1, m2 = 1, p = 1.         .0000000007     .0000000007     20.     0.
MG:       Jacobi, r = 1, m1 = 2, m2 = 2, p = 1.         .0000201252     .0000201329     8.      0.
MG:       Jacobi, r = 2, m1 = 2, m2 = 2, p = 1.         .0000000008     .0000000008     16.     0.
MG:       Jacobi, r = 3, m1 = 2, m2 = 2, p = 1.         .0000000000     .0000000000     24.     0.
MG:       Jacobi, r = 4, m1 = 2, m2 = 2, p = 1.         .0000000000     .0000000000     32.     0.
MG:       Jacobi, r = 1, m1 = 3, m2 = 3, p = 1.         .0000001291     .0000001295     12.     0.
MG:       Jacobi, r = 2, m1 = 3, m2 = 3, p = 1.         .0000000000     .0000000000     20.     0.
MG:       Jacobi, r = 3, m1 = 3, m2 = 3, p = 1.         .0000000000     .0000000000     36.     0.
MG:       Jacobi, r = 4, m1 = 3, m2 = 3, p = 1.         .0000000000     .0000000000     40.     0.
MG:       Jacobi, r = 1, m1 = 4, m2 = 4, p = 1.         .0000000008     .0000000008     16.     0.
MG:       Jacobi, r = 2, m1 = 4, m2 = 4, p = 1.         .0000000000     .0000000000     28.     0.
MG:       Jacobi, r = 3, m1 = 4, m2 = 4, p = 1.         .0000000000     .0000000000     40.     0.
MG:       Jacobi, r = 4, m1 = 4, m2 = 4, p = 1.         .0000000000     .0000000000     52.     0.
MG:       Jacobi, r = 1, m1 = 1, m2 = 1, p = 2.         .0000000032     .0000007273     16.     0.
MG:       Jacobi, r = 2, m1 = 1, m2 = 1, p = 2.         .0000000000     .0000000048     24.     0.
MG:       Jacobi, r = 3, m1 = 1, m2 = 1, p = 2.         .0000000000     .0000000000     40.     0.
MG:       Jacobi, r = 4, m1 = 1, m2 = 1, p = 2.         .0000000000     .0000000000     56.     0.
MG:       Jacobi, r = 1, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000084     20.     0.
MG:       Jacobi, r = 2, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     36.     0.
MG:       Jacobi, r = 3, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     60.     0.
MG:       Jacobi, r = 4, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     80.     0.
MG:       Jacobi, r = 1, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000001     24.     0.
MG:       Jacobi, r = 2, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     56.     0.
MG:       Jacobi, r = 3, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     80.     0.
MG:       Jacobi, r = 4, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     108.    0.
MG:       Jacobi, r = 1, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     36.     0.
MG:       Jacobi, r = 2, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     68.     0.
MG:       Jacobi, r = 3, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     100.    0.
MG:       Jacobi, r = 4, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     128.    100.
MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 1.         .0000005685     .0000005736     8.      0.
MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 1.         .0000000000     .0000000000     12.     0.
MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 1.         .0000000000     .0000000000     24.     0.
MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 1.         .0000000000     .0000000000     28.     0.
MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 1.         .0000000000     .0000000000     12.     0.
MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 1.         .0000000000     .0000000000     24.     0.
MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 1.         .0000000000     .0000000000     32.     0.
MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 1.         .0000000000     .0000000000     48.     0.
MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 1.         .0000000000     .0000000000     16.     0.
MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 1.         .0000000000     .0000000000     36.     0.
MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 1.         .0000000000     .0000000000     48.     0.
MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 1.         .0000000000     .0000000000     64.     0.
MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 1.         .0000000000     .0000000000     20.     0.
MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 1.         .0000000000     .0000000000     44.     0.
MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 1.         .0000000000     .0000000000     64.     0.
MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 1.         .0000000000     .0000000000     84.     0.
MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 2.         .0000000000     .0000000000     12.     0.
MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 2.         .0000000000     .0000000000     28.     0.
MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 2.         .0000000000     .0000000000     40.     0.
MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 2.         .0000000000     .0000000000     60.     0.
MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     20.     0.
MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     44.     0.
MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     64.     0.
MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     88.     0.
MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     28.     0.
MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     60.     0.
MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     88.     0.
MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     116.    0.
MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     36.     0.
MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     72.     0.
MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     108.    0.
MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     148.    104.

a=1+xy, c=1, u=x+y, Neumann (error = 0 with exact linear algebra)
MG:   Richardson, r = 1, m1 = 1, m2 = 1, p = 1.         .1682834371     .1701401832     8.      0.
MG:   Richardson, r = 2, m1 = 1, m2 = 1, p = 1.         .0291393933     .0295872349     8.      0.
MG:   Richardson, r = 3, m1 = 1, m2 = 1, p = 1.         .0050668851     .0051616101     16.     0.
MG:   Richardson, r = 4, m1 = 1, m2 = 1, p = 1.         .0008829806     .0009017995     20.     0.
MG:   Richardson, r = 1, m1 = 2, m2 = 2, p = 1.         .0293254461     .0297244992     8.      0.
MG:   Richardson, r = 2, m1 = 2, m2 = 2, p = 1.         .0008917406     .0009096386     20.     0.
MG:   Richardson, r = 3, m1 = 2, m2 = 2, p = 1.         .0000272760     .0000279526     28.     0.
MG:   Richardson, r = 4, m1 = 2, m2 = 2, p = 1.         .0000008368     .0000008607     32.     0.
MG:   Richardson, r = 1, m1 = 3, m2 = 3, p = 1.         .0051513316     .0052345193     12.     0.
MG:   Richardson, r = 2, m1 = 3, m2 = 3, p = 1.         .0000276823     .0000283919     24.     0.
MG:   Richardson, r = 3, m1 = 3, m2 = 3, p = 1.         .0000001499     .0000001549     32.     0.
MG:   Richardson, r = 4, m1 = 3, m2 = 3, p = 1.         .0000000007     .0000000008     48.     0.
MG:   Richardson, r = 1, m1 = 4, m2 = 4, p = 1.         .0009114261     .0009293160     16.     0.
MG:   Richardson, r = 2, m1 = 4, m2 = 4, p = 1.         .0000008720     .0000009034     28.     0.
MG:   Richardson, r = 3, m1 = 4, m2 = 4, p = 1.         .0000000008     .0000000008     44.     0.
MG:   Richardson, r = 4, m1 = 4, m2 = 4, p = 1.         .0000000000     .0000000000     56.     0.
MG:   Richardson, r = 1, m1 = 1, m2 = 1, p = 2.         .0000005046     .0000170471     16.     0.
MG:   Richardson, r = 2, m1 = 1, m2 = 1, p = 2.         .0000000083     .0000004412     24.     0.
MG:   Richardson, r = 3, m1 = 1, m2 = 1, p = 2.         .0000000002     .0000000189     44.     0.
MG:   Richardson, r = 4, m1 = 1, m2 = 1, p = 2.         .0000000000     .0000000010     56.     100.
MG:   Richardson, r = 1, m1 = 2, m2 = 2, p = 2.         .0000000206     .0000009825     20.     0.
MG:   Richardson, r = 2, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000025     40.     0.
MG:   Richardson, r = 3, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     60.     0.
MG:   Richardson, r = 4, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     84.     0.
MG:   Richardson, r = 1, m1 = 3, m2 = 3, p = 2.         .0000000026     .0000001491     28.     0.
MG:   Richardson, r = 2, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     52.     0.
MG:   Richardson, r = 3, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     84.     0.
MG:   Richardson, r = 4, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     108.    0.
MG:   Richardson, r = 1, m1 = 4, m2 = 4, p = 2.         .0000000005     .0000000352     32.     0.
MG:   Richardson, r = 2, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     68.     0.
MG:   Richardson, r = 3, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     108.    0.
MG:   Richardson, r = 4, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     140.    0.
MG:       Jacobi, r = 1, m1 = 1, m2 = 1, p = 1.         .0171879766     .0171983632     4.      0.
MG:       Jacobi, r = 2, m1 = 1, m2 = 1, p = 1.         .0002958620     .0002960412     12.     0.
MG:       Jacobi, r = 3, m1 = 1, m2 = 1, p = 1.         .0000050932     .0000050963     20.     0.
MG:       Jacobi, r = 4, m1 = 1, m2 = 1, p = 1.         .0000000876     .0000000876     24.     0.
MG:       Jacobi, r = 1, m1 = 2, m2 = 2, p = 1.         .0002997195     .0002999853     8.      0.
MG:       Jacobi, r = 2, m1 = 2, m2 = 2, p = 1.         .0000000904     .0000000905     16.     0.
MG:       Jacobi, r = 3, m1 = 2, m2 = 2, p = 1.         .0000000000     .0000000000     28.     0.
MG:       Jacobi, r = 4, m1 = 2, m2 = 2, p = 1.         .0000000000     .0000000000     36.     0.
MG:       Jacobi, r = 1, m1 = 3, m2 = 3, p = 1.         .0000052428     .0000052517     12.     0.
MG:       Jacobi, r = 2, m1 = 3, m2 = 3, p = 1.         .0000000000     .0000000000     24.     0.
MG:       Jacobi, r = 3, m1 = 3, m2 = 3, p = 1.         .0000000000     .0000000000     36.     0.
MG:       Jacobi, r = 4, m1 = 3, m2 = 3, p = 1.         .0000000000     .0000000000     48.     0.
MG:       Jacobi, r = 1, m1 = 4, m2 = 4, p = 1.         .0000000936     .0000000939     16.     0.
MG:       Jacobi, r = 2, m1 = 4, m2 = 4, p = 1.         .0000000000     .0000000000     32.     0.
MG:       Jacobi, r = 3, m1 = 4, m2 = 4, p = 1.         .0000000000     .0000000000     44.     0.
MG:       Jacobi, r = 4, m1 = 4, m2 = 4, p = 1.         .0000000000     .0000000000     60.     0.
MG:       Jacobi, r = 1, m1 = 1, m2 = 1, p = 2.         .0000000058     .0000012452     16.     0.
MG:       Jacobi, r = 2, m1 = 1, m2 = 1, p = 2.         .0000000000     .0000000121     28.     0.
MG:       Jacobi, r = 3, m1 = 1, m2 = 1, p = 2.         .0000000000     .0000000001     44.     0.
MG:       Jacobi, r = 4, m1 = 1, m2 = 1, p = 2.         .0000000000     .0000000000     56.     0.
MG:       Jacobi, r = 1, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000156     20.     0.
MG:       Jacobi, r = 2, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     40.     0.
MG:       Jacobi, r = 3, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     64.     96.
MG:       Jacobi, r = 4, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     80.     0.
MG:       Jacobi, r = 1, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000002     24.     0.
MG:       Jacobi, r = 2, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     52.     0.
MG:       Jacobi, r = 3, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     76.     0.
MG:       Jacobi, r = 4, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     104.    0.
MG:       Jacobi, r = 1, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     36.     0.
MG:       Jacobi, r = 2, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     64.     0.
MG:       Jacobi, r = 3, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     96.     0.
MG:       Jacobi, r = 4, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     128.    0.
MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 1.         .0000136452     .0000136939     8.      0.
MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 1.         .0000000001     .0000000001     16.     0.
MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 1.         .0000000000     .0000000000     20.     0.
MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 1.         .0000000000     .0000000000     28.     0.
MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 1.         .0000000005     .0000000005     12.     0.
MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 1.         .0000000000     .0000000000     24.     0.
MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 1.         .0000000000     .0000000000     36.     0.
MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 1.         .0000000000     .0000000000     44.     0.
MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 1.         .0000000000     .0000000000     20.     0.
MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 1.         .0000000000     .0000000000     32.     0.
MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 1.         .0000000000     .0000000000     48.     0.
MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 1.         .0000000000     .0000000000     68.     0.
MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 1.         .0000000000     .0000000000     20.     0.
MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 1.         .0000000000     .0000000000     44.     0.
MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 1.         .0000000000     .0000000000     60.     0.
MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 1.         .0000000000     .0000000000     84.     0.
MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 2.         .0000000000     .0000000000     12.     0.
MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 2.         .0000000000     .0000000000     32.     0.
MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 2.         .0000000000     .0000000000     44.     0.
MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 2.         .0000000000     .0000000000     56.     0.
MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     20.     0.
MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     44.     0.
MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     68.     0.
MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 2.         .0000000000     .0000000000     84.     0.
MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     32.     0.
MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     60.     0.
MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     88.     0.
MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 2.         .0000000000     .0000000000     116.    0.
MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     36.     0.
MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     72.     0.
MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     108.    0.
MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 2.         .0000000000     .0000000000     144.    0.

a=1, c=1, u=x^2-y^2, Neumann
MG:   Richardson, r = 1, m1 = 1, m2 = 1, p = 1.         .0004378261     .0072326137     8.      0.
MG:   Richardson, r = 2, m1 = 1, m2 = 1, p = 1.         .0000078173     .0055800823     12.     0.
MG:   Richardson, r = 3, m1 = 1, m2 = 1, p = 1.         .0000063913     .0055292304     16.     0.
MG:   Richardson, r = 4, m1 = 1, m2 = 1, p = 1.         .0000062459     .0055248819     20.     0.
MG:   Richardson, r = 1, m1 = 2, m2 = 2, p = 1.         .0000082085     .0055809129     12.     0.
MG:   Richardson, r = 2, m1 = 2, m2 = 2, p = 1.         .0000062513     .0055248850     16.     0.
MG:   Richardson, r = 3, m1 = 2, m2 = 2, p = 1.         .0000062299     .0055242937     24.     0.
MG:   Richardson, r = 4, m1 = 2, m2 = 2, p = 1.         .0000062293     .0055242760     32.     0.
MG:   Richardson, r = 1, m1 = 3, m2 = 3, p = 1.         .0000065191     .0055293172     12.     0.
MG:   Richardson, r = 2, m1 = 3, m2 = 3, p = 1.         .0000062304     .0055242940     24.     0.
MG:   Richardson, r = 3, m1 = 3, m2 = 3, p = 1.         .0000062293     .0055242754     36.     0.
MG:   Richardson, r = 4, m1 = 3, m2 = 3, p = 1.         .0000062293     .0055242752     44.     0.
MG:   Richardson, r = 1, m1 = 4, m2 = 4, p = 1.         .0000062747     .0055248968     16.     0.
MG:   Richardson, r = 2, m1 = 4, m2 = 4, p = 1.         .0000062294     .0055242761     28.     0.
MG:   Richardson, r = 3, m1 = 4, m2 = 4, p = 1.         .0000062293     .0055242752     44.     0.
MG:   Richardson, r = 4, m1 = 4, m2 = 4, p = 1.         .0000062293     .0055242752     60.     0.
MG:   Richardson, r = 1, m1 = 1, m2 = 1, p = 2.         .0000094709     .0062323143     16.     0.
MG:   Richardson, r = 2, m1 = 1, m2 = 1, p = 2.         .0000064934     .0055777389     32.     0.
MG:   Richardson, r = 3, m1 = 1, m2 = 1, p = 2.         .0000062558     .0055290945     44.     0.
MG:   Richardson, r = 4, m1 = 1, m2 = 1, p = 2.         .0000062330     .0055248680     60.     0.
MG:   Richardson, r = 1, m1 = 2, m2 = 2, p = 2.         .0000064969     .0055782897     20.     0.
MG:   Richardson, r = 2, m1 = 2, m2 = 2, p = 2.         .0000062335     .0055248691     40.     0.
MG:   Richardson, r = 3, m1 = 2, m2 = 2, p = 2.         .0000062295     .0055242934     64.     0.
MG:   Richardson, r = 4, m1 = 2, m2 = 2, p = 2.         .0000062293     .0055242760     84.     0.
MG:   Richardson, r = 1, m1 = 3, m2 = 3, p = 2.         .0000062591     .0055291851     28.     0.
MG:   Richardson, r = 2, m1 = 3, m2 = 3, p = 2.         .0000062296     .0055242934     56.     0.
MG:   Richardson, r = 3, m1 = 3, m2 = 3, p = 2.         .0000062293     .0055242754     88.     0.
MG:   Richardson, r = 4, m1 = 3, m2 = 3, p = 2.         .0000062293     .0055242752     112.    0.
MG:   Richardson, r = 1, m1 = 4, m2 = 4, p = 2.         .0000062357     .0055248829     36.     0.
MG:   Richardson, r = 2, m1 = 4, m2 = 4, p = 2.         .0000062293     .0055242760     68.     0.
MG:   Richardson, r = 3, m1 = 4, m2 = 4, p = 2.         .0000062293     .0055242752     108.    0.
MG:   Richardson, r = 4, m1 = 4, m2 = 4, p = 2.         .0000062293     .0055242752     140.    0.
MG:       Jacobi, r = 1, m1 = 1, m2 = 1, p = 1.         .0000062293     .0055242752     4.      0.
MG:       Jacobi, r = 2, m1 = 1, m2 = 1, p = 1.         .0000062293     .0055242752     12.     0.
MG:       Jacobi, r = 3, m1 = 1, m2 = 1, p = 1.         .0000062293     .0055242752     20.     0.
MG:       Jacobi, r = 4, m1 = 1, m2 = 1, p = 1.         .0000062293     .0055242752     24.     0.
MG:       Jacobi, r = 1, m1 = 2, m2 = 2, p = 1.         .0000062293     .0055242752     8.      0.
MG:       Jacobi, r = 2, m1 = 2, m2 = 2, p = 1.         .0000062293     .0055242752     16.     0.
MG:       Jacobi, r = 3, m1 = 2, m2 = 2, p = 1.         .0000062293     .0055242752     28.     0.
MG:       Jacobi, r = 4, m1 = 2, m2 = 2, p = 1.         .0000062293     .0055242752     36.     0.
MG:       Jacobi, r = 1, m1 = 3, m2 = 3, p = 1.         .0000062293     .0055242752     12.     0.
MG:       Jacobi, r = 2, m1 = 3, m2 = 3, p = 1.         .0000062293     .0055242752     24.     0.
MG:       Jacobi, r = 3, m1 = 3, m2 = 3, p = 1.         .0000062293     .0055242752     36.     0.
MG:       Jacobi, r = 4, m1 = 3, m2 = 3, p = 1.         .0000062293     .0055242752     40.     0.
MG:       Jacobi, r = 1, m1 = 4, m2 = 4, p = 1.         .0000062293     .0055242752     16.     0.
MG:       Jacobi, r = 2, m1 = 4, m2 = 4, p = 1.         .0000062293     .0055242752     28.     0.
MG:       Jacobi, r = 3, m1 = 4, m2 = 4, p = 1.         .0000062293     .0055242752     40.     0.
MG:       Jacobi, r = 4, m1 = 4, m2 = 4, p = 1.         .0000062293     .0055242752     56.     0.
MG:       Jacobi, r = 1, m1 = 1, m2 = 1, p = 2.         .0000062293     .0055242752     16.     0.
MG:       Jacobi, r = 2, m1 = 1, m2 = 1, p = 2.         .0000062293     .0055242752     28.     0.
MG:       Jacobi, r = 3, m1 = 1, m2 = 1, p = 2.         .0000062293     .0055242752     40.     0.
MG:       Jacobi, r = 4, m1 = 1, m2 = 1, p = 2.         .0000062293     .0055242752     52.     0.
MG:       Jacobi, r = 1, m1 = 2, m2 = 2, p = 2.         .0000062293     .0055242752     20.     0.
MG:       Jacobi, r = 2, m1 = 2, m2 = 2, p = 2.         .0000062293     .0055242752     40.     0.
MG:       Jacobi, r = 3, m1 = 2, m2 = 2, p = 2.         .0000062293     .0055242752     56.     0.
MG:       Jacobi, r = 4, m1 = 2, m2 = 2, p = 2.         .0000062293     .0055242752     80.     0.
MG:       Jacobi, r = 1, m1 = 3, m2 = 3, p = 2.         .0000062293     .0055242752     28.     0.
MG:       Jacobi, r = 2, m1 = 3, m2 = 3, p = 2.         .0000062293     .0055242752     52.     0.
MG:       Jacobi, r = 3, m1 = 3, m2 = 3, p = 2.         .0000062293     .0055242752     76.     0.
MG:       Jacobi, r = 4, m1 = 3, m2 = 3, p = 2.         .0000062293     .0055242752     104.    0.
MG:       Jacobi, r = 1, m1 = 4, m2 = 4, p = 2.         .0000062293     .0055242752     36.     0.
MG:       Jacobi, r = 2, m1 = 4, m2 = 4, p = 2.         .0000062293     .0055242752     64.     0.
MG:       Jacobi, r = 3, m1 = 4, m2 = 4, p = 2.         .0000062293     .0055242752     96.     0.
MG:       Jacobi, r = 4, m1 = 4, m2 = 4, p = 2.         .0000062293     .0055242752     128.    0.
MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 1.         .0000062601     .0055252170     8.      0.
MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 1.         .0000062295     .0055242783     16.     0.
MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 1.         .0000062293     .0055242752     20.     0.
MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 1.         .0000062293     .0055242752     32.     0.
MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 1.         .0000062304     .0055242879     12.     0.
MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 1.         .0000062293     .0055242752     24.     0.
MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 1.         .0000062293     .0055242752     32.     0.
MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 1.         .0000062293     .0055242752     48.     0.
MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 1.         .0000062295     .0055242764     16.     0.
MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 1.         .0000062293     .0055242752     36.     0.
MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 1.         .0000062293     .0055242752     48.     0.
MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 1.         .0000062293     .0055242752     68.     0.
MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 1.         .0000062294     .0055242754     20.     0.
MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 1.         .0000062293     .0055242752     44.     0.
MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 1.         .0000062293     .0055242752     60.     0.
MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 1.         .0000062293     .0055242752     88.     96.
MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 2.         .0000062518     .0055251896     16.     0.
MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 2.         .0000062294     .0055242780     28.     0.
MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 2.         .0000062293     .0055242752     44.     0.
MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 2.         .0000062293     .0055242752     56.     0.
MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 2.         .0000062302     .0055242874     24.     0.
MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 2.         .0000062293     .0055242752     44.     0.
MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 2.         .0000062293     .0055242752     64.     0.
MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 2.         .0000062293     .0055242752     88.     0.
MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 2.         .0000062295     .0055242763     28.     0.
MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 2.         .0000062293     .0055242752     60.     0.
MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 2.         .0000062293     .0055242752     88.     0.
MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 2.         .0000062293     .0055242752     116.    0.
MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 2.         .0000062294     .0055242754     36.     0.
MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 2.         .0000062293     .0055242752     76.     0.
MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 2.         .0000062293     .0055242752     108.    0.
MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 2.         .0000062293     .0055242752     148.    0.

a=1, c=1, u=\cos(\pi x), Neumann
MG:   Richardson, r = 1, m1 = 1, m2 = 1, p = 1.         .0054358964     .0390809425     4.      0.
MG:   Richardson, r = 2, m1 = 1, m2 = 1, p = 1.         .0004739971     .0142387214     12.     0.
MG:   Richardson, r = 3, m1 = 1, m2 = 1, p = 1.         .0000986245     .0136616518     16.     0.
MG:   Richardson, r = 4, m1 = 1, m2 = 1, p = 1.         .0000488106     .0136340656     24.     0.
MG:   Richardson, r = 1, m1 = 2, m2 = 2, p = 1.         .0006288170     .0145116072     8.      0.
MG:   Richardson, r = 2, m1 = 2, m2 = 2, p = 1.         .0000571341     .0136354312     16.     0.
MG:   Richardson, r = 3, m1 = 2, m2 = 2, p = 1.         .0000353715     .0136302541     28.     0.
MG:   Richardson, r = 4, m1 = 2, m2 = 2, p = 1.         .0000320064     .0136299067     32.     0.
MG:   Richardson, r = 1, m1 = 3, m2 = 3, p = 1.         .0001794943     .0136995623     12.     0.
MG:   Richardson, r = 2, m1 = 3, m2 = 3, p = 1.         .0000380616     .0136304427     24.     0.
MG:   Richardson, r = 3, m1 = 3, m2 = 3, p = 1.         .0000319522     .0136298885     32.     0.
MG:   Richardson, r = 4, m1 = 3, m2 = 3, p = 1.         .0000309306     .0136298285     48.     0.
MG:   Richardson, r = 1, m1 = 4, m2 = 4, p = 1.         .0000927813     .0136430435     12.     0.
MG:   Richardson, r = 2, m1 = 4, m2 = 4, p = 1.         .0000341310     .0136300105     28.     0.
MG:   Richardson, r = 3, m1 = 4, m2 = 4, p = 1.         .0000311711     .0136298380     48.     0.
MG:   Richardson, r = 4, m1 = 4, m2 = 4, p = 1.         .0000307036     .0136298158     64.     0.
MG:   Richardson, r = 1, m1 = 1, m2 = 1, p = 2.         .0000962161     .0154572028     16.     0.
MG:   Richardson, r = 2, m1 = 1, m2 = 1, p = 2.         .0000434194     .0137713190     32.     0.
MG:   Richardson, r = 3, m1 = 1, m2 = 1, p = 2.         .0000337677     .0136432949     44.     0.
MG:   Richardson, r = 4, m1 = 1, m2 = 1, p = 2.         .0000314236     .0136316147     60.     0.
MG:   Richardson, r = 1, m1 = 2, m2 = 2, p = 2.         .0000559057     .0137752384     24.     0.
MG:   Richardson, r = 2, m1 = 2, m2 = 2, p = 2.         .0000322004     .0136316833     44.     0.
MG:   Richardson, r = 3, m1 = 2, m2 = 2, p = 2.         .0000307052     .0136298958     64.     0.
MG:   Richardson, r = 4, m1 = 2, m2 = 2, p = 2.         .0000305904     .0136298200     88.     0.
MG:   Richardson, r = 1, m1 = 3, m2 = 3, p = 2.         .0000428556     .0136444298     32.     0.
MG:   Richardson, r = 2, m1 = 3, m2 = 3, p = 2.         .0000308208     .0136299127     56.     0.
MG:   Richardson, r = 3, m1 = 3, m2 = 3, p = 2.         .0000305853     .0136298169     88.     0.
MG:   Richardson, r = 4, m1 = 3, m2 = 3, p = 2.         .0000305779     .0136298112     112.    0.
MG:   Richardson, r = 1, m1 = 4, m2 = 4, p = 2.         .0000367499     .0136320551     36.     0.
MG:   Richardson, r = 2, m1 = 4, m2 = 4, p = 2.         .0000306219     .0136298316     68.     96.
MG:   Richardson, r = 3, m1 = 4, m2 = 4, p = 2.         .0000305782     .0136298121     104.    0.
MG:   Richardson, r = 4, m1 = 4, m2 = 4, p = 2.         .0000305775     .0136298103     136.    0.
MG:       Jacobi, r = 1, m1 = 1, m2 = 1, p = 1.         .0000986649     .0136465291     8.      0.
MG:       Jacobi, r = 2, m1 = 1, m2 = 1, p = 1.         .0000343911     .0136309715     8.      0.
MG:       Jacobi, r = 3, m1 = 1, m2 = 1, p = 1.         .0000312712     .0136305333     16.     0.
MG:       Jacobi, r = 4, m1 = 1, m2 = 1, p = 1.         .0000307536     .0136303998     20.     0.
MG:       Jacobi, r = 1, m1 = 2, m2 = 2, p = 1.         .0000509780     .0136307764     8.      0.
MG:       Jacobi, r = 2, m1 = 2, m2 = 2, p = 1.         .0000319269     .0136298910     16.     0.
MG:       Jacobi, r = 3, m1 = 2, m2 = 2, p = 1.         .0000307404     .0136298541     24.     0.
MG:       Jacobi, r = 4, m1 = 2, m2 = 2, p = 1.         .0000305968     .0136298431     32.     0.
MG:       Jacobi, r = 1, m1 = 3, m2 = 3, p = 1.         .0000430848     .0136308724     12.     0.
MG:       Jacobi, r = 2, m1 = 3, m2 = 3, p = 1.         .0000313215     .0136302129     24.     0.
MG:       Jacobi, r = 3, m1 = 3, m2 = 3, p = 1.         .0000306671     .0136301274     32.     0.
MG:       Jacobi, r = 4, m1 = 3, m2 = 3, p = 1.         .0000306063     .0136300818     44.     0.
MG:       Jacobi, r = 1, m1 = 4, m2 = 4, p = 1.         .0000390511     .0136301256     12.     0.
MG:       Jacobi, r = 2, m1 = 4, m2 = 4, p = 1.         .0000310072     .0136298752     28.     0.
MG:       Jacobi, r = 3, m1 = 4, m2 = 4, p = 1.         .0000306056     .0136298523     44.     0.
MG:       Jacobi, r = 4, m1 = 4, m2 = 4, p = 1.         .0000305757     .0136298427     52.     0.
MG:       Jacobi, r = 1, m1 = 1, m2 = 1, p = 2.         .0000308478     .0136338458     12.     0.
MG:       Jacobi, r = 2, m1 = 1, m2 = 1, p = 2.         .0000306125     .0136307996     28.     0.
MG:       Jacobi, r = 3, m1 = 1, m2 = 1, p = 2.         .0000306054     .0136305142     40.     0.
MG:       Jacobi, r = 4, m1 = 1, m2 = 1, p = 2.         .0000306042     .0136303989     52.     0.
MG:       Jacobi, r = 1, m1 = 2, m2 = 2, p = 2.         .0000306818     .0136299558     20.     0.
MG:       Jacobi, r = 2, m1 = 2, m2 = 2, p = 2.         .0000305757     .0136298669     40.     0.
MG:       Jacobi, r = 3, m1 = 2, m2 = 2, p = 2.         .0000305738     .0136298505     60.     0.
MG:       Jacobi, r = 4, m1 = 2, m2 = 2, p = 2.         .0000305735     .0136298419     80.     0.
MG:       Jacobi, r = 1, m1 = 3, m2 = 3, p = 2.         .0000307121     .0136304496     28.     0.
MG:       Jacobi, r = 2, m1 = 3, m2 = 3, p = 2.         .0000306019     .0136302000     52.     0.
MG:       Jacobi, r = 3, m1 = 3, m2 = 3, p = 2.         .0000306001     .0136301273     80.     0.
MG:       Jacobi, r = 4, m1 = 3, m2 = 3, p = 2.         .0000305999     .0136300827     108.    0.
MG:       Jacobi, r = 1, m1 = 4, m2 = 4, p = 2.         .0000306654     .0136299371     36.     0.
MG:       Jacobi, r = 2, m1 = 4, m2 = 4, p = 2.         .0000305746     .0136298687     68.     0.
MG:       Jacobi, r = 3, m1 = 4, m2 = 4, p = 2.         .0000305733     .0136298515     104.    0.
MG:       Jacobi, r = 4, m1 = 4, m2 = 4, p = 2.         .0000305730     .0136298423     136.    96.
MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 1.         .0000512079     .0136322464     8.      0.
MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 1.         .0000318329     .0136298421     12.     0.
MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 1.         .0000307050     .0136298126     24.     0.
MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 1.         .0000305914     .0136298102     28.     0.
MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 1.         .0000370979     .0136299923     12.     0.
MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 1.         .0000308494     .0136298150     24.     0.
MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 1.         .0000305923     .0136298102     36.     0.
MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 1.         .0000305783     .0136298099     48.     0.
MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 1.         .0000343745     .0136298904     16.     0.
MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 1.         .0000306928     .0136298120     32.     0.
MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 1.         .0000305817     .0136298100     52.     0.
MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 1.         .0000305776     .0136298099     64.     0.
MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 1.         .0000333035     .0136298603     24.     0.
MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 1.         .0000306429     .0136298110     40.     0.
MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 1.         .0000305793     .0136298100     64.     0.
MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 1.         .0000305775     .0136298099     84.     0.
MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 2.         .0000332675     .0136314300     16.     0.
MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 2.         .0000306015     .0136298197     28.     0.
MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 2.         .0000305780     .0136298102     40.     0.
MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 2.         .0000305774     .0136298099     60.     0.
MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 2.         .0000308897     .0136298604     20.     0.
MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 2.         .0000305786     .0136298107     48.     0.
MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 2.         .0000305774     .0136298099     64.     0.
MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 2.         .0000305774     .0136298099     88.     0.
MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 2.         .0000306843     .0136298308     28.     0.
MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 2.         .0000305779     .0136298103     60.     0.
MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 2.         .0000305774     .0136298099     88.     0.
MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 2.         .0000305774     .0136298099     116.    0.
MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 2.         .0000306505     .0136298241     40.     0.
MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 2.         .0000305778     .0136298101     72.     0.
MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 2.         .0000305774     .0136298099     112.    0.
MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 2.         .0000305774     .0136298099     144.    108.

a=1, c=1, u=\sin(\pi x)+\sin(\pi y), Neumann
MG:   Richardson, r = 1, m1 = 1, m2 = 1, p = 1.         .0695028503     .1050213908     4.      0.
MG:   Richardson, r = 2, m1 = 1, m2 = 1, p = 1.         .0042315813     .0204351142     12.     0.
MG:   Richardson, r = 3, m1 = 1, m2 = 1, p = 1.         .0002066308     .0193048663     16.     0.
MG:   Richardson, r = 4, m1 = 1, m2 = 1, p = 1.         .0001751082     .0192798781     24.     0.
MG:   Richardson, r = 1, m1 = 2, m2 = 2, p = 1.         .0037306023     .0215176674     8.      0.
MG:   Richardson, r = 2, m1 = 2, m2 = 2, p = 1.         .0002335703     .0192824401     16.     0.
MG:   Richardson, r = 3, m1 = 2, m2 = 2, p = 1.         .0000699977     .0192763915     28.     0.
MG:   Richardson, r = 4, m1 = 2, m2 = 2, p = 1.         .0000511290     .0192761127     32.     0.
MG:   Richardson, r = 1, m1 = 3, m2 = 3, p = 1.         .0006709051     .0194162422     12.     0.
MG:   Richardson, r = 2, m1 = 3, m2 = 3, p = 1.         .0000785120     .0192766178     24.     0.
MG:   Richardson, r = 3, m1 = 3, m2 = 3, p = 1.         .0000501485     .0192761084     36.     0.
MG:   Richardson, r = 4, m1 = 3, m2 = 3, p = 1.         .0000457822     .0192760629     44.     0.
MG:   Richardson, r = 1, m1 = 4, m2 = 4, p = 1.         .0003520284     .0192986948     16.     0.
MG:   Richardson, r = 2, m1 = 4, m2 = 4, p = 1.         .0000565372     .0192762118     28.     0.
MG:   Richardson, r = 3, m1 = 4, m2 = 4, p = 1.         .0000464192     .0192760707     44.     0.
MG:   Richardson, r = 4, m1 = 4, m2 = 4, p = 1.         .0000445657     .0192760536     60.     0.
MG:   Richardson, r = 1, m1 = 1, m2 = 1, p = 2.         .0001085239     .0216786406     12.     0.
MG:   Richardson, r = 2, m1 = 1, m2 = 1, p = 2.         .0000561681     .0194529958     28.     0.
MG:   Richardson, r = 3, m1 = 1, m2 = 1, p = 2.         .0000468198     .0192913672     40.     0.
MG:   Richardson, r = 4, m1 = 1, m2 = 1, p = 2.         .0000446456     .0192778485     56.     0.
MG:   Richardson, r = 1, m1 = 2, m2 = 2, p = 2.         .0000689349     .0194545713     24.     92.
MG:   Richardson, r = 2, m1 = 2, m2 = 2, p = 2.         .0000453642     .0192778626     40.     0.
MG:   Richardson, r = 3, m1 = 2, m2 = 2, p = 2.         .0000440173     .0192760999     64.     0.
MG:   Richardson, r = 4, m1 = 2, m2 = 2, p = 2.         .0000439286     .0192760507     80.     0.
MG:   Richardson, r = 1, m1 = 3, m2 = 3, p = 2.         .0000557968     .0192917440     28.     0.
MG:   Richardson, r = 2, m1 = 3, m2 = 3, p = 2.         .0000441123     .0192761012     56.     0.
MG:   Richardson, r = 3, m1 = 3, m2 = 3, p = 2.         .0000439253     .0192760491     80.     0.
MG:   Richardson, r = 4, m1 = 3, m2 = 3, p = 2.         .0000439215     .0192760487     108.    0.
MG:   Richardson, r = 1, m1 = 4, m2 = 4, p = 2.         .0000497369     .0192779604     36.     0.
MG:   Richardson, r = 2, m1 = 4, m2 = 4, p = 2.         .0000439496     .0192760510     68.     0.
MG:   Richardson, r = 3, m1 = 4, m2 = 4, p = 2.         .0000439215     .0192760487     100.    0.
MG:   Richardson, r = 4, m1 = 4, m2 = 4, p = 2.         .0000439213     .0192760487     136.    0.
MG:       Jacobi, r = 1, m1 = 1, m2 = 1, p = 1.         .0088256113     .0212346981     8.      0.
MG:       Jacobi, r = 2, m1 = 1, m2 = 1, p = 1.         .0000897298     .0192763387     8.      0.
MG:       Jacobi, r = 3, m1 = 1, m2 = 1, p = 1.         .0000435390     .0192760678     16.     0.
MG:       Jacobi, r = 4, m1 = 1, m2 = 1, p = 1.         .0000436320     .0192760567     24.     0.
MG:       Jacobi, r = 1, m1 = 2, m2 = 2, p = 1.         .0000910865     .0192772849     8.      0.
MG:       Jacobi, r = 2, m1 = 2, m2 = 2, p = 1.         .0000456536     .0192760733     16.     0.
MG:       Jacobi, r = 3, m1 = 2, m2 = 2, p = 1.         .0000441263     .0192760516     24.     0.
MG:       Jacobi, r = 4, m1 = 2, m2 = 2, p = 1.         .0000439415     .0192760492     36.     0.
MG:       Jacobi, r = 1, m1 = 3, m2 = 3, p = 1.         .0000613843     .0192765527     12.     0.
MG:       Jacobi, r = 2, m1 = 3, m2 = 3, p = 1.         .0000449599     .0192760669     24.     0.
MG:       Jacobi, r = 3, m1 = 3, m2 = 3, p = 1.         .0000440459     .0192760553     36.     0.
MG:       Jacobi, r = 4, m1 = 3, m2 = 3, p = 1.         .0000439545     .0192760542     48.     0.
MG:       Jacobi, r = 1, m1 = 4, m2 = 4, p = 1.         .0000552263     .0192762905     16.     0.
MG:       Jacobi, r = 2, m1 = 4, m2 = 4, p = 1.         .0000445416     .0192760555     28.     0.
MG:       Jacobi, r = 3, m1 = 4, m2 = 4, p = 1.         .0000439706     .0192760493     44.     0.
MG:       Jacobi, r = 4, m1 = 4, m2 = 4, p = 1.         .0000439228     .0192760488     60.     0.
MG:       Jacobi, r = 1, m1 = 1, m2 = 1, p = 2.         .0000440691     .0192768504     16.     0.
MG:       Jacobi, r = 2, m1 = 1, m2 = 1, p = 2.         .0000439442     .0192760567     28.     0.
MG:       Jacobi, r = 3, m1 = 1, m2 = 1, p = 2.         .0000439437     .0192760543     40.     0.
MG:       Jacobi, r = 4, m1 = 1, m2 = 1, p = 2.         .0000439436     .0192760542     60.     0.
MG:       Jacobi, r = 1, m1 = 2, m2 = 2, p = 2.         .0000440118     .0192760538     20.     0.
MG:       Jacobi, r = 2, m1 = 2, m2 = 2, p = 2.         .0000439184     .0192760489     44.     0.
MG:       Jacobi, r = 3, m1 = 2, m2 = 2, p = 2.         .0000439181     .0192760488     64.     0.
MG:       Jacobi, r = 4, m1 = 2, m2 = 2, p = 2.         .0000439181     .0192760488     84.     96.
MG:       Jacobi, r = 1, m1 = 3, m2 = 3, p = 2.         .0000440362     .0192760582     24.     0.
MG:       Jacobi, r = 2, m1 = 3, m2 = 3, p = 2.         .0000439438     .0192760542     52.     0.
MG:       Jacobi, r = 3, m1 = 3, m2 = 3, p = 2.         .0000439436     .0192760541     80.     0.
MG:       Jacobi, r = 4, m1 = 3, m2 = 3, p = 2.         .0000439436     .0192760541     104.    0.
MG:       Jacobi, r = 1, m1 = 4, m2 = 4, p = 2.         .0000440088     .0192760518     32.     0.
MG:       Jacobi, r = 2, m1 = 4, m2 = 4, p = 2.         .0000439183     .0192760488     68.     0.
MG:       Jacobi, r = 3, m1 = 4, m2 = 4, p = 2.         .0000439181     .0192760488     96.     0.
MG:       Jacobi, r = 4, m1 = 4, m2 = 4, p = 2.         .0000439181     .0192760488     132.    0.
MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 1.         .0000756504     .0192780611     8.      0.
MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 1.         .0000460658     .0192760782     12.     0.
MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 1.         .0000441990     .0192760514     24.     0.
MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 1.         .0000439606     .0192760491     28.     0.
MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 1.         .0000539447     .0192762387     12.     0.
MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 1.         .0000444445     .0192760537     24.     0.
MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 1.         .0000439587     .0192760491     36.     0.
MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 1.         .0000439240     .0192760487     44.     0.
MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 1.         .0000497231     .0192761340     20.     0.
MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 1.         .0000441486     .0192760509     32.     0.
MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 1.         .0000439322     .0192760488     52.     0.
MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 1.         .0000439218     .0192760487     64.     0.
MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 1.         .0000479810     .0192761019     24.     0.
MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 1.         .0000440469     .0192760499     44.     0.
MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 1.         .0000439258     .0192760488     64.     0.
MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 1.         .0000439215     .0192760487     84.     0.
MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 2.         .0000471408     .0192769312     12.     0.
MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 2.         .0000439425     .0192760534     32.     0.
MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 2.         .0000439216     .0192760488     44.     0.
MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 2.         .0000439213     .0192760487     56.     0.
MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 2.         .0000442054     .0192760609     24.     0.
MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 2.         .0000439217     .0192760488     44.     0.
MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 2.         .0000439213     .0192760487     64.     0.
MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 2.         .0000439213     .0192760487     88.     0.
MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 2.         .0000440112     .0192760515     32.     0.
MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 2.         .0000439215     .0192760487     56.     104.
MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 2.         .0000439213     .0192760487     88.     0.
MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 2.         .0000439213     .0192760487     116.    0.
MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 2.         .0000439933     .0192760505     36.     0.
MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 2.         .0000439214     .0192760487     76.     0.
MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 2.         .0000439213     .0192760487     108.    0.
MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 2.         .0000439213     .0192760487     148.    0.

a=11.0+10.0*\cos\pi x, c=1, u=xy, Neumann
MG:   Richardson, r = 1, m1 = 1, m2 = 1, p = 1.         .1876122832     .2279146610     4.      0.
MG:   Richardson, r = 2, m1 = 1, m2 = 1, p = 1.         .1585512167     .1728720726     12.     0.
MG:   Richardson, r = 3, m1 = 1, m2 = 1, p = 1.         .1351001050     .1422323846     16.     0.
MG:   Richardson, r = 4, m1 = 1, m2 = 1, p = 1.         .1153106056     .1194670210     24.     92.
MG:   Richardson, r = 1, m1 = 2, m2 = 2, p = 1.         .1614296035     .1795517778     8.      0.
MG:   Richardson, r = 2, m1 = 2, m2 = 2, p = 1.         .1172674639     .1224455920     16.     0.
MG:   Richardson, r = 3, m1 = 2, m2 = 2, p = 1.         .0855909162     .0878339215     24.     0.
MG:   Richardson, r = 4, m1 = 2, m2 = 2, p = 1.         .0625294072     .0637850535     36.     0.
MG:   Richardson, r = 1, m1 = 3, m2 = 3, p = 1.         .1391239894     .1500587771     12.     0.
MG:   Richardson, r = 2, m1 = 3, m2 = 3, p = 1.         .0866832755     .0893446082     24.     0.
MG:   Richardson, r = 3, m1 = 3, m2 = 3, p = 1.         .0541994692     .0553469584     32.     0.
MG:   Richardson, r = 4, m1 = 3, m2 = 3, p = 1.         .0339116915     .0346673935     48.     0.
MG:   Richardson, r = 1, m1 = 4, m2 = 4, p = 1.         .1197553077     .1271519449     16.     0.
MG:   Richardson, r = 2, m1 = 4, m2 = 4, p = 1.         .0639787023     .0656198076     28.     0.
MG:   Richardson, r = 3, m1 = 4, m2 = 4, p = 1.         .0342816165     .0350925869     44.     0.
MG:   Richardson, r = 4, m1 = 4, m2 = 4, p = 1.         .0183802288     .0191112979     60.     0.
MG:   Richardson, r = 1, m1 = 1, m2 = 1, p = 2.         .0023345911     .0169872913     12.     0.
MG:   Richardson, r = 2, m1 = 1, m2 = 1, p = 2.         .0002655194     .0088923772     28.     0.
MG:   Richardson, r = 3, m1 = 1, m2 = 1, p = 2.         .0001317085     .0065526935     44.     0.
MG:   Richardson, r = 4, m1 = 1, m2 = 1, p = 2.         .0000793071     .0055447479     56.     0.
MG:   Richardson, r = 1, m1 = 2, m2 = 2, p = 2.         .0002842696     .0090348207     20.     0.
MG:   Richardson, r = 2, m1 = 2, m2 = 2, p = 2.         .0000823650     .0055617285     44.     0.
MG:   Richardson, r = 3, m1 = 2, m2 = 2, p = 2.         .0000403830     .0047143020     64.     0.
MG:   Richardson, r = 4, m1 = 2, m2 = 2, p = 2.         .0000261580     .0043787749     88.     0.
MG:   Richardson, r = 1, m1 = 3, m2 = 3, p = 2.         .0001458822     .0066456793     28.     0.
MG:   Richardson, r = 2, m1 = 3, m2 = 3, p = 2.         .0000416069     .0047183202     56.     0.
MG:   Richardson, r = 3, m1 = 3, m2 = 3, p = 2.         .0000230491     .0042829364     88.     0.
MG:   Richardson, r = 4, m1 = 3, m2 = 3, p = 2.         .0000177457     .0041151168     112.    0.
MG:   Richardson, r = 1, m1 = 4, m2 = 4, p = 2.         .0000901827     .0056013172     36.     0.
MG:   Richardson, r = 2, m1 = 4, m2 = 4, p = 2.         .0000273896     .0043814525     68.     0.
MG:   Richardson, r = 3, m1 = 4, m2 = 4, p = 2.         .0000179582     .0041153801     104.    0.
MG:   Richardson, r = 4, m1 = 4, m2 = 4, p = 2.         .0000155688     .0040154803     144.    0.
MG:       Jacobi, r = 1, m1 = 1, m2 = 1, p = 1.         .1373365850     .1375621996     8.      0.
MG:       Jacobi, r = 2, m1 = 1, m2 = 1, p = 1.         .0847267020     .0850326358     12.     0.
MG:       Jacobi, r = 3, m1 = 1, m2 = 1, p = 1.         .0522137428     .0526724621     20.     0.
MG:       Jacobi, r = 4, m1 = 1, m2 = 1, p = 1.         .0320974572     .0328162737     24.     0.
MG:       Jacobi, r = 1, m1 = 2, m2 = 2, p = 1.         .0861963243     .0865034586     8.      0.
MG:       Jacobi, r = 2, m1 = 2, m2 = 2, p = 1.         .0336763634     .0343648729     16.     0.
MG:       Jacobi, r = 3, m1 = 2, m2 = 2, p = 1.         .0131573365     .0147989151     28.     0.
MG:       Jacobi, r = 4, m1 = 2, m2 = 2, p = 1.         .0051333198     .0084910564     36.     96.
MG:       Jacobi, r = 1, m1 = 3, m2 = 3, p = 1.         .0530468642     .0535046303     12.     0.
MG:       Jacobi, r = 2, m1 = 3, m2 = 3, p = 1.         .0124737046     .0141960308     20.     0.
MG:       Jacobi, r = 3, m1 = 3, m2 = 3, p = 1.         .0028130529     .0073272165     32.     0.
MG:       Jacobi, r = 4, m1 = 3, m2 = 3, p = 1.         .0005784984     .0067893981     44.     0.
MG:       Jacobi, r = 1, m1 = 4, m2 = 4, p = 1.         .0337957562     .0344833337     12.     0.
MG:       Jacobi, r = 2, m1 = 4, m2 = 4, p = 1.         .0051176375     .0084819693     24.     0.
MG:       Jacobi, r = 3, m1 = 4, m2 = 4, p = 1.         .0007531993     .0068036563     44.     0.
MG:       Jacobi, r = 4, m1 = 4, m2 = 4, p = 1.         .0001050809     .0067620969     52.     0.
MG:       Jacobi, r = 1, m1 = 1, m2 = 1, p = 2.         .0000126061     .0067795009     16.     0.
MG:       Jacobi, r = 2, m1 = 1, m2 = 1, p = 2.         .0000125125     .0067680571     28.     0.
MG:       Jacobi, r = 3, m1 = 1, m2 = 1, p = 2.         .0000125124     .0067672152     40.     0.
MG:       Jacobi, r = 4, m1 = 1, m2 = 1, p = 2.         .0000125119     .0067668035     52.     0.
MG:       Jacobi, r = 1, m1 = 2, m2 = 2, p = 2.         .0000333710     .0067639137     20.     0.
MG:       Jacobi, r = 2, m1 = 2, m2 = 2, p = 2.         .0000331558     .0067630569     40.     0.
MG:       Jacobi, r = 3, m1 = 2, m2 = 2, p = 2.         .0000331523     .0067625980     60.     0.
MG:       Jacobi, r = 4, m1 = 2, m2 = 2, p = 2.         .0000331513     .0067621798     80.     0.
MG:       Jacobi, r = 1, m1 = 3, m2 = 3, p = 2.         .0000125063     .0067663855     24.     0.
MG:       Jacobi, r = 2, m1 = 3, m2 = 3, p = 2.         .0000124969     .0067658216     52.     0.
MG:       Jacobi, r = 3, m1 = 3, m2 = 3, p = 2.         .0000124957     .0067652668     80.     0.
MG:       Jacobi, r = 4, m1 = 3, m2 = 3, p = 2.         .0000124946     .0067647380     104.    0.
MG:       Jacobi, r = 1, m1 = 4, m2 = 4, p = 2.         .0000332728     .0067632446     32.     0.
MG:       Jacobi, r = 2, m1 = 4, m2 = 4, p = 2.         .0000331261     .0067625323     68.     0.
MG:       Jacobi, r = 3, m1 = 4, m2 = 4, p = 2.         .0000331241     .0067619370     96.     0.
MG:       Jacobi, r = 4, m1 = 4, m2 = 4, p = 2.         .0000331229     .0067613649     132.    0.
MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 1.         .0659425953     .0660918607     4.      0.
MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 1.         .0166038366     .0170650194     16.     0.
MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 1.         .0041758966     .0057191896     24.     0.
MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 1.         .0010466401     .0040438771     28.     0.
MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 1.         .0176034345     .0180406651     12.     0.
MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 1.         .0011916193     .0040838316     24.     0.
MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 1.         .0000799977     .0039069005     36.     0.
MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 1.         .0000149053     .0039061203     48.     0.
MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 1.         .0047493706     .0061509589     16.     0.
MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 1.         .0000861475     .0039070335     32.     0.
MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 1.         .0000141011     .0039061180     52.     0.
MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 1.         .0000140062     .0039061176     64.     0.
MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 1.         .0012820405     .0041111973     24.     0.
MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 1.         .0000153151     .0039061239     40.     0.
MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 1.         .0000140078     .0039061177     64.     0.
MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 1.         .0000140204     .0039061177     84.     0.
MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 2.         .0000167739     .0039078986     12.     0.
MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 2.         .0000140734     .0039061299     32.     0.
MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 2.         .0000140336     .0039061180     44.     0.
MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 2.         .0000140322     .0039061177     56.     0.
MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 2.         .0000143601     .0039061783     24.     0.
MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 2.         .0000140334     .0039061184     44.     0.
MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 2.         .0000140321     .0039061177     64.     0.
MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 2.         .0000140320     .0039061177     88.     0.
MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 2.         .0000141185     .0039061384     28.     0.
MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 2.         .0000140324     .0039061179     60.     0.
MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 2.         .0000140321     .0039061177     88.     0.
MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 2.         .0000140320     .0039061177     116.    0.
MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 2.         .0000140837     .0039061313     36.     0.
MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 2.         .0000140322     .0039061178     76.     0.
MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 2.         .0000140320     .0039061177     108.    0.
MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 2.         .0000140320     .0039061177     148.    0.

a=1, c=1+xy, u=x^2-y^2, Neumann
MG:   Richardson, r = 1, m1 = 1, m2 = 1, p = 1.         .0004335031     .0072255096     4.      0.
MG:   Richardson, r = 2, m1 = 1, m2 = 1, p = 1.         .0000079018     .0055800862     12.     0.
MG:   Richardson, r = 3, m1 = 1, m2 = 1, p = 1.         .0000064103     .0055292311     16.     0.
MG:   Richardson, r = 4, m1 = 1, m2 = 1, p = 1.         .0000062553     .0055248823     24.     0.
MG:   Richardson, r = 1, m1 = 2, m2 = 2, p = 1.         .0000083455     .0055809214     8.      0.
MG:   Richardson, r = 2, m1 = 2, m2 = 2, p = 1.         .0000062616     .0055248855     16.     0.
MG:   Richardson, r = 3, m1 = 2, m2 = 2, p = 1.         .0000062356     .0055242939     28.     0.
MG:   Richardson, r = 4, m1 = 2, m2 = 2, p = 1.         .0000062343     .0055242761     32.     0.
MG:   Richardson, r = 1, m1 = 3, m2 = 3, p = 1.         .0000065427     .0055293180     12.     0.
MG:   Richardson, r = 2, m1 = 3, m2 = 3, p = 1.         .0000062364     .0055242942     24.     0.
MG:   Richardson, r = 3, m1 = 3, m2 = 3, p = 1.         .0000062343     .0055242755     36.     0.
MG:   Richardson, r = 4, m1 = 3, m2 = 3, p = 1.         .0000062343     .0055242753     44.     0.
MG:   Richardson, r = 1, m1 = 4, m2 = 4, p = 1.         .0000062880     .0055248974     16.     0.
MG:   Richardson, r = 2, m1 = 4, m2 = 4, p = 1.         .0000062347     .0055242762     28.     0.
MG:   Richardson, r = 3, m1 = 4, m2 = 4, p = 1.         .0000062343     .0055242753     48.     0.
MG:   Richardson, r = 4, m1 = 4, m2 = 4, p = 1.         .0000062343     .0055242753     64.     0.
MG:   Richardson, r = 1, m1 = 1, m2 = 1, p = 2.         .0000094840     .0062323128     16.     0.
MG:   Richardson, r = 2, m1 = 1, m2 = 1, p = 2.         .0000065014     .0055777384     32.     0.
MG:   Richardson, r = 3, m1 = 1, m2 = 1, p = 2.         .0000062620     .0055290945     44.     0.
MG:   Richardson, r = 4, m1 = 1, m2 = 1, p = 2.         .0000062384     .0055248681     60.     0.
MG:   Richardson, r = 1, m1 = 2, m2 = 2, p = 2.         .0000065049     .0055782886     20.     0.
MG:   Richardson, r = 2, m1 = 2, m2 = 2, p = 2.         .0000062392     .0055248692     44.     0.
MG:   Richardson, r = 3, m1 = 2, m2 = 2, p = 2.         .0000062345     .0055242935     64.     0.
MG:   Richardson, r = 4, m1 = 2, m2 = 2, p = 2.         .0000062343     .0055242761     88.     0.
MG:   Richardson, r = 1, m1 = 3, m2 = 3, p = 2.         .0000062659     .0055291848     28.     0.
MG:   Richardson, r = 2, m1 = 3, m2 = 3, p = 2.         .0000062347     .0055242935     56.     0.
MG:   Richardson, r = 3, m1 = 3, m2 = 3, p = 2.         .0000062343     .0055242755     88.     0.
MG:   Richardson, r = 4, m1 = 3, m2 = 3, p = 2.         .0000062343     .0055242753     112.    0.
MG:   Richardson, r = 1, m1 = 4, m2 = 4, p = 2.         .0000062420     .0055248830     36.     0.
MG:   Richardson, r = 2, m1 = 4, m2 = 4, p = 2.         .0000062343     .0055242761     68.     96.
MG:   Richardson, r = 3, m1 = 4, m2 = 4, p = 2.         .0000062343     .0055242753     100.    0.
MG:   Richardson, r = 4, m1 = 4, m2 = 4, p = 2.         .0000062343     .0055242753     136.    0.
MG:       Jacobi, r = 1, m1 = 1, m2 = 1, p = 1.         .0000062342     .0055242753     8.      0.
MG:       Jacobi, r = 2, m1 = 1, m2 = 1, p = 1.         .0000062343     .0055242753     8.      0.
MG:       Jacobi, r = 3, m1 = 1, m2 = 1, p = 1.         .0000062343     .0055242753     16.     0.
MG:       Jacobi, r = 4, m1 = 1, m2 = 1, p = 1.         .0000062343     .0055242753     24.     0.
MG:       Jacobi, r = 1, m1 = 2, m2 = 2, p = 1.         .0000062343     .0055242753     8.      0.
MG:       Jacobi, r = 2, m1 = 2, m2 = 2, p = 1.         .0000062343     .0055242753     16.     0.
MG:       Jacobi, r = 3, m1 = 2, m2 = 2, p = 1.         .0000062343     .0055242753     24.     0.
MG:       Jacobi, r = 4, m1 = 2, m2 = 2, p = 1.         .0000062343     .0055242753     32.     0.
MG:       Jacobi, r = 1, m1 = 3, m2 = 3, p = 1.         .0000062343     .0055242753     12.     0.
MG:       Jacobi, r = 2, m1 = 3, m2 = 3, p = 1.         .0000062343     .0055242753     20.     0.
MG:       Jacobi, r = 3, m1 = 3, m2 = 3, p = 1.         .0000062343     .0055242753     36.     0.
MG:       Jacobi, r = 4, m1 = 3, m2 = 3, p = 1.         .0000062343     .0055242753     40.     0.
MG:       Jacobi, r = 1, m1 = 4, m2 = 4, p = 1.         .0000062343     .0055242753     16.     0.
MG:       Jacobi, r = 2, m1 = 4, m2 = 4, p = 1.         .0000062343     .0055242753     28.     0.
MG:       Jacobi, r = 3, m1 = 4, m2 = 4, p = 1.         .0000062343     .0055242753     40.     0.
MG:       Jacobi, r = 4, m1 = 4, m2 = 4, p = 1.         .0000062343     .0055242753     56.     0.
MG:       Jacobi, r = 1, m1 = 1, m2 = 1, p = 2.         .0000062343     .0055242753     12.     0.
MG:       Jacobi, r = 2, m1 = 1, m2 = 1, p = 2.         .0000062343     .0055242753     28.     0.
MG:       Jacobi, r = 3, m1 = 1, m2 = 1, p = 2.         .0000062343     .0055242753     40.     0.
MG:       Jacobi, r = 4, m1 = 1, m2 = 1, p = 2.         .0000062343     .0055242753     52.     0.
MG:       Jacobi, r = 1, m1 = 2, m2 = 2, p = 2.         .0000062343     .0055242753     20.     0.
MG:       Jacobi, r = 2, m1 = 2, m2 = 2, p = 2.         .0000062343     .0055242753     40.     0.
MG:       Jacobi, r = 3, m1 = 2, m2 = 2, p = 2.         .0000062343     .0055242753     60.     0.
MG:       Jacobi, r = 4, m1 = 2, m2 = 2, p = 2.         .0000062343     .0055242753     80.     0.
MG:       Jacobi, r = 1, m1 = 3, m2 = 3, p = 2.         .0000062343     .0055242753     24.     0.
MG:       Jacobi, r = 2, m1 = 3, m2 = 3, p = 2.         .0000062343     .0055242753     56.     0.
MG:       Jacobi, r = 3, m1 = 3, m2 = 3, p = 2.         .0000062343     .0055242753     80.     0.
MG:       Jacobi, r = 4, m1 = 3, m2 = 3, p = 2.         .0000062343     .0055242753     108.    0.
MG:       Jacobi, r = 1, m1 = 4, m2 = 4, p = 2.         .0000062343     .0055242753     36.     0.
MG:       Jacobi, r = 2, m1 = 4, m2 = 4, p = 2.         .0000062343     .0055242753     68.     0.
MG:       Jacobi, r = 3, m1 = 4, m2 = 4, p = 2.         .0000062343     .0055242753     104.    0.
MG:       Jacobi, r = 4, m1 = 4, m2 = 4, p = 2.         .0000062343     .0055242753     136.    96.
MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 1.         .0000062624     .0055252170     8.      0.
MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 1.         .0000062342     .0055242783     12.     0.
MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 1.         .0000062342     .0055242753     24.     0.
MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 1.         .0000062343     .0055242753     28.     0.
MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 1.         .0000062345     .0055242879     12.     0.
MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 1.         .0000062342     .0055242753     24.     0.
MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 1.         .0000062343     .0055242753     36.     0.
MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 1.         .0000062343     .0055242753     44.     0.
MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 1.         .0000062340     .0055242764     20.     0.
MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 1.         .0000062343     .0055242753     32.     0.
MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 1.         .0000062343     .0055242753     48.     0.
MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 1.         .0000062343     .0055242753     68.     0.
MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 1.         .0000062341     .0055242755     20.     0.
MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 1.         .0000062343     .0055242753     44.     0.
MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 1.         .0000062343     .0055242753     64.     0.
MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 1.         .0000062343     .0055242753     84.     0.
MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 2.         .0000062559     .0055251897     12.     0.
MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 2.         .0000062343     .0055242781     32.     0.
MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 2.         .0000062343     .0055242753     40.     0.
MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 2.         .0000062343     .0055242753     60.     0.
MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 2.         .0000062346     .0055242874     20.     0.
MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 2.         .0000062343     .0055242753     44.     0.
MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 2.         .0000062343     .0055242753     68.     0.
MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 2.         .0000062343     .0055242753     88.     0.
MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 2.         .0000062342     .0055242764     32.     0.
MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 2.         .0000062343     .0055242753     60.     0.
MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 2.         .0000062343     .0055242753     88.     0.
MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 2.         .0000062343     .0055242753     116.    0.
MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 2.         .0000062342     .0055242755     36.     0.
MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 2.         .0000062343     .0055242753     72.     0.
MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 2.         .0000062343     .0055242753     112.    0.
MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 2.         .0000062343     .0055242753     148.    104.

a=1, c=2+\cos(\pi y), u=\cos(\pi x), Neumann
MG:   Richardson, r = 1, m1 = 1, m2 = 1, p = 1.         .0050343571     .0376165399     4.      0.
MG:   Richardson, r = 2, m1 = 1, m2 = 1, p = 1.         .0004439634     .0141948469     12.     0.
MG:   Richardson, r = 3, m1 = 1, m2 = 1, p = 1.         .0000918513     .0136603610     16.     0.
MG:   Richardson, r = 4, m1 = 1, m2 = 1, p = 1.         .0000453637     .0136339424     24.     0.
MG:   Richardson, r = 1, m1 = 2, m2 = 2, p = 1.         .0005794607     .0144391395     8.      0.
MG:   Richardson, r = 2, m1 = 2, m2 = 2, p = 1.         .0000528121     .0136352031     16.     0.
MG:   Richardson, r = 3, m1 = 2, m2 = 2, p = 1.         .0000329382     .0136302391     24.     0.
MG:   Richardson, r = 4, m1 = 2, m2 = 2, p = 1.         .0000298876     .0136299060     32.     0.
MG:   Richardson, r = 1, m1 = 3, m2 = 3, p = 1.         .0001636591     .0136942651     12.     0.
MG:   Richardson, r = 2, m1 = 3, m2 = 3, p = 1.         .0000352815     .0136304128     24.     0.
MG:   Richardson, r = 3, m1 = 3, m2 = 3, p = 1.         .0000298213     .0136298879     32.     0.
MG:   Richardson, r = 4, m1 = 3, m2 = 3, p = 1.         .0000289119     .0136298312     48.     0.
MG:   Richardson, r = 1, m1 = 4, m2 = 4, p = 1.         .0000841558     .0136420077     16.     0.
MG:   Richardson, r = 2, m1 = 4, m2 = 4, p = 1.         .0000317135     .0136300014     28.     0.
MG:   Richardson, r = 3, m1 = 4, m2 = 4, p = 1.         .0000291177     .0136298399     44.     0.
MG:   Richardson, r = 4, m1 = 4, m2 = 4, p = 1.         .0000287081     .0136298192     56.     0.
MG:   Richardson, r = 1, m1 = 1, m2 = 1, p = 2.         .0000891241     .0154566987     16.     0.
MG:   Richardson, r = 2, m1 = 1, m2 = 1, p = 2.         .0000402929     .0137712710     28.     0.
MG:   Richardson, r = 3, m1 = 1, m2 = 1, p = 2.         .0000314764     .0136432890     40.     0.
MG:   Richardson, r = 4, m1 = 1, m2 = 1, p = 2.         .0000293567     .0136316163     60.     92.
MG:   Richardson, r = 1, m1 = 2, m2 = 2, p = 2.         .0000517316     .0137751234     20.     0.
MG:   Richardson, r = 2, m1 = 2, m2 = 2, p = 2.         .0000300571     .0136316826     40.     0.
MG:   Richardson, r = 3, m1 = 2, m2 = 2, p = 2.         .0000287110     .0136298992     64.     0.
MG:   Richardson, r = 4, m1 = 2, m2 = 2, p = 2.         .0000286084     .0136298237     80.     0.
MG:   Richardson, r = 1, m1 = 3, m2 = 3, p = 2.         .0000397441     .0136443912     28.     0.
MG:   Richardson, r = 2, m1 = 3, m2 = 3, p = 2.         .0000288149     .0136299159     56.     0.
MG:   Richardson, r = 3, m1 = 3, m2 = 3, p = 2.         .0000286039     .0136298206     80.     0.
MG:   Richardson, r = 4, m1 = 3, m2 = 3, p = 2.         .0000285974     .0136298149     108.    0.
MG:   Richardson, r = 1, m1 = 4, m2 = 4, p = 2.         .0000341779     .0136320408     36.     0.
MG:   Richardson, r = 2, m1 = 4, m2 = 4, p = 2.         .0000286366     .0136298352     68.     0.
MG:   Richardson, r = 3, m1 = 4, m2 = 4, p = 2.         .0000285977     .0136298159     100.    0.
MG:   Richardson, r = 4, m1 = 4, m2 = 4, p = 2.         .0000285970     .0136298140     136.    0.
MG:       Jacobi, r = 1, m1 = 1, m2 = 1, p = 1.         .0000884709     .0136450228     4.      0.
MG:       Jacobi, r = 2, m1 = 1, m2 = 1, p = 1.         .0000318932     .0136309572     8.      0.
MG:       Jacobi, r = 3, m1 = 1, m2 = 1, p = 1.         .0000291981     .0136305338     16.     0.
MG:       Jacobi, r = 4, m1 = 1, m2 = 1, p = 1.         .0000287532     .0136304023     20.     0.
MG:       Jacobi, r = 1, m1 = 2, m2 = 2, p = 1.         .0000455261     .0136306388     8.      0.
MG:       Jacobi, r = 2, m1 = 2, m2 = 2, p = 1.         .0000296926     .0136298910     20.     0.
MG:       Jacobi, r = 3, m1 = 2, m2 = 2, p = 1.         .0000287260     .0136298574     24.     0.
MG:       Jacobi, r = 4, m1 = 2, m2 = 2, p = 1.         .0000286112     .0136298468     36.     0.
MG:       Jacobi, r = 1, m1 = 3, m2 = 3, p = 1.         .0000389202     .0136308103     12.     0.
MG:       Jacobi, r = 2, m1 = 3, m2 = 3, p = 1.         .0000291983     .0136302141     24.     0.
MG:       Jacobi, r = 3, m1 = 3, m2 = 3, p = 1.         .0000286733     .0136301304     36.     0.
MG:       Jacobi, r = 4, m1 = 3, m2 = 3, p = 1.         .0000286255     .0136300850     48.     0.
MG:       Jacobi, r = 1, m1 = 4, m2 = 4, p = 1.         .0000355074     .0136300942     16.     0.
MG:       Jacobi, r = 2, m1 = 4, m2 = 4, p = 1.         .0000289334     .0136298778     32.     0.
MG:       Jacobi, r = 3, m1 = 4, m2 = 4, p = 1.         .0000286172     .0136298560     44.     0.
MG:       Jacobi, r = 4, m1 = 4, m2 = 4, p = 1.         .0000285943     .0136298464     60.     0.
MG:       Jacobi, r = 1, m1 = 1, m2 = 1, p = 2.         .0000288704     .0136338567     16.     0.
MG:       Jacobi, r = 2, m1 = 1, m2 = 1, p = 2.         .0000286340     .0136308014     28.     0.
MG:       Jacobi, r = 3, m1 = 1, m2 = 1, p = 2.         .0000286267     .0136305170     40.     0.
MG:       Jacobi, r = 4, m1 = 1, m2 = 1, p = 2.         .0000286253     .0136304018     60.     0.
MG:       Jacobi, r = 1, m1 = 2, m2 = 2, p = 2.         .0000286989     .0136299592     20.     0.
MG:       Jacobi, r = 2, m1 = 2, m2 = 2, p = 2.         .0000285950     .0136298706     44.     0.
MG:       Jacobi, r = 3, m1 = 2, m2 = 2, p = 2.         .0000285931     .0136298542     64.     96.
MG:       Jacobi, r = 4, m1 = 2, m2 = 2, p = 2.         .0000285928     .0136298457     80.     0.
MG:       Jacobi, r = 1, m1 = 3, m2 = 3, p = 2.         .0000287312     .0136304520     24.     0.
MG:       Jacobi, r = 2, m1 = 3, m2 = 3, p = 2.         .0000286228     .0136302032     52.     0.
MG:       Jacobi, r = 3, m1 = 3, m2 = 3, p = 2.         .0000286210     .0136301304     80.     0.
MG:       Jacobi, r = 4, m1 = 3, m2 = 3, p = 2.         .0000286207     .0136300858     104.    0.
MG:       Jacobi, r = 1, m1 = 4, m2 = 4, p = 2.         .0000286823     .0136299405     32.     0.
MG:       Jacobi, r = 2, m1 = 4, m2 = 4, p = 2.         .0000285939     .0136298724     64.     0.
MG:       Jacobi, r = 3, m1 = 4, m2 = 4, p = 2.         .0000285926     .0136298552     100.    0.
MG:       Jacobi, r = 4, m1 = 4, m2 = 4, p = 2.         .0000285923     .0136298460     128.    0.
MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 1.         .0000454562     .0136321199     8.      0.
MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 1.         .0000295890     .0136298426     16.     0.
MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 1.         .0000286931     .0136298160     20.     0.
MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 1.         .0000286072     .0136298139     28.     0.
MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 1.         .0000339094     .0136299731     16.     0.
MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 1.         .0000288066     .0136298181     20.     0.
MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 1.         .0000286079     .0136298139     36.     0.
MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 1.         .0000285976     .0136298137     48.     0.
MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 1.         .0000317119     .0136298828     16.     0.
MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 1.         .0000286865     .0136298154     36.     0.
MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 1.         .0000286002     .0136298137     48.     0.
MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 1.         .0000285971     .0136298137     64.     0.
MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 1.         .0000308423     .0136298565     24.     0.
MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 1.         .0000286480     .0136298146     44.     0.
MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 1.         .0000285984     .0136298137     60.     0.
MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 1.         .0000285970     .0136298137     84.     0.
MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 2.         .0000310210     .0136314263     16.     0.
MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 2.         .0000286182     .0136298234     28.     0.
MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 2.         .0000285975     .0136298140     44.     0.
MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 2.         .0000285970     .0136298137     60.     0.
MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 2.         .0000288789     .0136298634     20.     0.
MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 2.         .0000285981     .0136298145     44.     0.
MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 2.         .0000285970     .0136298137     68.     0.
MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 2.         .0000285970     .0136298137     84.     0.
MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 2.         .0000286963     .0136298343     32.     0.
MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 2.         .0000285975     .0136298140     56.     104.
MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 2.         .0000285970     .0136298137     88.     0.
MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 2.         .0000285970     .0136298137     116.    0.
MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 2.         .0000286664     .0136298276     36.     0.
MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 2.         .0000285973     .0136298138     76.     0.
MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 2.         .0000285970     .0136298137     108.    0.
MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 2.         .0000285970     .0136298137     148.    0.

a=1, c=1, gamma=2+\cos(\pi x), u=\sin(\pi x)+\sin(\pi y), Robin
MG:   Richardson, r = 1, m1 = 1, m2 = 1, p = 1.         .0060856079     .0367856920     4.      0.
MG:   Richardson, r = 2, m1 = 1, m2 = 1, p = 1.         .0002786845     .0194991879     12.     0.
MG:   Richardson, r = 3, m1 = 1, m2 = 1, p = 1.         .0000711390     .0192923266     16.     96.
MG:   Richardson, r = 4, m1 = 1, m2 = 1, p = 1.         .0000444349     .0192778985     20.     0.
MG:   Richardson, r = 1, m1 = 2, m2 = 2, p = 1.         .0004243314     .0195468623     8.      0.
MG:   Richardson, r = 2, m1 = 2, m2 = 2, p = 1.         .0000501160     .0192779852     20.     0.
MG:   Richardson, r = 3, m1 = 2, m2 = 2, p = 1.         .0000358350     .0192760973     24.     0.
MG:   Richardson, r = 4, m1 = 2, m2 = 2, p = 1.         .0000333563     .0192760529     36.     0.
MG:   Richardson, r = 1, m1 = 3, m2 = 3, p = 1.         .0001248007     .0192962237     12.     0.
MG:   Richardson, r = 2, m1 = 3, m2 = 3, p = 1.         .0000375926     .0192761019     20.     0.
MG:   Richardson, r = 3, m1 = 3, m2 = 3, p = 1.         .0000333125     .0192760518     36.     0.
MG:   Richardson, r = 4, m1 = 3, m2 = 3, p = 1.         .0000326488     .0192760538     48.     0.
MG:   Richardson, r = 1, m1 = 4, m2 = 4, p = 1.         .0000714042     .0192787059     12.     0.
MG:   Richardson, r = 2, m1 = 4, m2 = 4, p = 1.         .0000347376     .0192760507     32.     0.
MG:   Richardson, r = 3, m1 = 4, m2 = 4, p = 1.         .0000327911     .0192760533     44.     0.
MG:   Richardson, r = 4, m1 = 4, m2 = 4, p = 1.         .0000325329     .0192760542     56.     0.
MG:   Richardson, r = 1, m1 = 1, m2 = 1, p = 2.         .0000941203     .0216732392     16.     0.
MG:   Richardson, r = 2, m1 = 1, m2 = 1, p = 2.         .0000435950     .0194524602     28.     0.
MG:   Richardson, r = 3, m1 = 1, m2 = 1, p = 2.         .0000350048     .0192912533     40.     0.
MG:   Richardson, r = 4, m1 = 1, m2 = 1, p = 2.         .0000331039     .0192778245     56.     0.
MG:   Richardson, r = 1, m1 = 2, m2 = 2, p = 2.         .0000555903     .0194535487     20.     0.
MG:   Richardson, r = 2, m1 = 2, m2 = 2, p = 2.         .0000337220     .0192778245     40.     0.
MG:   Richardson, r = 3, m1 = 2, m2 = 2, p = 2.         .0000325653     .0192761026     64.     0.
MG:   Richardson, r = 4, m1 = 2, m2 = 2, p = 2.         .0000324896     .0192760560     80.     0.
MG:   Richardson, r = 1, m1 = 3, m2 = 3, p = 2.         .0000431374     .0192914104     28.     0.
MG:   Richardson, r = 2, m1 = 3, m2 = 3, p = 2.         .0000326459     .0192761020     56.     0.
MG:   Richardson, r = 3, m1 = 3, m2 = 3, p = 2.         .0000324869     .0192760546     80.     0.
MG:   Richardson, r = 4, m1 = 3, m2 = 3, p = 2.         .0000324835     .0192760543     108.    0.
MG:   Richardson, r = 1, m1 = 4, m2 = 4, p = 2.         .0000375824     .0192778331     36.     0.
MG:   Richardson, r = 2, m1 = 4, m2 = 4, p = 2.         .0000325075     .0192760558     72.     0.
MG:   Richardson, r = 3, m1 = 4, m2 = 4, p = 2.         .0000324836     .0192760543     108.    0.
MG:   Richardson, r = 4, m1 = 4, m2 = 4, p = 2.         .0000324834     .0192760543     140.    0.
MG:       Jacobi, r = 1, m1 = 1, m2 = 1, p = 1.         .0001038624     .0192804431     8.      0.
MG:       Jacobi, r = 2, m1 = 1, m2 = 1, p = 1.         .0000359478     .0192760560     12.     0.
MG:       Jacobi, r = 3, m1 = 1, m2 = 1, p = 1.         .0000330890     .0192760547     20.     0.
MG:       Jacobi, r = 4, m1 = 1, m2 = 1, p = 1.         .0000326209     .0192760560     24.     0.
MG:       Jacobi, r = 1, m1 = 2, m2 = 2, p = 1.         .0000510544     .0192763218     8.      0.
MG:       Jacobi, r = 2, m1 = 2, m2 = 2, p = 1.         .0000336138     .0192760503     20.     0.
MG:       Jacobi, r = 3, m1 = 2, m2 = 2, p = 1.         .0000326242     .0192760541     28.     0.
MG:       Jacobi, r = 4, m1 = 2, m2 = 2, p = 1.         .0000325023     .0192760546     32.     0.
MG:       Jacobi, r = 1, m1 = 3, m2 = 3, p = 1.         .0000429981     .0192761228     12.     0.
MG:       Jacobi, r = 2, m1 = 3, m2 = 3, p = 1.         .0000330723     .0192760534     24.     0.
MG:       Jacobi, r = 3, m1 = 3, m2 = 3, p = 1.         .0000325461     .0192760558     32.     0.
MG:       Jacobi, r = 4, m1 = 3, m2 = 3, p = 1.         .0000324963     .0192760560     44.     0.
MG:       Jacobi, r = 1, m1 = 4, m2 = 4, p = 1.         .0000394822     .0192760683     12.     0.
MG:       Jacobi, r = 2, m1 = 4, m2 = 4, p = 1.         .0000328275     .0192760527     28.     0.
MG:       Jacobi, r = 3, m1 = 4, m2 = 4, p = 1.         .0000325082     .0192760544     44.     0.
MG:       Jacobi, r = 4, m1 = 4, m2 = 4, p = 1.         .0000324842     .0192760545     56.     0.
MG:       Jacobi, r = 1, m1 = 1, m2 = 1, p = 2.         .0000325835     .0192761285     16.     0.
MG:       Jacobi, r = 2, m1 = 1, m2 = 1, p = 2.         .0000324911     .0192760567     28.     0.
MG:       Jacobi, r = 3, m1 = 1, m2 = 1, p = 2.         .0000324906     .0192760565     40.     0.
MG:       Jacobi, r = 4, m1 = 1, m2 = 1, p = 2.         .0000324906     .0192760564     52.     0.
MG:       Jacobi, r = 1, m1 = 2, m2 = 2, p = 2.         .0000325414     .0192760543     20.     0.
MG:       Jacobi, r = 2, m1 = 2, m2 = 2, p = 2.         .0000324821     .0192760548     40.     0.
MG:       Jacobi, r = 3, m1 = 2, m2 = 2, p = 2.         .0000324819     .0192760547     60.     0.
MG:       Jacobi, r = 4, m1 = 2, m2 = 2, p = 2.         .0000324819     .0192760547     80.     0.
MG:       Jacobi, r = 1, m1 = 3, m2 = 3, p = 2.         .0000325492     .0192760555     28.     0.
MG:       Jacobi, r = 2, m1 = 3, m2 = 3, p = 2.         .0000324907     .0192760562     52.     0.
MG:       Jacobi, r = 3, m1 = 3, m2 = 3, p = 2.         .0000324905     .0192760561     76.     0.
MG:       Jacobi, r = 4, m1 = 3, m2 = 3, p = 2.         .0000324904     .0192760560     104.    0.
MG:       Jacobi, r = 1, m1 = 4, m2 = 4, p = 2.         .0000325389     .0192760535     32.     0.
MG:       Jacobi, r = 2, m1 = 4, m2 = 4, p = 2.         .0000324821     .0192760546     64.     0.
MG:       Jacobi, r = 3, m1 = 4, m2 = 4, p = 2.         .0000324820     .0192760546     96.     0.
MG:       Jacobi, r = 4, m1 = 4, m2 = 4, p = 2.         .0000324820     .0192760545     128.    0.
MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 1.         .0000488668     .0192770994     8.      0.
MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 1.         .0000334224     .0192760532     16.     0.
MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 1.         .0000325787     .0192760537     20.     0.
MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 1.         .0000324940     .0192760542     28.     0.
MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 1.         .0000377312     .0192760570     12.     0.
MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 1.         .0000326935     .0192760528     24.     0.
MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 1.         .0000324949     .0192760542     36.     0.
MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 1.         .0000324841     .0192760543     48.     0.
MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 1.         .0000355612     .0192760443     16.     0.
MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 1.         .0000325735     .0192760537     32.     0.
MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 1.         .0000324867     .0192760543     52.     0.
MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 1.         .0000324835     .0192760543     64.     0.
MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 1.         .0000346663     .0192760446     20.     0.
MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 1.         .0000325337     .0192760540     44.     104.
MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 1.         .0000324848     .0192760543     64.     0.
MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 1.         .0000324834     .0192760543     84.     0.
MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 2.         .0000343292     .0192768757     12.     0.
MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 2.         .0000324936     .0192760584     28.     0.
MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 2.         .0000324835     .0192760543     40.     0.
MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 2.         .0000324834     .0192760543     56.     0.
MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 2.         .0000326603     .0192760593     24.     0.
MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 2.         .0000324838     .0192760543     44.     0.
MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 2.         .0000324834     .0192760543     64.     0.
MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 2.         .0000324834     .0192760543     88.     0.
MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 2.         .0000325510     .0192760535     28.     0.
MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 2.         .0000324835     .0192760543     60.     0.
MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 2.         .0000324834     .0192760543     84.     0.
MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 2.         .0000324834     .0192760543     116.    0.
MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 2.         .0000325370     .0192760532     36.     0.
MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 2.         .0000324835     .0192760543     76.     0.
MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 2.         .0000324834     .0192760543     108.    0.
MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 2.         .0000324834     .0192760543     144.    0.

a=11.0+10.0*\cos\pi x, c=1+xy, gamma=2+\cos(\pi x), u=xy, Robin
MG:   Richardson, r = 1, m1 = 1, m2 = 1, p = 1.         .0554729263     .1245833796     8.      0.
MG:   Richardson, r = 2, m1 = 1, m2 = 1, p = 1.         .0128733208     .0567584012     12.     0.
MG:   Richardson, r = 3, m1 = 1, m2 = 1, p = 1.         .0039404586     .0345620318     16.     0.
MG:   Richardson, r = 4, m1 = 1, m2 = 1, p = 1.         .0022583834     .0237099987     20.     0.
MG:   Richardson, r = 1, m1 = 2, m2 = 2, p = 1.         .0144147737     .0644444683     8.      0.
MG:   Richardson, r = 2, m1 = 2, m2 = 2, p = 1.         .0027394669     .0264155247     20.     0.
MG:   Richardson, r = 3, m1 = 2, m2 = 2, p = 1.         .0012960781     .0144139033     24.     0.
MG:   Richardson, r = 4, m1 = 2, m2 = 2, p = 1.         .0006789477     .0092368615     32.     0.
MG:   Richardson, r = 1, m1 = 3, m2 = 3, p = 1.         .0057535858     .0432885026     12.     0.
MG:   Richardson, r = 2, m1 = 3, m2 = 3, p = 1.         .0014873859     .0156152543     24.     0.
MG:   Richardson, r = 3, m1 = 3, m2 = 3, p = 1.         .0005699442     .0081806653     36.     0.
MG:   Richardson, r = 4, m1 = 3, m2 = 3, p = 1.         .0002714460     .0056575480     44.     0.
MG:   Richardson, r = 1, m1 = 4, m2 = 4, p = 1.         .0037875242     .0317802685     16.     0.
MG:   Richardson, r = 2, m1 = 4, m2 = 4, p = 1.         .0008564355     .0103725952     28.     0.
MG:   Richardson, r = 3, m1 = 4, m2 = 4, p = 1.         .0002956523     .0057965379     44.     0.
MG:   Richardson, r = 4, m1 = 4, m2 = 4, p = 1.         .0001359881     .0045672381     56.     0.
MG:   Richardson, r = 1, m1 = 1, m2 = 1, p = 2.         .0007176820     .0153757374     16.     0.
MG:   Richardson, r = 2, m1 = 1, m2 = 1, p = 2.         .0002380084     .0084887242     28.     0.
MG:   Richardson, r = 3, m1 = 1, m2 = 1, p = 2.         .0001204251     .0063904538     40.     0.
MG:   Richardson, r = 4, m1 = 1, m2 = 1, p = 2.         .0000728873     .0054677306     56.     0.
MG:   Richardson, r = 1, m1 = 2, m2 = 2, p = 2.         .0002499523     .0085716178     20.     0.
MG:   Richardson, r = 2, m1 = 2, m2 = 2, p = 2.         .0000750194     .0054772877     40.     0.
MG:   Richardson, r = 3, m1 = 2, m2 = 2, p = 2.         .0000369480     .0046890661     64.     0.
MG:   Richardson, r = 4, m1 = 2, m2 = 2, p = 2.         .0000237718     .0043687385     80.     0.
MG:   Richardson, r = 1, m1 = 3, m2 = 3, p = 2.         .0001300452     .0064442160     28.     0.
MG:   Richardson, r = 2, m1 = 3, m2 = 3, p = 2.         .0000378859     .0046914205     52.     0.
MG:   Richardson, r = 3, m1 = 3, m2 = 3, p = 2.         .0000207994     .0042758250     84.     0.
MG:   Richardson, r = 4, m1 = 3, m2 = 3, p = 2.         .0000157704     .0041123959     108.    100.
MG:   Richardson, r = 1, m1 = 4, m2 = 4, p = 2.         .0000806208     .0055005085     32.     0.
MG:   Richardson, r = 2, m1 = 4, m2 = 4, p = 2.         .0000247754     .0043703907     68.     0.
MG:   Richardson, r = 3, m1 = 4, m2 = 4, p = 2.         .0000159540     .0041125722     104.    0.
MG:   Richardson, r = 4, m1 = 4, m2 = 4, p = 2.         .0000136861     .0040143796     132.    0.
MG:       Jacobi, r = 1, m1 = 1, m2 = 1, p = 1.         .0031868584     .0075890077     4.      0.
MG:       Jacobi, r = 2, m1 = 1, m2 = 1, p = 1.         .0000279198     .0067673084     8.      0.
MG:       Jacobi, r = 3, m1 = 1, m2 = 1, p = 1.         .0000192858     .0067660337     16.     0.
MG:       Jacobi, r = 4, m1 = 1, m2 = 1, p = 1.         .0000132277     .0067653774     20.     0.
MG:       Jacobi, r = 1, m1 = 2, m2 = 2, p = 1.         .0000619573     .0067650303     8.      0.
MG:       Jacobi, r = 2, m1 = 2, m2 = 2, p = 1.         .0000318314     .0067620783     20.     0.
MG:       Jacobi, r = 3, m1 = 2, m2 = 2, p = 1.         .0000292401     .0067611067     24.     0.
MG:       Jacobi, r = 4, m1 = 2, m2 = 2, p = 1.         .0000289388     .0067602402     32.     0.
MG:       Jacobi, r = 1, m1 = 3, m2 = 3, p = 1.         .0000225482     .0067654200     8.      0.
MG:       Jacobi, r = 2, m1 = 3, m2 = 3, p = 1.         .0000121854     .0067639773     24.     0.
MG:       Jacobi, r = 3, m1 = 3, m2 = 3, p = 1.         .0000119736     .0067627336     32.     0.
MG:       Jacobi, r = 4, m1 = 3, m2 = 3, p = 1.         .0000119436     .0067615394     44.     0.
MG:       Jacobi, r = 1, m1 = 4, m2 = 4, p = 1.         .0000378988     .0067626182     12.     0.
MG:       Jacobi, r = 2, m1 = 4, m2 = 4, p = 1.         .0000292863     .0067605921     28.     0.
MG:       Jacobi, r = 3, m1 = 4, m2 = 4, p = 1.         .0000288666     .0067590971     40.     0.
MG:       Jacobi, r = 4, m1 = 4, m2 = 4, p = 1.         .0000288303     .0067576526     56.     0.
MG:       Jacobi, r = 1, m1 = 1, m2 = 1, p = 2.         .0000120254     .0067783292     12.     0.
MG:       Jacobi, r = 2, m1 = 1, m2 = 1, p = 2.         .0000119596     .0067670579     28.     0.
MG:       Jacobi, r = 3, m1 = 1, m2 = 1, p = 2.         .0000119576     .0067660145     40.     0.
MG:       Jacobi, r = 4, m1 = 1, m2 = 1, p = 2.         .0000119565     .0067653801     52.     0.
MG:       Jacobi, r = 1, m1 = 2, m2 = 2, p = 2.         .0000290126     .0067632565     20.     0.
MG:       Jacobi, r = 2, m1 = 2, m2 = 2, p = 2.         .0000288646     .0067619902     44.     0.
MG:       Jacobi, r = 3, m1 = 2, m2 = 2, p = 2.         .0000288619     .0067610857     60.     0.
MG:       Jacobi, r = 4, m1 = 2, m2 = 2, p = 2.         .0000288609     .0067602234     84.     0.
MG:       Jacobi, r = 1, m1 = 3, m2 = 3, p = 2.         .0000119515     .0067652262     28.     0.
MG:       Jacobi, r = 2, m1 = 3, m2 = 3, p = 2.         .0000119425     .0067639841     56.     0.
MG:       Jacobi, r = 3, m1 = 3, m2 = 3, p = 2.         .0000119405     .0067627522     84.     0.
MG:       Jacobi, r = 4, m1 = 3, m2 = 3, p = 2.         .0000119386     .0067615533     108.    0.
MG:       Jacobi, r = 1, m1 = 4, m2 = 4, p = 2.         .0000289481     .0067621695     36.     0.
MG:       Jacobi, r = 2, m1 = 4, m2 = 4, p = 2.         .0000288356     .0067605755     68.     0.
MG:       Jacobi, r = 3, m1 = 4, m2 = 4, p = 2.         .0000288335     .0067590944     96.     0.
MG:       Jacobi, r = 4, m1 = 4, m2 = 4, p = 2.         .0000288319     .0067576476     128.    0.
MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 1.         .0000257780     .0039085377     8.      0.
MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 1.         .0000131985     .0039060725     16.     0.
MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 1.         .0000124344     .0039060321     20.     0.
MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 1.         .0000123382     .0039060286     28.     0.
MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 1.         .0000153647     .0039062226     12.     0.
MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 1.         .0000124662     .0039060345     24.     0.
MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 1.         .0000123291     .0039060285     36.     0.
MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 1.         .0000123202     .0039060281     48.     0.
MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 1.         .0000139703     .0039061171     16.     0.
MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 1.         .0000123770     .0039060306     32.     0.
MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 1.         .0000123219     .0039060282     48.     0.
MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 1.         .0000123196     .0039060281     68.     0.
MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 1.         .0000134716     .0039060864     20.     0.
MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 1.         .0000123508     .0039060294     40.     0.
MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 1.         .0000123205     .0039060281     60.     0.
MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 1.         .0000123196     .0039060281     84.     0.
MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 2.         .0000145082     .0039077319     16.     0.
MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 2.         .0000123512     .0039060389     28.     0.
MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 2.         .0000123207     .0039060284     44.     0.
MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 2.         .0000123196     .0039060281     56.     0.
MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 2.         .0000125699     .0039060769     20.     0.
MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 2.         .0000123205     .0039060286     44.     0.
MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 2.         .0000123196     .0039060281     68.     0.
MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 2.         .0000123195     .0039060281     84.     0.
MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 2.         .0000123854     .0039060436     32.     0.
MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 2.         .0000123198     .0039060283     56.     0.
MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 2.         .0000123195     .0039060281     88.     0.
MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 2.         .0000123195     .0039060281     116.    0.
MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 2.         .0000123599     .0039060379     36.     0.
MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 2.         .0000123197     .0039060282     72.     0.
MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 2.         .0000123195     .0039060281     112.    104.
MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 2.         .0000123195     .0039060281     144.    0.
> 

