On Tue, Mar 15, 2011 at 6:00 PM, Jed Brown <span dir="ltr"><<a href="mailto:jed@59a2.org">jed@59a2.org</a>></span> wrote:<br><div class="gmail_quote"><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 Tue, Mar 15, 2011 at 23:50, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">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">

That is not what I have here. I want to change a portion, but outside of that I want it to fall back to a global matrix.</blockquote></div><br></div><div>You were going to put your MatShell in the "pmat" slot, right? What are you going to use for preconditioning the other blocks? If it's the same as in "mat", you could use a MatNest for both operators and put the special sauce in that slot of the "pmat".</div>
</blockquote><div><br></div><div>This is the same objection Barry has. I do not want this tied to a particular matrix impl for the global matrix.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div>If you assemble the matrix entirely monolithically (using MatSetValues globally instead of MatSetValuesLocal on each physics separately), then we haven't written MatSetValues_Nest yet so maybe your MatShell is the right idea.</div>

<div><br></div><div><br></div><div>Related: What do people think about having MatGetSubMatrix() retain ownership of the submatrix?</div></blockquote><div><br></div><div>Does it matter with our interface? Should be transparent.</div>
<div><br></div><div>   Matt</div><div> </div></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>