[MOAB-dev] vertex to cell adjacencies in parallel

[Hong-Jun Kim]

I think you may use "filter_pstatus" as follows.
Otherwise, I think it is not a bad idea to just check all element status with ! and & operators and gather the right elements.

------------
Range all_ents, owned_ents, owned_non_ghost_ents;
moab.get_entities_by_handle( 0, all_ents );

pcomm.filter_pstatus(all_ents, PSTATUS_NOT_OWNED, PSTATUS_NOT, -1, &owned_ents);
pcomm.filter_pstatus(owned_ents, PSTATUS_GHOST, PSTATUS_NOT, -1, &owned_non_ghost_ents);
-------------

Hong-Jun


----- Original Message -----
From: "Lorenzo Alessio Botti" <ihabiamx at yahoo.it>
To: "Hong-Jun Kim" <hongjun at mcs.anl.gov>
Cc: "Tim Tautges" <tautges at mcs.anl.gov>, "MOAB" <moab-dev at mcs.anl.gov>
Sent: Tuesday, October 11, 2011 3:02:38 PM
Subject: Re: [MOAB-dev] vertex to cell adjacencies in parallel

This option is easier and more readable. May I also ask which is the best way to get all the owned entities, excluding the ghost ones.

Thanks a lot.
Lorenzo
Inviato da iPhone

Il giorno 11/ott/2011, alle ore 21:15, Hong-Jun Kim <hongjun at mcs.anl.gov> ha scritto:

> Sorry, it should be as follows.
> 
> rval = pcomm.get_shared_entities(-1, shared_owned_vols, 3, false, true);
> 
> -----------------------------
> Hong-Jun Kim
> Post-doc researcher
> MCS, Argonne National Laboratory
> 9700 S. Cass Ave. B240/R2147
> Argonne, IL 60439
> 630-252-4791
> hongjun at mcs.anl.gov
> -----------------------------
> 
> ----- Original Message -----
> From: "Tim Tautges" <tautges at mcs.anl.gov>
> To: "Hong-Jun Kim" <hongjun at mcs.anl.gov>
> Cc: "Lorenzo Alessio Botti" <ihabiamx at yahoo.it>, moab-dev at mcs.anl.gov
> Sent: Tuesday, October 11, 2011 1:55:33 PM
> Subject: Re: [MOAB-dev] vertex to cell adjacencies in parallel
> 
> I don't think that will return any region entities, as those aren't marked as interface.  You'll only get vertices, 
> faces, and edges from that.
> 
> - tim
> 
> On 10/11/2011 01:31 PM, Hong-Jun Kim wrote:
>> As Tim mentioned, you can reduce your code with the following 1 function to get shared and owned elements.
>> 
>> Range shared_owned_vols;
>> rval = pcomm.get_shared_entities(-1, shared_owned_vols, 3, true, true);
>> 
>> Thanks.
>> 
>> -----------------------------
>> Hong-Jun Kim
>> Post-doc researcher
>> MCS, Argonne National Laboratory
>> 9700 S. Cass Ave. B240/R2147
>> Argonne, IL 60439