-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathcmg.m
494 lines (493 loc) · 9.21 KB
/
cmg.m
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
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
% Matlab mesh
% untitled, Created by Gmsh
% ASCII
clear msh;
msh.nbNod = 156;
msh.POS = [
0 0 0;
50 0 0;
50 50 0;
0 50 0;
25 12.5 0;
25 37.5 0;
4.99999999998079 0 0;
9.999999999957065 0 0;
14.99999999993112 0 0;
19.99999999992002 0 0;
24.99999999991184 0 0;
29.99999999990366 0 0;
34.99999999991796 0 0;
39.99999999994531 0 0;
44.99999999997265 0 0;
50 4.99999999998079 0;
50 9.999999999957065 0;
50 14.99999999993112 0;
50 19.99999999992002 0;
50 24.99999999991184 0;
50 29.99999999990366 0;
50 34.99999999991796 0;
50 39.99999999994531 0;
50 44.99999999997265 0;
45.00000000003553 50 0;
40.00000000007105 50 0;
35.00000000010658 50 0;
30.00000000012907 50 0;
25.00000000012907 50 0;
20.00000000012907 50 0;
15.00000000010658 50 0;
10.00000000007105 50 0;
5.000000000035527 50 0;
0 45.00000000003553 0;
0 40.00000000007105 0;
0 35.00000000010658 0;
0 30.00000000012907 0;
0 25.00000000012907 0;
0 20.00000000012907 0;
0 15.00000000010658 0;
0 10.00000000007105 0;
0 5.000000000035527 0;
25 17.49999999997613 0;
25 22.49999999993005 0;
25 27.49999999988398 0;
25 32.49999999990396 0;
37.2946736212956 20.42247780841625 0;
12.59516289991639 25.02413034397097 0;
37.37499999999519 34.99999999995197 0;
12.38620336671171 12.38796633295073 0;
12.38620336668229 37.61203366689077 0;
37.45117723458146 9.942704549443722 0;
41.02768217262006 26.20681708521252 0;
31.98430367773825 42.59569757788206 0;
33.52065645399939 27.62352045815971 0;
41.93218085105658 41.93218085103228 0;
7.784431106749517 18.78337395530823 0;
7.677806154686131 31.49164090456035 0;
17.34181395284249 18.41041553066569 0;
17.34181395283336 31.58958446917754 0;
18.73473865028075 7.756104539895794 0;
18.73473865020767 42.24389546009613 0;
30.56579898795677 7.506730112635911 0;
6.890774708150047 42.3209511384917 0;
6.890774708117888 7.679048861445464 0;
32.37983960108629 14.66752401931072 0;
42.56017346997093 14.90300466898672 0;
43.7281162129909 7.260466960711899 0;
43.68749999999393 33.64509813674822 0;
43.64672044882065 22.45846899449088 0;
20.74450517299256 24.99999999991632 0;
30.39236114572125 34.83925510895128 0;
26.23977917942075 43.7499999999872 0;
30.21576868539011 21.46817031488806 0;
5.94074065862581 24.95717844729782 0;
6.661132040777021 36.83253545007161 0;
6.67065212571827 13.19201944808219 0;
24.60469360552855 5.745315957082457 0;
14.00108156061731 5.613259729206797 0;
14.00108156062677 44.38674027081271 0;
19.38673637524887 13.01167723868277 0;
19.42056023390442 37.0329939215696 0;
37.09040471870525 44.78404829764397 0;
37.58093879081318 3.550317783794356 0;
37.39575001453549 14.81316603162266 0;
45.78968934491832 37.49954815165543 0;
37.33417607157894 29.40272628837579 0;
29.42323218685812 29.89492061357919 0;
12.81609519004634 19.77361167639097 0;
12.77569305646387 30.2932536966563 0;
37.11761605902963 39.9967627979343 0;
37.45463721662492 24.89829585758075 0;
45.58307997903969 11.46836679549677 0;
45.27735147161695 45.27735147161347 0;
21.72127307856697 45.59904971399403 0;
29.13277538750872 25.31951038222275 0;
45.48003282576585 27.92856268685329 0;
32.55035386263415 46.63481582932901 0;
45.71839669597447 17.83541982918394 0;
41.18039725099873 37.61471798746444 0;
21.04593460124479 20.26499871658146 0;
21.04593460124512 29.73500128326727 0;
32.70657520933042 3.686235678138872 0;
3.63590096844852 32.32721866714011 0;
3.65145038725655 17.72164723047118 0;
17.6698967295092 3.488050307477434 0;
17.89141865792741 46.44593708898057 0;
29.30981852966352 39.22751370637736 0;
29.4685464416777 12.00002542500462 0;
33.60908517163058 31.35208449380359 0;
22.53176838987832 41.20003257995573 0;
28.18395760545304 3.841301103390393 0;
33.46029828472124 10.7181641399176 0;
4.220169857041123 45.64565454911664 0;
4.220169857031061 4.354345450862754 0;
3.840566646779452 28.31191698083295 0;
41.18985257689322 46.01967974575138 0;
3.802168703814964 21.63284274872984 0;
11.02608093441418 16.26259029415173 0;
11.02608093446585 33.73740970565959 0;
16.48122803320454 22.61008069111969 0;
16.55125619836578 27.42085195458329 0;
15.36530882990531 14.86025703240062 0;
15.36502343339244 35.1445819465842 0;
29.06172335537437 17.6282330835369 0;
15.06097657797002 39.91612910245519 0;
15.05698390394292 10.0897649617093 0;
41.22292612944455 11.10991039007585 0;
22.54372099788541 8.867057846177204 0;
9.10277239024915 27.38349674759589 0;
3.946035036296615 40.48768336712918 0;
3.946035036244857 9.512316632723936 0;
9.484209397842834 3.9347072109716 0;
9.484209397764548 46.06529278904343 0;
46.22775143953425 3.772248560454305 0;
9.197363332991181 22.72307540511682 0;
21.42916773879826 3.661082110677274 0;
10.65298860361849 41.13812985640286 0;
10.65298860355446 8.861870143540266 0;
41.03431462288697 18.58962993703214 0;
33.23581988242957 38.3318458382194 0;
33.86465528053914 18.29921446143479 0;
47.12172545703572 7.499999999968928 0;
42.07567904502429 3.634569152254107 0;
33.85964507798092 23.37935449515397 0;
46.6845608595421 41.84544253692431 0;
47.23493813916864 32.4999999999108 0;
29.15488734401026 46.59610268143965 0;
47.23906880304668 22.49999999991593 0;
40.98087821399079 30.43664083942836 0;
35.28034272238208 7.104755645415757 0;
26.72778627308357 8.410071740715097 0;
21.53271931631359 33.64979883815461 0;
21.53271931635309 16.35020116173494 0;
40.09160561644964 22.51513793654651 0;
39.73321792671085 7.161153277808234 0;
];
msh.MAX = max(msh.POS);
msh.MIN = min(msh.POS);
msh.LINES =[
1 7 0
7 8 0
8 9 0
9 10 0
10 11 0
11 12 0
12 13 0
13 14 0
14 15 0
15 2 0
2 16 0
16 17 0
17 18 0
18 19 0
19 20 0
20 21 0
21 22 0
22 23 0
23 24 0
24 3 0
3 25 0
25 26 0
26 27 0
27 28 0
28 29 0
29 30 0
30 31 0
31 32 0
32 33 0
33 4 0
4 34 0
34 35 0
35 36 0
36 37 0
37 38 0
38 39 0
39 40 0
40 41 0
41 42 0
42 1 0
5 43 0
43 44 0
44 45 0
45 46 0
46 6 0
];
msh.TRIANGLES =[
44 45 71 0
6 46 72 0
93 18 99 0
53 70 97 0
50 77 139 0
76 51 138 0
80 31 134 0
9 79 133 0
31 32 134 0
8 9 133 0
77 40 132 0
36 76 131 0
40 41 132 0
35 36 131 0
77 50 119 0
51 76 120 0
67 93 99 0
97 69 150 0
5 109 125 0
69 97 147 0
3 25 94 0
24 3 94 0
97 70 149 0
74 44 125 0
64 76 138 0
77 65 139 0
43 5 125 0
44 71 101 0
71 45 102 0
72 46 88 0
57 77 119 0
76 58 120 0
54 73 108 0
73 6 108 0
73 29 95 0
53 97 150 0
6 72 108 0
73 54 148 0
49 72 110 0
49 69 100 0
69 49 150 0
70 53 155 0
72 49 141 0
26 27 83 0
6 73 111 0
44 74 96 0
29 73 148 0
44 43 125 0
38 75 116 0
75 38 118 0
17 18 93 0
22 86 147 0
82 62 126 0
61 81 127 0
124 82 126 0
81 123 127 0
101 71 121 0
71 102 122 0
56 83 91 0
114 34 131 0
42 115 132 0
76 64 131 0
65 77 132 0
76 36 104 0
77 57 105 0
58 76 104 0
40 77 105 0
115 7 133 0
33 114 134 0
79 9 106 0
61 79 106 0
31 80 107 0
80 62 107 0
62 82 111 0
82 6 111 0
78 11 112 0
86 69 147 0
11 78 137 0
79 61 127 0
62 80 126 0
59 89 123 0
90 60 124 0
89 119 123 0
120 90 124 0
22 23 86 0
109 66 125 0
59 101 121 0
102 60 122 0
34 35 131 0
41 42 132 0
85 66 113 0
7 8 133 0
32 33 134 0
83 54 91 0
67 85 128 0
4 34 114 0
33 4 114 0
1 7 115 0
42 1 115 0
26 83 117 0
52 85 113 0
13 14 84 0
129 61 137 0
118 57 136 0
78 129 137 0
54 98 148 0
46 45 88 0
83 27 98 0
54 83 98 0
85 67 140 0
5 81 129 0
81 61 129 0
81 5 154 0
6 82 153 0
85 52 128 0
29 30 95 0
87 55 92 0
53 87 92 0
75 118 136 0
88 45 96 0
25 26 117 0
56 91 100 0
82 60 153 0
59 81 154 0
58 116 130 0
94 25 117 0
60 82 124 0
81 59 123 0
45 44 96 0
83 56 117 0
112 63 152 0
84 14 144 0
51 124 126 0
123 50 127 0
55 88 96 0
2 16 135 0
15 2 135 0
47 85 140 0
12 13 103 0
87 49 110 0
78 112 152 0
116 75 130 0
12 103 112 0
142 47 145 0
37 38 116 0
53 92 155 0
38 39 118 0
84 144 156 0
11 12 112 0
49 87 150 0
10 11 137 0
69 86 100 0
13 84 103 0
28 29 148 0
46 6 153 0
5 43 154 0
18 19 99 0
64 114 131 0
115 65 132 0
20 21 97 0
91 49 100 0
68 93 128 0
49 91 141 0
65 115 133 0
114 64 134 0
27 28 98 0
67 99 140 0
99 70 140 0
96 74 145 0
55 96 145 0
85 47 142 0
66 85 142 0
86 56 100 0
45 46 102 0
43 44 101 0
36 37 104 0
39 40 105 0
9 10 106 0
30 31 107 0
47 92 145 0
74 142 145 0
103 63 112 0
56 86 146 0
97 21 147 0
20 97 149 0
95 62 111 0
73 95 111 0
95 30 107 0
62 95 107 0
55 87 110 0
88 55 110 0
72 88 110 0
23 24 146 0
58 104 116 0
48 121 122 0
121 71 122 0
105 57 118 0
94 56 146 0
68 128 156 0
61 106 137 0
134 64 138 0
80 134 138 0
133 79 139 0
65 133 139 0
89 57 119 0
58 90 120 0
70 99 149 0
99 19 149 0
48 89 121 0
90 48 122 0
89 59 121 0
60 90 122 0
48 90 130 0
48 130 136 0
57 89 136 0
24 94 146 0
101 59 154 0
60 102 153 0
144 68 156 0
130 75 136 0
90 58 130 0
56 94 117 0
14 15 144 0
16 17 143 0
119 50 123 0
51 120 124 0
89 48 136 0
66 109 113 0
21 22 147 0
19 20 149 0
91 54 141 0
74 125 142 0
125 66 142 0
52 151 156 0
151 84 156 0
93 67 128 0
86 23 146 0
87 53 150 0
92 55 145 0
17 93 143 0
93 68 143 0
109 63 113 0
104 37 116 0
92 47 155 0
39 105 118 0
54 108 141 0
108 72 141 0
63 103 151 0
113 63 151 0
63 109 152 0
79 127 139 0
127 50 139 0
126 80 138 0
51 126 138 0
106 10 137 0
98 28 148 0
103 84 151 0
128 52 156 0
43 101 154 0
102 46 153 0
5 129 152 0
129 78 152 0
109 5 152 0
52 113 151 0
135 68 144 0
15 135 144 0
135 16 143 0
68 135 143 0
47 140 155 0
140 70 155 0
];
msh.PNT =[
1 0
2 0
3 0
4 0
5 0
6 0
];
save('compare_cmg.mat')