[petsc-users] submatrix times subvector

[Jed Brown]

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>