Inhomogeneous Dirichlet implemented for ellip. Closes gh-13.

jedbrown · May 7, 2009 · ac86dc9 · ac86dc9
1 parent 4a194b0
commit ac86dc9
Show file tree

Hide file tree

Showing 4 changed files with 227 additions and 51 deletions.
diff --git a/src/fs/tests/CMakeLists.txt b/src/fs/tests/CMakeLists.txt
@@ -10,3 +10,5 @@ add_test (fs-mf-op-proj-1  fs-ex1 -const_BDeg 7 -snes_mf_operator -ksp_converged
   -ksp_rtol 1e-10 -snes_monitor -frequency 2,2,2 -require_ptwise 2e-6,2e-6,9e-6 -require_grad 3e-6,7e-6,1.1e-4)
 
 dohp_add_test (bu-0 1 bu -dmesh_in "${Dohp_DATA_DIR}/dblock211.h5m")
+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)
+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)
diff --git a/src/fs/tests/ellip-e0-b2-p18.refout b/src/fs/tests/ellip-e0-b2-p18.refout
@@ -0,0 +1,17 @@
+  0 SNES Function norm 3.178176880421e+00 
+    0 KSP Residual norm 3.758510964659e+00 
+    1 KSP Residual norm 9.314056362651e-16 
+  1 SNES Function norm 1.521565617543e-01 
+    0 KSP Residual norm 1.994655602452e-01 
+    1 KSP Residual norm 4.814133587197e-17 
+  2 SNES Function norm 9.282240955142e-04 
+    0 KSP Residual norm 1.160261666396e-03 
+    1 KSP Residual norm 2.523918178477e-19 
+  3 SNES Function norm 4.248606159251e-08 
+    0 KSP Residual norm 5.059845889641e-08 
+    1 KSP Residual norm 1.468450792022e-23 
+  4 SNES Function norm 2.619432448725e-16 
+Algebraic residual        |x|_1 9.21e-16  |x|_2 2.62e-16  |x|_inf 1.67e-16
+Interpolation residual    |x|_1 2.43e-01  |x|_2 6.74e-02  |x|_inf 3.37e-02
+Pointwise solution error  |x|_1 3.80e-01  |x|_2 2.13e-01  |x|_inf 3.03e-01
+Pointwise gradient error  |x|_1 8.22e+00  |x|_2 7.19e+00  |x|_inf 1.59e+01
diff --git a/src/fs/tests/ellip-e0-b4-p18.refout b/src/fs/tests/ellip-e0-b4-p18.refout
@@ -0,0 +1,106 @@
+  0 SNES Function norm 4.934076170102e+00 
+    0 KSP Residual norm 5.170589949553e+01 
+    1 KSP Residual norm 6.792085968900e+00 
+    2 KSP Residual norm 2.786298241160e+00 
+    3 KSP Residual norm 1.506873739231e+00 
+    4 KSP Residual norm 7.477740709955e-01 
+    5 KSP Residual norm 3.748422051368e-01 
+    6 KSP Residual norm 1.741413878933e-01 
+    7 KSP Residual norm 7.840855423625e-02 
+    8 KSP Residual norm 3.894880717488e-02 
+    9 KSP Residual norm 1.938885736179e-02 
+   10 KSP Residual norm 1.005068812857e-02 
+   11 KSP Residual norm 6.108519793533e-03 
+   12 KSP Residual norm 3.430909619328e-03 
+   13 KSP Residual norm 1.902656749805e-03 
+   14 KSP Residual norm 1.027709993021e-03 
+   15 KSP Residual norm 5.396388657937e-04 
+   16 KSP Residual norm 2.928232176611e-04 
+  1 SNES Function norm 2.624947463465e+00 
+    0 KSP Residual norm 1.637128709669e+01 
+    1 KSP Residual norm 6.346260152425e+00 
+    2 KSP Residual norm 2.335832818640e+00 
+    3 KSP Residual norm 9.988470168762e-01 
+    4 KSP Residual norm 3.922788940362e-01 
+    5 KSP Residual norm 1.574277094450e-01 
+    6 KSP Residual norm 4.967489580574e-02 
+    7 KSP Residual norm 2.067656129127e-02 
+    8 KSP Residual norm 1.142549725147e-02 
+    9 KSP Residual norm 6.797256062821e-03 
+   10 KSP Residual norm 2.942851402326e-03 
+   11 KSP Residual norm 1.257297637599e-03 
+   12 KSP Residual norm 5.210047321609e-04 
+   13 KSP Residual norm 2.133245993830e-04 
+   14 KSP Residual norm 9.093411575835e-05 
+  2 SNES Function norm 1.087877916987e+00 
+    0 KSP Residual norm 5.390089094232e+00 
+    1 KSP Residual norm 1.802492093096e+00 
+    2 KSP Residual norm 6.283436802896e-01 
+    3 KSP Residual norm 2.654992552903e-01 
+    4 KSP Residual norm 1.009967053549e-01 
+    5 KSP Residual norm 3.571516348950e-02 
+    6 KSP Residual norm 1.247488955373e-02 
+    7 KSP Residual norm 4.797796170408e-03 
+    8 KSP Residual norm 2.291401856793e-03 
+    9 KSP Residual norm 1.253935366293e-03 
+   10 KSP Residual norm 6.224267349908e-04 
+   11 KSP Residual norm 2.854660315134e-04 
+   12 KSP Residual norm 1.425610617966e-04 
+   13 KSP Residual norm 6.540958392889e-05 
+   14 KSP Residual norm 2.795000878944e-05 
+  3 SNES Function norm 7.562338860336e-02 
+    0 KSP Residual norm 4.618062566565e-01 
+    1 KSP Residual norm 2.254873333412e-01 
+    2 KSP Residual norm 8.903748081300e-02 
+    3 KSP Residual norm 3.592112294645e-02 
+    4 KSP Residual norm 1.219649001148e-02 
+    5 KSP Residual norm 4.298024860438e-03 
+    6 KSP Residual norm 1.964384619717e-03 
+    7 KSP Residual norm 1.020155978482e-03 
+    8 KSP Residual norm 5.709073792501e-04 
+    9 KSP Residual norm 2.829194665885e-04 
+   10 KSP Residual norm 8.648123751309e-05 
+   11 KSP Residual norm 4.169914694815e-05 
+   12 KSP Residual norm 2.141240381749e-05 
+   13 KSP Residual norm 1.056903254502e-05 
+   14 KSP Residual norm 5.430687032717e-06 
+   15 KSP Residual norm 2.914213301318e-06 
+  4 SNES Function norm 4.624677605888e-04 
+    0 KSP Residual norm 2.470130476938e-03 
+    1 KSP Residual norm 1.245657318218e-03 
+    2 KSP Residual norm 4.712429489391e-04 
+    3 KSP Residual norm 1.741776319481e-04 
+    4 KSP Residual norm 6.615055847631e-05 
+    5 KSP Residual norm 2.408722536514e-05 
+    6 KSP Residual norm 1.017132004301e-05 
+    7 KSP Residual norm 4.014648325794e-06 
+    8 KSP Residual norm 1.863252697179e-06 
+    9 KSP Residual norm 1.037340802489e-06 
+   10 KSP Residual norm 5.391761676595e-07 
+   11 KSP Residual norm 2.783812663564e-07 
+   12 KSP Residual norm 1.112437504610e-07 
+   13 KSP Residual norm 4.428540930041e-08 
+   14 KSP Residual norm 2.045644517912e-08 
+  5 SNES Function norm 3.337073115677e-08 
+    0 KSP Residual norm 1.667420339594e-07 
+    1 KSP Residual norm 7.629137559881e-08 
+    2 KSP Residual norm 3.189920926387e-08 
+    3 KSP Residual norm 1.337282536996e-08 
+    4 KSP Residual norm 5.284115967178e-09 
+    5 KSP Residual norm 2.708090174828e-09 
+    6 KSP Residual norm 1.681878531801e-09 
+    7 KSP Residual norm 9.450272554581e-10 
+    8 KSP Residual norm 4.250031335518e-10 
+    9 KSP Residual norm 1.778078484098e-10 
+   10 KSP Residual norm 7.143399530159e-11 
+   11 KSP Residual norm 4.171459510235e-11 
+   12 KSP Residual norm 2.114921003974e-11 
+   13 KSP Residual norm 1.118465915671e-11 
+   14 KSP Residual norm 5.271654559672e-12 
+   15 KSP Residual norm 2.088132859954e-12 
+   16 KSP Residual norm 8.477479183515e-13 
+  6 SNES Function norm 2.592730097803e-13 
+Algebraic residual        |x|_1 6.18e-12  |x|_2 2.59e-13  |x|_inf 4.08e-14
+Interpolation residual    |x|_1 3.13e-02  |x|_2 1.57e-03  |x|_inf 2.55e-04
+Pointwise solution error  |x|_1 1.13e-03  |x|_2 7.50e-04  |x|_inf 2.66e-03
+Pointwise gradient error  |x|_1 4.85e+00  |x|_2 5.11e+00  |x|_inf 1.70e+01