<div dir="ltr"><br><div class="gmail_extra"><div class="gmail_quote">On Thu, Apr 18, 2013 at 8:44 PM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="im">Dharmendar Reddy <<a href="mailto:dharmareddy84@gmail.com">dharmareddy84@gmail.com</a>> writes:<br>

<br>
> Sorry Again. I always use to use reply and the mail went to petsc-users.<br>
> Only since last few days it is not going to petsc-users. I will make sure i<br>
> reply-all from now on.<br>
<br>
</div>We just changed the list default to make it easier for people to filter<br>
email and to have discussions on petsc-users without necessarily being<br>
overwhelmed by all list traffic.  For example, "subscribe without<br>
delivery" works properly now.  Discussed at length starting here:<br>
<br>
<a href="http://lists.mcs.anl.gov/pipermail/petsc-users/2013-April/017067.html" target="_blank">http://lists.mcs.anl.gov/pipermail/petsc-users/2013-April/017067.html</a><br>
<div class="im"><br>
>>    This procedure may not work if U has multiple fields per mesh node<br>
> right ? The field names are already set into Default section of the dm.<br>
<br>
</div>Yes, the name is a combination of the object name (the automatic name<br>
was "Vec_0xyyyyyyyy"), the field name you give to DMPlex, and the<br>
location (point/cell).<br>
<div class="im"><br></div></blockquote><div><br></div><div> it tired with <br> PetscObjectSetName(u,'myfield', ierr)<br></div><div>Now the data in vtu file is calle<br><br></div><div>myfieldPoint0    So the pattern is ..<fieldName><Point| Cell><field Id> ? where the field Name is the one set using PetscObjectSetName ?<br>
<br></div><div>But why is it not using the name set using PetscSectionSetFieldName ?<br><br></div><div>If i have a few fields per point then i can manage with the field Ids  but usually the aux data i look at has about 10 to 12 fields. It would make things easier if i can refer to them using Names in the paraview.   As the viewer is accessing the DM object, why not get the user set field names in the default section ?<br>
</div><div><br></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 class="im">
>   I set the field layout and boundary points using DMPlexCreateSection and<br>
> DMSetDefualtSection<br>
>   Now, my test problem has 567 nodes. Node 1 and 567 have Dirichlet BC.<br>
> When i do DMGetGlobalVector it has size 565, the solution vector of the<br>
> snes  also has size 565.  The solution vector written vtu file has 567<br>
> values with zero value at the boundary nodes. Am i missing a step here?<br>
><br>
> Usgin Vecsetvaluessection i apply bc  inside the subroutines passed to<br>
> dmsnessetfunction/jacobain. The bc values are propagated to the solver.<br>
<br>
</div>I think Matt envisions running DMGlobalToLocalBegin/End yourself,<br>
filling in those values, and then VecView on the local Vec.  I'm not<br>
thrilled with using local vectors to represent global state, but it is<br>
the first thing the implementation does when you give it a global<br>
vector, so it's the right thing for you to do.<br>
</blockquote></div><br><br clear="all"><br>-- <br>-----------------------------------------------------<br>Dharmendar Reddy Palle<br>Graduate Student<br>Microelectronics Research center,<br>University of Texas at Austin,<br>
10100 Burnet Road, Bldg. 160<br>MER 2.608F, TX 78758-4445<br>e-mail: <a href="mailto:dharmareddy84@gmail.com" target="_blank">dharmareddy84@gmail.com</a><br>Phone: +1-512-350-9082<br>United States of America.<br>Homepage: <a href="https://webspace.utexas.edu/~dpr342" target="_blank">https://webspace.utexas.edu/~dpr342</a><br>

</div></div>