<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<br></div><div><div><div><br></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<br>Reddy<br></div><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: +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></div></div>