<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, Jan 13, 2014 at 4:11 PM, Dharmendar Reddy <span dir="ltr"><<a href="mailto:dharmareddy84@gmail.com" target="_blank">dharmareddy84@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello,<br>
        Is there a petsc function which can help accomplish the following:<br>
I am solving a coupled system of equations:<br>
<br>
F1(x1,x2,x3) = 0<br>
<br>
F2(x1,x2,x3) = 0<br>
<br>
F3(x1,x2,x3) = 0<br>
<br>
<br>
I am using DMPlex to handle to mesh data. I would like to compute the<br>
norm of residual of each function.<br>
<br>
I assemble the function element wise. I am looking for some interface like,<br>
ComputNorm(dm, F, fieldIndex).<br></blockquote><div><br></div><div>Our thinking is generally that you would use</div><div><br></div><div>  DMCreateSubDM()</div><div><br></div><div>and then use VecScatterCreate() with the IS it returns. This will map F into</div>
<div>one of the F1/F2/F3.</div><div><br></div><div>Does this make sense?</div><div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

Thanks<br>
<span class="HOEnZb"><font color="#888888">Reddy<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</div></div>