Additional weird behavior: If I run with
-ksp_monitor -ksp_monitor_true_residual -ksp_converged_reason -pc_type
jacobi
I get:
0 KSP Residual norm 1.389113910577e-03
0 KSP preconditioned resid norm 1.389113910577e-03 true resid norm
-nan ||r(i)||/||b|| -nan
1 KSP Residual norm 2.214853823259e-19
1 KSP preconditioned resid norm 2.214853823259e-19 true resid norm
-nan ||r(i)||/||b|| -nan
Linear solve converged due to CONVERGED_RTOL iterations 1
