PS D:\Computational_geomechanics\geoalg_64_barry_final> mpiexec -n 1 defmod.exe -f test/test.inp -ksp_monitor -ksp_converged_reason -ksp_view Reading input ... Reading mesh data ... Forming [K] ... Memory used: 1268928.0000000000 10 % complete. 20 % complete. 30 % complete. 40 % complete. 50 % complete. 60 % complete. 70 % complete. 80 % complete. 90 % complete. 100 % complete. Applying constraints ... Memory used: 170643152.00000000 10 % complete 20 % complete 30 % complete 40 % complete 50 % complete 60 % complete 70 % complete 80 % complete 90 % complete 100 % complete Forming RHS ... Memory used: 171862896.00000000 Setting up solver ... Solving ... Solving ... Memory used: 176084896.00000000 Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0 PC_FAILED due to FACTOR_NUMERIC_ZEROPIVOT KSP Object: 1 MPI processes type: gmres restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement happy breakdown tolerance 1e-30 maximum iterations=100000, initial guess is zero tolerances: relative=1e-08, absolute=1e-50, divergence=10000. left preconditioning using PRECONDITIONED norm type for convergence test PC Object: 1 MPI processes type: sor type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. linear system matrix = precond matrix: Mat Object: 1 MPI processes type: seqaij rows=152466, cols=152466 total: nonzeros=11841674, allocated nonzeros=13907097 total number of mallocs used during MatSetValues calls=0 using I-node routines: found 93030 nodes, limit used is 5 Recovering stress ... Cleaning up ... Finished