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
|
import gmsh
if length(ARGS) < 1
println("Usage: julia explore.jl file.msh")
exit(0)
end
gmsh.initialize()
try
gmsh.open(ARGS[1])
catch e
exit(0)
end
# get all elementary entities in the model
entities = gmsh.model.getEntities()
for e in entities
# get the mesh nodes for each elementary entity
nodeTags, nodeCoords, nodeParams = gmsh.model.mesh.getNodes(e[1], e[2])
# get the mesh elements for each elementary entity
elemTypes, elemTags, elemNodeTags = gmsh.model.mesh.getElements(e[1], e[2])
# report some statistics
numElem = sum([length(t) for t in elemTags])
println(length(nodeTags), " mesh nodes and ", numElem, " mesh elements ",
"on entity ", e, " of type ", gmsh.model.getType(e[1], e[2]))
partitions = gmsh.model.getPartitions(e[1], e[2])
if length(partitions) > 0
println(" - Partition tag(s): ", partitions, " - parent entity",
gmsh.model.getParent(e[1], e[2]))
end
for t in elemTypes
name, dim, order, numv, parv = gmsh.model.mesh.getElementProperties(t)
println(" - Element type: ", name, ", order ", order)
println(" with ", numv, " nodes in param coord: ", parv)
end
end
gmsh.finalize()
|
