[MOAB-dev] Puzzled

Grindeanu, Iulian R. iulian at mcs.anl.gov
Mon Aug 11 18:12:34 CDT 2014


hello,
Yes, it does not add up. Your code seems fine.

 There are other issues I can't explain with this mesh file.
How was this mesh generated? 
It seems that elements have a "negative" orientation (negative volume).

moab uses the skinner to identify the boundary, a similar algorithm is used by paraview/visit to show the surface. if some elements are oriented differently, the skinner fails to identify the correct boundary, and it fails later in resolving the shared entities.

I can see only about 5x8x10~400  hexas, total, but mbsize reports 599; 

It seems that there are elements filling a cube 10x10x10
about half are hexas, but 5x10x10 is 500, (again, mbsize reports 599)

So is half of the volume filled with prisms, and the other half with hexas?

Iulian
________________________________________
From: moab-dev-bounces at mcs.anl.gov [moab-dev-bounces at mcs.anl.gov] on behalf of Gerd Heber [gheber at hdfgroup.org]
Sent: Monday, August 11, 2014 3:15 PM
To: moab-dev at mcs.anl.gov
Subject: [MOAB-dev] Puzzled

I'm probably doing something stupid, but here it is:

The attached mesh contains 599 hexahedra and 1,140 prisms.
I has 1,408 nodes. It was partitioned into two partitions with mbpart.
When I run the attached program with 'mpirun -np 2 feawdy feawd.h5m' the
output reads:

# hexes :       299     300
# prisms:       570     570
# vertexes:     781     781
# shared:       144     144
# owned:        781     637
# shared/owned: 144     0

The element counts add up nicely, but
I'm confused that the 'owned' counts don't add up to 1408.
What am I missing?

Thanks for your help.
G.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: feady0000.png
Type: image/png
Size: 118985 bytes
Desc: feady0000.png
URL: <http://lists.mcs.anl.gov/pipermail/moab-dev/attachments/20140811/5d77cd14/attachment-0001.png>


More information about the moab-dev mailing list