[MOAB-dev] In PyTAPS, how to create edge entities by pytaps?

[Paul Wilson]

Hi Huayi,

I think the issue is that MOAB (the underlying mesh database for this 
iMesh implementation) tries to conserve memory by not actually creating 
entities at every dimension.  When reading a VTK file, the vertices, 
faces and elements are explicitly generated, with the other entities 
defined implicitly by their connectivity.

It may be necessary to build these entities explicitly in order to get a 
count using some of the get_adjacency calls to determine the 
connectivity of the mesh.

Paul

On 06/17/2014 09:23 AM, Huayi Wei wrote:
> Hi, there,
>
> I loaded a tri mesh by PyTAPS. And I want to get the edge entities, 
> but edge entities doesn't exist in this mesh instance. So how to 
> create or get the edge entities ？
>
> ```
> >>> s = "4k-tri-plane.vtk"
> >>> mesh = iMesh.Mesh()
> >>> mesh.load(s)
> >>> mesh.getNumOfType(iBase.Type.all)
> 6085
> >>> mesh.getNumOfType(iBase.Type.vertex)
> 2085
> >>> mesh.getNumOfType(iBase.Type.face)
> 4000
> ```
>
> Best
>
> Huayi
>

-- 
-- ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ --
Paul Wilson ~ UW-Madison ~ 608-263-0807 ~ cal: http://bit.ly/pphw-cal
Professor, Engineering Physics. ~ http://cnerg.engr.wisc.edu
Faculty Director, Advanced Computing Infrastructure


-------------- next part --------------
A non-text attachment was scrubbed...
Name: smime.p7s
Type: application/pkcs7-signature
Size: 6244 bytes
Desc: S/MIME Cryptographic Signature
URL: <http://lists.mcs.anl.gov/pipermail/moab-dev/attachments/20140618/411fe8ed/attachment.bin>