Thanks, Jed. I used left Jacobi preconditioner (by setting M1 = diag(A)) with Matlab's QMR and BCGS. I am new to PETSc and I am not very fimilar with linear solvers. If I want to implement the QMR solver, what is the starting point?
Thanks in advance. Best Regards, Ce Qin
