[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