On Mon, Jan 13, 2014 at 4:11 PM, Dharmendar Reddy <[email protected]>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
