-
Notifications
You must be signed in to change notification settings - Fork 11
/
readMesh.py
133 lines (105 loc) · 2.48 KB
/
readMesh.py
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
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
'''
Author: Andy Wang
Description: Parses mesh files into Python lists
Important Variables:
Nodes (x,y)
Elements ( (node_index_1, node_index_2, node_index_3), attribute)
Faces (node_index_1, node_index_2)
Neighbors (element_index_1, element_index_2, element_index_3)
'''
from globalVars import meshPath, meshName, dimension
import re
filePath = meshPath + meshName + ".1"
'''
Nodes
(x, y)
(x, y, z)
'''
firstLine = True
numLines = 0
Nodes = []
for line in open(filePath + ".node"):
if firstLine:
n = int( line.split(" ")[0] )
firstLine = False
continue
line = re.split(" +", line.strip())
t = map(float, line[1 : dimension + 1])
Nodes.append( tuple(t) )
numLines += 1
if numLines >= n:
break
'''
Elements
((node_index_0, node_index_1, node_index_2), attribute)
((node_index_0, node_index_1, node_index_2, node_index_3), attribute)
'''
firstLine = True
numLines = 0
Elements = []
for line in open(filePath + ".ele"):
if firstLine:
n = int( line.split(" ")[0] )
firstLine = False
continue
line = re.split(" +", line.strip())
t = [int(node_index) - 1 for node_index in line[1:] ]
t, attr = tuple( sorted( t[ : 1 + dimension] ) ), t[1 + dimension]
Elements.append( (t, attr) )
numLines += 1
if numLines >= n:
break
'''
Faces
(node_index_0, node_index_1)
(node_index_0, node_index_1, node_index_2)
'''
firstLine = True
numLines = 0
Faces = []
for line in open(filePath + ".edge"):
if firstLine:
n = int( line.split(" ")[0] )
firstLine = False
continue
line = re.split(" +", line.strip())
t = [int(node_index) - 1 for node_index in line[1 : 1 + dimension] ]
Faces.append( tuple( sorted(t) ) )
numLines += 1
if numLines >= n:
break
'''
Neighbors
(element_index_0, element_index_1, element_index_2)
(element_index_0, element_index_1, element_index_2, element_index_3)
'''
firstLine = True
numLines = 0
Neighbors = []
for line in open(filePath + ".neigh"):
if firstLine:
n = int( line.split(" ")[0] )
firstLine = False
continue
line = re.split(" +", line.strip())
t = [int(element_index) - 1 for element_index in line[1:]]
t = sorted(element_index for element_index in t if element_index != -2)
Neighbors.append( tuple(t) )
numLines += 1
if numLines >= n:
break
'''
Poly
eh, this is basically a summary of all the information above
'''
firstLine = True
numLines = 0
for line in open(filePath + ".poly"):
if firstLine:
n = int( line.split(" ")[0] )
firstLine = False
continue
line = re.split(" +", line.strip())
numLines += 1
if numLines >= n:
break