<div dir="ltr">On Mon, Jan 14, 2013 at 3:27 PM, Tabrez Ali <span dir="ltr"><<a href="mailto:stali@geology.wisc.edu" target="_blank">stali@geology.wisc.edu</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello<br>
<br>
I am solving a system of equations of the form:<br>
<br>
|A C| |u1| = |f1|<br>
|C'B| |u2|   |f2|<br>
<br>
After each solve, I need to perform B*f2 before updating f. Should I use MatGetSubMatrix/<u></u>VecGetSubVector followed by MatMult or is there something simpler.<br></blockquote><div><br></div><div style>Yes, or let PCFIELDSPLIT do all the block solver stuff for you.</div>
</div></div></div>