<div dir="ltr"><div>Richardson damping factor does not seem to do anything.  I set it to zero and the solve was fine.</div><div><br></div><div>Am I missing something?</div><div><br></div><div>Thanks,</div><div>Mark</div><div><br></div><div>  Down solver (pre-smoother) on level 3 -------------------------------</div><div>    KSP Object:    (mg_levels_3_)     1024 MPI processes</div><div>      type: richardson</div><div>        Richardson: damping factor=0</div><div>      maximum iterations=4</div><div>      tolerances:  relative=1e-05, absolute=1e-50, divergence=10000</div><div>      left preconditioning</div><div>      using nonzero initial guess</div><div>      using NONE norm type for convergence test</div><div>    PC Object:    (mg_levels_3_)     1024 MPI processes</div><div>      type: sor</div><div>        SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1</div><div>      linear system matrix followed by preconditioner matrix:</div><div>      Mat Object:       1024 MPI processes</div><div>        type: shell</div><div>        rows=33427584, cols=33427584</div><div>      Mat Object:       1024 MPI processes</div><div>        type: mpiaij</div><div>        rows=33427584, cols=33427584</div><div>        total: nonzeros=1.7272e+09, allocated nonzeros=3.4544e+09</div><div>        total number of mallocs used during MatSetValues calls =0</div><div>          not using I-node (on process 0) routines</div><div><br></div></div>