Hi all,<br>I did a numerical experiment of ksp solve on a large sparse matrix.<br><br>For the first step, I load the matrix from petsc readable files. If I load the matrix as a MPIAIJ matrix and then apply the prometheus preconditioner with GMRES, it will not converge in 10000 iterations.<br>
<br>However, if I load the matrix as a MPIBAIJ matrix and then apply the prometheus preconditioner  with GEMRES, it will converge with in 800 iterations<br><br>Can anyone please provide some hints on how to explain this numerical result.<br>
<br>Thanks a lot,   <br><br>Yan<br>