<div dir="ltr"><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">I have a linear system which I can easily solve using the Pardiso library, which is a direct solver. But when I try Petsc, the solver does not converge. It reduces the norm a little but then gets stuck.</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">I tried several combinations of pc_type and ksp_type. It seems that pc=jacobi is better, since my linear system does not have the diagonal element in some rows.</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Is there some advice you can give me? Some way to analyse the matrix and find out why it does not work?</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br></div></div>