[petsc-users] DMPlex total point counts

[Matthew Knepley]

On Tue, Aug 2, 2022 at 4:00 PM Nicholas Arnold-Medabalimi <
narnoldm at umich.edu> wrote:

> Hi
>
> Currently, in my distributed DMPlex object, every so often, I want to look
> at total cell counts. To do this in my current implementation, I get the
> cell stratum from the DMplex and search through the PetscSF to remove any
> overlapped (not owned) cells and then reduce over all the ranks to get the
> total cell count. Is there a more graceful way of doing this already built
> into DM?
>

If all you want is the global cell count, this is probably the easiest way.
If you want a global numbering, you can use

  https://petsc.org/main/docs/manualpages/DMPLEX/DMPlexGetCellNumbering/

  Thanks,

     Matt


> Sincerely
> Nicholas
>
> --
> Nicholas Arnold-Medabalimi
>
> Ph.D. Candidate
> Computational Aeroscience Lab
> University of Michigan
>


-- 
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.cse.buffalo.edu/~knepley/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20220803/f3cb2012/attachment.html>