[petsc-users] Field wise residual norm

[Matthew Knepley]

On Mon, Jan 13, 2014 at 4:11 PM, Dharmendar Reddy
<dharmareddy84 at gmail.com>wrote:

> Hello,
>         Is there a petsc function which can help accomplish the following:
> I am solving a coupled system of equations:
>
> F1(x1,x2,x3) = 0
>
> F2(x1,x2,x3) = 0
>
> F3(x1,x2,x3) = 0
>
>
> I am using DMPlex to handle to mesh data. I would like to compute the
> norm of residual of each function.
>
> I assemble the function element wise. I am looking for some interface like,
> ComputNorm(dm, F, fieldIndex).
>

Our thinking is generally that you would use

  DMCreateSubDM()

and then use VecScatterCreate() with the IS it returns. This will map F into
one of the F1/F2/F3.

Does this make sense?

   Matt


> Thanks
> Reddy
>



-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140113/25dec901/attachment-0001.html>