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 1025 \times 1025 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 2049 \times 2049 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: model name : Intel(R) Xeon(R) CPU E3-1271 v3 @ 3.60GHz % uname -a Linux heine 4.2.0-30-generic #36-Ubuntu SMP Fri Feb 26 00:58:07 UTC 2016 x86_64 x86_64 x86_64 GNU/Linux % gsi -v v4.8.4 20160308184356 x86_64-unknown-linux-gnu "./configure 'CC=gcc -march=native' '--enable-single-host' '--enable-multiple-versions' '--enable-shared'" % gcc -v Using built-in specs. COLLECT_GCC=gcc COLLECT_LTO_WRAPPER=/usr/lib/gcc/x86_64-linux-gnu/5/lto-wrapper Target: x86_64-linux-gnu Configured with: ../src/configure -v --with-pkgversion='Ubuntu 5.2.1-22ubuntu2' --with-bugurl=file:///usr/share/doc/gcc-5/README.Bugs --enable-languages=c,ada,c++,java,go,d,fortran,objc,obj-c++ --prefix=/usr --program-suffix=-5 --enable-shared --enable-linker-build-id --libexecdir=/usr/lib --without-included-gettext --enable-threads=posix --libdir=/usr/lib --enable-nls --with-sysroot=/ --enable-clocale=gnu --enable-libstdcxx-debug --enable-libstdcxx-time=yes --with-default-libstdcxx-abi=new --enable-gnu-unique-object --disable-vtable-verify --enable-libmpx --enable-plugin --with-system-zlib --disable-browser-plugin --enable-java-awt=gtk --enable-gtk-cairo --with-java-home=/usr/lib/jvm/java-1.5.0-gcj-5-amd64/jre --enable-java-home --with-jvm-root-dir=/usr/lib/jvm/java-1.5.0-gcj-5-amd64 --with-jvm-jar-dir=/usr/lib/jvm-exports/java-1.5.0-gcj-5-amd64 --with-arch-directory=amd64 --with-ecj-jar=/usr/share/java/eclipse-ecj.jar --enable-objc-gc --enable-multiarch --disable-werror --with-arch-32=i686 --with-abi=m64 --with-multilib-list=m32,m64,mx32 --enable-multilib --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 5.2.1 20151010 (Ubuntu 5.2.1-22ubuntu2) I have not yet analyzed these results to see what the "best" combination of m1, m2, r, and p is. This version of final-tests.scm was modified to time various operations; here are the diffs: heine:~/WWW/615-2016/software> diff -up final-tests.scm final-tests.scm-new --- final-tests.scm 2016-03-21 16:25:41.214076520 -0400 +++ final-tests.scm-new 2016-03-11 12:15:28.206443217 -0500 @@ -7,12 +7,12 @@ (include "points.scm") (declare (core)) -(define number-of-refinements 7) +(define number-of-refinements 10) -(define triangulations (Triangulation->refined-triangulations coarse-square-triangulation number-of-refinements)) -(define spaces (triangulations->linear-finite-element-spaces triangulations)) -(define projectors (linear-finite-element-spaces->projectors spaces)) -(define injectors (projectors->injectors projectors)) +(define triangulations (time (Triangulation->refined-triangulations coarse-square-triangulation number-of-refinements))) +(define spaces (time (triangulations->linear-finite-element-spaces triangulations))) +(define projectors (time (linear-finite-element-spaces->projectors spaces))) +(define injectors (time (projectors->injectors projectors))) (define fes (vector-ref spaces number-of-refinements)) @@ -28,15 +28,15 @@ (define L2-inner-product-operator (->Linear-finite-element-operator fine-fes)) -(add-mass-terms! L2-inner-product-operator (lambda (p) 1.0)) +(time (add-mass-terms! L2-inner-product-operator (lambda (p) 1.0))) (define H1-inner-product-operator (->Linear-finite-element-operator fine-fes)) -(add-stiffness-terms! H1-inner-product-operator (lambda (p) 1.0)) +(time (add-stiffness-terms! H1-inner-product-operator (lambda (p) 1.0))) -(add-mass-terms! H1-inner-product-operator (lambda (p) 1.0)) +(time (add-mass-terms! H1-inner-product-operator (lambda (p) 1.0))) -(define coarse->fine-operator (Operator-adjoint (coarse<-fine fes fine-fes))) +(define coarse->fine-operator (time (Operator-adjoint (coarse<-fine fes fine-fes)))) (define dirichlet-tests (list @@ -512,11 +512,11 @@ (with-access problem-data (Problem-data a c f gamma g u) (let* ((check (interpolate-procedure->Linear-finite-element-vector u fine-fes)) - (operator (build-operator a c gamma fes)) - (rhs (build-rhs f g fes)) + (operator (time (build-operator a c gamma fes))) + (rhs (time (build-rhs f g fes))) (right-hand-sides (most-refined-rhs->all-refined-rhs rhs projectors)) (operators (most-refined-operator->all-refined-operators operator projectors injectors))) - + (for-each display (list #\newline (Problem-data-description problem-data) #\newline)) (for-each (lambda (smoother-maker smoother-name) heine:~/WWW/615-2016/software> gsc++ [ Meroon V3 Paques2001+1 $Revision: 1.2 $ ] Gambit v4.8.4 > (compile-file "all") "/home/lucier/WWW/615-2016/software/all.o2" > (load "all") "/home/lucier/WWW/615-2016/software/all.o2" > (compile-file "final-tests") "/home/lucier/WWW/615-2016/software/final-tests.o4" > (load "final-tests") (time (Triangulation->refined-triangulations coarse-square-triangulation number-of-refinements)) 10678 ms real time 10676 ms cpu time (9760 user, 916 system) 1 collection accounting for 4151 ms real time (3600 user, 552 system) 8868456768 bytes allocated 2772124 minor faults no major faults (time (triangulations->linear-finite-element-spaces triangulations)) 1524 ms real time 1520 ms cpu time (1364 user, 156 system) no collections 2728068512 bytes allocated 334345 minor faults no major faults (time (linear-finite-element-spaces->projectors spaces)) 4258 ms real time 4260 ms cpu time (3976 user, 284 system) 1 collection accounting for 3138 ms real time (3060 user, 76 system) 2711205664 bytes allocated 544476 minor faults no major faults (time (projectors->injectors projectors)) 556 ms real time 556 ms cpu time (532 user, 24 system) no collections 1979038592 bytes allocated 50512 minor faults no major faults (time (add-mass-terms! L2-inner-product-operator (lambda (p) 1.))) 12388 ms real time 12384 ms cpu time (12372 user, 12 system) 1 collection accounting for 3830 ms real time (3820 user, 8 system) 8187281472 bytes allocated 8068 minor faults no major faults (time (add-stiffness-terms! H1-inner-product-operator (lambda (p) 1.))) 13492 ms real time 13488 ms cpu time (13352 user, 136 system) 2 collections accounting for 4299 ms real time (4180 user, 120 system) 16106127424 bytes allocated 12924 minor faults no major faults (time (add-mass-terms! H1-inner-product-operator (lambda (p) 1.))) 12841 ms real time 12836 ms cpu time (12836 user, 0 system) 2 collections accounting for 4168 ms real time (4164 user, 0 system) 8187281472 bytes allocated 6 minor faults no major faults (time (Operator-adjoint (coarse<-fine fes fine-fes))) 7515 ms real time 7512 ms cpu time (7496 user, 16 system) 1 collection accounting for 2773 ms real time (2768 user, 4 system) 14098477248 bytes allocated 22803 minor faults no major faults "/home/lucier/WWW/615-2016/software/final-tests.o4" > (run-multigrid-tests symmetric-tests) (time (build-operator (instance-ref g1229 1) (instance-ref g1229 2) (instance-ref g1229 4) fes)) 7181 ms real time 7180 ms cpu time (7172 user, 8 system) 1 collection accounting for 2305 ms real time (2300 user, 4 system) 7162847824 bytes allocated 1460 minor faults no major faults (time (build-rhs (instance-ref g1229 5) (instance-ref g1229 6) fes)) 3593 ms real time 3592 ms cpu time (3592 user, 0 system) 1 collection accounting for 2285 ms real time (2284 user, 0 system) 1587233160 bytes allocated 4 minor faults no major faults 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. .0009048543 .0017644837 104. 0. MG: Richardson, r = 2, m1 = 1, m2 = 1, p = 1. .0000018274 .0000062713 200. 0. MG: Richardson, r = 3, m1 = 1, m2 = 1, p = 1. .0000000048 .0000000262 300. 0. MG: Richardson, r = 4, m1 = 1, m2 = 1, p = 1. .0000000008 .0000000009 396. 0. MG: Richardson, r = 1, m1 = 2, m2 = 2, p = 1. .0000028074 .0000122816 156. 0. MG: Richardson, r = 2, m1 = 2, m2 = 2, p = 1. .0000000008 .0000000010 308. 0. MG: Richardson, r = 3, m1 = 2, m2 = 2, p = 1. .0000000008 .0000000009 456. 0. MG: Richardson, r = 4, m1 = 2, m2 = 2, p = 1. .0000000007 .0000000008 604. 0. MG: Richardson, r = 1, m1 = 3, m2 = 3, p = 1. .0000000431 .0000002630 208. 0. MG: Richardson, r = 2, m1 = 3, m2 = 3, p = 1. .0000000007 .0000000008 416. 0. MG: Richardson, r = 3, m1 = 3, m2 = 3, p = 1. .0000000008 .0000000010 648. 0. MG: Richardson, r = 4, m1 = 3, m2 = 3, p = 1. .0000000007 .0000000008 812. 0. MG: Richardson, r = 1, m1 = 4, m2 = 4, p = 1. .0000000015 .0000000110 260. 0. MG: Richardson, r = 2, m1 = 4, m2 = 4, p = 1. .0000000007 .0000000008 516. 0. MG: Richardson, r = 3, m1 = 4, m2 = 4, p = 1. .0000000008 .0000000010 772. 0. MG: Richardson, r = 4, m1 = 4, m2 = 4, p = 1. .0000000007 .0000000008 1028. 0. MG: Richardson, r = 1, m1 = 1, m2 = 1, p = 2. .0000000008 .0000000023 168. 0. MG: Richardson, r = 2, m1 = 1, m2 = 1, p = 2. .0000000007 .0000000008 348. 0. MG: Richardson, r = 3, m1 = 1, m2 = 1, p = 2. .0000000009 .0000000010 544. 0. MG: Richardson, r = 4, m1 = 1, m2 = 1, p = 2. .0000000007 .0000000008 672. 2172. MG: Richardson, r = 1, m1 = 2, m2 = 2, p = 2. .0000000009 .0000000010 260. 0. MG: Richardson, r = 2, m1 = 2, m2 = 2, p = 2. .0000000007 .0000000008 516. 0. MG: Richardson, r = 3, m1 = 2, m2 = 2, p = 2. .0000000009 .0000000010 772. 0. MG: Richardson, r = 4, m1 = 2, m2 = 2, p = 2. .0000000007 .0000000008 1028. 0. MG: Richardson, r = 1, m1 = 3, m2 = 3, p = 2. .0000000009 .0000000010 348. 0. MG: Richardson, r = 2, m1 = 3, m2 = 3, p = 2. .0000000007 .0000000008 688. 0. MG: Richardson, r = 3, m1 = 3, m2 = 3, p = 2. .0000000009 .0000000010 1040. 2184. MG: Richardson, r = 4, m1 = 3, m2 = 3, p = 2. .0000000007 .0000000008 1384. 0. MG: Richardson, r = 1, m1 = 4, m2 = 4, p = 2. .0000000009 .0000000010 436. 0. MG: Richardson, r = 2, m1 = 4, m2 = 4, p = 2. .0000000007 .0000000008 864. 0. MG: Richardson, r = 3, m1 = 4, m2 = 4, p = 2. .0000000009 .0000000010 1300. 0. MG: Richardson, r = 4, m1 = 4, m2 = 4, p = 2. .0000000007 .0000000008 1768. 2156. MG: Jacobi, r = 1, m1 = 1, m2 = 1, p = 1. .0000000825 .0000002634 96. 0. MG: Jacobi, r = 2, m1 = 1, m2 = 1, p = 1. .0000000007 .0000000008 188. 0. MG: Jacobi, r = 3, m1 = 1, m2 = 1, p = 1. .0000000008 .0000000009 276. 0. MG: Jacobi, r = 4, m1 = 1, m2 = 1, p = 1. .0000000007 .0000000008 364. 0. MG: Jacobi, r = 1, m1 = 2, m2 = 2, p = 1. .0000000008 .0000000009 144. 0. MG: Jacobi, r = 2, m1 = 2, m2 = 2, p = 1. .0000000007 .0000000008 276. 0. MG: Jacobi, r = 3, m1 = 2, m2 = 2, p = 1. .0000000009 .0000000010 416. 0. MG: Jacobi, r = 4, m1 = 2, m2 = 2, p = 1. .0000000007 .0000000008 556. 0. MG: Jacobi, r = 1, m1 = 3, m2 = 3, p = 1. .0000000009 .0000000010 200. 2156. MG: Jacobi, r = 2, m1 = 3, m2 = 3, p = 1. .0000000007 .0000000008 368. 0. MG: Jacobi, r = 3, m1 = 3, m2 = 3, p = 1. .0000000009 .0000000010 552. 0. MG: Jacobi, r = 4, m1 = 3, m2 = 3, p = 1. .0000000007 .0000000008 736. 0. MG: Jacobi, r = 1, m1 = 4, m2 = 4, p = 1. .0000000009 .0000000010 236. 0. MG: Jacobi, r = 2, m1 = 4, m2 = 4, p = 1. .0000000007 .0000000008 464. 0. MG: Jacobi, r = 3, m1 = 4, m2 = 4, p = 1. .0000000009 .0000000010 696. 0. MG: Jacobi, r = 4, m1 = 4, m2 = 4, p = 1. .0000000007 .0000000008 924. 0. MG: Jacobi, r = 1, m1 = 1, m2 = 1, p = 2. .0000000008 .0000000047 172. 0. MG: Jacobi, r = 2, m1 = 1, m2 = 1, p = 2. .0000000007 .0000000008 316. 0. MG: Jacobi, r = 3, m1 = 1, m2 = 1, p = 2. .0000000009 .0000000010 472. 0. MG: Jacobi, r = 4, m1 = 1, m2 = 1, p = 2. .0000000007 .0000000008 628. 0. MG: Jacobi, r = 1, m1 = 2, m2 = 2, p = 2. .0000000009 .0000000010 236. 0. MG: Jacobi, r = 2, m1 = 2, m2 = 2, p = 2. .0000000007 .0000000008 464. 0. MG: Jacobi, r = 3, m1 = 2, m2 = 2, p = 2. .0000000009 .0000000010 696. 0. MG: Jacobi, r = 4, m1 = 2, m2 = 2, p = 2. .0000000007 .0000000008 932. 0. MG: Jacobi, r = 1, m1 = 3, m2 = 3, p = 2. .0000000009 .0000000010 320. 2192. MG: Jacobi, r = 2, m1 = 3, m2 = 3, p = 2. .0000000007 .0000000008 620. 0. MG: Jacobi, r = 3, m1 = 3, m2 = 3, p = 2. .0000000009 .0000000010 928. 0. MG: Jacobi, r = 4, m1 = 3, m2 = 3, p = 2. .0000000007 .0000000008 1232. 0. MG: Jacobi, r = 1, m1 = 4, m2 = 4, p = 2. .0000000009 .0000000010 392. 0. MG: Jacobi, r = 2, m1 = 4, m2 = 4, p = 2. .0000000007 .0000000008 776. 0. MG: Jacobi, r = 3, m1 = 4, m2 = 4, p = 2. .0000000009 .0000000010 1180. 2176. MG: Jacobi, r = 4, m1 = 4, m2 = 4, p = 2. .0000000007 .0000000008 1544. 0. MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 1. .0000000009 .0000000010 160. 0. MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 1. .0000000007 .0000000008 320. 0. MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 1. .0000000009 .0000000010 476. 0. MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 1. .0000000007 .0000000008 632. 0. MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 1. .0000000009 .0000000010 276. 0. MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 1. .0000000007 .0000000008 544. 0. MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 1. .0000000009 .0000000010 812. 0. MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 1. .0000000007 .0000000008 1076. 2140. MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 1. .0000000009 .0000000010 384. 0. MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 1. .0000000007 .0000000008 764. 0. MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 1. .0000000009 .0000000010 1148. 0. MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 1. .0000000007 .0000000008 1524. 0. MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 1. .0000000009 .0000000010 500. 0. MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 1. .0000000007 .0000000008 992. 0. MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 1. .0000000009 .0000000010 1484. 0. MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 1. .0000000007 .0000000008 1976. 0. MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 2. .0000000009 .0000000010 248. 0. MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 2. .0000000007 .0000000008 500. 0. MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 2. .0000000009 .0000000010 740. 2168. MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 2. .0000000007 .0000000008 984. 0. MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 2. .0000000009 .0000000010 416. 0. MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 2. .0000000007 .0000000008 832. 0. MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 2. .0000000009 .0000000010 1240. 0. MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 2. .0000000007 .0000000008 1660. 0. MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 2. .0000000009 .0000000010 588. 0. MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 2. .0000000007 .0000000008 1164. 0. MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 2. .0000000009 .0000000010 1744. 2168. MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 2. .0000000007 .0000000008 2324. 0. MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 2. .0000000009 .0000000010 752. 0. MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 2. .0000000007 .0000000008 1504. 0. MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 2. .0000000009 .0000000010 2248. 0. MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 2. .0000000007 .0000000008 2996. 0. (time (build-operator (instance-ref g1229 1) (instance-ref g1229 2) (instance-ref g1229 4) fes)) 7188 ms real time 7184 ms cpu time (7184 user, 0 system) 1 collection accounting for 2228 ms real time (2228 user, 0 system) 7160428816 bytes allocated 4 minor faults no major faults (time (build-rhs (instance-ref g1229 5) (instance-ref g1229 6) fes)) 3606 ms real time 3604 ms cpu time (3596 user, 8 system) 1 collection accounting for 2274 ms real time (2264 user, 8 system) 1788559752 bytes allocated no minor faults no major faults a=1, c=1, u=x, Neumann (error = 0 with exact linear algebra) MG: Richardson, r = 1, m1 = 1, m2 = 1, p = 1. .0228284973 .0232123035 104. 0. MG: Richardson, r = 2, m1 = 1, m2 = 1, p = 1. .0010829797 .0011038693 196. 0. MG: Richardson, r = 3, m1 = 1, m2 = 1, p = 1. .0000515504 .0000525744 296. 0. MG: Richardson, r = 4, m1 = 1, m2 = 1, p = 1. .0000024558 .0000025050 392. 0. MG: Richardson, r = 1, m1 = 2, m2 = 2, p = 1. .0011012850 .0011199854 156. 0. MG: Richardson, r = 2, m1 = 2, m2 = 2, p = 1. .0000025319 .0000025781 304. 0. MG: Richardson, r = 3, m1 = 2, m2 = 2, p = 1. .0000000047 .0000000049 452. 0. MG: Richardson, r = 4, m1 = 2, m2 = 2, p = 1. .0000000010 .0000000010 604. 0. MG: Richardson, r = 1, m1 = 3, m2 = 3, p = 1. .0000540970 .0000550186 208. 0. MG: Richardson, r = 2, m1 = 3, m2 = 3, p = 1. .0000000050 .0000000052 416. 0. MG: Richardson, r = 3, m1 = 3, m2 = 3, p = 1. .0000000010 .0000000010 644. 2128. MG: Richardson, r = 4, m1 = 3, m2 = 3, p = 1. .0000000010 .0000000010 812. 0. MG: Richardson, r = 1, m1 = 4, m2 = 4, p = 1. .0000027100 .0000027599 260. 0. MG: Richardson, r = 2, m1 = 4, m2 = 4, p = 1. .0000000010 .0000000010 516. 0. MG: Richardson, r = 3, m1 = 4, m2 = 4, p = 1. .0000000010 .0000000010 772. 0. MG: Richardson, r = 4, m1 = 4, m2 = 4, p = 1. .0000000010 .0000000010 1024. 0. MG: Richardson, r = 1, m1 = 1, m2 = 1, p = 2. .0000000010 .0000000029 172. 0. MG: Richardson, r = 2, m1 = 1, m2 = 1, p = 2. .0000000010 .0000000010 340. 0. MG: Richardson, r = 3, m1 = 1, m2 = 1, p = 2. .0000000010 .0000000010 524. 0. MG: Richardson, r = 4, m1 = 1, m2 = 1, p = 2. .0000000010 .0000000010 672. 0. MG: Richardson, r = 1, m1 = 2, m2 = 2, p = 2. .0000000010 .0000000010 264. 0. MG: Richardson, r = 2, m1 = 2, m2 = 2, p = 2. .0000000010 .0000000010 516. 0. MG: Richardson, r = 3, m1 = 2, m2 = 2, p = 2. .0000000010 .0000000010 772. 0. MG: Richardson, r = 4, m1 = 2, m2 = 2, p = 2. .0000000010 .0000000010 1024. 0. MG: Richardson, r = 1, m1 = 3, m2 = 3, p = 2. .0000000010 .0000000010 348. 0. MG: Richardson, r = 2, m1 = 3, m2 = 3, p = 2. .0000000010 .0000000010 688. 0. MG: Richardson, r = 3, m1 = 3, m2 = 3, p = 2. .0000000010 .0000000010 1040. 0. MG: Richardson, r = 4, m1 = 3, m2 = 3, p = 2. .0000000010 .0000000010 1388. 0. MG: Richardson, r = 1, m1 = 4, m2 = 4, p = 2. .0000000010 .0000000010 436. 0. MG: Richardson, r = 2, m1 = 4, m2 = 4, p = 2. .0000000010 .0000000010 868. 0. MG: Richardson, r = 3, m1 = 4, m2 = 4, p = 2. .0000000010 .0000000010 1320. 0. MG: Richardson, r = 4, m1 = 4, m2 = 4, p = 2. .0000000010 .0000000010 1748. 2160. MG: Jacobi, r = 1, m1 = 1, m2 = 1, p = 1. .0004352011 .0004356213 96. 0. MG: Jacobi, r = 2, m1 = 1, m2 = 1, p = 1. .0000003786 .0000003790 184. 0. MG: Jacobi, r = 3, m1 = 1, m2 = 1, p = 1. .0000000007 .0000000007 276. 0. MG: Jacobi, r = 4, m1 = 1, m2 = 1, p = 1. .0000000010 .0000000010 368. 0. MG: Jacobi, r = 1, m1 = 2, m2 = 2, p = 1. .0000003818 .0000003820 140. 0. MG: Jacobi, r = 2, m1 = 2, m2 = 2, p = 1. .0000000010 .0000000010 280. 0. MG: Jacobi, r = 3, m1 = 2, m2 = 2, p = 1. .0000000010 .0000000010 416. 0. MG: Jacobi, r = 4, m1 = 2, m2 = 2, p = 1. .0000000010 .0000000010 564. 0. MG: Jacobi, r = 1, m1 = 3, m2 = 3, p = 1. .0000000007 .0000000007 188. 0. MG: Jacobi, r = 2, m1 = 3, m2 = 3, p = 1. .0000000010 .0000000010 372. 0. MG: Jacobi, r = 3, m1 = 3, m2 = 3, p = 1. .0000000010 .0000000010 552. 0. MG: Jacobi, r = 4, m1 = 3, m2 = 3, p = 1. .0000000010 .0000000010 740. 0. MG: Jacobi, r = 1, m1 = 4, m2 = 4, p = 1. .0000000010 .0000000010 236. 0. MG: Jacobi, r = 2, m1 = 4, m2 = 4, p = 1. .0000000010 .0000000010 468. 0. MG: Jacobi, r = 3, m1 = 4, m2 = 4, p = 1. .0000000010 .0000000010 696. 0. MG: Jacobi, r = 4, m1 = 4, m2 = 4, p = 1. .0000000010 .0000000010 936. 0. MG: Jacobi, r = 1, m1 = 1, m2 = 1, p = 2. .0000000010 .0000000081 176. 2108. MG: Jacobi, r = 2, m1 = 1, m2 = 1, p = 2. .0000000010 .0000000010 320. 0. MG: Jacobi, r = 3, m1 = 1, m2 = 1, p = 2. .0000000010 .0000000010 476. 0. MG: Jacobi, r = 4, m1 = 1, m2 = 1, p = 2. .0000000010 .0000000010 620. 0. MG: Jacobi, r = 1, m1 = 2, m2 = 2, p = 2. .0000000010 .0000000010 236. 0. MG: Jacobi, r = 2, m1 = 2, m2 = 2, p = 2. .0000000010 .0000000010 472. 0. MG: Jacobi, r = 3, m1 = 2, m2 = 2, p = 2. .0000000010 .0000000010 696. 0. MG: Jacobi, r = 4, m1 = 2, m2 = 2, p = 2. .0000000010 .0000000010 928. 0. MG: Jacobi, r = 1, m1 = 3, m2 = 3, p = 2. .0000000010 .0000000010 316. 0. MG: Jacobi, r = 2, m1 = 3, m2 = 3, p = 2. .0000000010 .0000000010 624. 0. MG: Jacobi, r = 3, m1 = 3, m2 = 3, p = 2. .0000000010 .0000000010 936. 0. MG: Jacobi, r = 4, m1 = 3, m2 = 3, p = 2. .0000000010 .0000000010 1244. 0. MG: Jacobi, r = 1, m1 = 4, m2 = 4, p = 2. .0000000010 .0000000010 396. 0. MG: Jacobi, r = 2, m1 = 4, m2 = 4, p = 2. .0000000010 .0000000010 796. 0. MG: Jacobi, r = 3, m1 = 4, m2 = 4, p = 2. .0000000010 .0000000010 1184. 2152. MG: Jacobi, r = 4, m1 = 4, m2 = 4, p = 2. .0000000010 .0000000010 1552. 0. MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 1. .0000000013 .0000000014 164. 0. MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 1. .0000000010 .0000000010 324. 0. MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 1. .0000000010 .0000000010 476. 0. MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 1. .0000000010 .0000000010 628. 0. MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 1. .0000000010 .0000000010 276. 0. MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 1. .0000000010 .0000000010 544. 0. MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 1. .0000000010 .0000000010 812. 0. MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 1. .0000000010 .0000000010 1084. 0. MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 1. .0000000010 .0000000010 388. 0. MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 1. .0000000010 .0000000010 764. 0. MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 1. .0000000010 .0000000010 1144. 0. MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 1. .0000000010 .0000000010 1516. 0. MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 1. .0000000010 .0000000010 492. 0. MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 1. .0000000010 .0000000010 984. 0. MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 1. .0000000010 .0000000010 1476. 0. MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 1. .0000000010 .0000000010 1976. 0. MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 2. .0000000010 .0000000010 252. 0. MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 2. .0000000010 .0000000010 496. 0. MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 2. .0000000010 .0000000010 740. 0. MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 2. .0000000010 .0000000010 992. 0. MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 2. .0000000010 .0000000010 420. 0. MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 2. .0000000010 .0000000010 832. 0. MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 2. .0000000010 .0000000010 1248. 0. MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 2. .0000000010 .0000000010 1660. 0. MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 2. .0000000010 .0000000010 584. 0. MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 2. .0000000010 .0000000010 1164. 0. MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 2. .0000000010 .0000000010 1744. 2172. MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 2. .0000000010 .0000000010 2328. 0. MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 2. .0000000010 .0000000010 752. 0. MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 2. .0000000010 .0000000010 1500. 0. MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 2. .0000000010 .0000000010 2252. 0. MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 2. .0000000010 .0000000010 2992. 0. (time (build-operator (instance-ref g1229 1) (instance-ref g1229 2) (instance-ref g1229 4) fes)) 9374 ms real time 9372 ms cpu time (9368 user, 4 system) 2 collections accounting for 4427 ms real time (4420 user, 4 system) 7361755408 bytes allocated 2 minor faults no major faults (time (build-rhs (instance-ref g1229 5) (instance-ref g1229 6) fes)) 1318 ms real time 1316 ms cpu time (1316 user, 0 system) no collections 1587626376 bytes allocated no minor faults no major faults 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. .0856470377 .0868064770 104. 0. MG: Richardson, r = 2, m1 = 1, m2 = 1, p = 1. .0075870823 .0077350180 200. 0. MG: Richardson, r = 3, m1 = 1, m2 = 1, p = 1. .0006756698 .0006917431 296. 0. MG: Richardson, r = 4, m1 = 1, m2 = 1, p = 1. .0000603314 .0000619643 396. 0. MG: Richardson, r = 1, m1 = 2, m2 = 2, p = 1. .0076484846 .0077850532 156. 0. MG: Richardson, r = 2, m1 = 2, m2 = 2, p = 1. .0000611392 .0000627868 304. 0. MG: Richardson, r = 3, m1 = 2, m2 = 2, p = 1. .0000004911 .0000005078 452. 0. MG: Richardson, r = 4, m1 = 2, m2 = 2, p = 1. .0000000026 .0000000028 604. 0. MG: Richardson, r = 1, m1 = 3, m2 = 3, p = 1. .0006907017 .0007063809 208. 0. MG: Richardson, r = 2, m1 = 3, m2 = 3, p = 1. .0000005020 .0000005217 420. 0. MG: Richardson, r = 3, m1 = 3, m2 = 3, p = 1. .0000000009 .0000000009 624. 2136. MG: Richardson, r = 4, m1 = 3, m2 = 3, p = 1. .0000000013 .0000000013 820. 0. MG: Richardson, r = 1, m1 = 4, m2 = 4, p = 1. .0000630526 .0000650100 260. 0. MG: Richardson, r = 2, m1 = 4, m2 = 4, p = 1. .0000000029 .0000000033 520. 0. MG: Richardson, r = 3, m1 = 4, m2 = 4, p = 1. .0000000013 .0000000013 776. 0. MG: Richardson, r = 4, m1 = 4, m2 = 4, p = 1. .0000000013 .0000000013 1032. 0. MG: Richardson, r = 1, m1 = 1, m2 = 1, p = 2. .0000000013 .0000000809 172. 0. MG: Richardson, r = 2, m1 = 1, m2 = 1, p = 2. .0000000013 .0000000014 364. 0. MG: Richardson, r = 3, m1 = 1, m2 = 1, p = 2. .0000000013 .0000000013 520. 2132. MG: Richardson, r = 4, m1 = 1, m2 = 1, p = 2. .0000000013 .0000000013 672. 0. MG: Richardson, r = 1, m1 = 2, m2 = 2, p = 2. .0000000013 .0000000023 264. 0. MG: Richardson, r = 2, m1 = 2, m2 = 2, p = 2. .0000000013 .0000000013 516. 0. MG: Richardson, r = 3, m1 = 2, m2 = 2, p = 2. .0000000013 .0000000013 772. 0. MG: Richardson, r = 4, m1 = 2, m2 = 2, p = 2. .0000000013 .0000000013 1028. 0. MG: Richardson, r = 1, m1 = 3, m2 = 3, p = 2. .0000000013 .0000000013 348. 0. MG: Richardson, r = 2, m1 = 3, m2 = 3, p = 2. .0000000013 .0000000013 696. 2184. MG: Richardson, r = 3, m1 = 3, m2 = 3, p = 2. .0000000013 .0000000013 1032. 0. MG: Richardson, r = 4, m1 = 3, m2 = 3, p = 2. .0000000013 .0000000013 1380. 0. MG: Richardson, r = 1, m1 = 4, m2 = 4, p = 2. .0000000013 .0000000013 440. 0. MG: Richardson, r = 2, m1 = 4, m2 = 4, p = 2. .0000000013 .0000000013 868. 0. MG: Richardson, r = 3, m1 = 4, m2 = 4, p = 2. .0000000013 .0000000013 1304. 0. MG: Richardson, r = 4, m1 = 4, m2 = 4, p = 2. .0000000013 .0000000013 1732. 0. MG: Jacobi, r = 1, m1 = 1, m2 = 1, p = 1. .0035026489 .0035047628 96. 0. MG: Jacobi, r = 2, m1 = 1, m2 = 1, p = 1. .0000122852 .0000122926 188. 0. MG: Jacobi, r = 3, m1 = 1, m2 = 1, p = 1. .0000000417 .0000000417 280. 0. MG: Jacobi, r = 4, m1 = 1, m2 = 1, p = 1. .0000000011 .0000000012 368. 0. MG: Jacobi, r = 1, m1 = 2, m2 = 2, p = 1. .0000125223 .0000125335 140. 0. MG: Jacobi, r = 2, m1 = 2, m2 = 2, p = 1. .0000000011 .0000000011 276. 0. MG: Jacobi, r = 3, m1 = 2, m2 = 2, p = 1. .0000000013 .0000000013 416. 0. MG: Jacobi, r = 4, m1 = 2, m2 = 2, p = 1. .0000000013 .0000000013 552. 2180. MG: Jacobi, r = 1, m1 = 3, m2 = 3, p = 1. .0000000437 .0000000437 192. 0. MG: Jacobi, r = 2, m1 = 3, m2 = 3, p = 1. .0000000013 .0000000013 368. 0. MG: Jacobi, r = 3, m1 = 3, m2 = 3, p = 1. .0000000013 .0000000013 556. 0. MG: Jacobi, r = 4, m1 = 3, m2 = 3, p = 1. .0000000013 .0000000013 736. 0. MG: Jacobi, r = 1, m1 = 4, m2 = 4, p = 1. .0000000011 .0000000011 236. 0. MG: Jacobi, r = 2, m1 = 4, m2 = 4, p = 1. .0000000013 .0000000013 460. 0. MG: Jacobi, r = 3, m1 = 4, m2 = 4, p = 1. .0000000013 .0000000013 692. 0. MG: Jacobi, r = 4, m1 = 4, m2 = 4, p = 1. .0000000013 .0000000013 936. 0. MG: Jacobi, r = 1, m1 = 1, m2 = 1, p = 2. .0000000013 .0000000130 160. 0. MG: Jacobi, r = 2, m1 = 1, m2 = 1, p = 2. .0000000013 .0000000013 316. 0. MG: Jacobi, r = 3, m1 = 1, m2 = 1, p = 2. .0000000013 .0000000013 472. 0. MG: Jacobi, r = 4, m1 = 1, m2 = 1, p = 2. .0000000013 .0000000013 624. 0. MG: Jacobi, r = 1, m1 = 2, m2 = 2, p = 2. .0000000013 .0000000013 236. 0. MG: Jacobi, r = 2, m1 = 2, m2 = 2, p = 2. .0000000013 .0000000013 468. 0. MG: Jacobi, r = 3, m1 = 2, m2 = 2, p = 2. .0000000013 .0000000013 700. 0. MG: Jacobi, r = 4, m1 = 2, m2 = 2, p = 2. .0000000013 .0000000013 932. 2168. MG: Jacobi, r = 1, m1 = 3, m2 = 3, p = 2. .0000000013 .0000000013 316. 0. MG: Jacobi, r = 2, m1 = 3, m2 = 3, p = 2. .0000000013 .0000000013 620. 0. MG: Jacobi, r = 3, m1 = 3, m2 = 3, p = 2. .0000000013 .0000000013 928. 0. MG: Jacobi, r = 4, m1 = 3, m2 = 3, p = 2. .0000000013 .0000000013 1236. 0. MG: Jacobi, r = 1, m1 = 4, m2 = 4, p = 2. .0000000013 .0000000013 388. 0. MG: Jacobi, r = 2, m1 = 4, m2 = 4, p = 2. .0000000013 .0000000013 776. 0. MG: Jacobi, r = 3, m1 = 4, m2 = 4, p = 2. .0000000013 .0000000013 1160. 0. MG: Jacobi, r = 4, m1 = 4, m2 = 4, p = 2. .0000000013 .0000000013 1552. 0. MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 1. .0000001614 .0000001620 160. 0. MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 1. .0000000013 .0000000013 320. 0. MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 1. .0000000013 .0000000013 472. 0. MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 1. .0000000013 .0000000013 632. 0. MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 1. .0000000014 .0000000014 272. 0. MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 1. .0000000014 .0000000014 544. 0. MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 1. .0000000013 .0000000014 812. 2120. MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 1. .0000000013 .0000000014 1076. 0. MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 1. .0000000014 .0000000014 392. 0. MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 1. .0000000014 .0000000014 768. 0. MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 1. .0000000014 .0000000014 1148. 0. MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 1. .0000000014 .0000000014 1512. 0. MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 1. .0000000014 .0000000014 496. 0. MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 1. .0000000014 .0000000014 988. 0. MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 1. .0000000014 .0000000014 1480. 0. MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 1. .0000000014 .0000000014 1984. 0. MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 2. .0000000013 .0000000013 256. 0. MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 2. .0000000013 .0000000013 496. 2140. MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 2. .0000000013 .0000000013 744. 0. MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 2. .0000000013 .0000000013 988. 0. MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 2. .0000000014 .0000000014 420. 0. MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 2. .0000000014 .0000000014 832. 0. MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 2. .0000000013 .0000000014 1248. 0. MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 2. .0000000013 .0000000013 1660. 0. MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 2. .0000000014 .0000000014 588. 0. MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 2. .0000000014 .0000000014 1168. 0. MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 2. .0000000014 .0000000014 1740. 0. MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 2. .0000000014 .0000000014 2328. 0. MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 2. .0000000014 .0000000014 756. 0. MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 2. .0000000014 .0000000014 1504. 0. MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 2. .0000000014 .0000000014 2248. 0. MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 2. .0000000014 .0000000014 2996. 0. (time (build-operator (instance-ref g1229 1) (instance-ref g1229 2) (instance-ref g1229 4) fes)) 4990 ms real time 4988 ms cpu time (4988 user, 0 system) no collections 6254459152 bytes allocated no minor faults no major faults (time (build-rhs (instance-ref g1229 5) (instance-ref g1229 6) fes)) 3574 ms real time 3572 ms cpu time (3572 user, 0 system) 1 collection accounting for 2225 ms real time (2224 user, 0 system) 1788821896 bytes allocated no minor faults no major faults a=1, c=1, u=x^2-y^2, Neumann MG: Richardson, r = 1, m1 = 1, m2 = 1, p = 1. .0000528776 .0009273677 104. 0. MG: Richardson, r = 2, m1 = 1, m2 = 1, p = 1. .0000001346 .0006961121 200. 0. MG: Richardson, r = 3, m1 = 1, m2 = 1, p = 1. .0000000977 .0006909053 296. 0. MG: Richardson, r = 4, m1 = 1, m2 = 1, p = 1. .0000000973 .0006905616 396. 0. MG: Richardson, r = 1, m1 = 2, m2 = 2, p = 1. .0000002381 .0006961370 156. 0. MG: Richardson, r = 2, m1 = 2, m2 = 2, p = 1. .0000000973 .0006905616 308. 0. MG: Richardson, r = 3, m1 = 2, m2 = 2, p = 1. .0000000973 .0006905343 460. 0. MG: Richardson, r = 4, m1 = 2, m2 = 2, p = 1. .0000000973 .0006905339 604. 0. MG: Richardson, r = 1, m1 = 3, m2 = 3, p = 1. .0000000980 .0006909066 208. 0. MG: Richardson, r = 2, m1 = 3, m2 = 3, p = 1. .0000000973 .0006905343 416. 0. MG: Richardson, r = 3, m1 = 3, m2 = 3, p = 1. .0000000973 .0006905339 620. 2148. MG: Richardson, r = 4, m1 = 3, m2 = 3, p = 1. .0000000973 .0006905339 820. 0. MG: Richardson, r = 1, m1 = 4, m2 = 4, p = 1. .0000000974 .0006905618 264. 0. MG: Richardson, r = 2, m1 = 4, m2 = 4, p = 1. .0000000973 .0006905339 520. 0. MG: Richardson, r = 3, m1 = 4, m2 = 4, p = 1. .0000000973 .0006905339 780. 0. MG: Richardson, r = 4, m1 = 4, m2 = 4, p = 1. .0000000973 .0006905339 1036. 0. MG: Richardson, r = 1, m1 = 1, m2 = 1, p = 2. .0000001395 .0007732508 172. 0. MG: Richardson, r = 2, m1 = 1, m2 = 1, p = 2. .0000001004 .0006960696 336. 0. MG: Richardson, r = 3, m1 = 1, m2 = 1, p = 2. .0000000975 .0006909036 516. 0. MG: Richardson, r = 4, m1 = 1, m2 = 1, p = 2. .0000000973 .0006905614 672. 0. MG: Richardson, r = 1, m1 = 2, m2 = 2, p = 2. .0000001004 .0006960781 264. 0. MG: Richardson, r = 2, m1 = 2, m2 = 2, p = 2. .0000000973 .0006905614 520. 0. MG: Richardson, r = 3, m1 = 2, m2 = 2, p = 2. .0000000973 .0006905343 776. 0. MG: Richardson, r = 4, m1 = 2, m2 = 2, p = 2. .0000000973 .0006905339 1036. 0. MG: Richardson, r = 1, m1 = 3, m2 = 3, p = 2. .0000000975 .0006909049 352. 0. MG: Richardson, r = 2, m1 = 3, m2 = 3, p = 2. .0000000973 .0006905343 692. 0. MG: Richardson, r = 3, m1 = 3, m2 = 3, p = 2. .0000000973 .0006905339 1036. 0. MG: Richardson, r = 4, m1 = 3, m2 = 3, p = 2. .0000000973 .0006905339 1380. 0. MG: Richardson, r = 1, m1 = 4, m2 = 4, p = 2. .0000000973 .0006905616 440. 0. MG: Richardson, r = 2, m1 = 4, m2 = 4, p = 2. .0000000973 .0006905339 868. 0. MG: Richardson, r = 3, m1 = 4, m2 = 4, p = 2. .0000000973 .0006905339 1296. 0. MG: Richardson, r = 4, m1 = 4, m2 = 4, p = 2. .0000000973 .0006905339 1732. 2164. MG: Jacobi, r = 1, m1 = 1, m2 = 1, p = 1. .0000000973 .0006905339 96. 0. MG: Jacobi, r = 2, m1 = 1, m2 = 1, p = 1. .0000000973 .0006905339 184. 0. MG: Jacobi, r = 3, m1 = 1, m2 = 1, p = 1. .0000000973 .0006905339 276. 0. MG: Jacobi, r = 4, m1 = 1, m2 = 1, p = 1. .0000000973 .0006905339 368. 0. MG: Jacobi, r = 1, m1 = 2, m2 = 2, p = 1. .0000000973 .0006905339 140. 0. MG: Jacobi, r = 2, m1 = 2, m2 = 2, p = 1. .0000000973 .0006905339 276. 0. MG: Jacobi, r = 3, m1 = 2, m2 = 2, p = 1. .0000000973 .0006905339 416. 0. MG: Jacobi, r = 4, m1 = 2, m2 = 2, p = 1. .0000000973 .0006905339 552. 0. MG: Jacobi, r = 1, m1 = 3, m2 = 3, p = 1. .0000000973 .0006905339 192. 0. MG: Jacobi, r = 2, m1 = 3, m2 = 3, p = 1. .0000000973 .0006905339 372. 0. MG: Jacobi, r = 3, m1 = 3, m2 = 3, p = 1. .0000000973 .0006905339 556. 0. MG: Jacobi, r = 4, m1 = 3, m2 = 3, p = 1. .0000000973 .0006905339 740. 0. MG: Jacobi, r = 1, m1 = 4, m2 = 4, p = 1. .0000000973 .0006905339 232. 0. MG: Jacobi, r = 2, m1 = 4, m2 = 4, p = 1. .0000000973 .0006905339 460. 0. MG: Jacobi, r = 3, m1 = 4, m2 = 4, p = 1. .0000000973 .0006905339 692. 0. MG: Jacobi, r = 4, m1 = 4, m2 = 4, p = 1. .0000000973 .0006905339 928. 0. MG: Jacobi, r = 1, m1 = 1, m2 = 1, p = 2. .0000000973 .0006905339 172. 2172. MG: Jacobi, r = 2, m1 = 1, m2 = 1, p = 2. .0000000973 .0006905339 316. 0. MG: Jacobi, r = 3, m1 = 1, m2 = 1, p = 2. .0000000973 .0006905339 472. 0. MG: Jacobi, r = 4, m1 = 1, m2 = 1, p = 2. .0000000973 .0006905339 632. 0. MG: Jacobi, r = 1, m1 = 2, m2 = 2, p = 2. .0000000973 .0006905339 236. 0. MG: Jacobi, r = 2, m1 = 2, m2 = 2, p = 2. .0000000973 .0006905339 472. 0. MG: Jacobi, r = 3, m1 = 2, m2 = 2, p = 2. .0000000973 .0006905339 696. 0. MG: Jacobi, r = 4, m1 = 2, m2 = 2, p = 2. .0000000973 .0006905339 932. 0. MG: Jacobi, r = 1, m1 = 3, m2 = 3, p = 2. .0000000973 .0006905339 312. 0. MG: Jacobi, r = 2, m1 = 3, m2 = 3, p = 2. .0000000973 .0006905339 620. 0. MG: Jacobi, r = 3, m1 = 3, m2 = 3, p = 2. .0000000973 .0006905339 932. 0. MG: Jacobi, r = 4, m1 = 3, m2 = 3, p = 2. .0000000973 .0006905339 1240. 0. MG: Jacobi, r = 1, m1 = 4, m2 = 4, p = 2. .0000000973 .0006905339 392. 0. MG: Jacobi, r = 2, m1 = 4, m2 = 4, p = 2. .0000000973 .0006905339 776. 0. MG: Jacobi, r = 3, m1 = 4, m2 = 4, p = 2. .0000000973 .0006905339 1164. 2160. MG: Jacobi, r = 4, m1 = 4, m2 = 4, p = 2. .0000000973 .0006905339 1544. 0. MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 1. .0000000978 .0006906556 164. 0. MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 1. .0000000973 .0006905343 320. 0. MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 1. .0000000973 .0006905339 476. 0. MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 1. .0000000973 .0006905339 632. 0. MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 1. .0000000973 .0006905355 276. 0. MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 1. .0000000973 .0006905339 544. 0. MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 1. .0000000973 .0006905339 808. 0. MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 1. .0000000973 .0006905339 1080. 0. MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 1. .0000000973 .0006905341 384. 0. MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 1. .0000000973 .0006905339 764. 0. MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 1. .0000000973 .0006905339 1160. 0. MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 1. .0000000973 .0006905339 1528. 0. MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 1. .0000000973 .0006905340 500. 0. MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 1. .0000000973 .0006905339 988. 0. MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 1. .0000000973 .0006905339 1488. 0. MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 1. .0000000973 .0006905339 1972. 0. MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 2. .0000000977 .0006906520 252. 0. MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 2. .0000000973 .0006905343 500. 0. MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 2. .0000000973 .0006905339 740. 0. MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 2. .0000000973 .0006905339 988. 0. MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 2. .0000000973 .0006905355 420. 0. MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 2. .0000000973 .0006905339 832. 0. MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 2. .0000000973 .0006905339 1248. 0. MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 2. .0000000973 .0006905339 1664. 0. MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 2. .0000000973 .0006905341 588. 0. MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 2. .0000000973 .0006905339 1164. 0. MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 2. .0000000973 .0006905339 1744. 2152. MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 2. .0000000973 .0006905339 2328. 0. MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 2. .0000000973 .0006905340 760. 0. MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 2. .0000000973 .0006905339 1500. 0. MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 2. .0000000973 .0006905339 2248. 0. MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 2. .0000000973 .0006905339 3000. 0. (time (build-operator (instance-ref g1229 1) (instance-ref g1229 2) (instance-ref g1229 4) fes)) 7229 ms real time 7228 ms cpu time (7216 user, 12 system) 1 collection accounting for 2203 ms real time (2192 user, 12 system) 6254459152 bytes allocated no minor faults no major faults (time (build-rhs (instance-ref g1229 5) (instance-ref g1229 6) fes)) 3712 ms real time 3708 ms cpu time (3708 user, 0 system) 1 collection accounting for 2275 ms real time (2276 user, 0 system) 1786790216 bytes allocated no minor faults no major faults a=1, c=1, u=\cos(\pi x), Neumann MG: Richardson, r = 1, m1 = 1, m2 = 1, p = 1. .0010386820 .0081187908 104. 0. MG: Richardson, r = 2, m1 = 1, m2 = 1, p = 1. .0000235623 .0017302805 200. 0. MG: Richardson, r = 3, m1 = 1, m2 = 1, p = 1. .0000018507 .0017048451 296. 0. MG: Richardson, r = 4, m1 = 1, m2 = 1, p = 1. .0000007508 .0017039014 392. 0. MG: Richardson, r = 1, m1 = 2, m2 = 2, p = 1. .0000375563 .0017468472 156. 0. MG: Richardson, r = 2, m1 = 2, m2 = 2, p = 1. .0000009112 .0017039075 304. 0. MG: Richardson, r = 3, m1 = 2, m2 = 2, p = 1. .0000005498 .0017038243 456. 0. MG: Richardson, r = 4, m1 = 2, m2 = 2, p = 1. .0000005000 .0017038225 604. 0. MG: Richardson, r = 1, m1 = 3, m2 = 3, p = 1. .0000049978 .0017053328 212. 0. MG: Richardson, r = 2, m1 = 3, m2 = 3, p = 1. .0000005923 .0017038249 416. 0. MG: Richardson, r = 3, m1 = 3, m2 = 3, p = 1. .0000004991 .0017038224 648. 0. MG: Richardson, r = 4, m1 = 3, m2 = 3, p = 1. .0000004833 .0017038223 816. 0. MG: Richardson, r = 1, m1 = 4, m2 = 4, p = 1. .0000018163 .0017039597 264. 0. MG: Richardson, r = 2, m1 = 4, m2 = 4, p = 1. .0000005325 .0017038228 520. 0. MG: Richardson, r = 3, m1 = 4, m2 = 4, p = 1. .0000004870 .0017038223 776. 0. MG: Richardson, r = 4, m1 = 4, m2 = 4, p = 1. .0000004797 .0017038222 1032. 0. MG: Richardson, r = 1, m1 = 1, m2 = 1, p = 2. .0000014225 .0019088218 172. 0. MG: Richardson, r = 2, m1 = 1, m2 = 1, p = 2. .0000006651 .0017175962 340. 0. MG: Richardson, r = 3, m1 = 1, m2 = 1, p = 2. .0000005226 .0017047513 512. 0. MG: Richardson, r = 4, m1 = 1, m2 = 1, p = 2. .0000004889 .0017038930 680. 2124. MG: Richardson, r = 1, m1 = 2, m2 = 2, p = 2. .0000008580 .0017176338 260. 0. MG: Richardson, r = 2, m1 = 2, m2 = 2, p = 2. .0000005000 .0017038932 516. 0. MG: Richardson, r = 3, m1 = 2, m2 = 2, p = 2. .0000004791 .0017038232 776. 0. MG: Richardson, r = 4, m1 = 2, m2 = 2, p = 2. .0000004778 .0017038223 1032. 0. MG: Richardson, r = 1, m1 = 3, m2 = 3, p = 2. .0000006608 .0017047586 348. 0. MG: Richardson, r = 2, m1 = 3, m2 = 3, p = 2. .0000004805 .0017038232 696. 0. MG: Richardson, r = 3, m1 = 3, m2 = 3, p = 2. .0000004778 .0017038222 1036. 2152. MG: Richardson, r = 4, m1 = 3, m2 = 3, p = 2. .0000004777 .0017038222 1376. 0. MG: Richardson, r = 1, m1 = 4, m2 = 4, p = 2. .0000005676 .0017038948 436. 0. MG: Richardson, r = 2, m1 = 4, m2 = 4, p = 2. .0000004781 .0017038223 864. 0. MG: Richardson, r = 3, m1 = 4, m2 = 4, p = 2. .0000004777 .0017038222 1300. 0. MG: Richardson, r = 4, m1 = 4, m2 = 4, p = 2. .0000004777 .0017038222 1776. 2160. MG: Jacobi, r = 1, m1 = 1, m2 = 1, p = 1. .0000023565 .0017039610 96. 0. MG: Jacobi, r = 2, m1 = 1, m2 = 1, p = 1. .0000005368 .0017038246 188. 0. MG: Jacobi, r = 3, m1 = 1, m2 = 1, p = 1. .0000004881 .0017038236 280. 0. MG: Jacobi, r = 4, m1 = 1, m2 = 1, p = 1. .0000004801 .0017038233 368. 0. MG: Jacobi, r = 1, m1 = 2, m2 = 2, p = 1. .0000008335 .0017038256 140. 0. MG: Jacobi, r = 2, m1 = 2, m2 = 2, p = 1. .0000004992 .0017038224 280. 0. MG: Jacobi, r = 3, m1 = 2, m2 = 2, p = 1. .0000004804 .0017038223 416. 0. MG: Jacobi, r = 4, m1 = 2, m2 = 2, p = 1. .0000004781 .0017038223 560. 0. MG: Jacobi, r = 1, m1 = 3, m2 = 3, p = 1. .0000006777 .0017038248 192. 2136. MG: Jacobi, r = 2, m1 = 3, m2 = 3, p = 1. .0000004890 .0017038230 372. 0. MG: Jacobi, r = 3, m1 = 3, m2 = 3, p = 1. .0000004788 .0017038228 556. 0. MG: Jacobi, r = 4, m1 = 3, m2 = 3, p = 1. .0000004778 .0017038227 740. 0. MG: Jacobi, r = 1, m1 = 4, m2 = 4, p = 1. .0000006144 .0017038231 236. 0. MG: Jacobi, r = 2, m1 = 4, m2 = 4, p = 1. .0000004846 .0017038223 464. 0. MG: Jacobi, r = 3, m1 = 4, m2 = 4, p = 1. .0000004782 .0017038223 692. 0. MG: Jacobi, r = 4, m1 = 4, m2 = 4, p = 1. .0000004777 .0017038223 928. 0. MG: Jacobi, r = 1, m1 = 1, m2 = 1, p = 2. .0000004778 .0017038341 172. 0. MG: Jacobi, r = 2, m1 = 1, m2 = 1, p = 2. .0000004777 .0017038241 316. 0. MG: Jacobi, r = 3, m1 = 1, m2 = 1, p = 2. .0000004777 .0017038236 472. 0. MG: Jacobi, r = 4, m1 = 1, m2 = 1, p = 2. .0000004777 .0017038233 628. 0. MG: Jacobi, r = 1, m1 = 2, m2 = 2, p = 2. .0000004777 .0017038225 236. 0. MG: Jacobi, r = 2, m1 = 2, m2 = 2, p = 2. .0000004777 .0017038223 468. 0. MG: Jacobi, r = 3, m1 = 2, m2 = 2, p = 2. .0000004777 .0017038223 696. 0. MG: Jacobi, r = 4, m1 = 2, m2 = 2, p = 2. .0000004777 .0017038223 928. 0. MG: Jacobi, r = 1, m1 = 3, m2 = 3, p = 2. .0000004778 .0017038234 316. 2160. MG: Jacobi, r = 2, m1 = 3, m2 = 3, p = 2. .0000004777 .0017038229 620. 0. MG: Jacobi, r = 3, m1 = 3, m2 = 3, p = 2. .0000004777 .0017038228 932. 0. MG: Jacobi, r = 4, m1 = 3, m2 = 3, p = 2. .0000004777 .0017038227 1244. 0. MG: Jacobi, r = 1, m1 = 4, m2 = 4, p = 2. .0000004777 .0017038225 392. 0. MG: Jacobi, r = 2, m1 = 4, m2 = 4, p = 2. .0000004777 .0017038223 776. 0. MG: Jacobi, r = 3, m1 = 4, m2 = 4, p = 2. .0000004777 .0017038223 1176. 2140. MG: Jacobi, r = 4, m1 = 4, m2 = 4, p = 2. .0000004777 .0017038223 1568. 0. MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 1. .0000008079 .0017040281 164. 0. MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 1. .0000004974 .0017038232 324. 0. MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 1. .0000004797 .0017038222 476. 0. MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 1. .0000004779 .0017038222 632. 0. MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 1. .0000005808 .0017038252 276. 0. MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 1. .0000004820 .0017038222 548. 0. MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 1. .0000004779 .0017038222 816. 0. MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 1. .0000004777 .0017038222 1084. 2108. MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 1. .0000005377 .0017038227 384. 0. MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 1. .0000004795 .0017038222 764. 0. MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 1. .0000004778 .0017038222 1140. 0. MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 1. .0000004777 .0017038222 1520. 0. MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 1. .0000005208 .0017038224 496. 0. MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 1. .0000004787 .0017038222 992. 0. MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 1. .0000004777 .0017038222 1480. 0. MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 1. .0000004777 .0017038222 1976. 0. MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 2. .0000005173 .0017040181 252. 0. MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 2. .0000004780 .0017038230 500. 0. MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 2. .0000004777 .0017038222 740. 2136. MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 2. .0000004777 .0017038222 984. 0. MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 2. .0000004814 .0017038246 420. 0. MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 2. .0000004777 .0017038222 832. 0. MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 2. .0000004777 .0017038222 1252. 0. MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 2. .0000004777 .0017038222 1668. 0. MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 2. .0000004784 .0017038224 584. 0. MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 2. .0000004777 .0017038222 1168. 0. MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 2. .0000004777 .0017038222 1748. 2160. MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 2. .0000004777 .0017038222 2320. 0. MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 2. .0000004779 .0017038223 760. 0. MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 2. .0000004777 .0017038222 1508. 0. MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 2. .0000004777 .0017038222 2260. 0. MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 2. .0000004777 .0017038222 3004. 0. (time (build-operator (instance-ref g1229 1) (instance-ref g1229 2) (instance-ref g1229 4) fes)) 7162 ms real time 7156 ms cpu time (7156 user, 0 system) 1 collection accounting for 2189 ms real time (2188 user, 0 system) 7160428816 bytes allocated no minor faults no major faults (time (build-rhs (instance-ref g1229 5) (instance-ref g1229 6) fes)) 3747 ms real time 3748 ms cpu time (3748 user, 0 system) 1 collection accounting for 2244 ms real time (2244 user, 0 system) 1788297608 bytes allocated no minor faults no major faults a=1, c=1, u=\sin(\pi x)+\sin(\pi y), Neumann MG: Richardson, r = 1, m1 = 1, m2 = 1, p = 1. .0282498208 .0320317159 100. 2144. MG: Richardson, r = 2, m1 = 1, m2 = 1, p = 1. .0007175888 .0025415635 196. 0. MG: Richardson, r = 3, m1 = 1, m2 = 1, p = 1. .0000137503 .0024109195 296. 0. MG: Richardson, r = 4, m1 = 1, m2 = 1, p = 1. .0000060576 .0024096730 392. 0. MG: Richardson, r = 1, m1 = 2, m2 = 2, p = 1. .0005845087 .0025595527 156. 0. MG: Richardson, r = 2, m1 = 2, m2 = 2, p = 1. .0000080577 .0024096835 304. 0. MG: Richardson, r = 3, m1 = 2, m2 = 2, p = 1. .0000011128 .0024095719 456. 0. MG: Richardson, r = 4, m1 = 2, m2 = 2, p = 1. .0000007989 .0024095705 604. 0. MG: Richardson, r = 1, m1 = 3, m2 = 3, p = 1. .0000448832 .0024123359 208. 0. MG: Richardson, r = 2, m1 = 3, m2 = 3, p = 1. .0000012567 .0024095724 424. 0. MG: Richardson, r = 3, m1 = 3, m2 = 3, p = 1. .0000007843 .0024095704 640. 2140. MG: Richardson, r = 4, m1 = 3, m2 = 3, p = 1. .0000007163 .0024095703 824. 0. MG: Richardson, r = 1, m1 = 4, m2 = 4, p = 1. .0000123912 .0024097720 264. 0. MG: Richardson, r = 2, m1 = 4, m2 = 4, p = 1. .0000008841 .0024095707 520. 0. MG: Richardson, r = 3, m1 = 4, m2 = 4, p = 1. .0000007263 .0024095704 780. 0. MG: Richardson, r = 4, m1 = 4, m2 = 4, p = 1. .0000006973 .0024095703 1040. 0. MG: Richardson, r = 1, m1 = 1, m2 = 1, p = 2. .0000016752 .0026970169 172. 0. MG: Richardson, r = 2, m1 = 1, m2 = 1, p = 2. .0000008754 .0024287318 336. 0. MG: Richardson, r = 3, m1 = 1, m2 = 1, p = 2. .0000007314 .0024108375 536. 2136. MG: Richardson, r = 4, m1 = 1, m2 = 1, p = 2. .0000006981 .0024096626 668. 0. MG: Richardson, r = 1, m1 = 2, m2 = 2, p = 2. .0000010735 .0024287451 256. 0. MG: Richardson, r = 2, m1 = 2, m2 = 2, p = 2. .0000007090 .0024096626 516. 0. MG: Richardson, r = 3, m1 = 2, m2 = 2, p = 2. .0000006886 .0024095713 776. 0. MG: Richardson, r = 4, m1 = 2, m2 = 2, p = 2. .0000006873 .0024095703 1024. 0. MG: Richardson, r = 1, m1 = 3, m2 = 3, p = 2. .0000008702 .0024108395 352. 0. MG: Richardson, r = 2, m1 = 3, m2 = 3, p = 2. .0000006900 .0024095714 696. 0. MG: Richardson, r = 3, m1 = 3, m2 = 3, p = 2. .0000006873 .0024095703 1036. 0. MG: Richardson, r = 4, m1 = 3, m2 = 3, p = 2. .0000006873 .0024095703 1380. 0. MG: Richardson, r = 1, m1 = 4, m2 = 4, p = 2. .0000007762 .0024096629 432. 0. MG: Richardson, r = 2, m1 = 4, m2 = 4, p = 2. .0000006876 .0024095703 868. 0. MG: Richardson, r = 3, m1 = 4, m2 = 4, p = 2. .0000006873 .0024095703 1324. 0. MG: Richardson, r = 4, m1 = 4, m2 = 4, p = 2. .0000006873 .0024095703 1740. 2124. MG: Jacobi, r = 1, m1 = 1, m2 = 1, p = 1. .0012224480 .0027030208 96. 0. MG: Jacobi, r = 2, m1 = 1, m2 = 1, p = 1. .0000019064 .0024095712 188. 0. MG: Jacobi, r = 3, m1 = 1, m2 = 1, p = 1. .0000006817 .0024095703 276. 0. MG: Jacobi, r = 4, m1 = 1, m2 = 1, p = 1. .0000006824 .0024095703 364. 0. MG: Jacobi, r = 1, m1 = 2, m2 = 2, p = 1. .0000017595 .0024095736 144. 0. MG: Jacobi, r = 2, m1 = 2, m2 = 2, p = 1. .0000007144 .0024095704 276. 0. MG: Jacobi, r = 3, m1 = 2, m2 = 2, p = 1. .0000006905 .0024095703 416. 0. MG: Jacobi, r = 4, m1 = 2, m2 = 2, p = 1. .0000006876 .0024095703 572. 0. MG: Jacobi, r = 1, m1 = 3, m2 = 3, p = 1. .0000009662 .0024095713 188. 0. MG: Jacobi, r = 2, m1 = 3, m2 = 3, p = 1. .0000007031 .0024095703 376. 0. MG: Jacobi, r = 3, m1 = 3, m2 = 3, p = 1. .0000006889 .0024095703 560. 0. MG: Jacobi, r = 4, m1 = 3, m2 = 3, p = 1. .0000006874 .0024095703 744. 0. MG: Jacobi, r = 1, m1 = 4, m2 = 4, p = 1. .0000008690 .0024095708 236. 0. MG: Jacobi, r = 2, m1 = 4, m2 = 4, p = 1. .0000006970 .0024095703 464. 0. MG: Jacobi, r = 3, m1 = 4, m2 = 4, p = 1. .0000006881 .0024095703 692. 0. MG: Jacobi, r = 4, m1 = 4, m2 = 4, p = 1. .0000006873 .0024095703 928. 0. MG: Jacobi, r = 1, m1 = 1, m2 = 1, p = 2. .0000006873 .0024095713 164. 0. MG: Jacobi, r = 2, m1 = 1, m2 = 1, p = 2. .0000006873 .0024095703 320. 0. MG: Jacobi, r = 3, m1 = 1, m2 = 1, p = 2. .0000006873 .0024095703 476. 0. MG: Jacobi, r = 4, m1 = 1, m2 = 1, p = 2. .0000006873 .0024095703 628. 0. MG: Jacobi, r = 1, m1 = 2, m2 = 2, p = 2. .0000006873 .0024095703 236. 0. MG: Jacobi, r = 2, m1 = 2, m2 = 2, p = 2. .0000006873 .0024095703 468. 0. MG: Jacobi, r = 3, m1 = 2, m2 = 2, p = 2. .0000006873 .0024095703 700. 0. MG: Jacobi, r = 4, m1 = 2, m2 = 2, p = 2. .0000006873 .0024095703 984. 0. MG: Jacobi, r = 1, m1 = 3, m2 = 3, p = 2. .0000006873 .0024095703 312. 0. MG: Jacobi, r = 2, m1 = 3, m2 = 3, p = 2. .0000006873 .0024095703 624. 0. MG: Jacobi, r = 3, m1 = 3, m2 = 3, p = 2. .0000006873 .0024095703 940. 0. MG: Jacobi, r = 4, m1 = 3, m2 = 3, p = 2. .0000006873 .0024095703 1240. 0. MG: Jacobi, r = 1, m1 = 4, m2 = 4, p = 2. .0000006873 .0024095703 388. 0. MG: Jacobi, r = 2, m1 = 4, m2 = 4, p = 2. .0000006873 .0024095703 780. 0. MG: Jacobi, r = 3, m1 = 4, m2 = 4, p = 2. .0000006873 .0024095703 1172. 2152. MG: Jacobi, r = 4, m1 = 4, m2 = 4, p = 2. .0000006873 .0024095703 1564. 0. MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 1. .0000012253 .0024096845 160. 0. MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 1. .0000007209 .0024095710 316. 0. MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 1. .0000006916 .0024095703 472. 0. MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 1. .0000006879 .0024095703 632. 0. MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 1. .0000008458 .0024095717 272. 0. MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 1. .0000006955 .0024095703 540. 0. MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 1. .0000006878 .0024095703 824. 0. MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 1. .0000006873 .0024095703 1080. 0. MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 1. .0000007786 .0024095706 384. 0. MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 1. .0000006908 .0024095703 768. 0. MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 1. .0000006874 .0024095703 1152. 0. MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 1. .0000006873 .0024095703 1532. 0. MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 1. .0000007511 .0024095704 500. 0. MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 1. .0000006892 .0024095703 992. 0. MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 1. .0000006873 .0024095703 1472. 0. MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 1. .0000006873 .0024095703 1968. 0. MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 2. .0000007369 .0024096769 248. 0. MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 2. .0000006876 .0024095709 496. 0. MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 2. .0000006873 .0024095703 736. 0. MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 2. .0000006873 .0024095703 984. 0. MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 2. .0000006907 .0024095713 416. 0. MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 2. .0000006873 .0024095703 832. 0. MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 2. .0000006873 .0024095703 1252. 0. MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 2. .0000006873 .0024095703 1668. 0. MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 2. .0000006876 .0024095704 588. 0. MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 2. .0000006873 .0024095703 1168. 0. MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 2. .0000006873 .0024095703 1752. 0. MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 2. .0000006873 .0024095703 2328. 0. MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 2. .0000006873 .0024095703 752. 0. MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 2. .0000006873 .0024095703 1504. 0. MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 2. .0000006873 .0024095703 2252. 0. MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 2. .0000006873 .0024095703 3008. 0. (time (build-operator (instance-ref g1229 1) (instance-ref g1229 2) (instance-ref g1229 4) fes)) 7000 ms real time 6996 ms cpu time (6996 user, 0 system) 1 collection accounting for 1983 ms real time (1980 user, 0 system) 7361755408 bytes allocated 239 minor faults no major faults (time (build-rhs (instance-ref g1229 5) (instance-ref g1229 6) fes)) 1406 ms real time 1404 ms cpu time (1404 user, 0 system) no collections 1789018504 bytes allocated no minor faults no major faults a=11.0+10.0*\cos\pi x, c=1, u=xy, Neumann MG: Richardson, r = 1, m1 = 1, m2 = 1, p = 1. .1753278103 .2003961697 100. 0. MG: Richardson, r = 2, m1 = 1, m2 = 1, p = 1. .1400530584 .1473350360 200. 0. MG: Richardson, r = 3, m1 = 1, m2 = 1, p = 1. .1125019824 .1155522973 296. 0. MG: Richardson, r = 4, m1 = 1, m2 = 1, p = 1. .0904611189 .0920072722 396. 0. MG: Richardson, r = 1, m1 = 2, m2 = 2, p = 1. .1425015101 .1518651694 156. 0. MG: Richardson, r = 2, m1 = 2, m2 = 2, p = 1. .0919996570 .0939336932 304. 0. MG: Richardson, r = 3, m1 = 2, m2 = 2, p = 1. .0595813487 .0602603487 456. 0. MG: Richardson, r = 4, m1 = 2, m2 = 2, p = 1. .0386089661 .0389336504 608. 0. MG: Richardson, r = 1, m1 = 3, m2 = 3, p = 1. .1158048260 .1206038211 208. 0. MG: Richardson, r = 2, m1 = 3, m2 = 3, p = 1. .0603570447 .0611578200 424. 2112. MG: Richardson, r = 3, m1 = 3, m2 = 3, p = 1. .0315331308 .0318029752 616. 0. MG: Richardson, r = 4, m1 = 3, m2 = 3, p = 1. .0164810150 .0166112857 820. 0. MG: Richardson, r = 1, m1 = 4, m2 = 4, p = 1. .0939475785 .0967688303 260. 0. MG: Richardson, r = 2, m1 = 4, m2 = 4, p = 1. .0395323612 .0399472649 516. 0. MG: Richardson, r = 3, m1 = 4, m2 = 4, p = 1. .0166693906 .0168091863 776. 0. MG: Richardson, r = 4, m1 = 4, m2 = 4, p = 1. .0070311329 .0071058690 1036. 0. MG: Richardson, r = 1, m1 = 1, m2 = 1, p = 2. .0000166111 .0016817835 188. 0. MG: Richardson, r = 2, m1 = 1, m2 = 1, p = 2. .0000034155 .0009467560 364. 0. MG: Richardson, r = 3, m1 = 1, m2 = 1, p = 2. .0000013849 .0007422770 500. 0. MG: Richardson, r = 4, m1 = 1, m2 = 1, p = 2. .0000007785 .0006521391 672. 0. MG: Richardson, r = 1, m1 = 2, m2 = 2, p = 2. .0000035446 .0009487689 260. 0. MG: Richardson, r = 2, m1 = 2, m2 = 2, p = 2. .0000008061 .0006523043 516. 0. MG: Richardson, r = 3, m1 = 2, m2 = 2, p = 2. .0000004284 .0005734856 772. 0. MG: Richardson, r = 4, m1 = 2, m2 = 2, p = 2. .0000003205 .0005399575 1024. 0. MG: Richardson, r = 1, m1 = 3, m2 = 3, p = 2. .0000014964 .0007435323 348. 0. MG: Richardson, r = 2, m1 = 3, m2 = 3, p = 2. .0000004430 .0005735328 696. 0. MG: Richardson, r = 3, m1 = 3, m2 = 3, p = 2. .0000003003 .0005299086 1040. 0. MG: Richardson, r = 4, m1 = 3, m2 = 3, p = 2. .0000002583 .0005119099 1388. 0. MG: Richardson, r = 1, m1 = 4, m2 = 4, p = 2. .0000008813 .0006529257 436. 0. MG: Richardson, r = 2, m1 = 4, m2 = 4, p = 2. .0000003368 .0005399877 872. 0. MG: Richardson, r = 3, m1 = 4, m2 = 4, p = 2. .0000002613 .0005119122 1364. 2156. MG: Richardson, r = 4, m1 = 4, m2 = 4, p = 2. .0000002387 .0005008008 1736. 0. MG: Jacobi, r = 1, m1 = 1, m2 = 1, p = 1. .1141930700 .1142452197 96. 0. MG: Jacobi, r = 2, m1 = 1, m2 = 1, p = 1. .0585787704 .0586100167 188. 0. MG: Jacobi, r = 3, m1 = 1, m2 = 1, p = 1. .0300166193 .0300414089 276. 0. MG: Jacobi, r = 4, m1 = 1, m2 = 1, p = 1. .0153425626 .0153724280 368. 0. MG: Jacobi, r = 1, m1 = 2, m2 = 2, p = 1. .0597805600 .0598168807 144. 0. MG: Jacobi, r = 2, m1 = 2, m2 = 2, p = 1. .0161995442 .0162298134 276. 0. MG: Jacobi, r = 3, m1 = 2, m2 = 2, p = 1. .0043888445 .0044717663 412. 0. MG: Jacobi, r = 4, m1 = 2, m2 = 2, p = 1. .0011867488 .0014577460 552. 0. MG: Jacobi, r = 1, m1 = 3, m2 = 3, p = 1. .0307017644 .0307299754 188. 0. MG: Jacobi, r = 2, m1 = 3, m2 = 3, p = 1. .0041775989 .0042645506 372. 0. MG: Jacobi, r = 3, m1 = 3, m2 = 3, p = 1. .0005446632 .0010060973 556. 0. MG: Jacobi, r = 4, m1 = 3, m2 = 3, p = 1. .0000644132 .0008481809 744. 0. MG: Jacobi, r = 1, m1 = 4, m2 = 4, p = 1. .0163196338 .0163505761 236. 0. MG: Jacobi, r = 2, m1 = 4, m2 = 4, p = 1. .0011924705 .0014624586 464. 0. MG: Jacobi, r = 3, m1 = 4, m2 = 4, p = 1. .0000840710 .0008498906 736. 0. MG: Jacobi, r = 4, m1 = 4, m2 = 4, p = 1. .0000050683 .0008457326 956. 2156. MG: Jacobi, r = 1, m1 = 1, m2 = 1, p = 2. .0000001955 .0008458622 160. 0. MG: Jacobi, r = 2, m1 = 1, m2 = 1, p = 2. .0000001954 .0008457383 316. 0. MG: Jacobi, r = 3, m1 = 1, m2 = 1, p = 2. .0000001954 .0008457343 472. 0. MG: Jacobi, r = 4, m1 = 1, m2 = 1, p = 2. .0000001954 .0008457334 628. 0. MG: Jacobi, r = 1, m1 = 2, m2 = 2, p = 2. .0000005185 .0008457240 240. 0. MG: Jacobi, r = 2, m1 = 2, m2 = 2, p = 2. .0000005185 .0008457209 472. 0. MG: Jacobi, r = 3, m1 = 2, m2 = 2, p = 2. .0000005185 .0008457200 704. 0. MG: Jacobi, r = 4, m1 = 2, m2 = 2, p = 2. .0000005185 .0008457192 932. 2168. MG: Jacobi, r = 1, m1 = 3, m2 = 3, p = 2. .0000001954 .0008457327 312. 0. MG: Jacobi, r = 2, m1 = 3, m2 = 3, p = 2. .0000001954 .0008457315 628. 0. MG: Jacobi, r = 3, m1 = 3, m2 = 3, p = 2. .0000001954 .0008457304 936. 0. MG: Jacobi, r = 4, m1 = 3, m2 = 3, p = 2. .0000001954 .0008457293 1240. 0. MG: Jacobi, r = 1, m1 = 4, m2 = 4, p = 2. .0000005185 .0008457213 400. 0. MG: Jacobi, r = 2, m1 = 4, m2 = 4, p = 2. .0000005185 .0008457199 816. 2124. MG: Jacobi, r = 3, m1 = 4, m2 = 4, p = 2. .0000005185 .0008457187 1168. 0. MG: Jacobi, r = 4, m1 = 4, m2 = 4, p = 2. .0000005185 .0008457176 1552. 0. MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 1. .0383371530 .0383604348 160. 0. MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 1. .0056120375 .0056361738 316. 0. MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 1. .0008203165 .0009550031 476. 0. MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 1. .0001193443 .0005026667 628. 0. MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 1. .0061596029 .0061824211 276. 0. MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 1. .0001456883 .0005095728 540. 0. MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 1. .0000033013 .0004882917 808. 0. MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 1. .0000002252 .0004882807 1088. 0. MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 1. .0009998515 .0011132244 388. 0. MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 1. .0000036950 .0004882945 768. 0. MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 1. .0000002186 .0004882807 1152. 0. MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 1. .0000002187 .0004882807 1524. 0. MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 1. .0001622675 .0005145643 496. 0. MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 1. .0000002319 .0004882807 984. 0. MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 1. .0000002187 .0004882807 1476. 0. MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 1. .0000002190 .0004882807 1964. 0. MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 2. .0000002615 .0004884743 248. 0. MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 2. .0000002198 .0004882815 496. 0. MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 2. .0000002192 .0004882807 740. 0. MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 2. .0000002192 .0004882807 984. 0. MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 2. .0000002237 .0004882831 416. 0. MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 2. .0000002192 .0004882807 832. 0. MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 2. .0000002192 .0004882807 1244. 0. MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 2. .0000002192 .0004882807 1652. 0. MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 2. .0000002200 .0004882809 580. 0. MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 2. .0000002192 .0004882807 1168. 2128. MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 2. .0000002192 .0004882807 1748. 0. MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 2. .0000002192 .0004882807 2336. 0. MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 2. .0000002194 .0004882807 756. 0. MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 2. .0000002192 .0004882807 1496. 0. MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 2. .0000002192 .0004882807 2252. 0. MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 2. .0000002192 .0004882807 3008. 2160. (time (build-operator (instance-ref g1229 1) (instance-ref g1229 2) (instance-ref g1229 4) fes)) 4242 ms real time 4240 ms cpu time (4240 user, 0 system) no collections 5851805968 bytes allocated no minor faults no major faults (time (build-rhs (instance-ref g1229 5) (instance-ref g1229 6) fes)) 3589 ms real time 3588 ms cpu time (3588 user, 0 system) 1 collection accounting for 2250 ms real time (2252 user, 0 system) 1788821896 bytes allocated no minor faults no major faults a=1, c=1+xy, u=x^2-y^2, Neumann MG: Richardson, r = 1, m1 = 1, m2 = 1, p = 1. .0000533967 .0009296052 104. 0. MG: Richardson, r = 2, m1 = 1, m2 = 1, p = 1. .0000001385 .0006961122 200. 0. MG: Richardson, r = 3, m1 = 1, m2 = 1, p = 1. .0000000977 .0006909053 296. 0. MG: Richardson, r = 4, m1 = 1, m2 = 1, p = 1. .0000000974 .0006905616 392. 0. MG: Richardson, r = 1, m1 = 2, m2 = 2, p = 1. .0000002481 .0006961373 156. 0. MG: Richardson, r = 2, m1 = 2, m2 = 2, p = 1. .0000000974 .0006905616 304. 0. MG: Richardson, r = 3, m1 = 2, m2 = 2, p = 1. .0000000974 .0006905343 456. 0. MG: Richardson, r = 4, m1 = 2, m2 = 2, p = 1. .0000000974 .0006905339 608. 0. MG: Richardson, r = 1, m1 = 3, m2 = 3, p = 1. .0000000982 .0006909066 212. 0. MG: Richardson, r = 2, m1 = 3, m2 = 3, p = 1. .0000000974 .0006905343 416. 0. MG: Richardson, r = 3, m1 = 3, m2 = 3, p = 1. .0000000974 .0006905339 648. 2164. MG: Richardson, r = 4, m1 = 3, m2 = 3, p = 1. .0000000974 .0006905339 820. 0. MG: Richardson, r = 1, m1 = 4, m2 = 4, p = 1. .0000000975 .0006905618 268. 0. MG: Richardson, r = 2, m1 = 4, m2 = 4, p = 1. .0000000974 .0006905339 520. 0. MG: Richardson, r = 3, m1 = 4, m2 = 4, p = 1. .0000000974 .0006905339 776. 0. MG: Richardson, r = 4, m1 = 4, m2 = 4, p = 1. .0000000974 .0006905339 1032. 0. MG: Richardson, r = 1, m1 = 1, m2 = 1, p = 2. .0000001396 .0007732508 172. 0. MG: Richardson, r = 2, m1 = 1, m2 = 1, p = 2. .0000001005 .0006960696 340. 0. MG: Richardson, r = 3, m1 = 1, m2 = 1, p = 2. .0000000976 .0006909036 520. 0. MG: Richardson, r = 4, m1 = 1, m2 = 1, p = 2. .0000000974 .0006905614 672. 0. MG: Richardson, r = 1, m1 = 2, m2 = 2, p = 2. .0000001005 .0006960781 264. 0. MG: Richardson, r = 2, m1 = 2, m2 = 2, p = 2. .0000000974 .0006905614 516. 0. MG: Richardson, r = 3, m1 = 2, m2 = 2, p = 2. .0000000974 .0006905343 776. 0. MG: Richardson, r = 4, m1 = 2, m2 = 2, p = 2. .0000000974 .0006905339 1024. 0. MG: Richardson, r = 1, m1 = 3, m2 = 3, p = 2. .0000000976 .0006909049 352. 0. MG: Richardson, r = 2, m1 = 3, m2 = 3, p = 2. .0000000974 .0006905343 692. 0. MG: Richardson, r = 3, m1 = 3, m2 = 3, p = 2. .0000000974 .0006905339 1036. 0. MG: Richardson, r = 4, m1 = 3, m2 = 3, p = 2. .0000000974 .0006905339 1384. 0. MG: Richardson, r = 1, m1 = 4, m2 = 4, p = 2. .0000000974 .0006905616 436. 0. MG: Richardson, r = 2, m1 = 4, m2 = 4, p = 2. .0000000974 .0006905339 864. 0. MG: Richardson, r = 3, m1 = 4, m2 = 4, p = 2. .0000000974 .0006905339 1320. 0. MG: Richardson, r = 4, m1 = 4, m2 = 4, p = 2. .0000000974 .0006905339 1748. 2112. MG: Jacobi, r = 1, m1 = 1, m2 = 1, p = 1. .0000000974 .0006905339 96. 0. MG: Jacobi, r = 2, m1 = 1, m2 = 1, p = 1. .0000000974 .0006905339 184. 0. MG: Jacobi, r = 3, m1 = 1, m2 = 1, p = 1. .0000000974 .0006905339 276. 0. MG: Jacobi, r = 4, m1 = 1, m2 = 1, p = 1. .0000000974 .0006905339 368. 0. MG: Jacobi, r = 1, m1 = 2, m2 = 2, p = 1. .0000000974 .0006905339 144. 0. MG: Jacobi, r = 2, m1 = 2, m2 = 2, p = 1. .0000000974 .0006905339 276. 0. MG: Jacobi, r = 3, m1 = 2, m2 = 2, p = 1. .0000000974 .0006905339 416. 0. MG: Jacobi, r = 4, m1 = 2, m2 = 2, p = 1. .0000000974 .0006905339 564. 0. MG: Jacobi, r = 1, m1 = 3, m2 = 3, p = 1. .0000000974 .0006905339 188. 0. MG: Jacobi, r = 2, m1 = 3, m2 = 3, p = 1. .0000000974 .0006905339 372. 0. MG: Jacobi, r = 3, m1 = 3, m2 = 3, p = 1. .0000000974 .0006905339 556. 0. MG: Jacobi, r = 4, m1 = 3, m2 = 3, p = 1. .0000000974 .0006905339 740. 0. MG: Jacobi, r = 1, m1 = 4, m2 = 4, p = 1. .0000000974 .0006905339 236. 0. MG: Jacobi, r = 2, m1 = 4, m2 = 4, p = 1. .0000000974 .0006905339 464. 0. MG: Jacobi, r = 3, m1 = 4, m2 = 4, p = 1. .0000000974 .0006905339 696. 0. MG: Jacobi, r = 4, m1 = 4, m2 = 4, p = 1. .0000000974 .0006905339 932. 0. MG: Jacobi, r = 1, m1 = 1, m2 = 1, p = 2. .0000000974 .0006905339 172. 2080. MG: Jacobi, r = 2, m1 = 1, m2 = 1, p = 2. .0000000974 .0006905339 312. 0. MG: Jacobi, r = 3, m1 = 1, m2 = 1, p = 2. .0000000974 .0006905339 472. 0. MG: Jacobi, r = 4, m1 = 1, m2 = 1, p = 2. .0000000974 .0006905339 632. 0. MG: Jacobi, r = 1, m1 = 2, m2 = 2, p = 2. .0000000974 .0006905339 240. 0. MG: Jacobi, r = 2, m1 = 2, m2 = 2, p = 2. .0000000974 .0006905339 468. 0. MG: Jacobi, r = 3, m1 = 2, m2 = 2, p = 2. .0000000974 .0006905339 704. 0. MG: Jacobi, r = 4, m1 = 2, m2 = 2, p = 2. .0000000974 .0006905339 932. 0. MG: Jacobi, r = 1, m1 = 3, m2 = 3, p = 2. .0000000974 .0006905339 316. 0. MG: Jacobi, r = 2, m1 = 3, m2 = 3, p = 2. .0000000974 .0006905339 628. 0. MG: Jacobi, r = 3, m1 = 3, m2 = 3, p = 2. .0000000974 .0006905339 932. 0. MG: Jacobi, r = 4, m1 = 3, m2 = 3, p = 2. .0000000974 .0006905339 1244. 0. MG: Jacobi, r = 1, m1 = 4, m2 = 4, p = 2. .0000000974 .0006905339 392. 0. MG: Jacobi, r = 2, m1 = 4, m2 = 4, p = 2. .0000000974 .0006905339 804. 0. MG: Jacobi, r = 3, m1 = 4, m2 = 4, p = 2. .0000000974 .0006905339 1196. 2128. MG: Jacobi, r = 4, m1 = 4, m2 = 4, p = 2. .0000000974 .0006905339 1552. 0. MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 1. .0000000979 .0006906556 164. 0. MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 1. .0000000974 .0006905343 320. 0. MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 1. .0000000974 .0006905339 476. 0. MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 1. .0000000974 .0006905339 628. 0. MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 1. .0000000974 .0006905355 272. 0. MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 1. .0000000974 .0006905339 544. 0. MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 1. .0000000974 .0006905339 812. 0. MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 1. .0000000974 .0006905339 1080. 0. MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 1. .0000000974 .0006905341 384. 0. MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 1. .0000000974 .0006905339 764. 0. MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 1. .0000000974 .0006905339 1140. 0. MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 1. .0000000974 .0006905339 1532. 0. MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 1. .0000000974 .0006905340 496. 0. MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 1. .0000000974 .0006905339 988. 0. MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 1. .0000000974 .0006905339 1488. 0. MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 1. .0000000974 .0006905339 1980. 0. MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 2. .0000000977 .0006906520 252. 0. MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 2. .0000000974 .0006905343 496. 0. MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 2. .0000000974 .0006905339 740. 0. MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 2. .0000000974 .0006905339 988. 0. MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 2. .0000000974 .0006905355 420. 0. MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 2. .0000000974 .0006905339 832. 0. MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 2. .0000000974 .0006905339 1244. 0. MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 2. .0000000974 .0006905339 1660. 0. MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 2. .0000000974 .0006905341 588. 0. MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 2. .0000000974 .0006905339 1168. 0. MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 2. .0000000974 .0006905339 1756. 2184. MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 2. .0000000974 .0006905339 2324. 0. MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 2. .0000000974 .0006905340 756. 0. MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 2. .0000000974 .0006905339 1500. 0. MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 2. .0000000974 .0006905339 2248. 0. MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 2. .0000000974 .0006905339 2996. 0. (time (build-operator (instance-ref g1229 1) (instance-ref g1229 2) (instance-ref g1229 4) fes)) 9511 ms real time 9508 ms cpu time (9508 user, 0 system) 2 collections accounting for 4456 ms real time (4456 user, 0 system) 7361755408 bytes allocated no minor faults no major faults (time (build-rhs (instance-ref g1229 5) (instance-ref g1229 6) fes)) 1538 ms real time 1536 ms cpu time (1536 user, 0 system) no collections 1786790216 bytes allocated no minor faults no major faults a=1, c=2+\cos(\pi y), u=\cos(\pi x), Neumann MG: Richardson, r = 1, m1 = 1, m2 = 1, p = 1. .0009729238 .0077032496 104. 0. MG: Richardson, r = 2, m1 = 1, m2 = 1, p = 1. .0000225402 .0017288199 200. 0. MG: Richardson, r = 3, m1 = 1, m2 = 1, p = 1. .0000017413 .0017048403 296. 0. MG: Richardson, r = 4, m1 = 1, m2 = 1, p = 1. .0000007001 .0017039012 396. 0. MG: Richardson, r = 1, m1 = 2, m2 = 2, p = 1. .0000350702 .0017435099 156. 0. MG: Richardson, r = 2, m1 = 2, m2 = 2, p = 1. .0000008464 .0017039070 308. 0. MG: Richardson, r = 3, m1 = 2, m2 = 2, p = 1. .0000005124 .0017038243 452. 0. MG: Richardson, r = 4, m1 = 2, m2 = 2, p = 1. .0000004670 .0017038225 604. 0. MG: Richardson, r = 1, m1 = 3, m2 = 3, p = 1. .0000046311 .0017052842 208. 0. MG: Richardson, r = 2, m1 = 3, m2 = 3, p = 1. .0000005497 .0017038249 428. 0. MG: Richardson, r = 3, m1 = 3, m2 = 3, p = 1. .0000004660 .0017038224 612. 0. MG: Richardson, r = 4, m1 = 3, m2 = 3, p = 1. .0000004519 .0017038223 816. 0. MG: Richardson, r = 1, m1 = 4, m2 = 4, p = 1. .0000016638 .0017039552 264. 0. MG: Richardson, r = 2, m1 = 4, m2 = 4, p = 1. .0000004951 .0017038228 520. 0. MG: Richardson, r = 3, m1 = 4, m2 = 4, p = 1. .0000004550 .0017038223 776. 0. MG: Richardson, r = 4, m1 = 4, m2 = 4, p = 1. .0000004486 .0017038222 1028. 0. MG: Richardson, r = 1, m1 = 1, m2 = 1, p = 2. .0000013180 .0019088210 176. 0. MG: Richardson, r = 2, m1 = 1, m2 = 1, p = 2. .0000006176 .0017175961 368. 0. MG: Richardson, r = 3, m1 = 1, m2 = 1, p = 2. .0000004874 .0017047513 504. 0. MG: Richardson, r = 4, m1 = 1, m2 = 1, p = 2. .0000004569 .0017038930 672. 0. MG: Richardson, r = 1, m1 = 2, m2 = 2, p = 2. .0000007943 .0017176336 264. 0. MG: Richardson, r = 2, m1 = 2, m2 = 2, p = 2. .0000004669 .0017038932 516. 0. MG: Richardson, r = 3, m1 = 2, m2 = 2, p = 2. .0000004481 .0017038232 776. 0. MG: Richardson, r = 4, m1 = 2, m2 = 2, p = 2. .0000004469 .0017038223 1024. 0. MG: Richardson, r = 1, m1 = 3, m2 = 3, p = 2. .0000006132 .0017047585 364. 0. MG: Richardson, r = 2, m1 = 3, m2 = 3, p = 2. .0000004494 .0017038232 700. 2172. MG: Richardson, r = 3, m1 = 3, m2 = 3, p = 2. .0000004469 .0017038222 1036. 0. MG: Richardson, r = 4, m1 = 3, m2 = 3, p = 2. .0000004469 .0017038222 1380. 0. MG: Richardson, r = 1, m1 = 4, m2 = 4, p = 2. .0000005282 .0017038948 436. 0. MG: Richardson, r = 2, m1 = 4, m2 = 4, p = 2. .0000004472 .0017038223 868. 0. MG: Richardson, r = 3, m1 = 4, m2 = 4, p = 2. .0000004469 .0017038222 1308. 2180. MG: Richardson, r = 4, m1 = 4, m2 = 4, p = 2. .0000004469 .0017038222 1728. 0. MG: Jacobi, r = 1, m1 = 1, m2 = 1, p = 1. .0000021033 .0017039480 96. 0. MG: Jacobi, r = 2, m1 = 1, m2 = 1, p = 1. .0000004978 .0017038246 184. 0. MG: Jacobi, r = 3, m1 = 1, m2 = 1, p = 1. .0000004558 .0017038236 272. 0. MG: Jacobi, r = 4, m1 = 1, m2 = 1, p = 1. .0000004489 .0017038233 368. 0. MG: Jacobi, r = 1, m1 = 2, m2 = 2, p = 1. .0000007430 .0017038252 144. 0. MG: Jacobi, r = 2, m1 = 2, m2 = 2, p = 1. .0000004643 .0017038224 276. 0. MG: Jacobi, r = 3, m1 = 2, m2 = 2, p = 1. .0000004490 .0017038223 412. 0. MG: Jacobi, r = 4, m1 = 2, m2 = 2, p = 1. .0000004471 .0017038223 548. 2172. MG: Jacobi, r = 1, m1 = 3, m2 = 3, p = 1. .0000006119 .0017038246 188. 0. MG: Jacobi, r = 2, m1 = 3, m2 = 3, p = 1. .0000004559 .0017038230 372. 0. MG: Jacobi, r = 3, m1 = 3, m2 = 3, p = 1. .0000004477 .0017038228 556. 0. MG: Jacobi, r = 4, m1 = 3, m2 = 3, p = 1. .0000004469 .0017038227 736. 0. MG: Jacobi, r = 1, m1 = 4, m2 = 4, p = 1. .0000005586 .0017038230 232. 0. MG: Jacobi, r = 2, m1 = 4, m2 = 4, p = 1. .0000004523 .0017038223 464. 0. MG: Jacobi, r = 3, m1 = 4, m2 = 4, p = 1. .0000004472 .0017038223 712. 0. MG: Jacobi, r = 4, m1 = 4, m2 = 4, p = 1. .0000004469 .0017038223 964. 2156. MG: Jacobi, r = 1, m1 = 1, m2 = 1, p = 2. .0000004469 .0017038341 160. 0. MG: Jacobi, r = 2, m1 = 1, m2 = 1, p = 2. .0000004469 .0017038241 316. 0. MG: Jacobi, r = 3, m1 = 1, m2 = 1, p = 2. .0000004469 .0017038236 472. 0. MG: Jacobi, r = 4, m1 = 1, m2 = 1, p = 2. .0000004469 .0017038233 628. 0. MG: Jacobi, r = 1, m1 = 2, m2 = 2, p = 2. .0000004469 .0017038225 236. 0. MG: Jacobi, r = 2, m1 = 2, m2 = 2, p = 2. .0000004469 .0017038223 468. 0. MG: Jacobi, r = 3, m1 = 2, m2 = 2, p = 2. .0000004469 .0017038223 696. 0. MG: Jacobi, r = 4, m1 = 2, m2 = 2, p = 2. .0000004468 .0017038223 928. 2172. MG: Jacobi, r = 1, m1 = 3, m2 = 3, p = 2. .0000004469 .0017038234 312. 0. MG: Jacobi, r = 2, m1 = 3, m2 = 3, p = 2. .0000004469 .0017038230 624. 0. MG: Jacobi, r = 3, m1 = 3, m2 = 3, p = 2. .0000004469 .0017038228 932. 0. MG: Jacobi, r = 4, m1 = 3, m2 = 3, p = 2. .0000004469 .0017038227 1244. 0. MG: Jacobi, r = 1, m1 = 4, m2 = 4, p = 2. .0000004469 .0017038225 392. 0. MG: Jacobi, r = 2, m1 = 4, m2 = 4, p = 2. .0000004469 .0017038223 788. 2148. MG: Jacobi, r = 3, m1 = 4, m2 = 4, p = 2. .0000004468 .0017038223 1164. 0. MG: Jacobi, r = 4, m1 = 4, m2 = 4, p = 2. .0000004468 .0017038223 1552. 0. MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 1. .0000007166 .0017040279 164. 0. MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 1. .0000004624 .0017038232 316. 0. MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 1. .0000004484 .0017038222 476. 0. MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 1. .0000004470 .0017038222 632. 0. MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 1. .0000005309 .0017038252 272. 0. MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 1. .0000004501 .0017038222 540. 0. MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 1. .0000004470 .0017038222 812. 0. MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 1. .0000004469 .0017038222 1080. 0. MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 1. .0000004961 .0017038227 384. 0. MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 1. .0000004483 .0017038222 764. 0. MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 1. .0000004469 .0017038222 1140. 0. MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 1. .0000004469 .0017038222 1528. 0. MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 1. .0000004824 .0017038224 496. 0. MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 1. .0000004477 .0017038222 984. 0. MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 1. .0000004469 .0017038222 1484. 0. MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 1. .0000004469 .0017038222 1976. 0. MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 2. .0000004824 .0017040181 252. 0. MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 2. .0000004471 .0017038231 496. 0. MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 2. .0000004469 .0017038222 740. 0. MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 2. .0000004469 .0017038222 988. 0. MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 2. .0000004501 .0017038247 424. 0. MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 2. .0000004469 .0017038222 832. 0. MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 2. .0000004469 .0017038222 1244. 0. MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 2. .0000004469 .0017038222 1660. 0. MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 2. .0000004474 .0017038224 584. 0. MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 2. .0000004469 .0017038222 1168. 2160. MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 2. .0000004469 .0017038222 1744. 0. MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 2. .0000004469 .0017038222 2328. 0. MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 2. .0000004470 .0017038223 756. 0. MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 2. .0000004469 .0017038222 1500. 0. MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 2. .0000004469 .0017038222 2252. 0. MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 2. .0000004469 .0017038222 3000. 2184. (time (build-operator (instance-ref g1229 1) (instance-ref g1229 2) (instance-ref g1229 4) fes)) 7220 ms real time 7216 ms cpu time (7212 user, 4 system) 1 collection accounting for 2225 ms real time (2220 user, 4 system) 7162853712 bytes allocated no minor faults no major faults (time (build-rhs (instance-ref g1229 5) (instance-ref g1229 6) fes)) 1519 ms real time 1516 ms cpu time (1516 user, 0 system) no collections 1788559752 bytes allocated 1 minor fault no major faults 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. .0007421081 .0051777775 100. 0. MG: Richardson, r = 2, m1 = 1, m2 = 1, p = 1. .0000072749 .0024292940 196. 0. MG: Richardson, r = 3, m1 = 1, m2 = 1, p = 1. .0000011423 .0024108409 300. 0. MG: Richardson, r = 4, m1 = 1, m2 = 1, p = 1. .0000006935 .0024096627 392. 0. MG: Richardson, r = 1, m1 = 2, m2 = 2, p = 1. .0000140203 .0024301381 156. 0. MG: Richardson, r = 2, m1 = 2, m2 = 2, p = 1. .0000007867 .0024096629 304. 0. MG: Richardson, r = 3, m1 = 2, m2 = 2, p = 1. .0000005596 .0024095713 456. 0. MG: Richardson, r = 4, m1 = 2, m2 = 2, p = 1. .0000005210 .0024095704 604. 0. MG: Richardson, r = 1, m1 = 3, m2 = 3, p = 1. .0000024483 .0024108602 208. 0. MG: Richardson, r = 2, m1 = 3, m2 = 3, p = 1. .0000005874 .0024095714 416. 0. MG: Richardson, r = 3, m1 = 3, m2 = 3, p = 1. .0000005203 .0024095703 616. 0. MG: Richardson, r = 4, m1 = 3, m2 = 3, p = 1. .0000005100 .0024095703 816. 0. MG: Richardson, r = 1, m1 = 4, m2 = 4, p = 1. .0000011991 .0024096651 264. 0. MG: Richardson, r = 2, m1 = 4, m2 = 4, p = 1. .0000005427 .0024095703 520. 0. MG: Richardson, r = 3, m1 = 4, m2 = 4, p = 1. .0000005122 .0024095703 772. 0. MG: Richardson, r = 4, m1 = 4, m2 = 4, p = 1. .0000005082 .0024095703 1060. 0. MG: Richardson, r = 1, m1 = 1, m2 = 1, p = 2. .0000014483 .0026970079 184. 0. MG: Richardson, r = 2, m1 = 1, m2 = 1, p = 2. .0000006774 .0024287308 348. 2120. MG: Richardson, r = 3, m1 = 1, m2 = 1, p = 2. .0000005455 .0024108373 504. 0. MG: Richardson, r = 4, m1 = 1, m2 = 1, p = 2. .0000005166 .0024096625 672. 0. MG: Richardson, r = 1, m1 = 2, m2 = 2, p = 2. .0000008634 .0024287432 260. 0. MG: Richardson, r = 2, m1 = 2, m2 = 2, p = 2. .0000005260 .0024096625 516. 0. MG: Richardson, r = 3, m1 = 2, m2 = 2, p = 2. .0000005085 .0024095714 776. 0. MG: Richardson, r = 4, m1 = 2, m2 = 2, p = 2. .0000005074 .0024095704 1028. 0. MG: Richardson, r = 1, m1 = 3, m2 = 3, p = 2. .0000006711 .0024108389 348. 2108. MG: Richardson, r = 2, m1 = 3, m2 = 3, p = 2. .0000005097 .0024095714 696. 0. MG: Richardson, r = 3, m1 = 3, m2 = 3, p = 2. .0000005074 .0024095703 1036. 0. MG: Richardson, r = 4, m1 = 3, m2 = 3, p = 2. .0000005074 .0024095703 1380. 0. MG: Richardson, r = 1, m1 = 4, m2 = 4, p = 2. .0000005852 .0024096627 436. 0. MG: Richardson, r = 2, m1 = 4, m2 = 4, p = 2. .0000005077 .0024095704 864. 0. MG: Richardson, r = 3, m1 = 4, m2 = 4, p = 2. .0000005074 .0024095703 1328. 2128. MG: Richardson, r = 4, m1 = 4, m2 = 4, p = 2. .0000005074 .0024095703 1724. 0. MG: Jacobi, r = 1, m1 = 1, m2 = 1, p = 1. .0000023723 .0024096019 100. 0. MG: Jacobi, r = 2, m1 = 1, m2 = 1, p = 1. .0000005610 .0024095703 188. 0. MG: Jacobi, r = 3, m1 = 1, m2 = 1, p = 1. .0000005167 .0024095703 276. 0. MG: Jacobi, r = 4, m1 = 1, m2 = 1, p = 1. .0000005094 .0024095703 364. 0. MG: Jacobi, r = 1, m1 = 2, m2 = 2, p = 1. .0000008289 .0024095710 140. 0. MG: Jacobi, r = 2, m1 = 2, m2 = 2, p = 1. .0000005252 .0024095703 280. 0. MG: Jacobi, r = 3, m1 = 2, m2 = 2, p = 1. .0000005096 .0024095703 432. 2148. MG: Jacobi, r = 4, m1 = 2, m2 = 2, p = 1. .0000005077 .0024095703 556. 0. MG: Jacobi, r = 1, m1 = 3, m2 = 3, p = 1. .0000006755 .0024095705 188. 0. MG: Jacobi, r = 2, m1 = 3, m2 = 3, p = 1. .0000005165 .0024095703 372. 0. MG: Jacobi, r = 3, m1 = 3, m2 = 3, p = 1. .0000005082 .0024095703 556. 0. MG: Jacobi, r = 4, m1 = 3, m2 = 3, p = 1. .0000005075 .0024095703 740. 0. MG: Jacobi, r = 1, m1 = 4, m2 = 4, p = 1. .0000006191 .0024095704 236. 0. MG: Jacobi, r = 2, m1 = 4, m2 = 4, p = 1. .0000005128 .0024095703 468. 0. MG: Jacobi, r = 3, m1 = 4, m2 = 4, p = 1. .0000005078 .0024095703 716. 0. MG: Jacobi, r = 4, m1 = 4, m2 = 4, p = 1. .0000005074 .0024095703 924. 0. MG: Jacobi, r = 1, m1 = 1, m2 = 1, p = 2. .0000005074 .0024095704 160. 0. MG: Jacobi, r = 2, m1 = 1, m2 = 1, p = 2. .0000005074 .0024095703 316. 0. MG: Jacobi, r = 3, m1 = 1, m2 = 1, p = 2. .0000005074 .0024095703 472. 0. MG: Jacobi, r = 4, m1 = 1, m2 = 1, p = 2. .0000005074 .0024095703 628. 0. MG: Jacobi, r = 1, m1 = 2, m2 = 2, p = 2. .0000005074 .0024095703 236. 0. MG: Jacobi, r = 2, m1 = 2, m2 = 2, p = 2. .0000005074 .0024095703 468. 0. MG: Jacobi, r = 3, m1 = 2, m2 = 2, p = 2. .0000005074 .0024095703 752. 0. MG: Jacobi, r = 4, m1 = 2, m2 = 2, p = 2. .0000005074 .0024095703 924. 0. MG: Jacobi, r = 1, m1 = 3, m2 = 3, p = 2. .0000005074 .0024095703 312. 0. MG: Jacobi, r = 2, m1 = 3, m2 = 3, p = 2. .0000005074 .0024095703 620. 0. MG: Jacobi, r = 3, m1 = 3, m2 = 3, p = 2. .0000005074 .0024095703 924. 0. MG: Jacobi, r = 4, m1 = 3, m2 = 3, p = 2. .0000005074 .0024095703 1232. 0. MG: Jacobi, r = 1, m1 = 4, m2 = 4, p = 2. .0000005074 .0024095703 416. 0. MG: Jacobi, r = 2, m1 = 4, m2 = 4, p = 2. .0000005074 .0024095703 784. 2124. MG: Jacobi, r = 3, m1 = 4, m2 = 4, p = 2. .0000005074 .0024095703 1156. 0. MG: Jacobi, r = 4, m1 = 4, m2 = 4, p = 2. .0000005074 .0024095703 1544. 0. MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 1. .0000007736 .0024096824 160. 0. MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 1. .0000005222 .0024095709 316. 0. MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 1. .0000005089 .0024095703 472. 0. MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 1. .0000005075 .0024095703 628. 0. MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 1. .0000005902 .0024095713 272. 0. MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 1. .0000005107 .0024095703 540. 0. MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 1. .0000005076 .0024095703 804. 0. MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 1. .0000005074 .0024095703 1080. 0. MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 1. .0000005558 .0024095704 384. 0. MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 1. .0000005088 .0024095703 764. 0. MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 1. .0000005074 .0024095703 1148. 0. MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 1. .0000005074 .0024095703 1524. 0. MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 1. .0000005418 .0024095703 496. 0. MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 1. .0000005082 .0024095703 992. 0. MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 1. .0000005074 .0024095703 1488. 0. MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 1. .0000005074 .0024095703 1972. 0. MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 2. .0000005356 .0024096768 252. 0. MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 2. .0000005075 .0024095709 500. 0. MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 2. .0000005074 .0024095703 740. 0. MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 2. .0000005074 .0024095703 988. 0. MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 2. .0000005094 .0024095713 416. 0. MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 2. .0000005074 .0024095703 828. 0. MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 2. .0000005074 .0024095703 1252. 0. MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 2. .0000005074 .0024095703 1672. 0. MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 2. .0000005077 .0024095704 588. 0. MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 2. .0000005074 .0024095703 1168. 0. MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 2. .0000005074 .0024095703 1748. 0. MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 2. .0000005074 .0024095703 2332. 0. MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 2. .0000005075 .0024095703 752. 0. MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 2. .0000005074 .0024095703 1504. 0. MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 2. .0000005074 .0024095703 2248. 0. MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 2. .0000005074 .0024095703 3008. 2128. (time (build-operator (instance-ref g1229 1) (instance-ref g1229 2) (instance-ref g1229 4) fes)) 7280 ms real time 7276 ms cpu time (7276 user, 0 system) 1 collection accounting for 2245 ms real time (2244 user, 0 system) 7364180304 bytes allocated no minor faults no major faults (time (build-rhs (instance-ref g1229 5) (instance-ref g1229 6) fes)) 1397 ms real time 1396 ms cpu time (1392 user, 4 system) no collections 1789280648 bytes allocated no minor faults no major faults 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. .0316518694 .0855727315 100. 0. MG: Richardson, r = 2, m1 = 1, m2 = 1, p = 1. .0047757699 .0346159366 200. 0. MG: Richardson, r = 3, m1 = 1, m2 = 1, p = 1. .0017430696 .0189542431 296. 0. MG: Richardson, r = 4, m1 = 1, m2 = 1, p = 1. .0010095611 .0116471522 392. 0. MG: Richardson, r = 1, m1 = 2, m2 = 2, p = 1. .0056426647 .0395919136 156. 0. MG: Richardson, r = 2, m1 = 2, m2 = 2, p = 1. .0011924499 .0130363203 308. 0. MG: Richardson, r = 3, m1 = 2, m2 = 2, p = 1. .0004368371 .0058088106 456. 0. MG: Richardson, r = 4, m1 = 2, m2 = 2, p = 1. .0001932900 .0030630199 604. 0. MG: Richardson, r = 1, m1 = 3, m2 = 3, p = 1. .0025668134 .0239542834 208. 0. MG: Richardson, r = 2, m1 = 3, m2 = 3, p = 1. .0004961397 .0063265476 412. 0. MG: Richardson, r = 3, m1 = 3, m2 = 3, p = 1. .0001493700 .0024786954 620. 0. MG: Richardson, r = 4, m1 = 3, m2 = 3, p = 1. .0000601814 .0012428185 816. 0. MG: Richardson, r = 1, m1 = 4, m2 = 4, p = 1. .0016133546 .0158151174 264. 0. MG: Richardson, r = 2, m1 = 4, m2 = 4, p = 1. .0002365424 .0034986893 516. 0. MG: Richardson, r = 3, m1 = 4, m2 = 4, p = 1. .0000641439 .0012915059 772. 0. MG: Richardson, r = 4, m1 = 4, m2 = 4, p = 1. .0000234644 .0007071436 1088. 0. MG: Richardson, r = 1, m1 = 1, m2 = 1, p = 2. .0000155725 .0016582046 192. 2168. MG: Richardson, r = 2, m1 = 1, m2 = 1, p = 2. .0000033327 .0009441643 340. 0. MG: Richardson, r = 3, m1 = 1, m2 = 1, p = 2. .0000013685 .0007416840 504. 0. MG: Richardson, r = 4, m1 = 1, m2 = 1, p = 2. .0000007649 .0006519458 668. 0. MG: Richardson, r = 1, m1 = 2, m2 = 2, p = 2. .0000034363 .0009458762 260. 0. MG: Richardson, r = 2, m1 = 2, m2 = 2, p = 2. .0000007890 .0006520953 520. 0. MG: Richardson, r = 3, m1 = 2, m2 = 2, p = 2. .0000004041 .0005734469 772. 0. MG: Richardson, r = 4, m1 = 2, m2 = 2, p = 2. .0000002919 .0005399459 1032. 0. MG: Richardson, r = 1, m1 = 3, m2 = 3, p = 2. .0000014631 .0007428193 348. 0. MG: Richardson, r = 2, m1 = 3, m2 = 3, p = 2. .0000004170 .0005734917 696. 0. MG: Richardson, r = 3, m1 = 3, m2 = 3, p = 2. .0000002703 .0005299011 1040. 0. MG: Richardson, r = 4, m1 = 3, m2 = 3, p = 2. .0000002281 .0005119073 1380. 0. MG: Richardson, r = 1, m1 = 4, m2 = 4, p = 2. .0000008535 .0006526798 440. 0. MG: Richardson, r = 2, m1 = 4, m2 = 4, p = 2. .0000003064 .0005399750 892. 0. MG: Richardson, r = 3, m1 = 4, m2 = 4, p = 2. .0000002307 .0005119096 1304. 2144. MG: Richardson, r = 4, m1 = 4, m2 = 4, p = 2. .0000002092 .0005007998 1736. 0. MG: Jacobi, r = 1, m1 = 1, m2 = 1, p = 1. .0006777496 .0011133256 96. 0. MG: Jacobi, r = 2, m1 = 1, m2 = 1, p = 1. .0000014328 .0008457385 184. 0. MG: Jacobi, r = 3, m1 = 1, m2 = 1, p = 1. .0000003364 .0008457321 276. 0. MG: Jacobi, r = 4, m1 = 1, m2 = 1, p = 1. .0000002054 .0008457305 368. 0. MG: Jacobi, r = 1, m1 = 2, m2 = 2, p = 1. .0000031664 .0008457349 140. 0. MG: Jacobi, r = 2, m1 = 2, m2 = 2, p = 1. .0000005035 .0008457190 300. 0. MG: Jacobi, r = 3, m1 = 2, m2 = 2, p = 1. .0000004575 .0008457170 416. 0. MG: Jacobi, r = 4, m1 = 2, m2 = 2, p = 1. .0000004527 .0008457152 552. 0. MG: Jacobi, r = 1, m1 = 3, m2 = 3, p = 1. .0000004070 .0008457322 188. 0. MG: Jacobi, r = 2, m1 = 3, m2 = 3, p = 1. .0000001905 .0008457278 376. 0. MG: Jacobi, r = 3, m1 = 3, m2 = 3, p = 1. .0000001873 .0008457252 556. 0. MG: Jacobi, r = 4, m1 = 3, m2 = 3, p = 1. .0000001869 .0008457228 736. 0. MG: Jacobi, r = 1, m1 = 4, m2 = 4, p = 1. .0000006088 .0008457204 236. 0. MG: Jacobi, r = 2, m1 = 4, m2 = 4, p = 1. .0000004587 .0008457159 464. 0. MG: Jacobi, r = 3, m1 = 4, m2 = 4, p = 1. .0000004520 .0008457129 724. 2136. MG: Jacobi, r = 4, m1 = 4, m2 = 4, p = 1. .0000004515 .0008457099 924. 0. MG: Jacobi, r = 1, m1 = 1, m2 = 1, p = 2. .0000001869 .0008458594 160. 0. MG: Jacobi, r = 2, m1 = 1, m2 = 1, p = 2. .0000001868 .0008457363 320. 0. MG: Jacobi, r = 3, m1 = 1, m2 = 1, p = 2. .0000001868 .0008457319 472. 0. MG: Jacobi, r = 4, m1 = 1, m2 = 1, p = 2. .0000001868 .0008457305 624. 0. MG: Jacobi, r = 1, m1 = 2, m2 = 2, p = 2. .0000004514 .0008457228 236. 0. MG: Jacobi, r = 2, m1 = 2, m2 = 2, p = 2. .0000004514 .0008457188 468. 0. MG: Jacobi, r = 3, m1 = 2, m2 = 2, p = 2. .0000004514 .0008457169 704. 2140. MG: Jacobi, r = 4, m1 = 2, m2 = 2, p = 2. .0000004514 .0008457152 928. 0. MG: Jacobi, r = 1, m1 = 3, m2 = 3, p = 2. .0000001868 .0008457303 316. 0. MG: Jacobi, r = 2, m1 = 3, m2 = 3, p = 2. .0000001868 .0008457277 624. 0. MG: Jacobi, r = 3, m1 = 3, m2 = 3, p = 2. .0000001868 .0008457252 924. 0. MG: Jacobi, r = 4, m1 = 3, m2 = 3, p = 2. .0000001868 .0008457228 1244. 0. MG: Jacobi, r = 1, m1 = 4, m2 = 4, p = 2. .0000004514 .0008457191 416. 0. MG: Jacobi, r = 2, m1 = 4, m2 = 4, p = 2. .0000004514 .0008457159 768. 0. MG: Jacobi, r = 3, m1 = 4, m2 = 4, p = 2. .0000004514 .0008457129 1152. 0. MG: Jacobi, r = 4, m1 = 4, m2 = 4, p = 2. .0000004514 .0008457099 1536. 0. MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 1. .0000005202 .0004884854 160. 0. MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 1. .0000002063 .0004882815 316. 0. MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 1. .0000001943 .0004882805 476. 0. MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 1. .0000001928 .0004882805 632. 0. MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 1. .0000002408 .0004882835 272. 0. MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 1. .0000001948 .0004882805 536. 0. MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 1. .0000001926 .0004882805 804. 0. MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 1. .0000001925 .0004882805 1068. 0. MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 1. .0000002185 .0004882810 388. 0. MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 1. .0000001934 .0004882805 768. 0. MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 1. .0000001925 .0004882805 1140. 0. MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 1. .0000001925 .0004882805 1524. 0. MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 1. .0000002107 .0004882807 496. 0. MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 1. .0000001930 .0004882805 988. 0. MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 1. .0000001925 .0004882805 1476. 0. MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 1. .0000001925 .0004882805 1972. 2124. MG: Gauss-Seidel, r = 1, m1 = 1, m2 = 1, p = 2. .0000002262 .0004884740 252. 0. MG: Gauss-Seidel, r = 2, m1 = 1, m2 = 1, p = 2. .0000001929 .0004882813 496. 0. MG: Gauss-Seidel, r = 3, m1 = 1, m2 = 1, p = 2. .0000001925 .0004882805 744. 0. MG: Gauss-Seidel, r = 4, m1 = 1, m2 = 1, p = 2. .0000001925 .0004882805 984. 0. MG: Gauss-Seidel, r = 1, m1 = 2, m2 = 2, p = 2. .0000001959 .0004882829 416. 0. MG: Gauss-Seidel, r = 2, m1 = 2, m2 = 2, p = 2. .0000001925 .0004882805 832. 0. MG: Gauss-Seidel, r = 3, m1 = 2, m2 = 2, p = 2. .0000001925 .0004882805 1244. 0. MG: Gauss-Seidel, r = 4, m1 = 2, m2 = 2, p = 2. .0000001925 .0004882805 1656. 0. MG: Gauss-Seidel, r = 1, m1 = 3, m2 = 3, p = 2. .0000001930 .0004882807 584. 0. MG: Gauss-Seidel, r = 2, m1 = 3, m2 = 3, p = 2. .0000001925 .0004882805 1164. 0. MG: Gauss-Seidel, r = 3, m1 = 3, m2 = 3, p = 2. .0000001925 .0004882805 1748. 0. MG: Gauss-Seidel, r = 4, m1 = 3, m2 = 3, p = 2. .0000001925 .0004882805 2328. 0. MG: Gauss-Seidel, r = 1, m1 = 4, m2 = 4, p = 2. .0000001926 .0004882806 748. 0. MG: Gauss-Seidel, r = 2, m1 = 4, m2 = 4, p = 2. .0000001925 .0004882805 1492. 0. MG: Gauss-Seidel, r = 3, m1 = 4, m2 = 4, p = 2. .0000001925 .0004882805 2260. 0. MG: Gauss-Seidel, r = 4, m1 = 4, m2 = 4, p = 2. .0000001925 .0004882805 2996. 0. >