Poor residuals from PETSc iterative solver for 3D elasticity equation

General Discussion
2

Please consider:

  • using a lower -ksp_rtol, e.g., -ksp_rtol 1e-8 (default is 1e-5);
  • setting bs = 3 in your set(AFEA); command
  • switching to -ksp_pc_side right (default is left).

Does that give you a good enough solution?