<div dir="ltr"><div dir="ltr"><br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div style="word-wrap:break-word;line-break:after-white-space"><div><br><div><br></div>   For a diagonal system with this modest range of values Jacobi should converge in a single iteration. </div></div></blockquote><div><br></div><div>This is what I wanted to confirm (and my expectation also). There could be a bug in the way we are setting up the linear operators in the preconditioner and the matrix-free solver. We need to do some debugging.  </div><div><br></div><div>(with regard to the diagonal).<br></div><div><br></div><div>We have printed the matrix and viewed it in Matlab. It is a diagonal matrix. </div></div></div>