<div dir="ltr">On Thu, Apr 25, 2013 at 7:03 PM, Dharmendar Reddy <span dir="ltr"><<a href="mailto:dharmareddy84@gmail.com" target="_blank">dharmareddy84@gmail.com</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote">
<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><div>Hello,<br></div>         I need to mark the nodes on the boundary of a mesh to impose Dirichlet BC.  I am doing the following:<br>
<br></div><div>     call DMPlexGetHeightStratum(subdm,1,facetIdStart,facetIdend,ierr)<br>
     do facetId=facetIdStart,facetIdEnd<br>       ! check if a facet is internal, i.e., numSharedCell=2 because a facet is<br>       ! shared by two cells<br>       call DMPlexGetSupportSize(subdm,facetId,numSharedCell,ierr)<br>

       if(numSharedCell==1) then ! facet is a boundary face<br>         ! mark the nodes of this face as boundary<br>         <br>       end if<br>     end do</div></div></blockquote><div><br></div><div style>In SNES ex12, I do this. You could start with</div>
<div style><br></div><div style>      ierr = DMPlexMarkBoundaryFaces(dm, label);CHKERRQ(ierr);</div><div style><br></div><div style>but what you have is fine. Then I use</div><div style><br></div><div>      ierr  = DMPlexLabelComplete(dm, label);CHKERRQ(ierr);<br>
</div><div><br></div><div style>which does what you want. Underneath it does exactly what you proposed, namely</div><div style>mark the closure of each point in the label.</div><div style><br></div><div style>  Thanks,</div>
<div style><br></div><div style>     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><div><div>How do i get the 0-cells (ie nodes) forming given dim-cell ?<br>
<br></div><div>Am i using the right approach ?<br></div><div>Also, is there DMPlex functionality to mark boundary nodes ?<br><br></div><div>Thanks<span class=""><font color="#888888"><br>Reddy<br></font></span></div><span class=""><font color="#888888"><div>
-- <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: <a href="tel:%2B1-512-350-9082" value="+15123509082" target="_blank">+1-512-350-9082</a><br>United States of America.<br>Homepage: <a href="https://webspace.utexas.edu/~dpr342" target="_blank">https://webspace.utexas.edu/~dpr342</a><br>

</div></font></span></div></div></div>
</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>