1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
|
import gmsh
import sys
if len(sys.argv) < 2:
print("Usage: " + sys.argv[0] + " file.msh")
exit(0)
gmsh.initialize()
try:
gmsh.open(sys.argv[1])
except:
exit(0)
print("Model name: " + gmsh.model.getCurrent())
# get all elementary entities in the model
entities = gmsh.model.getEntities()
for e in entities:
print("Entity " + str(e) + " of type " + gmsh.model.getType(e[0], e[1]))
# get the mesh nodes for each elementary entity
nodeTags, nodeCoords, nodeParams = gmsh.model.mesh.getNodes(e[0], e[1])
# get the mesh elements for each elementary entity
elemTypes, elemTags, elemNodeTags = gmsh.model.mesh.getElements(e[0], e[1])
# count number of elements
numElem = sum(len(i) for i in elemTags)
print(" - mesh has " + str(len(nodeTags)) + " nodes and " + str(numElem) +
" elements")
boundary = gmsh.model.getBoundary([e])
print(" - boundary entities " + str(boundary))
partitions = gmsh.model.getPartitions(e[0], e[1])
if len(partitions):
print(" - Partition tag(s): " + str(partitions) + " - parent entity " +
str(gmsh.model.getParent(e[0], e[1])))
for t in elemTypes:
name, dim, order, numv, parv, _ = gmsh.model.mesh.getElementProperties(
t)
print(" - Element type: " + name + ", order " + str(order) + " (" +
str(numv) + " nodes in param coord: " + str(parv) + ")")
# all mesh node coordinates
nodeTags, nodeCoords, _ = gmsh.model.mesh.getNodes()
x = dict(zip(nodeTags, nodeCoords[0::3]))
y = dict(zip(nodeTags, nodeCoords[1::3]))
z = dict(zip(nodeTags, nodeCoords[2::3]))
gmsh.finalize()
|
