On Mon, Jan 14, 2013 at 3:27 PM, Tabrez Ali <stali at geology.wisc.edu> wrote:
> Hello > > I am solving a system of equations of the form: > > |A C| |u1| = |f1| > |C'B| |u2| |f2| > > After each solve, I need to perform B*f2 before updating f. Should I use > MatGetSubMatrix/**VecGetSubVector followed by MatMult or is there > something simpler. > Yes, or let PCFIELDSPLIT do all the block solver stuff for you. -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130114/bf660934/attachment.html>
