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
|
// jvx.cmd
// Create jvx file for JavaView.
/*
Assumptions: 3D soapfilm model, linear model, not torus or symmetry group
(use the "detorus" command if necessary to convert torus or symmetry
to unwrapped surface, but remember that detorus alters the surface)
Does facets only, not edges.
Does facets satisfying "show" criterion.
Facet color is frontcolor on both sides.
*/
// Usage: jvx >>> "filename.jvx"
define vertex attribute jvx_number integer; // for implicit vertex order
define facet attribute fjvxnum integer; // for implicit facet order
define rgb_jvx integer[16][3]; // color definitions
rgb_jvx[1][1] := 0; rgb_jvx[1][2] := 0; rgb_jvx[1][3] := 0;
rgb_jvx[2][1] := 0; rgb_jvx[2][2] := 0; rgb_jvx[2][3] := 255;
rgb_jvx[3][1] := 0; rgb_jvx[3][2] := 255; rgb_jvx[3][3] := 0;
rgb_jvx[4][1] := 0; rgb_jvx[4][2] := 255; rgb_jvx[4][3] := 255;
rgb_jvx[5][1] := 255; rgb_jvx[5][2] := 0; rgb_jvx[5][3] := 0;
rgb_jvx[6][1] := 255; rgb_jvx[6][2] := 0; rgb_jvx[6][3] := 255;
rgb_jvx[7][1] := 255; rgb_jvx[7][2] := 127; rgb_jvx[7][3] := 0;
rgb_jvx[8][1] := 160; rgb_jvx[8][2] := 160; rgb_jvx[8][3] := 160;
rgb_jvx[9][1] := 80; rgb_jvx[9][2] := 80; rgb_jvx[9][3] := 80;
rgb_jvx[10][1] := 80; rgb_jvx[10][2] := 200; rgb_jvx[10][3] := 255;
rgb_jvx[11][1] := 127; rgb_jvx[11][2] := 255; rgb_jvx[11][3] := 127;
rgb_jvx[12][1] := 127; rgb_jvx[12][2] := 255; rgb_jvx[12][3] := 255;
rgb_jvx[13][1] := 255; rgb_jvx[13][2] := 127; rgb_jvx[13][3] := 127;
rgb_jvx[14][1] := 255; rgb_jvx[13][2] := 127; rgb_jvx[13][3] := 255;
rgb_jvx[15][1] := 255; rgb_jvx[15][2] := 255; rgb_jvx[15][3] := 0;
rgb_jvx[16][1] := 255; rgb_jvx[16][2] := 255; rgb_jvx[16][3] := 255;
jvx := {
local maxx,minx,maxy,miny,maxz,minz,jvxnum,facetnum;
local colornum,numcolors;
local facetcolorcount,edgecolorcount;
define facetcolorcount integer[16]; // to see which color facets are present
define edgecolorcount integer[16]; // to see which color edges are present
if torus then
{ errprintf "Cannot run 'jvx' command in torus mode. Do 'detorus' before 'jvx'.\n";
errprintf " WARNING: 'detorus' alters the surface, so save it first!\n";
abort;
};
if symmetry_group then
{ errprintf "Cannot run 'jvx' command in symmetry group mode. Do 'detorus' before 'jvx'.\n";
errprintf " WARNING: 'detorus' alters the surface, so save it first!\n";
abort;
};
if space_dimension != 3 then
{ errprintf "The 'jvx' command must be run in three-dimensional space.\n";
abort;
};
if surface_dimension == 1 then
{ errprintf "The 'jvx' command is not meant for the string model.\n";
abort;
};
if simplex_representation then
{ errprintf "The 'jvx' command is not meant for the simplex model.\n";
abort;
};
if lagrange_order >= 2 then
{ errprintf "The 'jvx' command is meant for the linear model, not quadratic or Lagrange.\n";
abort;
};
if rgb_colors then
{ errprintf "The 'jvx' command does not do RGB colors; do rgb_colors off.\n";
abort;
};
// First, bounding box.
maxx := max(vertex,x); minx := min(vertex,x);
maxy := max(vertex,y); miny := min(vertex,y);
maxz := max(vertex,z); minz := min(vertex,z);
printf "<?xml version=\"1.0\" encoding=\"ISO-8859-1\" standalone=\"yes\"?>\n";
printf "<jvx-model>\n";
printf "<meta generator=\"Surface Evolver 2.16; jvx.cmd\"/>\n";
// printf "<meta date=\"%s\"/>\n",date_and_time;
printf "<title> %s </title>\n",datafilename;
printf "<version>1.0</version>\n";
/*
printf "<authors>\n";
printf " <author>\n";
printf " <firstname>Kenneth</firstname>\n";
printf " <lastname>Brakke</lastname>\n";
printf " <affiliation>\n";
printf " <organization>Susquehanna University</organization>\n";
printf " <address>\n";
printf " <line>Mathematics Department</line>\n";
printf " <line>Susquehanna University</line>\n";
printf " <line>Selinsgrove, PA 17870</line>\n";
printf " <line>USA</line>\n";
printf " </address>\n";
printf " </affiliation>\n";
printf " <email>brakke@susqu.edu</email>\n";
printf " <url>http://www.susqu.edu/brakke</url>\n";
printf " </author>\n";
printf "</authors>\n";
printf "<description>\n";
printf " <abstract>One-line description of surface.</abstract>\n";
printf " <detail>\n";
printf " Paragraph of description here.\n";
printf " </detail>\n";
printf " <keywords>\n";
printf " <keyword>minimal surface</keyword>\n";
printf " </keywords>\n";
printf " <msc2000>\n";
printf " <primary>49Q20</primary>\n";
printf " <secondary> optional </secondary>\n";
printf " <secondary> optional </secondary>\n";
printf " </msc2000>\n";
printf " <software>Surface Evolver 2.16</software>\n";
printf "</description>\n";
*/
printf "<geometries>\n";
printf " <geometry name=\"surface1\">\n";
printf " <pointSet dim=\"3\" point=\"hide\">\n";
printf " <points num=\"%d\">\n",vertex_count;
jvxnum := 0;
foreach vertex vv do
{ printf " <p> %f %f %f </p>\n",vv.x,vv.y,vv.z;
vv.jvx_number := jvxnum;
jvxnum += 1;
};
printf " </points>\n";
printf " </pointSet>\n";
facetnum := 0;
foreach facet ff where ff.show and color != clear do
{ ff.fjvxnum := facetnum; facetnum += 1; };
colornum := 1;
while colornum <= 16 do { facetcolorcount[colornum] := 0; colornum += 1; };
foreach facet ff where ff.show and color != clear do
facetcolorcount[ff.color+1] += 1;
colornum := 1; numcolors := 0;
while colornum <= 16 do
{ if facetcolorcount[colornum] > 0 then numcolors += 1; colornum += 1; };
if numcolors > 1 then
printf " <faceSet face=\"show\" edge=\"hide\" color=\"show\">\n"
else printf " <faceSet face=\"show\" edge=\"hide\">\n";
printf " <faces num=\"%d\">\n",sum(facets,show);
foreach facet ff where ff.show and color != clear do
{ printf " <f> %d %d %d </f>\n",ff.vertex[1].jvx_number,
ff.vertex[2].jvx_number, ff.vertex[3].jvx_number;
};
if numcolors <= 1 then
{ colornum := 1;
while colornum <= 16 do
{ if facetcolorcount[colornum] > 0 then
printf " <color> %d %d %d </color>\n",rgb_jvx[colornum][1],
rgb_jvx[colornum][2],rgb_jvx[colornum][3];
colornum += 1;
};
};
printf " </faces>\n";
if numcolors > 1 then
{
printf " <colors>\n";
foreach facet ff where ff.show and color != clear do
{ printf " <c> %d %d %d </c>\n",rgb_jvx[ff.color+1][1],
rgb_jvx[ff.color+1][2],rgb_jvx[ff.color+1][3];
};
printf " </colors>\n";
};
printf " <neighbours num=\"%d\">\n",sum(facets,show);
foreach facet ff where ff.show and color != clear do
printf " <nb> %d %d %d </nb>\n",
ff.edge[2].valence == 2 ?
max(ff.edge[2].facet where id != ff.id,fjvxnum) : -1,
ff.edge[3].valence == 2 ?
max(ff.edge[3].facet where id != ff.id,fjvxnum) : -1,
ff.edge[1].valence == 2 ?
max(ff.edge[1].facet where id != ff.id,fjvxnum) : -1;
printf " </neighbours>\n";
// printf " <edges num=\"%d\">\n",sum(edges,show);
// foreach edge ee where ee.show do
// { printf " <e> %d %d </e>\n",ee.vertex[1].jvx_number,
// ee.vertex[2].jvx_number;
// };
// printf " <color> 0 0 0 </color>\n";
// printf " </edges>\n";
printf " <edges>\n";
printf " <thickness> 1.0 </thickness>\n";
printf " <color type=\"rgb\">0 0 0</color>\n";
printf " </edges>\n";
printf " </faceSet>\n";
printf " <bndbox visible=\"hide\">\n";
printf " <p> %f %f %f </p>\n",minx,miny,minz;
printf " <p> %f %f %f </p>\n",maxx,maxy,maxz;
printf " </bndbox>\n";
printf " <center visible=\"hide\">\n";
printf " <p> %f %f %f </p>\n",(maxx+minx)/2,(maxy+miny)/2,
(maxz+minz)/2;
printf " </center>\n";
printf " </geometry>\n";
colornum := 1;
while colornum <= 16 do { edgecolorcount[colornum] := 0; colornum += 1; };
foreach edge ee where ee.show do
edgecolorcount[ee.color+1] += 1;
colornum := 1; numcolors := 0;
while colornum <= 16 do
{ if edgecolorcount[colornum] > 0 then numcolors += 1; colornum += 1; };
printf " <geometry name=\"surface1 edges\">\n";
printf " <pointSet dim=\"3\" point=\"hide\">\n";
printf " <points num=\"%d\">\n",2*sum(edges,show);
foreach edge ee where ee.show do
{ printf " <p> %f %f %f </p>\n",ee.vertex[1].x,ee.vertex[1].y,ee.vertex[1].z;
printf " <p> %f %f %f </p>\n",ee.vertex[2].x,ee.vertex[2].y,ee.vertex[2].z;
};
printf " </points>\n";
printf " </pointSet>\n";
if numcolors > 1 then
printf " <lineSet color=\"show\">\n"
else printf " <lineSet>\n";
printf " <lines num=\"%d\">\n",sum(edges,show);
jvxnum := 0;
foreach edge ee where ee.show do
{ printf " <l> %d %d </l>\n",jvxnum,jvxnum+1;
jvxnum += 2;
};
printf " <thickness> 2 </thickness>\n";
if numcolors <= 1 then
{ colornum := 1;
while colornum <= 16 do
{ if edgecolorcount[colornum] > 0 then
printf " <color> %d %d %d </color>\n",rgb_jvx[colornum][1],
rgb_jvx[colornum][2],rgb_jvx[colornum][3];
colornum += 1;
};
};
printf " </lines>\n";
if ( numcolors > 1 ) then
{
printf " <colors>\n";
foreach edge ee where ee.show do
{ printf " <c> %d %d %d </c>\n", rgb_jvx[ee.color+1][1],
rgb_jvx[ee.color+1][2],rgb_jvx[ee.color+1][3];
};
printf " </colors>\n";
};
printf " </lineSet>\n";
printf " <bndbox visible=\"hide\">\n";
printf " <p> %f %f %f </p>\n",minx,miny,minz;
printf " <p> %f %f %f </p>\n",maxx,maxy,maxz;
printf " </bndbox>\n";
printf " <center visible=\"hide\">\n";
printf " <p> %f %f %f </p>\n",(maxx+minx)/2,(maxy+miny)/2,
(maxz+minz)/2;
printf " </center>\n";
printf " </geometry>\n";
printf "</geometries>\n";
printf "</jvx-model>\n";
}
// End jvx.cmd
/* Usage: Set "show" criterion for edges and facets, if you want, and do
jvx >>> "filename.jvx"
*/
|
