<div class="gmail_quote">On Tue, Mar 15, 2011 at 23:50, 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;">
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>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>
<div><br></div><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>