<div dir="ltr"><div dir="ltr">On Thu, Aug 8, 2019 at 11:59 AM Daniel Mckinnell via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div style="color:black;font:10pt arial">
<div style="font-family:arial,helvetica;font-size:10pt;color:black">
<div id="gmail-m_-3115044098188067955yiv2247701322">

<div style="color:black;font:10pt arial"><font size="2">Hi,</font></div>

<div style="color:black;font:10pt arial"><font size="2">I have a DMPlex mesh stored in a dm object with an 
associated PetscSection containing the information on the layout of 
field0 and field1, I can create a local Vec with the layout from the 
PetscSection and edit that data. I have a couple of questions about what I can do from this point:<br>
</font></div>

<div style="color:black;font:10pt arial"><font size="2"><br>
</font></div>

<div style="color:black;font:10pt arial"><font size="2">1. Is there a way of attaching my field data to the dm object so that everything associated with my mesh is contained in a single object?</font></div></div></div></div></blockquote><div><br></div><div>You can use</div><div><br></div><div>  <a href="https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/DM/DMGetNamedLocalVector.html">https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/DM/DMGetNamedLocalVector.html</a></div><div><br></div><div>We tend to discourage this organization, but it might work for you.  If you use DMCreateLocalVector(), then</div><div>VecGetDM() will get you back the DM.</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="color:black;font:10pt arial"><div style="font-family:arial,helvetica;font-size:10pt;color:black"><div id="gmail-m_-3115044098188067955yiv2247701322">

<div style="color:black;font:10pt arial"><font size="2">2. Is there a way of creating a Vec object with the data from just field0?</font></div></div></div></div></blockquote><div><br></div><div>Yes, use DMCreateSubDM() selecting only field0, and then DMCreateLocalVector()</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="color:black;font:10pt arial"><div style="font-family:arial,helvetica;font-size:10pt;color:black"><div id="gmail-m_-3115044098188067955yiv2247701322">

<div style="color:black;font:10pt arial"><font size="2">3. What is the best way to export my data to a VTK file, everything I have tried so far has just created meaningless vector fields on the vertices?</font></div></div></div></div></blockquote><div><br></div><div>VTK files only understand vertex fields and cell fields. Is that what you have?</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:1px solid rgb(204,204,204);padding-left:1ex"><div style="color:black;font:10pt arial"><div style="font-family:arial,helvetica;font-size:10pt;color:black"><div id="gmail-m_-3115044098188067955yiv2247701322">

<div style="color:black;font:10pt arial"><font size="2">Any help would be greatly appreciated.</font></div>

<div style="color:black;font:10pt arial"><font size="2">Thanks,</font></div>

<div style="color:black;font:10pt arial"><font size="2">Daniel Mckinnell<br>
</font></div>

</div>
</div>
</div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>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><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>