Add elasticity test. This is a nonlinear vector problem and it is large

enough that it takes a couple seconds to run. It uses the (now default) RCM ordering. Closes gh-3.
jedbrown · May 11, 2009 · 806d0f9 · 806d0f9
1 parent 0c089ca
commit 806d0f9
Show file tree

Hide file tree

Showing 2 changed files with 173 additions and 0 deletions.
diff --git a/src/fs/tests/CMakeLists.txt b/src/fs/tests/CMakeLists.txt
@@ -13,3 +13,4 @@ dohp_add_test (bu-0 1 bu -dmesh_in "${Dohp_DATA_DIR}/dblock211.h5m" -dfs_orderin
 dohp_add_test (ellip-e0-b2-p18 1 ellip -dmesh_in "${Dohp_DATA_DIR}/dblock4.h5m" -snes_monitor -ksp_monitor -pc_type lu -exact_a 1.6 -exact_b 1.7 -exact_c 1.8 -exact 0 -const_bdeg 2 -ellip_p 1.8 -dfs_ordering_type natural)
 dohp_add_test (ellip-e0-b4-p18 1 ellip -dmesh_in "${Dohp_DATA_DIR}/dblock4.h5m" -snes_monitor -ksp_monitor -pc_type lu -exact_a 1.6 -exact_b 1.7 -exact_c 1.8 -exact 0 -const_bdeg 4 -ellip_p 1.8 -djac_tensor_no_unroll -dfs_ordering_type natural)
 dohp_add_test (ellip-e0-b4-p17-sse 1 ellip -dmesh_in "${Dohp_DATA_DIR}/dblock4.h5m" -snes_monitor -ksp_monitor -pc_type lu -exact_a 1.6 -exact_b 1.7 -exact_c 1.8 -exact 0 -const_bdeg 4 -ellip_p 1.7 -dfs_ordering_type natural)
+dohp_add_test (elast-e0-b4-g01 1 elast -dmesh_in "${Dohp_DATA_DIR}/dblock4.h5m" -snes_monitor_short -ksp_monitor_short -pc_type ilu -exact_a 0.98 -exact_b 1.01 -exact_c 1.01 -exact 0 -const_bdeg 4 -elast_gamma 0.1)
diff --git a/src/fs/tests/refout/elast-e0-b4-g01.refout b/src/fs/tests/refout/elast-e0-b4-g01.refout
@@ -0,0 +1,172 @@
+  0 SNES Function norm 2068.06 
+    0 KSP Residual norm 26.2867 
+    1 KSP Residual norm 17.0493 
+    2 KSP Residual norm 10.1709 
+    3 KSP Residual norm 6.00881 
+    4 KSP Residual norm 3.73965 
+    5 KSP Residual norm 2.16856 
+    6 KSP Residual norm 1.28789 
+    7 KSP Residual norm 0.804729 
+    8 KSP Residual norm 0.501676 
+    9 KSP Residual norm 0.338871 
+   10 KSP Residual norm 0.226614 
+   11 KSP Residual norm 0.149809 
+   12 KSP Residual norm 0.0954723 
+   13 KSP Residual norm 0.0570307 
+   14 KSP Residual norm 0.0347478 
+   15 KSP Residual norm 0.0209891 
+   16 KSP Residual norm 0.011769 
+   17 KSP Residual norm 0.00627491 
+   18 KSP Residual norm 0.00363572 
+   19 KSP Residual norm 0.00219991 
+   20 KSP Residual norm 0.00137035 
+   21 KSP Residual norm 0.000861783 
+   22 KSP Residual norm 0.000538043 
+   23 KSP Residual norm 0.00034024 
+   24 KSP Residual norm 0.000209286 
+  1 SNES Function norm 1830.73 
+    0 KSP Residual norm 22.1006 
+    1 KSP Residual norm 12.5316 
+    2 KSP Residual norm 6.66328 
+    3 KSP Residual norm 4.0004 
+    4 KSP Residual norm 2.51057 
+    5 KSP Residual norm 1.5194 
+    6 KSP Residual norm 0.827732 
+    7 KSP Residual norm 0.488091 
+    8 KSP Residual norm 0.291075 
+    9 KSP Residual norm 0.180007 
+   10 KSP Residual norm 0.118989 
+   11 KSP Residual norm 0.0787991 
+   12 KSP Residual norm 0.0507496 
+   13 KSP Residual norm 0.0329098 
+   14 KSP Residual norm 0.0211721 
+   15 KSP Residual norm 0.0131877 
+   16 KSP Residual norm 0.0077054 
+   17 KSP Residual norm 0.00434589 
+   18 KSP Residual norm 0.00265317 
+   19 KSP Residual norm 0.00152053 
+   20 KSP Residual norm 0.00089518 
+   21 KSP Residual norm 0.000535795 
+   22 KSP Residual norm 0.000320573 
+   23 KSP Residual norm 0.000201162 
+  2 SNES Function norm 937.526 
+    0 KSP Residual norm 9.48704 
+    1 KSP Residual norm 4.64943 
+    2 KSP Residual norm 2.37771 
+    3 KSP Residual norm 1.42404 
+    4 KSP Residual norm 0.743332 
+    5 KSP Residual norm 0.465683 
+    6 KSP Residual norm 0.266948 
+    7 KSP Residual norm 0.185781 
+    8 KSP Residual norm 0.12421 
+    9 KSP Residual norm 0.0764301 
+   10 KSP Residual norm 0.0454382 
+   11 KSP Residual norm 0.0253821 
+   12 KSP Residual norm 0.0125302 
+   13 KSP Residual norm 0.00696667 
+   14 KSP Residual norm 0.00408337 
+   15 KSP Residual norm 0.00247332 
+   16 KSP Residual norm 0.00157885 
+   17 KSP Residual norm 0.00098719 
+   18 KSP Residual norm 0.000632389 
+   19 KSP Residual norm 0.000393504 
+   20 KSP Residual norm 0.000228992 
+   21 KSP Residual norm 0.00012962 
+   22 KSP Residual norm 7.58914e-05 
+  3 SNES Function norm 126.903 
+    0 KSP Residual norm 1.59157 
+    1 KSP Residual norm 0.726655 
+    2 KSP Residual norm 0.352525 
+    3 KSP Residual norm 0.210034 
+    4 KSP Residual norm 0.118838 
+    5 KSP Residual norm 0.0739782 
+    6 KSP Residual norm 0.0449744 
+    7 KSP Residual norm 0.0302518 
+    8 KSP Residual norm 0.018787 
+    9 KSP Residual norm 0.0114696 
+   10 KSP Residual norm 0.00651417 
+   11 KSP Residual norm 0.00346173 
+   12 KSP Residual norm 0.00158909 
+   13 KSP Residual norm 0.000889244 
+   14 KSP Residual norm 0.000544063 
+   15 KSP Residual norm 0.000316031 
+   16 KSP Residual norm 0.000197128 
+   17 KSP Residual norm 0.000123053 
+   18 KSP Residual norm 7.37743e-05 
+   19 KSP Residual norm 4.54079e-05 
+   20 KSP Residual norm 2.80028e-05 
+   21 KSP Residual norm 1.54982e-05 
+  4 SNES Function norm 6.30591 
+    0 KSP Residual norm 0.0873773 
+    1 KSP Residual norm 0.0383946 
+    2 KSP Residual norm 0.0175348 
+    3 KSP Residual norm 0.0101397 
+    4 KSP Residual norm 0.00594617 
+    5 KSP Residual norm 0.00359432 
+    6 KSP Residual norm 0.00216183 
+    7 KSP Residual norm 0.00144617 
+    8 KSP Residual norm 0.000930069 
+    9 KSP Residual norm 0.0005881 
+   10 KSP Residual norm 0.000331411 
+   11 KSP Residual norm 0.00018192 
+   12 KSP Residual norm 8.70036e-05 
+   13 KSP Residual norm 4.72141e-05 
+   14 KSP Residual norm 2.80055e-05 
+   15 KSP Residual norm 1.37371e-05 
+   16 KSP Residual norm 7.25828e-06 
+   17 KSP Residual norm 4.02833e-06 
+   18 KSP Residual norm 2.28195e-06 
+   19 KSP Residual norm 1.39259e-06 
+   20 KSP Residual norm 9.10399e-07 
+   21 KSP Residual norm 5.78269e-07 
+  5 SNES Function norm 0.0261414 
+    0 KSP Residual norm 0.000376452 
+    1 KSP Residual norm 0.00016355 
+    2 KSP Residual norm 7.25982e-05 
+    3 KSP Residual norm 4.10521e-05 
+    4 KSP Residual norm 2.38329e-05 
+    5 KSP Residual norm 1.42496e-05 
+    6 KSP Residual norm 8.5265e-06 
+    7 KSP Residual norm 5.6316e-06 
+    8 KSP Residual norm 3.7122e-06 
+    9 KSP Residual norm 2.38538e-06 
+   10 KSP Residual norm 1.3687e-06 
+   11 KSP Residual norm 7.77879e-07 
+   12 KSP Residual norm 3.91131e-07 
+   13 KSP Residual norm 2.13147e-07 
+   14 KSP Residual norm 1.25473e-07 
+   15 KSP Residual norm 5.80952e-08 
+   16 KSP Residual norm 2.90278e-08 
+   17 KSP Residual norm 1.54276e-08 
+   18 KSP Residual norm 8.37984e-09 
+   19 KSP Residual norm 4.57782e-09 
+   20 KSP Residual norm 2.61206e-09 
+  6 SNES Function norm 5.86076e-07 
+    0 KSP Residual norm 8.8121e-09 
+    1 KSP Residual norm 4.08074e-09 
+    2 KSP Residual norm 2.10258e-09 
+    3 KSP Residual norm 1.35597e-09 
+    4 KSP Residual norm 8.448e-10 
+    5 KSP Residual norm 5.167e-10 
+    6 KSP Residual norm 2.985e-10 
+    7 KSP Residual norm 1.822e-10 
+    8 KSP Residual norm 1.149e-10 
+    9 KSP Residual norm 7.366e-11 
+   10 KSP Residual norm 4.618e-11 
+   11 KSP Residual norm 2.423e-11 
+   12 KSP Residual norm 1.310e-11 
+   13 KSP Residual norm < 1.e-11
+   14 KSP Residual norm < 1.e-11
+   15 KSP Residual norm < 1.e-11
+   16 KSP Residual norm < 1.e-11
+   17 KSP Residual norm < 1.e-11
+   18 KSP Residual norm < 1.e-11
+   19 KSP Residual norm < 1.e-11
+   20 KSP Residual norm < 1.e-11
+   21 KSP Residual norm < 1.e-11
+   22 KSP Residual norm < 1.e-11
+  7 SNES Function norm < 1.e-11
+Algebraic residual        |x|_1 1.52e-10  |x|_2 3.80e-12  |x|_inf 4.48e-13
+Interpolation residual    |x|_1 2.11e+01  |x|_2 4.75e-01  |x|_inf 2.55e-02
+Pointwise solution error  |x|_1 1.08e-01  |x|_2 1.50e-01  |x|_inf 2.56e-02
+Pointwise gradient error  |x|_1 5.20e-01  |x|_2 1.50e-01  |x|_inf 1.86e-01