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
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
|
package require vtk
package require vtkinteraction
# Create 2D and 3D quadratic cells and extract their edges
# Quadratic triangle
vtkPoints triPoints
triPoints SetNumberOfPoints 6
triPoints InsertPoint 0 2.0 0.0 0.0
triPoints InsertPoint 1 3.0 0.0 0.0
triPoints InsertPoint 2 2.5 0.8 0.0
triPoints InsertPoint 3 2.5 0.0 0.0
triPoints InsertPoint 4 2.75 0.4 0.0
triPoints InsertPoint 5 2.25 0.4 0.0
vtkFloatArray triScalars
triScalars SetNumberOfTuples 6
triScalars InsertValue 0 0.0
triScalars InsertValue 1 0.0
triScalars InsertValue 2 0.0
triScalars InsertValue 3 1.0
triScalars InsertValue 4 1.0
triScalars InsertValue 5 0.0
vtkQuadraticTriangle aTri
[aTri GetPointIds] SetId 0 0
[aTri GetPointIds] SetId 1 1
[aTri GetPointIds] SetId 2 2
[aTri GetPointIds] SetId 3 3
[aTri GetPointIds] SetId 4 4
[aTri GetPointIds] SetId 5 5
vtkUnstructuredGrid aTriGrid
aTriGrid Allocate 1 1
aTriGrid InsertNextCell [aTri GetCellType] [aTri GetPointIds]
aTriGrid SetPoints triPoints
[aTriGrid GetPointData] SetScalars triScalars
# Quadratic quadrilateral
vtkPoints quadPoints
quadPoints SetNumberOfPoints 8
quadPoints InsertPoint 0 4.0 0.0 0.0
quadPoints InsertPoint 1 5.0 0.0 0.0
quadPoints InsertPoint 2 5.0 1.0 0.0
quadPoints InsertPoint 3 4.0 1.0 0.0
quadPoints InsertPoint 4 4.5 0.0 0.0
quadPoints InsertPoint 5 5.0 0.5 0.0
quadPoints InsertPoint 6 4.5 1.0 0.0
quadPoints InsertPoint 7 4.0 0.5 0.0
vtkFloatArray quadScalars
quadScalars SetNumberOfTuples 8
quadScalars InsertValue 0 0.0
quadScalars InsertValue 1 0.0
quadScalars InsertValue 2 1.0
quadScalars InsertValue 3 1.0
quadScalars InsertValue 4 1.0
quadScalars InsertValue 5 0.0
quadScalars InsertValue 6 0.0
quadScalars InsertValue 7 0.0
vtkQuadraticQuad aQuad
[aQuad GetPointIds] SetId 0 0
[aQuad GetPointIds] SetId 1 1
[aQuad GetPointIds] SetId 2 2
[aQuad GetPointIds] SetId 3 3
[aQuad GetPointIds] SetId 4 4
[aQuad GetPointIds] SetId 5 5
[aQuad GetPointIds] SetId 6 6
[aQuad GetPointIds] SetId 7 7
vtkUnstructuredGrid aQuadGrid
aQuadGrid Allocate 1 1
aQuadGrid InsertNextCell [aQuad GetCellType] [aQuad GetPointIds]
aQuadGrid SetPoints quadPoints
[aQuadGrid GetPointData] SetScalars quadScalars
# Quadratic tetrahedron
vtkPoints tetPoints
tetPoints SetNumberOfPoints 10
tetPoints InsertPoint 0 6.0 0.0 0.0
tetPoints InsertPoint 1 7.0 0.0 0.0
tetPoints InsertPoint 2 6.5 0.8 0.0
tetPoints InsertPoint 3 6.5 0.4 1.0
tetPoints InsertPoint 4 6.5 0.0 0.0
tetPoints InsertPoint 5 6.75 0.4 0.0
tetPoints InsertPoint 6 6.25 0.4 0.0
tetPoints InsertPoint 7 6.25 0.2 0.5
tetPoints InsertPoint 8 6.75 0.2 0.5
tetPoints InsertPoint 9 6.50 0.6 0.5
vtkFloatArray tetScalars
tetScalars SetNumberOfTuples 10
tetScalars InsertValue 0 1.0
tetScalars InsertValue 1 1.0
tetScalars InsertValue 2 1.0
tetScalars InsertValue 3 1.0
tetScalars InsertValue 4 0.0
tetScalars InsertValue 5 0.0
tetScalars InsertValue 6 0.0
tetScalars InsertValue 7 0.0
tetScalars InsertValue 8 0.0
tetScalars InsertValue 9 0.0
vtkQuadraticTetra aTet
[aTet GetPointIds] SetId 0 0
[aTet GetPointIds] SetId 1 1
[aTet GetPointIds] SetId 2 2
[aTet GetPointIds] SetId 3 3
[aTet GetPointIds] SetId 4 4
[aTet GetPointIds] SetId 5 5
[aTet GetPointIds] SetId 6 6
[aTet GetPointIds] SetId 7 7
[aTet GetPointIds] SetId 8 8
[aTet GetPointIds] SetId 9 9
vtkUnstructuredGrid aTetGrid
aTetGrid Allocate 1 1
aTetGrid InsertNextCell [aTet GetCellType] [aTet GetPointIds]
aTetGrid SetPoints tetPoints
[aTetGrid GetPointData] SetScalars tetScalars
# Quadratic hexahedron
vtkPoints hexPoints
hexPoints SetNumberOfPoints 20
hexPoints InsertPoint 0 8 0 0
hexPoints InsertPoint 1 9 0 0
hexPoints InsertPoint 2 9 1 0
hexPoints InsertPoint 3 8 1 0
hexPoints InsertPoint 4 8 0 1
hexPoints InsertPoint 5 9 0 1
hexPoints InsertPoint 6 9 1 1
hexPoints InsertPoint 7 8 1 1
hexPoints InsertPoint 8 8.5 0 0
hexPoints InsertPoint 9 9 0.5 0
hexPoints InsertPoint 10 8.5 1 0
hexPoints InsertPoint 11 8 0.5 0
hexPoints InsertPoint 12 8.5 0 1
hexPoints InsertPoint 13 9 0.5 1
hexPoints InsertPoint 14 8.5 1 1
hexPoints InsertPoint 15 8 0.5 1
hexPoints InsertPoint 16 8 0 0.5
hexPoints InsertPoint 17 9 0 0.5
hexPoints InsertPoint 18 9 1 0.5
hexPoints InsertPoint 19 8 1 0.5
vtkFloatArray hexScalars
hexScalars SetNumberOfTuples 20
hexScalars InsertValue 0 1.0
hexScalars InsertValue 1 1.0
hexScalars InsertValue 2 1.0
hexScalars InsertValue 3 1.0
hexScalars InsertValue 4 1.0
hexScalars InsertValue 5 1.0
hexScalars InsertValue 6 1.0
hexScalars InsertValue 7 1.0
hexScalars InsertValue 8 0.0
hexScalars InsertValue 9 0.0
hexScalars InsertValue 10 0.0
hexScalars InsertValue 11 0.0
hexScalars InsertValue 12 0.0
hexScalars InsertValue 13 0.0
hexScalars InsertValue 14 0.0
hexScalars InsertValue 15 0.0
hexScalars InsertValue 16 0.0
hexScalars InsertValue 17 0.0
hexScalars InsertValue 18 0.0
hexScalars InsertValue 19 0.0
vtkQuadraticHexahedron aHex
[aHex GetPointIds] SetId 0 0
[aHex GetPointIds] SetId 1 1
[aHex GetPointIds] SetId 2 2
[aHex GetPointIds] SetId 3 3
[aHex GetPointIds] SetId 4 4
[aHex GetPointIds] SetId 5 5
[aHex GetPointIds] SetId 6 6
[aHex GetPointIds] SetId 7 7
[aHex GetPointIds] SetId 8 8
[aHex GetPointIds] SetId 9 9
[aHex GetPointIds] SetId 10 10
[aHex GetPointIds] SetId 11 11
[aHex GetPointIds] SetId 12 12
[aHex GetPointIds] SetId 13 13
[aHex GetPointIds] SetId 14 14
[aHex GetPointIds] SetId 15 15
[aHex GetPointIds] SetId 16 16
[aHex GetPointIds] SetId 17 17
[aHex GetPointIds] SetId 18 18
[aHex GetPointIds] SetId 19 19
vtkUnstructuredGrid aHexGrid
aHexGrid Allocate 1 1
aHexGrid InsertNextCell [aHex GetCellType] [aHex GetPointIds]
aHexGrid SetPoints hexPoints
[aHexGrid GetPointData] SetScalars hexScalars
# Quadratic wedge
vtkPoints wedgePoints
wedgePoints SetNumberOfPoints 15
wedgePoints InsertPoint 0 10 0 0
wedgePoints InsertPoint 1 11 0 0
wedgePoints InsertPoint 2 10 1 0
wedgePoints InsertPoint 3 10 0 1
wedgePoints InsertPoint 4 11 0 1
wedgePoints InsertPoint 5 10 1 1
wedgePoints InsertPoint 6 10.5 0 0
wedgePoints InsertPoint 7 10.5 0.5 0
wedgePoints InsertPoint 8 10 0.5 0
wedgePoints InsertPoint 9 10.5 0 1
wedgePoints InsertPoint 10 10.5 0.5 1
wedgePoints InsertPoint 11 10 0.5 1
wedgePoints InsertPoint 12 10 0 0.5
wedgePoints InsertPoint 13 11 0 0.5
wedgePoints InsertPoint 14 10 1 0.5
vtkFloatArray wedgeScalars
wedgeScalars SetNumberOfTuples 15
wedgeScalars InsertValue 0 1.0
wedgeScalars InsertValue 1 1.0
wedgeScalars InsertValue 2 1.0
wedgeScalars InsertValue 3 1.0
wedgeScalars InsertValue 4 1.0
wedgeScalars InsertValue 5 1.0
wedgeScalars InsertValue 6 0.0
wedgeScalars InsertValue 7 0.0
wedgeScalars InsertValue 8 0.0
wedgeScalars InsertValue 9 0.0
wedgeScalars InsertValue 10 0.0
wedgeScalars InsertValue 11 0.0
wedgeScalars InsertValue 12 0.0
wedgeScalars InsertValue 13 0.0
wedgeScalars InsertValue 14 0.0
vtkQuadraticWedge aWedge
[aWedge GetPointIds] SetId 0 0
[aWedge GetPointIds] SetId 1 1
[aWedge GetPointIds] SetId 2 2
[aWedge GetPointIds] SetId 3 3
[aWedge GetPointIds] SetId 4 4
[aWedge GetPointIds] SetId 5 5
[aWedge GetPointIds] SetId 6 6
[aWedge GetPointIds] SetId 7 7
[aWedge GetPointIds] SetId 8 8
[aWedge GetPointIds] SetId 9 9
[aWedge GetPointIds] SetId 10 10
[aWedge GetPointIds] SetId 11 11
[aWedge GetPointIds] SetId 12 12
[aWedge GetPointIds] SetId 13 13
[aWedge GetPointIds] SetId 14 14
#aWedge DebugOn
#puts aWedge
vtkUnstructuredGrid aWedgeGrid
aWedgeGrid Allocate 1 1
aWedgeGrid InsertNextCell [aWedge GetCellType] [aWedge GetPointIds]
aWedgeGrid SetPoints wedgePoints
[aWedgeGrid GetPointData] SetScalars wedgeScalars
# Quadratic pyramid
vtkPoints pyraPoints
pyraPoints SetNumberOfPoints 13
pyraPoints InsertPoint 0 12 0 0
pyraPoints InsertPoint 1 13 0 0
pyraPoints InsertPoint 2 13 1 0
pyraPoints InsertPoint 3 12 1 0
pyraPoints InsertPoint 4 12 0 1
pyraPoints InsertPoint 5 12.5 0 0
pyraPoints InsertPoint 6 13 0.5 0
pyraPoints InsertPoint 7 12.5 1 0
pyraPoints InsertPoint 8 12 0.5 0
pyraPoints InsertPoint 9 12 0 0.5
pyraPoints InsertPoint 10 12.5 0 0.5
pyraPoints InsertPoint 11 12.5 0.5 0.5
pyraPoints InsertPoint 12 12 0.5 0.5
vtkFloatArray pyraScalars
pyraScalars SetNumberOfTuples 13
pyraScalars InsertValue 0 1.0
pyraScalars InsertValue 1 1.0
pyraScalars InsertValue 2 1.0
pyraScalars InsertValue 3 1.0
pyraScalars InsertValue 4 1.0
pyraScalars InsertValue 5 1.0
pyraScalars InsertValue 6 0.0
pyraScalars InsertValue 7 0.0
pyraScalars InsertValue 8 0.0
pyraScalars InsertValue 9 0.0
pyraScalars InsertValue 10 0.0
pyraScalars InsertValue 11 0.0
pyraScalars InsertValue 12 0.0
vtkQuadraticPyramid aPyramid
[aPyramid GetPointIds] SetId 0 0
[aPyramid GetPointIds] SetId 1 1
[aPyramid GetPointIds] SetId 2 2
[aPyramid GetPointIds] SetId 3 3
[aPyramid GetPointIds] SetId 4 4
[aPyramid GetPointIds] SetId 5 5
[aPyramid GetPointIds] SetId 6 6
[aPyramid GetPointIds] SetId 7 7
[aPyramid GetPointIds] SetId 8 8
[aPyramid GetPointIds] SetId 9 9
[aPyramid GetPointIds] SetId 10 10
[aPyramid GetPointIds] SetId 11 11
[aPyramid GetPointIds] SetId 12 12
#aPyramid DebugOn
#puts aPyramid
vtkUnstructuredGrid aPyraGrid
aPyraGrid Allocate 1 1
aPyraGrid InsertNextCell [aPyramid GetCellType] [aPyramid GetPointIds]
aPyraGrid SetPoints pyraPoints
[aPyraGrid GetPointData] SetScalars pyraScalars
#vtkUnstructuredGridWriter writer
# writer SetInput aWedgeGrid
# writer SetFileName foo.vtk
# writer Write
# Append the quadratic cells together
vtkAppendFilter appendF
appendF AddInput aTriGrid
appendF AddInput aQuadGrid
appendF AddInput aTetGrid
appendF AddInput aHexGrid
appendF AddInput aWedgeGrid
appendF AddInput aPyraGrid
# Extract the edges
vtkExtractEdges extract
extract SetInputConnection [appendF GetOutputPort]
vtkShrinkPolyData shrink
shrink SetInputConnection [extract GetOutputPort]
shrink SetShrinkFactor 0.90
vtkDataSetMapper aMapper
aMapper SetInputConnection [shrink GetOutputPort]
#aMapper ScalarVisibilityOff
vtkActor aActor
aActor SetMapper aMapper
[aActor GetProperty] SetRepresentationToWireframe
[aActor GetProperty] SetAmbient 1.0
# Create the rendering related stuff.
vtkRenderer ren1
vtkRenderWindow renWin
renWin AddRenderer ren1
vtkRenderWindowInteractor iren
iren SetRenderWindow renWin
ren1 SetBackground .1 .2 .3
renWin SetSize 400 150
# specify properties
ren1 AddActor aActor
renWin Render
[ren1 GetActiveCamera] Dolly 3.0
ren1 ResetCameraClippingRange
# render the image
#
iren AddObserver UserEvent {wm deiconify .vtkInteract}
iren Initialize
wm withdraw .
|
