[petsc-users] Boundary Nodes of DMPlex object

[Dharmendar Reddy]

Should i be using :

DMPlexGetTransitiveClosure ?




On Thu, Apr 25, 2013 at 6:03 PM, Dharmendar Reddy
<dharmareddy84 at gmail.com>wrote:

> 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
>



-- 
-----------------------------------------------------
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/88628118/attachment.html>