<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Mar 26, 2015 at 10:36 AM, Ataollah Mesgarnejad <span dir="ltr"><<a href="mailto:amesga1@tigers.lsu.edu" target="_blank">amesga1@tigers.lsu.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr">Dear all,<div><br></div><div>I was wondering if someone can tell me how you get a field's Vec (subVec of the global Vec) from a distributed DMPlex. It seems that <span style="font-family:Helvetica;font-size:12px">DMCreateFieldIS works for sequential DMPlex but gives incorrect values for a distributed DMPlex. <br></span></div></div></blockquote><div><br></div><div>It seems like what you want is:</div><div><br></div><div> <a href="https://bitbucket.org/petsc/petsc/src/fdac1df83b6d73782209c263cf9d6b0b33600136/src/vec/vec/utils/vsection.c?at=master#cl-229">https://bitbucket.org/petsc/petsc/src/fdac1df83b6d73782209c263cf9d6b0b33600136/src/vec/vec/utils/vsection.c?at=master#cl-229</a></div><div><br></div><div>and also it seems like DMCreateFieldIS() is broken somehow. We should fix it. It currently passes all the tests for the FieldSplit</div><div>solvers, so I am not sure what is going wrong here.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr"><div><span style="font-family:Helvetica;font-size:12px"></span>Many thanks,
</div><div>Ata</div></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">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></div>