[petsc-users] *****SPAM*****Re: how to check if cell is local owned in DMPlex

Matthew Knepley knepley at gmail.com
Mon Feb 19 17:30:00 CST 2018


On Mon, Feb 19, 2018 at 3:11 PM, Danyang Su <danyang.su at gmail.com> wrote:

> Hi Matt,
>
> Would you please let me know how to check if a cell is local owned? When
> overlap is 0 in DMPlexDistribute, all the cells are local owned. How about
> overlap > 0? It sounds like impossible to check by node because a cell can
> be local owned even if none of the nodes in this cell is local owned.
>

If a cell is in the PetscSF, then it is not locally owned. The local nodes
in the SF are sorted, so I use
PetscFindInt (
http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Sys/PetscFindInt.html
).

  Thanks,

    Matt


> Thanks,
>
> Danyang
>
-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener

https://www.cse.buffalo.edu/~knepley/ <http://www.caam.rice.edu/~mk51/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20180219/dbfb1028/attachment.html>


More information about the petsc-users mailing list