<div class="gmail_quote">On Thu, Feb 9, 2012 at 11:56 PM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov">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 class="im"><div class="gmail_quote">On Thu, Feb 9, 2012 at 23:51, 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">


Good point.  It is still strange, however, that SNES insists on using the same Mat obtained from DMGetMatrix for J and P. </blockquote></div><br></div><div>Oh, you're concerned about when SNES creates a Mat automatically (getting it out of the DM)?</div>

</blockquote><div> </div><div>Yes. And, by the way, preallocation for J and P could be different (e.g., mass matrix vs lumped mass matrix for projection problems).  As could be the preferred MatType.</div><div><br></div>
<div>
Dmitry.</div></div><br>