Thanks, Jed.<div><br></div><div>I used left Jacobi preconditioner (by setting M1 = diag(A)) with Matlab's QMR and BCGS. </div><div>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? </div>
<div><br></div><div>Thanks in advance.</div><div><br></div><div>Best Regards,</div><div>Ce Qin</div>