<div dir="ltr"><div>We fixed the issue. The linear operator in the matrix-free solver was set up correctly, whereas the linear operator in the preconditioner was not set up correctly (it was lagging in time). After setting the linear operators correctly, we see that the diagonal system is solved in 1 iteration. Thank you.<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Sep 30, 2021 at 5:49 PM Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> wrote:<br></div><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 dir="ltr"><div dir="ltr">On Thu, Sep 30, 2021 at 6:58 PM Amneet Bhalla <<a href="mailto:mail2amneet@gmail.com" target="_blank">mail2amneet@gmail.com</a>> wrote:<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 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><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>
</blockquote></div><br clear="all"><div>Can you send us the matrix? This definitely should converge in 1 iterate now, so something I do not understand is going on.</div><div>I will take any format you've got :)</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div><br></div>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>
</blockquote></div>