[petsc-users] Boundary Nodes of DMPlex object

[Dharmendar Reddy]

Hello,
         I need to mark the nodes on the boundary of a mesh to impose
Dirichlet BC.  I am doing the following:

     call DMPlexGetHeightStratum(subdm,1,facetIdStart,facetIdend,ierr)
     do facetId=facetIdStart,facetIdEnd
       ! check if a facet is internal, i.e., numSharedCell=2 because a
facet is
       ! shared by two cells
       call DMPlexGetSupportSize(subdm,facetId,numSharedCell,ierr)
       if(numSharedCell==1) then ! facet is a boundary face
         ! mark the nodes of this face as boundary

       end if
     end do

How do i get the 0-cells (ie nodes) forming given dim-cell ?

Am i using the right approach ?
Also, is there DMPlex functionality to mark boundary nodes ?

Thanks
Reddy
-- 
-----------------------------------------------------
Dharmendar Reddy Palle
Graduate Student
Microelectronics Research center,
University of Texas at Austin,
10100 Burnet Road, Bldg. 160
MER 2.608F, TX 78758-4445
e-mail: dharmareddy84 at gmail.com
Phone: +1-512-350-9082
United States of America.
Homepage: https://webspace.utexas.edu/~dpr342
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130425/eb41966e/attachment-0001.html>