[MOAB-dev] get uniquely owned vertices

[Grindeanu, Iulian R.]

Hi Nico,
"shared" entities means entities that are shared with other tasks


________________________________
From: moab-dev-bounces at mcs.anl.gov [moab-dev-bounces at mcs.anl.gov] on behalf of Nico Schlömer [nico.schloemer at gmail.com]
Sent: Monday, December 07, 2015 7:45 AM
To: Vijay S. Mahadevan
Cc: moab-dev at mcs.anl.gov
Subject: Re: [MOAB-dev] get uniquely owned vertices

Thanks Vijay for the explanation.
So both `get_entities_by_dimension` and `get_shared_entities`  return owned and ghosted entities?

Cheers,
Nico

On Mon, Dec 7, 2015 at 2:31 PM Vijay S. Mahadevan <vijay.m at gmail.com<mailto:vijay.m at gmail.com>> wrote:
Hi Nico,

You can use pstatus tag to filter out [1] your desired entities once
you have a superset of entities you care about. All calls on MOAB
Core/Interface return matching entities that are local (owned + shared
+ ghosted) only. So if you want only owned entities, or only ghosted
ones, you can do something like this:
```
  moab::Range allverts, ownverts, ghostverts;
  ierr = mb->get_entities_by_dimension(0, 0, allverts);
  ierr = pcomm->filter_pstatus(allverts, PSTATUS_NOT_OWNED,
PSTATUS_NOT, -1, &ownverts);
  ierr = pcomm->filter_pstatus(allverts, PSTATUS_GHOST, PSTATUS_OR,
-1, &ghostverts);
```

Vijay

[1] ftp://ftp.mcs.anl.gov/pub/fathom/moab-docs/classmoab_1_1ParallelComm.html#ad151b3b024cfd0ab7adc00f77e2551c2

On Mon, Dec 7, 2015 at 6:25 AM, Nico Schlömer <nico.schloemer at gmail.com<mailto:nico.schloemer at gmail.com>> wrote:
> Hi everyone,
>
> When reading a mesh in parallel, I would like to get the vertices such that
> every vertex belong to exactly one process. When reading with
> ```
>   moab::Range verts;
>   ierr = mb->get_entities_by_dimension(0, 0, verts);
> ```
> the vertices along the process interfaces will be shared, i.e., they appear
> on two processes.
>
> I used to believe the above does the right thing and I have to use get the
> shared entities separately
> ```
>   moab::Range shared;
>   ierr = mcomm_->get_shared_entities(-1, shared, 0);
>   all.merge(shared);
> ```
> but it seems I'm on the wrong track.
>
> Any hints?
>
> Cheers,
> Nico
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/moab-dev/attachments/20151207/04d23d74/attachment-0001.html>