<br><br><div class="gmail_quote">On Fri, Feb 10, 2012 at 3:59 PM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

<div class="HOEnZb"><div class="h5">On Fri, Feb 10, 2012 at 3:58 PM, Dmitry Karpeev <span dir="ltr"><<a href="mailto:karpeev@mcs.anl.gov" target="_blank">karpeev@mcs.anl.gov</a>></span> wrote:<br></div></div><div class="gmail_quote">

<div><div class="h5"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div><div><br><br><div class="gmail_quote">On Fri, Feb 10, 2012 at 3:53 PM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@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">

<div><div class="gmail_quote">On Fri, Feb 10, 2012 at 15:50, Dmitry Karpeev <span dir="ltr"><<a href="mailto:karpeev@mcs.anl.gov" target="_blank">karpeev@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">





Why not simply allow DM to return two matrices as SNES and KSP would want anyway? </blockquote></div><br></div><div>How do you specify the reduced stencil to use for the preconditioning matrix? Maybe it wouldn't add too many more if statements to the current code which is already polluted with them.</div>





</blockquote></div><br></div></div><div>It seems to me that the sparsity pattern (and even the type of the matrices) is up to the particular DM type,</div><div>so it is, in principle, opaque.  DMDA can control it by manipulating stencil width or widths, while FEM-oriented DM </div>




<div>types (e.g., libMesh) could control it via the element type, etc. The generic DM interface doesn't have a way to specify </div><div>the sparsity pattern. As it should be, in my opinion.</div></blockquote><div><br>


</div></div></div><div>This makes no sense to me unless the two matrices being returned have different sparsity patterns, and we have no</div><div>good way for specifying that, or even a good rationale for doing it.</div>

</div></blockquote><div><br></div><div>The rationale is that the preconditioning matrix can have a different sparsity pattern. </div><div>L2 projection leads to a mass matrix, and the lumped mass matrix is a good preconditioner.</div>

<div>It is not DM's job to specify the sparsity pattern -- it is up to a particular implementation.</div><div><br></div><div>Dmitry. </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

<div class="gmail_quote"><div><br></div>
<div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span><font color="#888888"><div>
Dmitry.</div>
</font></span></blockquote></div><div class="HOEnZb"><div class="h5"><br><br clear="all"><div><br></div>-- <br>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<br>
</div></div></blockquote></div><br>