Hi all, I did a numerical experiment of ksp solve on a large sparse matrix.
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. 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 Can anyone please provide some hints on how to explain this numerical result. Thanks a lot, Yan -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20090521/ff85e998/attachment.htm>
