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
|
# ------------------------------------------------------------------------------
#
# Gmsh Julia tutorial 8
#
# Post-processing, image export and animations
#
# ------------------------------------------------------------------------------
import gmsh
# In addition to creating geometries and meshes, the Julia API can also be used
# to manipulate post-processing datasets (called "views" in Gmsh).
gmsh.initialize()
# We first create a simple geometry
lc = 1e-2
gmsh.model.geo.addPoint(0, 0, 0, lc, 1)
gmsh.model.geo.addPoint(.1, 0, 0, lc, 2)
gmsh.model.geo.addPoint(.1, .3, 0, lc, 3)
gmsh.model.geo.addPoint(0, .3, 0, lc, 4)
gmsh.model.geo.addLine(1, 2, 1)
gmsh.model.geo.addLine(3, 2, 2)
gmsh.model.geo.addLine(3, 4, 3)
gmsh.model.geo.addLine(4, 1, 4)
gmsh.model.geo.addCurveLoop([4, 1, -2, 3], 1)
gmsh.model.geo.addPlaneSurface([1], 1)
gmsh.model.geo.synchronize()
# We merge some post-processing views to work on
gmsh.merge(abspath(joinpath(@__DIR__, "..", "view1.pos")))
gmsh.merge(abspath(joinpath(@__DIR__, "..", "view1.pos")))
gmsh.merge(abspath(joinpath(@__DIR__, "..", "view4.pos"))) # contains 2 views inside
# Gmsh can read post-processing views in various formats. Here the `view1.pos'
# and `view4.pos' files are in the Gmsh "parsed" format, which is interpreted by
# the GEO script parser. The parsed format should only be used for relatively
# small datasets of course: for larger datasets using e.g. MSH files is much
# more efficient. Post-processing views can also be created directly from the
# Julia API.
# We then set some general options:
gmsh.option.setNumber("General.Trackball", 0)
gmsh.option.setNumber("General.RotationX", 0)
gmsh.option.setNumber("General.RotationY", 0)
gmsh.option.setNumber("General.RotationZ", 0)
white = (255, 255, 255)
black = (0, 0, 0)
# Color options are special
# Setting a color option of "X.Y" actually sets the option "X.Color.Y"
# Sets "General.Color.Background", etc.
gmsh.option.setColor("General.Background", white[1], white[2], white[3])
gmsh.option.setColor("General.Foreground", black[1], black[2], black[3])
gmsh.option.setColor("General.Text", black[1], black[2], black[3])
gmsh.option.setNumber("General.Orthographic", 0)
gmsh.option.setNumber("General.Axes", 0)
gmsh.option.setNumber("General.SmallAxes", 0)
# Show the GUI:
if !("-nopopup" in ARGS)
gmsh.fltk.initialize()
end
# We also set some options for each post-processing view:
# If we were to follow the geo example blindly, we would read the number of
# views from the relevant option value, and use the gmsh.option.setNumber() and
# gmsh.option.setString() functions. A nicer way is to use gmsh.view.getTags()
# and to use the gmsh.view.option.setNumber() and gmsh.view.option.setString()
# functions.
v = gmsh.view.getTags()
if length(v) != 4
gmsh.logger.write("Wrong number of views!", "error")
gmsh.finalize()
exit()
end
# We set some options for each post-processing view:
gmsh.view.option.setNumber(v[1], "IntervalsType", 2)
gmsh.view.option.setNumber(v[1], "OffsetZ", 0.05)
gmsh.view.option.setNumber(v[1], "RaiseZ", 0)
gmsh.view.option.setNumber(v[1], "Light", 1)
gmsh.view.option.setNumber(v[1], "ShowScale", 0)
gmsh.view.option.setNumber(v[1], "SmoothNormals", 1)
gmsh.view.option.setNumber(v[2], "IntervalsType", 1)
# Note that we can't yet set the ColorTable in API
gmsh.view.option.setNumber(v[2], "NbIso", 10)
gmsh.view.option.setNumber(v[2], "ShowScale", 0)
gmsh.view.option.setString(v[3], "Name", "Test...")
gmsh.view.option.setNumber(v[3], "Axes", 1)
gmsh.view.option.setNumber(v[3], "IntervalsType", 2)
gmsh.view.option.setNumber(v[3], "Type", 2)
gmsh.view.option.setNumber(v[3], "AutoPosition", 0)
gmsh.view.option.setNumber(v[3], "PositionX", 85)
gmsh.view.option.setNumber(v[3], "PositionY", 50)
gmsh.view.option.setNumber(v[3], "Width", 200)
gmsh.view.option.setNumber(v[3], "Height", 130)
gmsh.view.option.setNumber(v[4], "Visible", 0)
# You can save an MPEG movie directly by selecting `File->Export' in the
# GUI. Several predefined animations are setup, for looping on all the time
# steps in views, or for looping between views.
# But the API can be used to build much more complex animations, by changing
# options at run-time and re-rendering the graphics. Each frame can then be
# saved to disk as an image, and multiple frames can be encoded to form a
# movie. Below is an example of such a custom animation.
t = 0 # Initial step
for num in 1:3
global t
# Set time step
for vv in v
gmsh.view.option.setNumber(vv, "TimeStep", t)
end
current_step = gmsh.view.option.getNumber(v[1], "TimeStep")
max_step = gmsh.view.option.getNumber(v[1], "NbTimeStep") - 1
if current_step < max_step
t = t + 1
else
t = 0
end
gmsh.view.option.setNumber(v[1], "RaiseZ",
gmsh.view.option.getNumber(v[1], "RaiseZ") +
0.01 / gmsh.view.option.getNumber(v[1], "Max") * t)
if num == 3
# Resize the graphics when num == 3, to create 640x480 frames
gmsh.option.setNumber("General.GraphicsWidth",
gmsh.option.getNumber("General.MenuWidth") + 640)
gmsh.option.setNumber("General.GraphicsHeight", 480)
end
frames = 50
for num2 in 0:frames-1
# Incrementally rotate the scene
gmsh.option.setNumber("General.RotationX",
gmsh.option.getNumber("General.RotationX") + 10)
gmsh.option.setNumber("General.RotationY",
gmsh.option.getNumber("General.RotationX") / 3)
gmsh.option.setNumber("General.RotationZ",
gmsh.option.getNumber("General.RotationZ") + 0.1)
# Draw the scene
gmsh.graphics.draw()
if num == 3
# Uncomment the following lines to save each frame to an image file
# gmsh.write("t8-{}.gif".format(num2))
# gmsh.write("t8-{}.ppm".format(num2))
# gmsh.write("t8-{}.jpg".format(num2))
end
end
if num == 3
# Here we could make a system call to generate a movie. For example,
# with ffmeg:
# run("ffmpeg -i t8-%d.jpg t8.mpg")
end
end
if !("-nopopup" in ARGS)
gmsh.fltk.run()
end
gmsh.finalize()
|
