<div class="gmail_quote">On Wed, Oct 27, 2010 at 12:27, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Right, but how do you know what values you should put on the diagonal?  </blockquote><div><br></div><div>You have to look at the mesh and/or compute coefficients on the neighboring elements.  This is one reason I never use MatZeroRows, but it's possible for someone to have the information available.</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Should there be an option to just leave the value that is currently there?</blockquote><div><br></div><div>
Maybe, I was just pointing out the case that I assumed Lisandro was concerned with when he suggested it.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
 Should we change it so we always just leave the value that is already there?</blockquote></div><br><div>That would seem quite surprising to me.  They need the matrix scaling to be consistent with however they evaluate residuals, and it seems like an awkward bit of work to compute the exact matrix diagonal during residual evaluation.</div>
<div><br></div><div>Jed</div>