Add elasticity example

* The shell Jacobian and preconditioner (for linear elements) is correct even with the nonlinear strain tensor. * Spectral convergence to manufactured solution is observed for a single element. * There is some error with high order and multiple elements. This needs to be fixed before implementing an indefinite vector problem (Stokes). References gh-3
jedbrown · May 9, 2009 · fd8b403 · fd8b403
1 parent bb227fc
commit fd8b403
Show file tree

Hide file tree

Showing 2 changed files with 592 additions and 0 deletions.
diff --git a/src/fs/tests/CMakeLists.txt b/src/fs/tests/CMakeLists.txt
@@ -2,6 +2,7 @@ dohp_link_executable (fs-ex1 ex1.c)
 dohp_link_executable (ellip ellip.c)
 # dohp_link_executable (cellip cellip.c)
 dohp_link_executable (bu bu.c)
+dohp_link_executable (elast elast.c)
 
 dohp_add_test (fs-mf-projection 1 fs-ex1 -snes_mf -ksp_type minres -const_BDeg 6 -snes_monitor -ksp_converged_reason)
 dohp_add_test (fs-mf-op-proj-0  1 fs-ex1 -snes_mf_operator -ksp_type minres -pc_type jacobi -const_BDeg 10 -snes_monitor -ksp_converged_reason -ksp_rtol 1e-10 -frequency 5,4,3 -require_ptwise 2e-6,3e-6,2.3e-5 -require_grad 7e-7,3e-6,2e-4)