Update test output for elasticity.

This is still using a St. Venant-Kirchoff constitutive relation which goes haywire for moderate compression. One effect is that the resulting linear systems lose positive definiteness. The case simulated here is slightly beyond the physical regime and indeed, the system becomes indefinite partway through the nonlinear solve.
jedbrown · Apr 21, 2011 · cef7e4d · cef7e4d
1 parent 4787838
commit cef7e4d
Show file tree

Hide file tree

Showing 2 changed files with 345 additions and 168 deletions.
diff --git a/src/fs/tests/CMakeLists.txt b/src/fs/tests/CMakeLists.txt
@@ -34,4 +34,4 @@ dohp_add_test (ellip-e0-b3-p16 1 ellip -dmesh_in "${Dohp_DATA_DIR}/dblock4.h5m"
 dohp_add_test (ellip-e0-b3-p14-sse 1 ellip -dmesh_in "${Dohp_DATA_DIR}/dblock4.h5m" -snes_monitor_short -ksp_monitor_short -pc_type lu -exact_a 1.6 -exact_b 1.7 -exact_c 1.8 -exact 0 -const_bdeg 3 -ellip_p 1.4)
 dohp_add_test (ellip-e0-b3-p17-lam1-sse 1 ellip -dmesh_in "${Dohp_DATA_DIR}/dblock4.h5m" -snes_monitor_short -ksp_monitor_short -pc_type lu -exact 0 -const_bdeg 3 -ellip_p 1.7 -ellip_lam 1)
 dohp_add_test (ellip-morph 1 ellip -dmesh_in "${Dohp_DATA_DIR}/dblock4.h5m" -morph 1 -twist 0.4 -ksp_monitor_short -snes_max_it 1 -ksp_rtol 1e-8 -ksp_type cg -pc_type lu)
-dohp_add_test (elast-e0-b3-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 3 -elast_gamma 0.1)
+dohp_add_test (elast-e0-b3-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 3 -elast_scale 0.1 -dfs_ordering_type natural)