├── README.md
├── build
    ├── demo.cpp
    ├── demo.msh
    ├── elements.txt
    └── nodes.txt
└── demo
    ├── README
    ├── demo.cpp
    ├── oneDExample.cpp
    ├── threeDDemo.cpp
    └── twoDExample.cpp


/README.md:
--------------------------------------------------------------------------------
1 | # LearnGmsh
2 | Some demo programs of learning Gmsh.
3 | 
4 | 
5 | 
6 | demo      ：  storage all example programs
7 | build     ：  Compile and run directory
8 | 


--------------------------------------------------------------------------------
/build/demo.cpp:
--------------------------------------------------------------------------------
  1 | #include <iostream>
  2 | #include <gmsh.h>
  3 | #include <vector> 
  4 | #include <fstream>
  5 | #include <iomanip>
  6 | using namespace std;
  7 | 
  8 | 
  9 | void writeNodes(vector<double>& nodes, string s)
 10 | {
 11 | 	string filename = s + ".txt";
 12 | 	ofstream outfileb(filename.c_str(), ios::out);
 13 | 	for(int i = 0; i < nodes.size(); i++)
 14 | 	{
 15 | 		outfileb<<(nodes[i])<<endl;
 16 | 	}
 17 | 	outfileb.close();
 18 | }
 19 | 
 20 | 
 21 | void writeElements(vector<std::vector<std::size_t> >& elements, string s)
 22 | {
 23 | 	string filename = s + ".txt";
 24 | 	ofstream outfileb(filename.c_str(), ios::out);
 25 | 	for(int i = 0; i < elements[1].size(); i++)
 26 | 	{
 27 | 		outfileb<<(elements[1][i])<<endl;
 28 | 	}
 29 | 	outfileb.close();
 30 | }
 31 | 
 32 | 
 33 | int main(int argc, char **argv)
 34 | {
 35 | 	// Before using any functions in the C++ API, Gmsh must be initialized.
 36 | 	gmsh::initialize();
 37 | 	// By default Gmsh will not print out any messages: in order to output
 38 | 	// messages on the terminal, just set the standard Gmsh option
 39 | 	// "General.Terminal" (same format and meaning as in .geo files) using
 40 | 	// gmsh::option::setNumber():
 41 | 	//gmsh::option::setNumber("General.Terminal", 1);
 42 | 
 43 | 	// This adds a new model, named "t1". If gmsh::model::add() is not called, a
 44 | 	// new default (unnamed) model will be created on the fly, if necessary.
 45 | 	gmsh::model::add("Test1");
 46 | 	// The C++ API provides direct access to the internal CAD kernels. The
 47 |     // built-in CAD kernel was used in t1.geo: the corresponding API functions
 48 |     // live in the "gmsh::model::geo" namespace. To create geometrical points with
 49 |     // the built-in CAD kernel, one thus uses gmsh::model::geo::addPoint():
 50 |     //
 51 |     // - the first 3 arguments are the point coordinates (x, y, z)
 52 |     //
 53 |     // - the next (optional) argument is the target mesh size close to the point
 54 |     //
 55 |     // - the last (optional) argument is the point tag
 56 | 
 57 | 
 58 | 
 59 |   	double lcc = 0.05;//2e-1;
 60 |   	gmsh::model::geo::addPoint(1.0/3, 1.0/3,  0, lcc, 100);
 61 |   	gmsh::model::geo::addPoint(2.0/3, 1.0/3,  0, lcc, 101);
 62 |   	gmsh::model::geo::addPoint(2.0/3, 2.0/3,  0, lcc, 102);
 63 |   	gmsh::model::geo::addPoint(1.0/3, 2.0/3,  0, lcc, 103);
 64 |   	gmsh::model::geo::addLine(100, 101, 100);
 65 |   	gmsh::model::geo::addLine(101, 102, 101);
 66 |   	gmsh::model::geo::addLine(102, 103, 102);
 67 |   	gmsh::model::geo::addLine(103, 100, 103);
 68 |   	gmsh::model::geo::addCurveLoop({100, 101, 102, 103}, 100);
 69 |   	//gmsh::model::geo::addPlaneSurface({100}, 100);
 70 | 
 71 | 
 72 | 
 73 |   	double lc = 0.2;//2e-1;
 74 |   	gmsh::model::geo::addPoint(-1, 0,  0, lc, 1);
 75 |   	gmsh::model::geo::addPoint(1,  0,  0, lc, 2);
 76 |   	gmsh::model::geo::addPoint(1,  1,  0, lc, 3);
 77 |   	gmsh::model::geo::addPoint(-1, 1,  0, lc, 4);
 78 |   	gmsh::model::geo::addPoint(0,  0,  0, lc, 555);
 79 | 
 80 | 
 81 |   	// The API to create lines with the built-in kernel follows the same
 82 |   	// conventions: the first 2 arguments are point tags, the last (optional one)
 83 |   	// is the line tag.
 84 |   	gmsh::model::geo::addLine(1, 555, 1);
 85 |   	gmsh::model::geo::addLine(555, 2, 555);
 86 |   	gmsh::model::geo::addLine(2, 3, 2);
 87 |   	gmsh::model::geo::addLine(3, 4, 3);
 88 |   	gmsh::model::geo::addLine(4, 1, 4);
 89 | 
 90 |   	// The philosophy to construct curve loops and surfaces is similar: the first
 91 |   	// argument is now a vector of integers.
 92 |   	gmsh::model::geo::addCurveLoop({1, 555, 2, 3, 4}, 1);
 93 |   	gmsh::model::geo::addPlaneSurface({1, -100}, 1);
 94 | 
 95 | 
 96 |   	gmsh::model::geo::addPoint(1,  2,  0, lc, 5);
 97 |   	gmsh::model::geo::addPoint(-1, 2,  0, lc, 6);
 98 |   	gmsh::model::geo::addLine(3, 5, 5);
 99 |   	gmsh::model::geo::addLine(5, 6, 6);
100 |   	gmsh::model::geo::addLine(6, 4, 7);
101 |   	gmsh::model::geo::addCurveLoop({-3, 5, 6, 7}, 2);
102 | 	gmsh::model::geo::addPlaneSurface({2}, 2);
103 | 
104 | 
105 | 
106 | 
107 |   	// Physical groups are defined by providing the dimension of the group (0 for
108 |   	// physical points, 1 for physical curves, 2 for physical surfaces and 3 for
109 |   	// phsyical volumes) followed by a vector of entity tags. The last (optional)
110 |   	// argument is the tag of the new group to create.
111 | 
112 | 	gmsh::model::setPhysicalName(2, 6, "My surface");
113 | 
114 | 	// Before it can be meshed, the internal CAD representation must be
115 | 	// synchronized with the Gmsh model, which will create the relevant Gmsh data
116 | 	// structures. This is achieved by the gmsh::model::geo::synchronize() API
117 | 	// call for the built-in CAD kernel. Synchronizations can be called at any
118 | 	// time, but they involve a non trivial amount of processing; so while you
119 | 	// could synchronize the internal CAD data after every CAD command, it is
120 | 	// usually better to minimize the number of synchronization points.
121 | 	gmsh::model::geo::synchronize();
122 | 
123 | 	// We can then generate a 2D mesh...
124 | 	gmsh::model::mesh::generate(2);
125 | 	//gmsh::model::mesh::refine();
126 | 
127 | 	// get mesh information
128 | 	// Nodes
129 | 	vector<double> nodes,y;
130 | 	vector<std::size_t> nodeTags;
131 | 	gmsh::model::mesh::getNodes(nodeTags, nodes, y, -2, -1, false, true);
132 | 	cout<<"The number of nodes: "<<nodes.size()/3<<endl;
133 | #if 0
134 | 	cout<<endl;
135 | 	for(int i = 0; i < nodes.size(); i++)
136 | 	{
137 | 		cout<<nodes[i]<<endl;
138 | 	}
139 | #endif
140 | 	writeNodes(nodes, "nodes");
141 | 
142 | 	//Elements
143 | 	vector<int> elementTypes;
144 | 	vector<std::vector<std::size_t> > elementTags, nodeTags2;
145 | 	gmsh::model::mesh::getElements(elementTypes, elementTags, nodeTags2, -2, -1);
146 | 	cout<<"The number of elements: "<<nodeTags2[1].size()/3<<endl;
147 | #if 0
148 | 	cout<<endl;
149 | 	for(int j = 0; j < nodeTags2[1].size(); j++)
150 | 	{
151 | 		cout<<nodeTags2[1][j]<<endl;
152 | 		
153 | 	}
154 | #endif
155 | 	writeElements(nodeTags2, "elements");
156 | 
157 | 	// ... and save it to disk
158 | 	gmsh::write("demo.msh");
159 | 
160 | 	// Remember that by default, if physical groups are defined, Gmsh will export
161 | 	// in the output mesh file only those elements that belong to at least one
162 | 	// physical group. To force Gmsh to save all elements, you can use
163 | 	//
164 | 	// gmsh::option::setNumber("Mesh.SaveAll", 1);
165 | 
166 | 	// This should be called at the end:
167 | 	gmsh::finalize();
168 | 	return 0;
169 | }
170 | 


--------------------------------------------------------------------------------
/build/demo.msh:
--------------------------------------------------------------------------------
   1 | $MeshFormat
   2 | 4.1 0 8
   3 | $EndMeshFormat
   4 | $PhysicalNames
   5 | 1
   6 | 2 6 "My surface"
   7 | $EndPhysicalNames
   8 | $Entities
   9 | 11 12 2 0
  10 | 1 -1 0 0 0 
  11 | 2 1 0 0 0 
  12 | 3 1 1 0 0 
  13 | 4 -1 1 0 0 
  14 | 5 1 2 0 0 
  15 | 6 -1 2 0 0 
  16 | 100 0.3333333333333333 0.3333333333333333 0 0 
  17 | 101 0.6666666666666666 0.3333333333333333 0 0 
  18 | 102 0.6666666666666666 0.6666666666666666 0 0 
  19 | 103 0.3333333333333333 0.6666666666666666 0 0 
  20 | 555 0 0 0 0 
  21 | 1 -1 0 0 0 0 0 0 2 1 -555 
  22 | 2 1 0 0 1 1 0 0 2 2 -3 
  23 | 3 -1 1 0 1 1 0 0 2 3 -4 
  24 | 4 -1 0 0 -1 1 0 0 2 4 -1 
  25 | 5 1 1 0 1 2 0 0 2 3 -5 
  26 | 6 -1 2 0 1 2 0 0 2 5 -6 
  27 | 7 -1 1 0 -1 2 0 0 2 6 -4 
  28 | 100 0.3333333333333333 0.3333333333333333 0 0.6666666666666666 0.3333333333333333 0 0 2 100 -101 
  29 | 101 0.6666666666666666 0.3333333333333333 0 0.6666666666666666 0.6666666666666666 0 0 2 101 -102 
  30 | 102 0.3333333333333333 0.6666666666666666 0 0.6666666666666666 0.6666666666666666 0 0 2 102 -103 
  31 | 103 0.3333333333333333 0.3333333333333333 0 0.3333333333333333 0.6666666666666666 0 0 2 103 -100 
  32 | 555 0 0 0 1 0 0 0 2 555 -2 
  33 | 1 -1 0 0 1 1 0 0 9 1 555 2 3 4 100 101 102 103 
  34 | 2 -1 1 0 1 2 0 0 4 -3 5 6 7 
  35 | $EndEntities
  36 | $Nodes
  37 | 25 239 1 239
  38 | 0 1 0 1
  39 | 1
  40 | -1 0 0
  41 | 0 2 0 1
  42 | 2
  43 | 1 0 0
  44 | 0 3 0 1
  45 | 3
  46 | 1 1 0
  47 | 0 4 0 1
  48 | 4
  49 | -1 1 0
  50 | 0 5 0 1
  51 | 5
  52 | 1 2 0
  53 | 0 6 0 1
  54 | 6
  55 | -1 2 0
  56 | 0 100 0 1
  57 | 7
  58 | 0.3333333333333333 0.3333333333333333 0
  59 | 0 101 0 1
  60 | 8
  61 | 0.6666666666666666 0.3333333333333333 0
  62 | 0 102 0 1
  63 | 9
  64 | 0.6666666666666666 0.6666666666666666 0
  65 | 0 103 0 1
  66 | 10
  67 | 0.3333333333333333 0.6666666666666666 0
  68 | 0 555 0 1
  69 | 11
  70 | 0 0 0
  71 | 1 1 0 4
  72 | 12
  73 | 13
  74 | 14
  75 | 15
  76 | -0.7999999999999998 0 0
  77 | -0.6000000000013869 0 0
  78 | -0.4000000000016644 0 0
  79 | -0.2000000000008322 0 0
  80 | 1 2 0 4
  81 | 16
  82 | 17
  83 | 18
  84 | 19
  85 | 1 0.1999999999995579 0
  86 | 1 0.3999999999989749 0
  87 | 1 0.5999999999989468 0
  88 | 1 0.7999999999994734 0
  89 | 1 3 0 9
  90 | 20
  91 | 21
  92 | 22
  93 | 23
  94 | 24
  95 | 25
  96 | 26
  97 | 27
  98 | 28
  99 | 0.8000000000005548 1 0
 100 | 0.6000000000011096 1 0
 101 | 0.4000000000016644 1 0
 102 | 0.2000000000022192 1 0
 103 | 2.752242878045763e-12 1 0
 104 | -0.1999999999977808 1 0
 105 | -0.3999999999983355 1 0
 106 | -0.5999999999988903 1 0
 107 | -0.7999999999994454 1 0
 108 | 1 4 0 4
 109 | 29
 110 | 30
 111 | 31
 112 | 32
 113 | -1 0.7999999999999998 0
 114 | -1 0.6000000000013869 0
 115 | -1 0.4000000000016644 0
 116 | -1 0.2000000000008322 0
 117 | 1 5 0 4
 118 | 33
 119 | 34
 120 | 35
 121 | 36
 122 | 1 1.2 0
 123 | 1 1.4 0
 124 | 1 1.6 0
 125 | 1 1.8 0
 126 | 1 6 0 9
 127 | 37
 128 | 38
 129 | 39
 130 | 40
 131 | 41
 132 | 42
 133 | 43
 134 | 44
 135 | 45
 136 | 0.8000000000005548 2 0
 137 | 0.6000000000011096 2 0
 138 | 0.4000000000016644 2 0
 139 | 0.2000000000022192 2 0
 140 | 2.752242878045763e-12 2 0
 141 | -0.1999999999977808 2 0
 142 | -0.3999999999983355 2 0
 143 | -0.5999999999988903 2 0
 144 | -0.7999999999994454 2 0
 145 | 1 7 0 4
 146 | 46
 147 | 47
 148 | 48
 149 | 49
 150 | -1 1.8 0
 151 | -1 1.6 0
 152 | -1 1.4 0
 153 | -1 1.2 0
 154 | 1 100 0 6
 155 | 50
 156 | 51
 157 | 52
 158 | 53
 159 | 54
 160 | 55
 161 | 0.3809523809522818 0.3333333333333333 0
 162 | 0.428571428571204 0.3333333333333333 0
 163 | 0.4761904761901122 0.3333333333333333 0
 164 | 0.5238095238091398 0.3333333333333333 0
 165 | 0.5714285714283083 0.3333333333333333 0
 166 | 0.6190476190474922 0.3333333333333333 0
 167 | 1 101 0 6
 168 | 56
 169 | 57
 170 | 58
 171 | 59
 172 | 60
 173 | 61
 174 | 0.6666666666666666 0.3809523809522818 0
 175 | 0.6666666666666666 0.428571428571204 0
 176 | 0.6666666666666666 0.4761904761901122 0
 177 | 0.6666666666666666 0.5238095238091398 0
 178 | 0.6666666666666666 0.5714285714283083 0
 179 | 0.6666666666666666 0.6190476190474922 0
 180 | 1 102 0 6
 181 | 62
 182 | 63
 183 | 64
 184 | 65
 185 | 66
 186 | 67
 187 | 0.619047619047585 0.6666666666666666 0
 188 | 0.5714285714288503 0.6666666666666666 0
 189 | 0.5238095238101156 0.6666666666666666 0
 190 | 0.4761904761911269 0.6666666666666666 0
 191 | 0.4285714285718624 0.6666666666666666 0
 192 | 0.3809523809525978 0.6666666666666666 0
 193 | 1 103 0 6
 194 | 68
 195 | 69
 196 | 70
 197 | 71
 198 | 72
 199 | 73
 200 | 0.3333333333333333 0.619047619047585 0
 201 | 0.3333333333333333 0.5714285714288503 0
 202 | 0.3333333333333333 0.5238095238101156 0
 203 | 0.3333333333333333 0.4761904761911269 0
 204 | 0.3333333333333333 0.4285714285718624 0
 205 | 0.3333333333333333 0.3809523809525978 0
 206 | 1 555 0 4
 207 | 74
 208 | 75
 209 | 76
 210 | 77
 211 | 0.1999999999995579 0 0
 212 | 0.3999999999989749 0 0
 213 | 0.5999999999989468 0 0
 214 | 0.7999999999994734 0 0
 215 | 2 1 0 112
 216 | 78
 217 | 79
 218 | 80
 219 | 81
 220 | 82
 221 | 83
 222 | 84
 223 | 85
 224 | 86
 225 | 87
 226 | 88
 227 | 89
 228 | 90
 229 | 91
 230 | 92
 231 | 93
 232 | 94
 233 | 95
 234 | 96
 235 | 97
 236 | 98
 237 | 99
 238 | 100
 239 | 101
 240 | 102
 241 | 103
 242 | 104
 243 | 105
 244 | 106
 245 | 107
 246 | 108
 247 | 109
 248 | 110
 249 | 111
 250 | 112
 251 | 113
 252 | 114
 253 | 115
 254 | 116
 255 | 117
 256 | 118
 257 | 119
 258 | 120
 259 | 121
 260 | 122
 261 | 123
 262 | 124
 263 | 125
 264 | 126
 265 | 127
 266 | 128
 267 | 129
 268 | 130
 269 | 131
 270 | 132
 271 | 133
 272 | 134
 273 | 135
 274 | 136
 275 | 137
 276 | 138
 277 | 139
 278 | 140
 279 | 141
 280 | 142
 281 | 143
 282 | 144
 283 | 145
 284 | 146
 285 | 147
 286 | 148
 287 | 149
 288 | 150
 289 | 151
 290 | 152
 291 | 153
 292 | 154
 293 | 155
 294 | 156
 295 | 157
 296 | 158
 297 | 159
 298 | 160
 299 | 161
 300 | 162
 301 | 163
 302 | 164
 303 | 165
 304 | 166
 305 | 167
 306 | 168
 307 | 169
 308 | 170
 309 | 171
 310 | 172
 311 | 173
 312 | 174
 313 | 175
 314 | 176
 315 | 177
 316 | 178
 317 | 179
 318 | 180
 319 | 181
 320 | 182
 321 | 183
 322 | 184
 323 | 185
 324 | 186
 325 | 187
 326 | 188
 327 | 189
 328 | 0.2958339746560507 0.5003511884053465 0
 329 | -0.1081633915501419 0.1592072083925976 0
 330 | -0.1108437584743492 0.8294705067528795 0
 331 | -0.4873317348560665 0.1537537459264083 0
 332 | 0.2978735985523417 0.4018276742117401 0
 333 | 0.2963081020478399 0.598169519149394 0
 334 | -0.4917481652168882 0.8047642467942955 0
 335 | 0.7126542208250464 0.4960317460313927 0
 336 | 0.7103681632491897 0.5997809164845965 0
 337 | 0.7192051839042464 0.4057267509221013 0
 338 | 0.4999999999996596 0.2880384773622469 0
 339 | 0.5912698412696455 0.2873457791747926 0
 340 | 0.4087301587299836 0.2873277251368429 0
 341 | 0.5000000000005556 0.7124740222122581 0
 342 | 0.4087301587305098 0.7126890462638911 0
 343 | 0.5912698412700127 0.7126542208248416 0
 344 | -0.8211992221817153 0.5003818391998363 0
 345 | 0.7000000000008323 0.8522839056540795 0
 346 | 0.5315766701891277 0.8516958021763784 0
 347 | 0.8615507300186402 0.3061351350669481 0
 348 | 0.8648950053732909 0.4585171844733877 0
 349 | 0.2999999999992664 0.1477160943457985 0
 350 | 0.1381569764327036 0.1450194384994088 0
 351 | 0.2216464980975118 0.2596951904518771 0
 352 | 0.09702896965992366 0.2755410174344728 0
 353 | 0.1733014375709231 0.3656206376991714 0
 354 | 0.06035485297121138 0.3859558951408046 0
 355 | 0.1410337865163363 0.476675579146753 0
 356 | 0.01552899940124144 0.5065986829960127 0
 357 | 0.09768120477965672 0.5905936309841155 0
 358 | -0.07437179610427359 0.4080446702290406 0
 359 | -0.1243167235152927 0.5408789556239633 0
 360 | -0.2380261500256744 0.4254897838248473 0
 361 | -0.280182845601839 0.5744214910552736 0
 362 | -0.4148115954935033 0.4503506990048431 0
 363 | 0.4544063625667717 0.1445384716924746 0
 364 | 0.2775164009539346 0.8607772267332072 0
 365 | 0.1345361949963018 0.8565084240633661 0
 366 | 0.2199202455177193 0.7468417045904041 0
 367 | 0.8522839056549192 0.6999999999992101 0
 368 | 0.7016082155818444 0.148335107831305 0
 369 | 0.1957545890086183 0.5566252281667632 0
 370 | 0.1014425066936531 0.7200374667401355 0
 371 | -0.4589750768651193 0.6096896064020738 0
 372 | -0.600239178814179 0.4765343557597732 0
 373 | -0.6507992133330598 0.6549118413406578 0
 374 | -0.7279842749120053 0.3366882702900913 0
 375 | -0.285786942619385 0.835764962470251 0
 376 | -0.2909222210256258 0.1544900762522539 0
 377 | 0.3219913095396041 0.7671370328189047 0
 378 | 0.3979485348897127 0.8265880515275813 0
 379 | 0.8572858816741598 0.1459424957319209 0
 380 | 0.8550355706230863 0.8550355706214505 0
 381 | 0.3099321787795289 0.2479378763192695 0
 382 | -0.8326336369212892 0.7000000000006933 0
 383 | 0.1805850811343698 0.6502247863008124 0
 384 | -0.02446531775749456 0.6692847761762384 0
 385 | 0.2135038104782744 0.4396958670213075 0
 386 | -0.1872487677799896 0.6890195973592328 0
 387 | 0.2462874310717214 0.6706268136537157 0
 388 | -0.699999999999168 0.8381009998872744 0
 389 | 0.2430948782156112 0.3356612520601484 0
 390 | 0.2893445950755124 0.3579370359274245 0
 391 | 0.2772437813567042 0.4483547612050965 0
 392 | 0.2733128845517848 0.5517193125398439 0
 393 | 0.4562104770239036 0.7632627154987139 0
 394 | 0.4539405081035917 0.7043518234616393 0
 395 | 0.3533977439845711 0.7136734434503321 0
 396 | 0.6387155063280899 0.2953106068832634 0
 397 | 0.6917291831636807 0.2698934122142613 0
 398 | 0.6190476190476051 0.2345281928367354 0
 399 | 0.7298209010440151 0.3431829404499504 0
 400 | 0.5442757894054403 0.2376356490352939 0
 401 | 0.5838675973201216 0.1530074842791618 0
 402 | 0.6387155063281595 0.7046893931165181 0
 403 | 0.6915377289279495 0.7301802798670602 0
 404 | 0.6167864657225938 0.7688019738124764 0
 405 | 0.7327446458372554 0.6641613466125691 0
 406 | 0.7771389876551131 0.6015326348956914 0
 407 | 0.5405770836376127 0.7527991176888631 0
 408 | 0.5461567451824386 0.2959373144478 0
 409 | 0.7164944187042894 0.5527321527430318 0
 410 | 0.7814073945461081 0.5182136448334282 0
 411 | 0.7646652763020477 0.451880005468722 0
 412 | 0.3506365984088644 0.2903527227550221 0
 413 | 0.3793477084775123 0.2260928787017893 0
 414 | 0.4553399115801914 0.2438773508365683 0
 415 | 0.4537663950142302 0.297182044000465 0
 416 | 0.7059716028729348 0.4516800814367065 0
 417 | 0.5454170040294294 0.7022521388118592 0
 418 | -0.6667791058363611 0.1716608381814571 0
 419 | 0.8551450586458863 0.5756526928401329 0
 420 | 0.8042661295233491 0.3828207993913124 0
 421 | 0.2474526988904948 0.6117038004446091 0
 422 | -0.1861363021329183 0.2900747194009808 0
 423 | -0.3601238844296241 0.2996091928791109 0
 424 | 0.2299424534984453 0.4995017075941209 0
 425 | 0.2918569665125996 0.6454162136170729 0
 426 | 0.7743536086549332 0.2321607961698178 0
 427 | 0.7732930922841463 0.769385141997859 0
 428 | -0.5433494848118832 0.3212236554622404 0
 429 | -0.03234350262720026 0.2842368334998377 0
 430 | 0.2906396624401005 0.3049384958596506 0
 431 | -0.8426453339958503 0.1779636978615194 0
 432 | 0.01722665139299713 0.1865974904011301 0
 433 | -0.854502525391576 0.8545025253917271 0
 434 | 0.2944460864164176 0.7039584323284174 0
 435 | 0.01722506778019735 0.8134221155189779 0
 436 | -0.3432787920055922 0.6980208585058402 0
 437 | 0.2460468985237876 0.3911089324938525 0
 438 | -0.8734351539347203 0.3334142806782022 0
 439 | 0.3857404265819021 0.758363865807174 0
 440 | 2 2 0 50
 441 | 190
 442 | 191
 443 | 192
 444 | 193
 445 | 194
 446 | 195
 447 | 196
 448 | 197
 449 | 198
 450 | 199
 451 | 200
 452 | 201
 453 | 202
 454 | 203
 455 | 204
 456 | 205
 457 | 206
 458 | 207
 459 | 208
 460 | 209
 461 | 210
 462 | 211
 463 | 212
 464 | 213
 465 | 214
 466 | 215
 467 | 216
 468 | 217
 469 | 218
 470 | 219
 471 | 220
 472 | 221
 473 | 222
 474 | 223
 475 | 224
 476 | 225
 477 | 226
 478 | 227
 479 | 228
 480 | 229
 481 | 230
 482 | 231
 483 | 232
 484 | 233
 485 | 234
 486 | 235
 487 | 236
 488 | 237
 489 | 238
 490 | 239
 491 | 0.3000000000020203 1.173205080756751 0
 492 | -0.09999999999761203 1.173205080757016 -0
 493 | -0.837244843819962 1.509667905494991 -0
 494 | 0.4896860819409486 1.836450076080011 0
 495 | 0.1000000000024857 1.826794919243583 0
 496 | -0.2999999999980582 1.826794919242641 -0
 497 | 0.8362933879519492 1.511334914258793 0
 498 | -0.4999999999986053 1.173205080757443 -0
 499 | 0.6991525261446764 1.169354409010789 0
 500 | -0.68611490773137 1.808097945376386 -0
 501 | -0.2999999999981168 1.173205080757334 -0
 502 | -0.1999999999974288 1.346410161514085 -0
 503 | 2.56615099540151e-12 1.346410161513837 0
 504 | -0.09999999999711794 1.519615242270834 -0
 505 | 0.1000000000028235 1.519615242270536 0
 506 | -0.2999999999971742 1.519615242271185 -0
 507 | 0.2000000000023782 1.346410161513627 0
 508 | 0.3000000000025264 1.519615242270207 0
 509 | 0.4000000000020144 1.346410161513375 0
 510 | 0.5152045956149124 1.518583659286369 0
 511 | -0.0999999999975143 1.826794919243266 -0
 512 | -0.3999999999979094 1.346410161514659 -0
 513 | -0.5152153411280574 1.518553193817405 -0
 514 | -0.6062371528807424 1.32988406126764 -0
 515 | -0.6996652585422465 1.169737621024827 -0
 516 | 0.1000000000023873 1.173205080756872 0
 517 | -0.4890113227668438 1.832387037178297 -0
 518 | 0.8306178565233185 1.699886423395244 0
 519 | 0.6060852408238202 1.32926046601995 0
 520 | 0.3000000000019418 1.826794919243593 0
 521 | 0.5008729611625509 1.169705019550144 0
 522 | 0.8121168520916187 1.316365973602896 0
 523 | 0.6843755678592193 1.810881348927808 0
 524 | -0.8409589541417488 1.695206648784086 -0
 525 | -0.8123553311347743 1.316260425560569 -0
 526 | 0.6840187486595389 1.614188552469801 0
 527 | -0.6862062166523103 1.612645533464884 -0
 528 | 0.3921357983894884 1.676815347028361 0
 529 | 0.2000000000029278 1.67320508075671 0
 530 | 2.871856580966955e-12 1.673205080756968 0
 531 | -0.1999999999972526 1.67320508075692 -0
 532 | -0.392029948794212 1.676129958493429 -0
 533 | -0.8535898384856323 1.146410161514219 -0
 534 | 0.8535898384862793 1.85358983848613 0
 535 | -0.8535898384857241 1.853589838485873 -0
 536 | 0.853589838486725 1.146410161513423 0
 537 | -0.5494530597339548 1.687093305175372 -0
 538 | 0.5530841584928216 1.69138379675847 0
 539 | 0.6907437650283679 1.457946713127562 0
 540 | -0.6897502242148376 1.45812322962231 -0
 541 | $EndNodes
 542 | $Elements
 543 | 25 499 1 499
 544 | 0 1 15 1
 545 | 1 1 
 546 | 0 2 15 1
 547 | 2 2 
 548 | 0 3 15 1
 549 | 3 3 
 550 | 0 4 15 1
 551 | 4 4 
 552 | 0 5 15 1
 553 | 5 5 
 554 | 0 6 15 1
 555 | 6 6 
 556 | 0 100 15 1
 557 | 7 7 
 558 | 0 101 15 1
 559 | 8 8 
 560 | 0 102 15 1
 561 | 9 9 
 562 | 0 103 15 1
 563 | 10 10 
 564 | 0 555 15 1
 565 | 11 11 
 566 | 1 1 1 5
 567 | 12 1 12 
 568 | 13 12 13 
 569 | 14 13 14 
 570 | 15 14 15 
 571 | 16 15 11 
 572 | 1 2 1 5
 573 | 17 2 16 
 574 | 18 16 17 
 575 | 19 17 18 
 576 | 20 18 19 
 577 | 21 19 3 
 578 | 1 3 1 10
 579 | 22 3 20 
 580 | 23 20 21 
 581 | 24 21 22 
 582 | 25 22 23 
 583 | 26 23 24 
 584 | 27 24 25 
 585 | 28 25 26 
 586 | 29 26 27 
 587 | 30 27 28 
 588 | 31 28 4 
 589 | 1 4 1 5
 590 | 32 4 29 
 591 | 33 29 30 
 592 | 34 30 31 
 593 | 35 31 32 
 594 | 36 32 1 
 595 | 1 5 1 5
 596 | 37 3 33 
 597 | 38 33 34 
 598 | 39 34 35 
 599 | 40 35 36 
 600 | 41 36 5 
 601 | 1 6 1 10
 602 | 42 5 37 
 603 | 43 37 38 
 604 | 44 38 39 
 605 | 45 39 40 
 606 | 46 40 41 
 607 | 47 41 42 
 608 | 48 42 43 
 609 | 49 43 44 
 610 | 50 44 45 
 611 | 51 45 6 
 612 | 1 7 1 5
 613 | 52 6 46 
 614 | 53 46 47 
 615 | 54 47 48 
 616 | 55 48 49 
 617 | 56 49 4 
 618 | 1 100 1 7
 619 | 57 7 50 
 620 | 58 50 51 
 621 | 59 51 52 
 622 | 60 52 53 
 623 | 61 53 54 
 624 | 62 54 55 
 625 | 63 55 8 
 626 | 1 101 1 7
 627 | 64 8 56 
 628 | 65 56 57 
 629 | 66 57 58 
 630 | 67 58 59 
 631 | 68 59 60 
 632 | 69 60 61 
 633 | 70 61 9 
 634 | 1 102 1 7
 635 | 71 9 62 
 636 | 72 62 63 
 637 | 73 63 64 
 638 | 74 64 65 
 639 | 75 65 66 
 640 | 76 66 67 
 641 | 77 67 10 
 642 | 1 103 1 7
 643 | 78 10 68 
 644 | 79 68 69 
 645 | 80 69 70 
 646 | 81 70 71 
 647 | 82 71 72 
 648 | 83 72 73 
 649 | 84 73 7 
 650 | 1 555 1 5
 651 | 85 11 74 
 652 | 86 74 75 
 653 | 87 75 76 
 654 | 88 76 77 
 655 | 89 77 2 
 656 | 2 1 2 282
 657 | 90 97 98 170 
 658 | 91 96 22 128 
 659 | 92 98 18 169 
 660 | 93 153 95 154 
 661 | 94 118 147 148 
 662 | 95 1 12 181 
 663 | 96 114 115 116 
 664 | 97 142 83 171 
 665 | 98 100 99 101 
 666 | 99 94 122 123 
 667 | 100 32 1 181 
 668 | 101 113 76 151 
 669 | 102 149 147 176 
 670 | 103 153 155 177 
 671 | 104 122 94 124 
 672 | 105 172 110 173 
 673 | 106 95 96 154 
 674 | 107 99 113 163 
 675 | 108 119 142 171 
 676 | 109 102 101 103 
 677 | 110 96 128 143 
 678 | 111 142 119 174 
 679 | 112 82 141 187 
 680 | 113 114 116 127 
 681 | 114 22 114 128 
 682 | 115 113 150 164 
 683 | 116 126 172 173 
 684 | 117 101 99 131 
 685 | 118 94 123 132 
 686 | 119 117 156 169 
 687 | 120 80 25 125 
 688 | 121 117 155 156 
 689 | 122 103 101 139 
 690 | 123 74 99 100 
 691 | 124 114 23 115 
 692 | 125 56 8 149 
 693 | 126 9 61 155 
 694 | 127 107 105 119 
 695 | 128 96 143 157 
 696 | 129 97 149 176 
 697 | 130 155 117 177 
 698 | 131 121 84 123 
 699 | 132 17 18 98 
 700 | 133 75 76 113 
 701 | 134 21 22 96 
 702 | 135 141 78 174 
 703 | 136 152 153 154 
 704 | 137 147 146 148 
 705 | 138 80 125 136 
 706 | 139 11 100 182 
 707 | 140 147 118 176 
 708 | 141 95 153 177 
 709 | 142 80 134 185 
 710 | 143 95 21 96 
 711 | 144 97 17 98 
 712 | 145 99 75 113 
 713 | 146 134 80 136 
 714 | 147 115 24 185 
 715 | 148 20 21 95 
 716 | 149 16 17 97 
 717 | 150 74 75 99 
 718 | 151 76 77 118 
 719 | 152 18 19 117 
 720 | 153 22 23 114 
 721 | 154 24 25 80 
 722 | 155 15 11 79 
 723 | 156 100 101 102 
 724 | 157 23 24 115 
 725 | 158 11 74 100 
 726 | 159 18 117 169 
 727 | 160 168 81 178 
 728 | 161 81 14 126 
 729 | 162 124 168 178 
 730 | 163 15 79 126 
 731 | 164 156 86 159 
 732 | 165 106 107 134 
 733 | 166 78 142 174 
 734 | 167 105 103 135 
 735 | 168 123 84 138 
 736 | 169 76 118 151 
 737 | 170 26 84 125 
 738 | 171 25 26 125 
 739 | 172 135 141 174 
 740 | 173 87 149 170 
 741 | 174 110 112 173 
 742 | 175 128 127 189 
 743 | 176 109 106 134 
 744 | 177 116 115 120 
 745 | 178 149 97 170 
 746 | 179 104 103 105 
 747 | 180 150 113 151 
 748 | 181 161 87 170 
 749 | 182 119 105 174 
 750 | 183 79 11 182 
 751 | 184 14 15 126 
 752 | 185 116 120 133 
 753 | 186 141 135 187 
 754 | 187 24 80 185 
 755 | 188 150 88 164 
 756 | 189 143 91 157 
 757 | 190 13 14 81 
 758 | 191 26 27 84 
 759 | 192 30 31 94 
 760 | 193 84 121 186 
 761 | 194 106 105 107 
 762 | 195 30 94 132 
 763 | 196 3 20 130 
 764 | 197 2 16 129 
 765 | 198 19 3 130 
 766 | 199 77 2 129 
 767 | 200 84 27 138 
 768 | 201 29 30 132 
 769 | 202 107 119 133 
 770 | 203 126 79 172 
 771 | 204 120 107 133 
 772 | 205 109 134 136 
 773 | 206 122 121 123 
 774 | 207 131 99 163 
 775 | 208 87 56 149 
 776 | 209 61 86 155 
 777 | 210 125 84 186 
 778 | 211 102 103 104 
 779 | 212 112 121 122 
 780 | 213 20 95 130 
 781 | 214 16 97 129 
 782 | 215 118 77 129 
 783 | 216 117 19 130 
 784 | 217 27 28 138 
 785 | 218 13 81 168 
 786 | 219 134 120 185 
 787 | 220 132 123 138 
 788 | 221 107 120 134 
 789 | 222 136 125 186 
 790 | 223 12 168 181 
 791 | 224 105 135 174 
 792 | 225 81 126 173 
 793 | 226 112 111 121 
 794 | 227 31 32 188 
 795 | 228 4 29 183 
 796 | 229 28 4 183 
 797 | 230 154 96 157 
 798 | 231 133 119 171 
 799 | 232 110 111 112 
 800 | 233 163 113 164 
 801 | 234 118 148 151 
 802 | 235 122 124 178 
 803 | 236 110 109 111 
 804 | 237 156 159 160 
 805 | 238 116 133 137 
 806 | 239 112 122 178 
 807 | 240 101 131 180 
 808 | 241 162 7 180 
 809 | 242 111 109 136 
 810 | 243 98 160 161 
 811 | 244 108 109 110 
 812 | 245 171 83 175 
 813 | 246 12 13 168 
 814 | 247 72 71 141 
 815 | 248 73 72 82 
 816 | 249 71 70 78 
 817 | 250 69 68 83 
 818 | 251 70 69 142 
 819 | 252 104 105 106 
 820 | 253 127 116 184 
 821 | 254 129 97 176 
 822 | 255 117 130 177 
 823 | 256 130 95 177 
 824 | 257 118 129 176 
 825 | 258 108 106 109 
 826 | 259 32 181 188 
 827 | 260 114 127 128 
 828 | 261 172 79 179 
 829 | 262 124 94 188 
 830 | 263 175 10 184 
 831 | 264 7 73 140 
 832 | 265 68 10 175 
 833 | 266 173 112 178 
 834 | 267 104 106 108 
 835 | 268 111 136 186 
 836 | 269 10 145 184 
 837 | 270 116 137 184 
 838 | 271 66 65 144 
 839 | 272 65 64 91 
 840 | 273 62 9 152 
 841 | 274 8 55 146 
 842 | 275 67 66 92 
 843 | 276 60 59 159 
 844 | 277 54 53 158 
 845 | 278 53 52 88 
 846 | 279 59 58 85 
 847 | 280 61 60 86 
 848 | 281 55 54 89 
 849 | 282 51 50 90 
 850 | 283 57 56 87 
 851 | 284 63 62 93 
 852 | 285 64 63 167 
 853 | 286 58 57 166 
 854 | 287 52 51 165 
 855 | 288 10 67 145 
 856 | 289 50 7 162 
 857 | 290 131 162 180 
 858 | 291 179 79 182 
 859 | 292 139 101 180 
 860 | 293 108 110 172 
 861 | 294 94 31 188 
 862 | 295 29 132 183 
 863 | 296 82 72 141 
 864 | 297 71 78 141 
 865 | 298 78 70 142 
 866 | 299 69 83 142 
 867 | 300 135 103 187 
 868 | 301 108 172 179 
 869 | 302 132 138 183 
 870 | 303 100 102 182 
 871 | 304 91 143 144 
 872 | 305 143 92 144 
 873 | 306 138 28 183 
 874 | 307 137 171 175 
 875 | 308 121 111 186 
 876 | 309 160 98 169 
 877 | 310 155 86 156 
 878 | 311 8 147 149 
 879 | 312 153 9 155 
 880 | 313 73 82 140 
 881 | 314 83 68 175 
 882 | 315 140 82 187 
 883 | 316 103 139 187 
 884 | 317 137 133 171 
 885 | 318 88 150 158 
 886 | 319 150 89 158 
 887 | 320 87 161 166 
 888 | 321 161 85 166 
 889 | 322 93 157 167 
 890 | 323 157 91 167 
 891 | 324 90 164 165 
 892 | 325 164 88 165 
 893 | 326 102 104 179 
 894 | 327 102 179 182 
 895 | 328 156 160 169 
 896 | 329 92 143 189 
 897 | 330 104 108 179 
 898 | 331 168 124 181 
 899 | 332 148 89 150 
 900 | 333 160 85 161 
 901 | 334 93 154 157 
 902 | 335 90 163 164 
 903 | 336 120 115 185 
 904 | 337 7 140 180 
 905 | 338 81 173 178 
 906 | 339 145 127 184 
 907 | 340 148 150 151 
 908 | 341 162 131 163 
 909 | 342 98 161 170 
 910 | 343 8 146 147 
 911 | 344 152 9 153 
 912 | 345 65 91 144 
 913 | 346 92 66 144 
 914 | 347 91 64 167 
 915 | 348 93 62 152 
 916 | 349 59 85 159 
 917 | 350 53 88 158 
 918 | 351 86 60 159 
 919 | 352 89 54 158 
 920 | 353 85 58 166 
 921 | 354 88 52 165 
 922 | 355 55 89 146 
 923 | 356 51 90 165 
 924 | 357 57 87 166 
 925 | 358 63 93 167 
 926 | 359 127 145 189 
 927 | 360 67 92 145 
 928 | 361 90 50 162 
 929 | 362 137 175 184 
 930 | 363 159 85 160 
 931 | 364 93 152 154 
 932 | 365 146 89 148 
 933 | 366 90 162 163 
 934 | 367 139 140 187 
 935 | 368 181 124 188 
 936 | 369 145 92 189 
 937 | 370 143 128 189 
 938 | 371 140 139 180 
 939 | 2 2 2 128
 940 | 372 213 224 214 
 941 | 373 198 221 218 
 942 | 374 222 225 217 
 943 | 375 223 226 199 
 944 | 376 221 238 218 
 945 | 377 213 239 224 
 946 | 378 222 237 225 
 947 | 379 226 236 199 
 948 | 380 27 214 28 
 949 | 381 190 208 206 
 950 | 382 213 214 197 
 951 | 383 197 214 27 
 952 | 384 190 220 208 
 953 | 385 206 215 190 
 954 | 386 208 220 218 
 955 | 387 211 213 197 
 956 | 388 24 191 25 
 957 | 389 201 211 200 
 958 | 390 191 201 200 
 959 | 391 205 211 201 
 960 | 392 190 215 23 
 961 | 393 23 215 24 
 962 | 394 191 200 25 
 963 | 395 42 210 41 
 964 | 396 218 220 198 
 965 | 397 202 215 206 
 966 | 398 191 202 201 
 967 | 399 202 206 204 
 968 | 400 41 210 194 
 969 | 401 24 215 191 
 970 | 402 26 200 197 
 971 | 403 206 207 204 
 972 | 404 35 217 196 
 973 | 405 202 204 203 
 974 | 406 25 200 26 
 975 | 407 199 216 44 
 976 | 408 200 211 197 
 977 | 409 206 208 207 
 978 | 410 202 203 201 
 979 | 411 205 212 211 
 980 | 412 36 217 35 
 981 | 413 35 196 34 
 982 | 414 26 197 27 
 983 | 415 43 216 195 
 984 | 416 43 195 42 
 985 | 417 195 210 42 
 986 | 418 45 199 44 
 987 | 419 44 216 43 
 988 | 420 191 215 202 
 989 | 421 212 213 211 
 990 | 422 39 193 38 
 991 | 423 22 190 23 
 992 | 424 208 209 207 
 993 | 425 40 219 39 
 994 | 426 39 219 193 
 995 | 427 22 220 190 
 996 | 428 21 220 22 
 997 | 429 208 218 209 
 998 | 430 41 194 40 
 999 | 431 20 198 21 
1000 | 432 194 219 40 
1001 | 433 198 220 21 
1002 | 434 203 205 201 
1003 | 435 48 192 47 
1004 | 436 34 221 33 
1005 | 437 38 222 37 
1006 | 438 47 223 46 
1007 | 439 49 224 48 
1008 | 440 196 221 34 
1009 | 441 193 222 38 
1010 | 442 192 223 47 
1011 | 443 48 224 192 
1012 | 444 217 225 196 
1013 | 445 28 232 4 
1014 | 446 4 232 49 
1015 | 447 5 233 36 
1016 | 448 6 234 45 
1017 | 449 37 233 5 
1018 | 450 46 234 6 
1019 | 451 33 235 3 
1020 | 452 3 235 20 
1021 | 453 199 236 216 
1022 | 454 209 227 207 
1023 | 455 194 229 228 
1024 | 456 194 228 219 
1025 | 457 210 229 194 
1026 | 458 219 227 193 
1027 | 459 207 228 204 
1028 | 460 228 229 204 
1029 | 461 204 229 203 
1030 | 462 210 230 229 
1031 | 463 195 230 210 
1032 | 464 229 230 203 
1033 | 465 219 228 227 
1034 | 466 227 228 207 
1035 | 467 216 231 195 
1036 | 468 203 230 205 
1037 | 469 195 231 230 
1038 | 470 205 231 212 
1039 | 471 230 231 205 
1040 | 472 192 226 223 
1041 | 473 212 239 213 
1042 | 474 218 238 209 
1043 | 475 193 237 222 
1044 | 476 196 238 221 
1045 | 477 224 239 192 
1046 | 478 226 239 212 
1047 | 479 209 238 225 
1048 | 480 221 235 33 
1049 | 481 49 232 224 
1050 | 482 222 233 37 
1051 | 483 223 234 46 
1052 | 484 214 232 28 
1053 | 485 36 233 217 
1054 | 486 45 234 199 
1055 | 487 20 235 198 
1056 | 488 198 235 221 
1057 | 489 217 233 222 
1058 | 490 199 234 223 
1059 | 491 224 232 214 
1060 | 492 225 237 209 
1061 | 493 212 236 226 
1062 | 494 231 236 212 
1063 | 495 209 237 227 
1064 | 496 225 238 196 
1065 | 497 192 239 226 
1066 | 498 216 236 231 
1067 | 499 227 237 193 
1068 | $EndElements
1069 | 


--------------------------------------------------------------------------------
/build/elements.txt:
--------------------------------------------------------------------------------
   1 | 97
   2 | 98
   3 | 170
   4 | 96
   5 | 22
   6 | 128
   7 | 98
   8 | 18
   9 | 169
  10 | 153
  11 | 95
  12 | 154
  13 | 118
  14 | 147
  15 | 148
  16 | 1
  17 | 12
  18 | 181
  19 | 114
  20 | 115
  21 | 116
  22 | 142
  23 | 83
  24 | 171
  25 | 100
  26 | 99
  27 | 101
  28 | 94
  29 | 122
  30 | 123
  31 | 32
  32 | 1
  33 | 181
  34 | 113
  35 | 76
  36 | 151
  37 | 149
  38 | 147
  39 | 176
  40 | 153
  41 | 155
  42 | 177
  43 | 122
  44 | 94
  45 | 124
  46 | 172
  47 | 110
  48 | 173
  49 | 95
  50 | 96
  51 | 154
  52 | 99
  53 | 113
  54 | 163
  55 | 119
  56 | 142
  57 | 171
  58 | 102
  59 | 101
  60 | 103
  61 | 96
  62 | 128
  63 | 143
  64 | 142
  65 | 119
  66 | 174
  67 | 82
  68 | 141
  69 | 187
  70 | 114
  71 | 116
  72 | 127
  73 | 22
  74 | 114
  75 | 128
  76 | 113
  77 | 150
  78 | 164
  79 | 126
  80 | 172
  81 | 173
  82 | 101
  83 | 99
  84 | 131
  85 | 94
  86 | 123
  87 | 132
  88 | 117
  89 | 156
  90 | 169
  91 | 80
  92 | 25
  93 | 125
  94 | 117
  95 | 155
  96 | 156
  97 | 103
  98 | 101
  99 | 139
 100 | 74
 101 | 99
 102 | 100
 103 | 114
 104 | 23
 105 | 115
 106 | 56
 107 | 8
 108 | 149
 109 | 9
 110 | 61
 111 | 155
 112 | 107
 113 | 105
 114 | 119
 115 | 96
 116 | 143
 117 | 157
 118 | 97
 119 | 149
 120 | 176
 121 | 155
 122 | 117
 123 | 177
 124 | 121
 125 | 84
 126 | 123
 127 | 17
 128 | 18
 129 | 98
 130 | 75
 131 | 76
 132 | 113
 133 | 21
 134 | 22
 135 | 96
 136 | 141
 137 | 78
 138 | 174
 139 | 152
 140 | 153
 141 | 154
 142 | 147
 143 | 146
 144 | 148
 145 | 80
 146 | 125
 147 | 136
 148 | 11
 149 | 100
 150 | 182
 151 | 147
 152 | 118
 153 | 176
 154 | 95
 155 | 153
 156 | 177
 157 | 80
 158 | 134
 159 | 185
 160 | 95
 161 | 21
 162 | 96
 163 | 97
 164 | 17
 165 | 98
 166 | 99
 167 | 75
 168 | 113
 169 | 134
 170 | 80
 171 | 136
 172 | 115
 173 | 24
 174 | 185
 175 | 20
 176 | 21
 177 | 95
 178 | 16
 179 | 17
 180 | 97
 181 | 74
 182 | 75
 183 | 99
 184 | 76
 185 | 77
 186 | 118
 187 | 18
 188 | 19
 189 | 117
 190 | 22
 191 | 23
 192 | 114
 193 | 24
 194 | 25
 195 | 80
 196 | 15
 197 | 11
 198 | 79
 199 | 100
 200 | 101
 201 | 102
 202 | 23
 203 | 24
 204 | 115
 205 | 11
 206 | 74
 207 | 100
 208 | 18
 209 | 117
 210 | 169
 211 | 168
 212 | 81
 213 | 178
 214 | 81
 215 | 14
 216 | 126
 217 | 124
 218 | 168
 219 | 178
 220 | 15
 221 | 79
 222 | 126
 223 | 156
 224 | 86
 225 | 159
 226 | 106
 227 | 107
 228 | 134
 229 | 78
 230 | 142
 231 | 174
 232 | 105
 233 | 103
 234 | 135
 235 | 123
 236 | 84
 237 | 138
 238 | 76
 239 | 118
 240 | 151
 241 | 26
 242 | 84
 243 | 125
 244 | 25
 245 | 26
 246 | 125
 247 | 135
 248 | 141
 249 | 174
 250 | 87
 251 | 149
 252 | 170
 253 | 110
 254 | 112
 255 | 173
 256 | 128
 257 | 127
 258 | 189
 259 | 109
 260 | 106
 261 | 134
 262 | 116
 263 | 115
 264 | 120
 265 | 149
 266 | 97
 267 | 170
 268 | 104
 269 | 103
 270 | 105
 271 | 150
 272 | 113
 273 | 151
 274 | 161
 275 | 87
 276 | 170
 277 | 119
 278 | 105
 279 | 174
 280 | 79
 281 | 11
 282 | 182
 283 | 14
 284 | 15
 285 | 126
 286 | 116
 287 | 120
 288 | 133
 289 | 141
 290 | 135
 291 | 187
 292 | 24
 293 | 80
 294 | 185
 295 | 150
 296 | 88
 297 | 164
 298 | 143
 299 | 91
 300 | 157
 301 | 13
 302 | 14
 303 | 81
 304 | 26
 305 | 27
 306 | 84
 307 | 30
 308 | 31
 309 | 94
 310 | 84
 311 | 121
 312 | 186
 313 | 106
 314 | 105
 315 | 107
 316 | 30
 317 | 94
 318 | 132
 319 | 3
 320 | 20
 321 | 130
 322 | 2
 323 | 16
 324 | 129
 325 | 19
 326 | 3
 327 | 130
 328 | 77
 329 | 2
 330 | 129
 331 | 84
 332 | 27
 333 | 138
 334 | 29
 335 | 30
 336 | 132
 337 | 107
 338 | 119
 339 | 133
 340 | 126
 341 | 79
 342 | 172
 343 | 120
 344 | 107
 345 | 133
 346 | 109
 347 | 134
 348 | 136
 349 | 122
 350 | 121
 351 | 123
 352 | 131
 353 | 99
 354 | 163
 355 | 87
 356 | 56
 357 | 149
 358 | 61
 359 | 86
 360 | 155
 361 | 125
 362 | 84
 363 | 186
 364 | 102
 365 | 103
 366 | 104
 367 | 112
 368 | 121
 369 | 122
 370 | 20
 371 | 95
 372 | 130
 373 | 16
 374 | 97
 375 | 129
 376 | 118
 377 | 77
 378 | 129
 379 | 117
 380 | 19
 381 | 130
 382 | 27
 383 | 28
 384 | 138
 385 | 13
 386 | 81
 387 | 168
 388 | 134
 389 | 120
 390 | 185
 391 | 132
 392 | 123
 393 | 138
 394 | 107
 395 | 120
 396 | 134
 397 | 136
 398 | 125
 399 | 186
 400 | 12
 401 | 168
 402 | 181
 403 | 105
 404 | 135
 405 | 174
 406 | 81
 407 | 126
 408 | 173
 409 | 112
 410 | 111
 411 | 121
 412 | 31
 413 | 32
 414 | 188
 415 | 4
 416 | 29
 417 | 183
 418 | 28
 419 | 4
 420 | 183
 421 | 154
 422 | 96
 423 | 157
 424 | 133
 425 | 119
 426 | 171
 427 | 110
 428 | 111
 429 | 112
 430 | 163
 431 | 113
 432 | 164
 433 | 118
 434 | 148
 435 | 151
 436 | 122
 437 | 124
 438 | 178
 439 | 110
 440 | 109
 441 | 111
 442 | 156
 443 | 159
 444 | 160
 445 | 116
 446 | 133
 447 | 137
 448 | 112
 449 | 122
 450 | 178
 451 | 101
 452 | 131
 453 | 180
 454 | 162
 455 | 7
 456 | 180
 457 | 111
 458 | 109
 459 | 136
 460 | 98
 461 | 160
 462 | 161
 463 | 108
 464 | 109
 465 | 110
 466 | 171
 467 | 83
 468 | 175
 469 | 12
 470 | 13
 471 | 168
 472 | 72
 473 | 71
 474 | 141
 475 | 73
 476 | 72
 477 | 82
 478 | 71
 479 | 70
 480 | 78
 481 | 69
 482 | 68
 483 | 83
 484 | 70
 485 | 69
 486 | 142
 487 | 104
 488 | 105
 489 | 106
 490 | 127
 491 | 116
 492 | 184
 493 | 129
 494 | 97
 495 | 176
 496 | 117
 497 | 130
 498 | 177
 499 | 130
 500 | 95
 501 | 177
 502 | 118
 503 | 129
 504 | 176
 505 | 108
 506 | 106
 507 | 109
 508 | 32
 509 | 181
 510 | 188
 511 | 114
 512 | 127
 513 | 128
 514 | 172
 515 | 79
 516 | 179
 517 | 124
 518 | 94
 519 | 188
 520 | 175
 521 | 10
 522 | 184
 523 | 7
 524 | 73
 525 | 140
 526 | 68
 527 | 10
 528 | 175
 529 | 173
 530 | 112
 531 | 178
 532 | 104
 533 | 106
 534 | 108
 535 | 111
 536 | 136
 537 | 186
 538 | 10
 539 | 145
 540 | 184
 541 | 116
 542 | 137
 543 | 184
 544 | 66
 545 | 65
 546 | 144
 547 | 65
 548 | 64
 549 | 91
 550 | 62
 551 | 9
 552 | 152
 553 | 8
 554 | 55
 555 | 146
 556 | 67
 557 | 66
 558 | 92
 559 | 60
 560 | 59
 561 | 159
 562 | 54
 563 | 53
 564 | 158
 565 | 53
 566 | 52
 567 | 88
 568 | 59
 569 | 58
 570 | 85
 571 | 61
 572 | 60
 573 | 86
 574 | 55
 575 | 54
 576 | 89
 577 | 51
 578 | 50
 579 | 90
 580 | 57
 581 | 56
 582 | 87
 583 | 63
 584 | 62
 585 | 93
 586 | 64
 587 | 63
 588 | 167
 589 | 58
 590 | 57
 591 | 166
 592 | 52
 593 | 51
 594 | 165
 595 | 10
 596 | 67
 597 | 145
 598 | 50
 599 | 7
 600 | 162
 601 | 131
 602 | 162
 603 | 180
 604 | 179
 605 | 79
 606 | 182
 607 | 139
 608 | 101
 609 | 180
 610 | 108
 611 | 110
 612 | 172
 613 | 94
 614 | 31
 615 | 188
 616 | 29
 617 | 132
 618 | 183
 619 | 82
 620 | 72
 621 | 141
 622 | 71
 623 | 78
 624 | 141
 625 | 78
 626 | 70
 627 | 142
 628 | 69
 629 | 83
 630 | 142
 631 | 135
 632 | 103
 633 | 187
 634 | 108
 635 | 172
 636 | 179
 637 | 132
 638 | 138
 639 | 183
 640 | 100
 641 | 102
 642 | 182
 643 | 91
 644 | 143
 645 | 144
 646 | 143
 647 | 92
 648 | 144
 649 | 138
 650 | 28
 651 | 183
 652 | 137
 653 | 171
 654 | 175
 655 | 121
 656 | 111
 657 | 186
 658 | 160
 659 | 98
 660 | 169
 661 | 155
 662 | 86
 663 | 156
 664 | 8
 665 | 147
 666 | 149
 667 | 153
 668 | 9
 669 | 155
 670 | 73
 671 | 82
 672 | 140
 673 | 83
 674 | 68
 675 | 175
 676 | 140
 677 | 82
 678 | 187
 679 | 103
 680 | 139
 681 | 187
 682 | 137
 683 | 133
 684 | 171
 685 | 88
 686 | 150
 687 | 158
 688 | 150
 689 | 89
 690 | 158
 691 | 87
 692 | 161
 693 | 166
 694 | 161
 695 | 85
 696 | 166
 697 | 93
 698 | 157
 699 | 167
 700 | 157
 701 | 91
 702 | 167
 703 | 90
 704 | 164
 705 | 165
 706 | 164
 707 | 88
 708 | 165
 709 | 102
 710 | 104
 711 | 179
 712 | 102
 713 | 179
 714 | 182
 715 | 156
 716 | 160
 717 | 169
 718 | 92
 719 | 143
 720 | 189
 721 | 104
 722 | 108
 723 | 179
 724 | 168
 725 | 124
 726 | 181
 727 | 148
 728 | 89
 729 | 150
 730 | 160
 731 | 85
 732 | 161
 733 | 93
 734 | 154
 735 | 157
 736 | 90
 737 | 163
 738 | 164
 739 | 120
 740 | 115
 741 | 185
 742 | 7
 743 | 140
 744 | 180
 745 | 81
 746 | 173
 747 | 178
 748 | 145
 749 | 127
 750 | 184
 751 | 148
 752 | 150
 753 | 151
 754 | 162
 755 | 131
 756 | 163
 757 | 98
 758 | 161
 759 | 170
 760 | 8
 761 | 146
 762 | 147
 763 | 152
 764 | 9
 765 | 153
 766 | 65
 767 | 91
 768 | 144
 769 | 92
 770 | 66
 771 | 144
 772 | 91
 773 | 64
 774 | 167
 775 | 93
 776 | 62
 777 | 152
 778 | 59
 779 | 85
 780 | 159
 781 | 53
 782 | 88
 783 | 158
 784 | 86
 785 | 60
 786 | 159
 787 | 89
 788 | 54
 789 | 158
 790 | 85
 791 | 58
 792 | 166
 793 | 88
 794 | 52
 795 | 165
 796 | 55
 797 | 89
 798 | 146
 799 | 51
 800 | 90
 801 | 165
 802 | 57
 803 | 87
 804 | 166
 805 | 63
 806 | 93
 807 | 167
 808 | 127
 809 | 145
 810 | 189
 811 | 67
 812 | 92
 813 | 145
 814 | 90
 815 | 50
 816 | 162
 817 | 137
 818 | 175
 819 | 184
 820 | 159
 821 | 85
 822 | 160
 823 | 93
 824 | 152
 825 | 154
 826 | 146
 827 | 89
 828 | 148
 829 | 90
 830 | 162
 831 | 163
 832 | 139
 833 | 140
 834 | 187
 835 | 181
 836 | 124
 837 | 188
 838 | 145
 839 | 92
 840 | 189
 841 | 143
 842 | 128
 843 | 189
 844 | 140
 845 | 139
 846 | 180
 847 | 213
 848 | 224
 849 | 214
 850 | 198
 851 | 221
 852 | 218
 853 | 222
 854 | 225
 855 | 217
 856 | 223
 857 | 226
 858 | 199
 859 | 221
 860 | 238
 861 | 218
 862 | 213
 863 | 239
 864 | 224
 865 | 222
 866 | 237
 867 | 225
 868 | 226
 869 | 236
 870 | 199
 871 | 27
 872 | 214
 873 | 28
 874 | 190
 875 | 208
 876 | 206
 877 | 213
 878 | 214
 879 | 197
 880 | 197
 881 | 214
 882 | 27
 883 | 190
 884 | 220
 885 | 208
 886 | 206
 887 | 215
 888 | 190
 889 | 208
 890 | 220
 891 | 218
 892 | 211
 893 | 213
 894 | 197
 895 | 24
 896 | 191
 897 | 25
 898 | 201
 899 | 211
 900 | 200
 901 | 191
 902 | 201
 903 | 200
 904 | 205
 905 | 211
 906 | 201
 907 | 190
 908 | 215
 909 | 23
 910 | 23
 911 | 215
 912 | 24
 913 | 191
 914 | 200
 915 | 25
 916 | 42
 917 | 210
 918 | 41
 919 | 218
 920 | 220
 921 | 198
 922 | 202
 923 | 215
 924 | 206
 925 | 191
 926 | 202
 927 | 201
 928 | 202
 929 | 206
 930 | 204
 931 | 41
 932 | 210
 933 | 194
 934 | 24
 935 | 215
 936 | 191
 937 | 26
 938 | 200
 939 | 197
 940 | 206
 941 | 207
 942 | 204
 943 | 35
 944 | 217
 945 | 196
 946 | 202
 947 | 204
 948 | 203
 949 | 25
 950 | 200
 951 | 26
 952 | 199
 953 | 216
 954 | 44
 955 | 200
 956 | 211
 957 | 197
 958 | 206
 959 | 208
 960 | 207
 961 | 202
 962 | 203
 963 | 201
 964 | 205
 965 | 212
 966 | 211
 967 | 36
 968 | 217
 969 | 35
 970 | 35
 971 | 196
 972 | 34
 973 | 26
 974 | 197
 975 | 27
 976 | 43
 977 | 216
 978 | 195
 979 | 43
 980 | 195
 981 | 42
 982 | 195
 983 | 210
 984 | 42
 985 | 45
 986 | 199
 987 | 44
 988 | 44
 989 | 216
 990 | 43
 991 | 191
 992 | 215
 993 | 202
 994 | 212
 995 | 213
 996 | 211
 997 | 39
 998 | 193
 999 | 38
1000 | 22
1001 | 190
1002 | 23
1003 | 208
1004 | 209
1005 | 207
1006 | 40
1007 | 219
1008 | 39
1009 | 39
1010 | 219
1011 | 193
1012 | 22
1013 | 220
1014 | 190
1015 | 21
1016 | 220
1017 | 22
1018 | 208
1019 | 218
1020 | 209
1021 | 41
1022 | 194
1023 | 40
1024 | 20
1025 | 198
1026 | 21
1027 | 194
1028 | 219
1029 | 40
1030 | 198
1031 | 220
1032 | 21
1033 | 203
1034 | 205
1035 | 201
1036 | 48
1037 | 192
1038 | 47
1039 | 34
1040 | 221
1041 | 33
1042 | 38
1043 | 222
1044 | 37
1045 | 47
1046 | 223
1047 | 46
1048 | 49
1049 | 224
1050 | 48
1051 | 196
1052 | 221
1053 | 34
1054 | 193
1055 | 222
1056 | 38
1057 | 192
1058 | 223
1059 | 47
1060 | 48
1061 | 224
1062 | 192
1063 | 217
1064 | 225
1065 | 196
1066 | 28
1067 | 232
1068 | 4
1069 | 4
1070 | 232
1071 | 49
1072 | 5
1073 | 233
1074 | 36
1075 | 6
1076 | 234
1077 | 45
1078 | 37
1079 | 233
1080 | 5
1081 | 46
1082 | 234
1083 | 6
1084 | 33
1085 | 235
1086 | 3
1087 | 3
1088 | 235
1089 | 20
1090 | 199
1091 | 236
1092 | 216
1093 | 209
1094 | 227
1095 | 207
1096 | 194
1097 | 229
1098 | 228
1099 | 194
1100 | 228
1101 | 219
1102 | 210
1103 | 229
1104 | 194
1105 | 219
1106 | 227
1107 | 193
1108 | 207
1109 | 228
1110 | 204
1111 | 228
1112 | 229
1113 | 204
1114 | 204
1115 | 229
1116 | 203
1117 | 210
1118 | 230
1119 | 229
1120 | 195
1121 | 230
1122 | 210
1123 | 229
1124 | 230
1125 | 203
1126 | 219
1127 | 228
1128 | 227
1129 | 227
1130 | 228
1131 | 207
1132 | 216
1133 | 231
1134 | 195
1135 | 203
1136 | 230
1137 | 205
1138 | 195
1139 | 231
1140 | 230
1141 | 205
1142 | 231
1143 | 212
1144 | 230
1145 | 231
1146 | 205
1147 | 192
1148 | 226
1149 | 223
1150 | 212
1151 | 239
1152 | 213
1153 | 218
1154 | 238
1155 | 209
1156 | 193
1157 | 237
1158 | 222
1159 | 196
1160 | 238
1161 | 221
1162 | 224
1163 | 239
1164 | 192
1165 | 226
1166 | 239
1167 | 212
1168 | 209
1169 | 238
1170 | 225
1171 | 221
1172 | 235
1173 | 33
1174 | 49
1175 | 232
1176 | 224
1177 | 222
1178 | 233
1179 | 37
1180 | 223
1181 | 234
1182 | 46
1183 | 214
1184 | 232
1185 | 28
1186 | 36
1187 | 233
1188 | 217
1189 | 45
1190 | 234
1191 | 199
1192 | 20
1193 | 235
1194 | 198
1195 | 198
1196 | 235
1197 | 221
1198 | 217
1199 | 233
1200 | 222
1201 | 199
1202 | 234
1203 | 223
1204 | 224
1205 | 232
1206 | 214
1207 | 225
1208 | 237
1209 | 209
1210 | 212
1211 | 236
1212 | 226
1213 | 231
1214 | 236
1215 | 212
1216 | 209
1217 | 237
1218 | 227
1219 | 225
1220 | 238
1221 | 196
1222 | 192
1223 | 239
1224 | 226
1225 | 216
1226 | 236
1227 | 231
1228 | 227
1229 | 237
1230 | 193
1231 | 


--------------------------------------------------------------------------------
/build/nodes.txt:
--------------------------------------------------------------------------------
  1 | -1
  2 | 0
  3 | 0
  4 | 1
  5 | 0
  6 | 0
  7 | 1
  8 | 1
  9 | 0
 10 | -1
 11 | 1
 12 | 0
 13 | 1
 14 | 2
 15 | 0
 16 | -1
 17 | 2
 18 | 0
 19 | 0.333333
 20 | 0.333333
 21 | 0
 22 | 0.666667
 23 | 0.333333
 24 | 0
 25 | 0.666667
 26 | 0.666667
 27 | 0
 28 | 0.333333
 29 | 0.666667
 30 | 0
 31 | 0
 32 | 0
 33 | 0
 34 | -0.8
 35 | 0
 36 | 0
 37 | -0.6
 38 | 0
 39 | 0
 40 | -0.4
 41 | 0
 42 | 0
 43 | -0.2
 44 | 0
 45 | 0
 46 | 1
 47 | 0.2
 48 | 0
 49 | 1
 50 | 0.4
 51 | 0
 52 | 1
 53 | 0.6
 54 | 0
 55 | 1
 56 | 0.8
 57 | 0
 58 | 0.8
 59 | 1
 60 | 0
 61 | 0.6
 62 | 1
 63 | 0
 64 | 0.4
 65 | 1
 66 | 0
 67 | 0.2
 68 | 1
 69 | 0
 70 | 2.75224e-12
 71 | 1
 72 | 0
 73 | -0.2
 74 | 1
 75 | 0
 76 | -0.4
 77 | 1
 78 | 0
 79 | -0.6
 80 | 1
 81 | 0
 82 | -0.8
 83 | 1
 84 | 0
 85 | -1
 86 | 0.8
 87 | 0
 88 | -1
 89 | 0.6
 90 | 0
 91 | -1
 92 | 0.4
 93 | 0
 94 | -1
 95 | 0.2
 96 | 0
 97 | 1
 98 | 1.2
 99 | 0
100 | 1
101 | 1.4
102 | 0
103 | 1
104 | 1.6
105 | 0
106 | 1
107 | 1.8
108 | 0
109 | 0.8
110 | 2
111 | 0
112 | 0.6
113 | 2
114 | 0
115 | 0.4
116 | 2
117 | 0
118 | 0.2
119 | 2
120 | 0
121 | 2.75224e-12
122 | 2
123 | 0
124 | -0.2
125 | 2
126 | 0
127 | -0.4
128 | 2
129 | 0
130 | -0.6
131 | 2
132 | 0
133 | -0.8
134 | 2
135 | 0
136 | -1
137 | 1.8
138 | 0
139 | -1
140 | 1.6
141 | 0
142 | -1
143 | 1.4
144 | 0
145 | -1
146 | 1.2
147 | 0
148 | 0.380952
149 | 0.333333
150 | 0
151 | 0.428571
152 | 0.333333
153 | 0
154 | 0.47619
155 | 0.333333
156 | 0
157 | 0.52381
158 | 0.333333
159 | 0
160 | 0.571429
161 | 0.333333
162 | 0
163 | 0.619048
164 | 0.333333
165 | 0
166 | 0.666667
167 | 0.380952
168 | 0
169 | 0.666667
170 | 0.428571
171 | 0
172 | 0.666667
173 | 0.47619
174 | 0
175 | 0.666667
176 | 0.52381
177 | 0
178 | 0.666667
179 | 0.571429
180 | 0
181 | 0.666667
182 | 0.619048
183 | 0
184 | 0.619048
185 | 0.666667
186 | 0
187 | 0.571429
188 | 0.666667
189 | 0
190 | 0.52381
191 | 0.666667
192 | 0
193 | 0.47619
194 | 0.666667
195 | 0
196 | 0.428571
197 | 0.666667
198 | 0
199 | 0.380952
200 | 0.666667
201 | 0
202 | 0.333333
203 | 0.619048
204 | 0
205 | 0.333333
206 | 0.571429
207 | 0
208 | 0.333333
209 | 0.52381
210 | 0
211 | 0.333333
212 | 0.47619
213 | 0
214 | 0.333333
215 | 0.428571
216 | 0
217 | 0.333333
218 | 0.380952
219 | 0
220 | 0.2
221 | 0
222 | 0
223 | 0.4
224 | 0
225 | 0
226 | 0.6
227 | 0
228 | 0
229 | 0.8
230 | 0
231 | 0
232 | 0.295834
233 | 0.500351
234 | 0
235 | -0.108163
236 | 0.159207
237 | 0
238 | -0.110844
239 | 0.829471
240 | 0
241 | -0.487332
242 | 0.153754
243 | 0
244 | 0.297874
245 | 0.401828
246 | 0
247 | 0.296308
248 | 0.59817
249 | 0
250 | -0.491748
251 | 0.804764
252 | 0
253 | 0.712654
254 | 0.496032
255 | 0
256 | 0.710368
257 | 0.599781
258 | 0
259 | 0.719205
260 | 0.405727
261 | 0
262 | 0.5
263 | 0.288038
264 | 0
265 | 0.59127
266 | 0.287346
267 | 0
268 | 0.40873
269 | 0.287328
270 | 0
271 | 0.5
272 | 0.712474
273 | 0
274 | 0.40873
275 | 0.712689
276 | 0
277 | 0.59127
278 | 0.712654
279 | 0
280 | -0.821199
281 | 0.500382
282 | 0
283 | 0.7
284 | 0.852284
285 | 0
286 | 0.531577
287 | 0.851696
288 | 0
289 | 0.861551
290 | 0.306135
291 | 0
292 | 0.864895
293 | 0.458517
294 | 0
295 | 0.3
296 | 0.147716
297 | 0
298 | 0.138157
299 | 0.145019
300 | 0
301 | 0.221646
302 | 0.259695
303 | 0
304 | 0.097029
305 | 0.275541
306 | 0
307 | 0.173301
308 | 0.365621
309 | 0
310 | 0.0603549
311 | 0.385956
312 | 0
313 | 0.141034
314 | 0.476676
315 | 0
316 | 0.015529
317 | 0.506599
318 | 0
319 | 0.0976812
320 | 0.590594
321 | 0
322 | -0.0743718
323 | 0.408045
324 | 0
325 | -0.124317
326 | 0.540879
327 | 0
328 | -0.238026
329 | 0.42549
330 | 0
331 | -0.280183
332 | 0.574421
333 | 0
334 | -0.414812
335 | 0.450351
336 | 0
337 | 0.454406
338 | 0.144538
339 | 0
340 | 0.277516
341 | 0.860777
342 | 0
343 | 0.134536
344 | 0.856508
345 | 0
346 | 0.21992
347 | 0.746842
348 | 0
349 | 0.852284
350 | 0.7
351 | 0
352 | 0.701608
353 | 0.148335
354 | 0
355 | 0.195755
356 | 0.556625
357 | 0
358 | 0.101443
359 | 0.720037
360 | 0
361 | -0.458975
362 | 0.60969
363 | 0
364 | -0.600239
365 | 0.476534
366 | 0
367 | -0.650799
368 | 0.654912
369 | 0
370 | -0.727984
371 | 0.336688
372 | 0
373 | -0.285787
374 | 0.835765
375 | 0
376 | -0.290922
377 | 0.15449
378 | 0
379 | 0.321991
380 | 0.767137
381 | 0
382 | 0.397949
383 | 0.826588
384 | 0
385 | 0.857286
386 | 0.145942
387 | 0
388 | 0.855036
389 | 0.855036
390 | 0
391 | 0.309932
392 | 0.247938
393 | 0
394 | -0.832634
395 | 0.7
396 | 0
397 | 0.180585
398 | 0.650225
399 | 0
400 | -0.0244653
401 | 0.669285
402 | 0
403 | 0.213504
404 | 0.439696
405 | 0
406 | -0.187249
407 | 0.68902
408 | 0
409 | 0.246287
410 | 0.670627
411 | 0
412 | -0.7
413 | 0.838101
414 | 0
415 | 0.243095
416 | 0.335661
417 | 0
418 | 0.289345
419 | 0.357937
420 | 0
421 | 0.277244
422 | 0.448355
423 | 0
424 | 0.273313
425 | 0.551719
426 | 0
427 | 0.45621
428 | 0.763263
429 | 0
430 | 0.453941
431 | 0.704352
432 | 0
433 | 0.353398
434 | 0.713673
435 | 0
436 | 0.638716
437 | 0.295311
438 | 0
439 | 0.691729
440 | 0.269893
441 | 0
442 | 0.619048
443 | 0.234528
444 | 0
445 | 0.729821
446 | 0.343183
447 | 0
448 | 0.544276
449 | 0.237636
450 | 0
451 | 0.583868
452 | 0.153007
453 | 0
454 | 0.638716
455 | 0.704689
456 | 0
457 | 0.691538
458 | 0.73018
459 | 0
460 | 0.616786
461 | 0.768802
462 | 0
463 | 0.732745
464 | 0.664161
465 | 0
466 | 0.777139
467 | 0.601533
468 | 0
469 | 0.540577
470 | 0.752799
471 | 0
472 | 0.546157
473 | 0.295937
474 | 0
475 | 0.716494
476 | 0.552732
477 | 0
478 | 0.781407
479 | 0.518214
480 | 0
481 | 0.764665
482 | 0.45188
483 | 0
484 | 0.350637
485 | 0.290353
486 | 0
487 | 0.379348
488 | 0.226093
489 | 0
490 | 0.45534
491 | 0.243877
492 | 0
493 | 0.453766
494 | 0.297182
495 | 0
496 | 0.705972
497 | 0.45168
498 | 0
499 | 0.545417
500 | 0.702252
501 | 0
502 | -0.666779
503 | 0.171661
504 | 0
505 | 0.855145
506 | 0.575653
507 | 0
508 | 0.804266
509 | 0.382821
510 | 0
511 | 0.247453
512 | 0.611704
513 | 0
514 | -0.186136
515 | 0.290075
516 | 0
517 | -0.360124
518 | 0.299609
519 | 0
520 | 0.229942
521 | 0.499502
522 | 0
523 | 0.291857
524 | 0.645416
525 | 0
526 | 0.774354
527 | 0.232161
528 | 0
529 | 0.773293
530 | 0.769385
531 | 0
532 | -0.543349
533 | 0.321224
534 | 0
535 | -0.0323435
536 | 0.284237
537 | 0
538 | 0.29064
539 | 0.304938
540 | 0
541 | -0.842645
542 | 0.177964
543 | 0
544 | 0.0172267
545 | 0.186597
546 | 0
547 | -0.854503
548 | 0.854503
549 | 0
550 | 0.294446
551 | 0.703958
552 | 0
553 | 0.0172251
554 | 0.813422
555 | 0
556 | -0.343279
557 | 0.698021
558 | 0
559 | 0.246047
560 | 0.391109
561 | 0
562 | -0.873435
563 | 0.333414
564 | 0
565 | 0.38574
566 | 0.758364
567 | 0
568 | 0.3
569 | 1.17321
570 | 0
571 | -0.1
572 | 1.17321
573 | -0
574 | -0.837245
575 | 1.50967
576 | -0
577 | 0.489686
578 | 1.83645
579 | 0
580 | 0.1
581 | 1.82679
582 | 0
583 | -0.3
584 | 1.82679
585 | -0
586 | 0.836293
587 | 1.51133
588 | 0
589 | -0.5
590 | 1.17321
591 | -0
592 | 0.699153
593 | 1.16935
594 | 0
595 | -0.686115
596 | 1.8081
597 | -0
598 | -0.3
599 | 1.17321
600 | -0
601 | -0.2
602 | 1.34641
603 | -0
604 | 2.56615e-12
605 | 1.34641
606 | 0
607 | -0.1
608 | 1.51962
609 | -0
610 | 0.1
611 | 1.51962
612 | 0
613 | -0.3
614 | 1.51962
615 | -0
616 | 0.2
617 | 1.34641
618 | 0
619 | 0.3
620 | 1.51962
621 | 0
622 | 0.4
623 | 1.34641
624 | 0
625 | 0.515205
626 | 1.51858
627 | 0
628 | -0.1
629 | 1.82679
630 | -0
631 | -0.4
632 | 1.34641
633 | -0
634 | -0.515215
635 | 1.51855
636 | -0
637 | -0.606237
638 | 1.32988
639 | -0
640 | -0.699665
641 | 1.16974
642 | -0
643 | 0.1
644 | 1.17321
645 | 0
646 | -0.489011
647 | 1.83239
648 | -0
649 | 0.830618
650 | 1.69989
651 | 0
652 | 0.606085
653 | 1.32926
654 | 0
655 | 0.3
656 | 1.82679
657 | 0
658 | 0.500873
659 | 1.16971
660 | 0
661 | 0.812117
662 | 1.31637
663 | 0
664 | 0.684376
665 | 1.81088
666 | 0
667 | -0.840959
668 | 1.69521
669 | -0
670 | -0.812355
671 | 1.31626
672 | -0
673 | 0.684019
674 | 1.61419
675 | 0
676 | -0.686206
677 | 1.61265
678 | -0
679 | 0.392136
680 | 1.67682
681 | 0
682 | 0.2
683 | 1.67321
684 | 0
685 | 2.87186e-12
686 | 1.67321
687 | 0
688 | -0.2
689 | 1.67321
690 | -0
691 | -0.39203
692 | 1.67613
693 | -0
694 | -0.85359
695 | 1.14641
696 | -0
697 | 0.85359
698 | 1.85359
699 | 0
700 | -0.85359
701 | 1.85359
702 | -0
703 | 0.85359
704 | 1.14641
705 | 0
706 | -0.549453
707 | 1.68709
708 | -0
709 | 0.553084
710 | 1.69138
711 | 0
712 | 0.690744
713 | 1.45795
714 | 0
715 | -0.68975
716 | 1.45812
717 | -0
718 | 


--------------------------------------------------------------------------------
/demo/README:
--------------------------------------------------------------------------------
1 | 
2 | g++ *.cpp -std=c++11 -Iinclude -lgmsh -o exe
3 | 
4 | ./exe
5 | 
6 | demo.cpp        : a test case
7 | oneDExample.cpp : get a one-dimensional grid(Segment)
8 | twoDExample.cpp : get a two-dimensional grid(Triangle)
9 | 


--------------------------------------------------------------------------------
/demo/demo.cpp:
--------------------------------------------------------------------------------
  1 | #include <iostream>
  2 | #include <gmsh.h>
  3 | #include <vector> 
  4 | #include <fstream>
  5 | #include <iomanip>
  6 | using namespace std;
  7 | 
  8 | 
  9 | void writeNodes(vector<double>& nodes, string s)
 10 | {
 11 | 	string filename = s + ".txt";
 12 | 	ofstream outfileb(filename.c_str(), ios::out);
 13 | 	for(int i = 0; i < nodes.size(); i++)
 14 | 	{
 15 | 		outfileb<<(nodes[i])<<endl;
 16 | 	}
 17 | 	outfileb.close();
 18 | }
 19 | 
 20 | 
 21 | void writeElements(vector<std::vector<std::size_t> >& elements, string s)
 22 | {
 23 | 	string filename = s + ".txt";
 24 | 	ofstream outfileb(filename.c_str(), ios::out);
 25 | 	for(int i = 0; i < elements[1].size(); i++)
 26 | 	{
 27 | 		outfileb<<(elements[1][i])<<endl;
 28 | 	}
 29 | 	outfileb.close();
 30 | }
 31 | 
 32 | 
 33 | int main(int argc, char **argv)
 34 | {
 35 | 	// Before using any functions in the C++ API, Gmsh must be initialized.
 36 | 	gmsh::initialize();
 37 | 	// By default Gmsh will not print out any messages: in order to output
 38 | 	// messages on the terminal, just set the standard Gmsh option
 39 | 	// "General.Terminal" (same format and meaning as in .geo files) using
 40 | 	// gmsh::option::setNumber():
 41 | 	//gmsh::option::setNumber("General.Terminal", 1);
 42 | 
 43 | 	// This adds a new model, named "t1". If gmsh::model::add() is not called, a
 44 | 	// new default (unnamed) model will be created on the fly, if necessary.
 45 | 	gmsh::model::add("Test1");
 46 | 	// The C++ API provides direct access to the internal CAD kernels. The
 47 |     // built-in CAD kernel was used in t1.geo: the corresponding API functions
 48 |     // live in the "gmsh::model::geo" namespace. To create geometrical points with
 49 |     // the built-in CAD kernel, one thus uses gmsh::model::geo::addPoint():
 50 |     //
 51 |     // - the first 3 arguments are the point coordinates (x, y, z)
 52 |     //
 53 |     // - the next (optional) argument is the target mesh size close to the point
 54 |     //
 55 |     // - the last (optional) argument is the point tag
 56 | 
 57 | 
 58 | 
 59 |   	double lcc = 0.05;//2e-1;
 60 |   	gmsh::model::geo::addPoint(1.0/3, 1.0/3,  0, lcc, 100);
 61 |   	gmsh::model::geo::addPoint(2.0/3, 1.0/3,  0, lcc, 101);
 62 |   	gmsh::model::geo::addPoint(2.0/3, 2.0/3,  0, lcc, 102);
 63 |   	gmsh::model::geo::addPoint(1.0/3, 2.0/3,  0, lcc, 103);
 64 |   	gmsh::model::geo::addLine(100, 101, 100);
 65 |   	gmsh::model::geo::addLine(101, 102, 101);
 66 |   	gmsh::model::geo::addLine(102, 103, 102);
 67 |   	gmsh::model::geo::addLine(103, 100, 103);
 68 |   	gmsh::model::geo::addCurveLoop({100, 101, 102, 103}, 100);
 69 |   	//gmsh::model::geo::addPlaneSurface({100}, 100);
 70 | 
 71 | 
 72 | 
 73 |   	double lc = 0.2;//2e-1;
 74 |   	gmsh::model::geo::addPoint(-1, 0,  0, lc, 1);
 75 |   	gmsh::model::geo::addPoint(1,  0,  0, lc, 2);
 76 |   	gmsh::model::geo::addPoint(1,  1,  0, lc, 3);
 77 |   	gmsh::model::geo::addPoint(-1, 1,  0, lc, 4);
 78 |   	gmsh::model::geo::addPoint(0,  0,  0, lc, 555);
 79 | 
 80 | 
 81 |   	// The API to create lines with the built-in kernel follows the same
 82 |   	// conventions: the first 2 arguments are point tags, the last (optional one)
 83 |   	// is the line tag.
 84 |   	gmsh::model::geo::addLine(1, 555, 1);
 85 |   	gmsh::model::geo::addLine(555, 2, 555);
 86 |   	gmsh::model::geo::addLine(2, 3, 2);
 87 |   	gmsh::model::geo::addLine(3, 4, 3);
 88 |   	gmsh::model::geo::addLine(4, 1, 4);
 89 | 
 90 |   	// The philosophy to construct curve loops and surfaces is similar: the first
 91 |   	// argument is now a vector of integers.
 92 |   	gmsh::model::geo::addCurveLoop({1, 555, 2, 3, 4}, 1);
 93 |   	gmsh::model::geo::addPlaneSurface({1, -100}, 1);
 94 | 
 95 | 
 96 |   	gmsh::model::geo::addPoint(1,  2,  0, lc, 5);
 97 |   	gmsh::model::geo::addPoint(-1, 2,  0, lc, 6);
 98 |   	gmsh::model::geo::addLine(3, 5, 5);
 99 |   	gmsh::model::geo::addLine(5, 6, 6);
100 |   	gmsh::model::geo::addLine(6, 4, 7);
101 |   	gmsh::model::geo::addCurveLoop({-3, 5, 6, 7}, 2);
102 | 	gmsh::model::geo::addPlaneSurface({2}, 2);
103 | 
104 | 
105 | 
106 | 
107 |   	// Physical groups are defined by providing the dimension of the group (0 for
108 |   	// physical points, 1 for physical curves, 2 for physical surfaces and 3 for
109 |   	// phsyical volumes) followed by a vector of entity tags. The last (optional)
110 |   	// argument is the tag of the new group to create.
111 | 
112 | 	gmsh::model::setPhysicalName(2, 6, "My surface");
113 | 
114 | 	// Before it can be meshed, the internal CAD representation must be
115 | 	// synchronized with the Gmsh model, which will create the relevant Gmsh data
116 | 	// structures. This is achieved by the gmsh::model::geo::synchronize() API
117 | 	// call for the built-in CAD kernel. Synchronizations can be called at any
118 | 	// time, but they involve a non trivial amount of processing; so while you
119 | 	// could synchronize the internal CAD data after every CAD command, it is
120 | 	// usually better to minimize the number of synchronization points.
121 | 	gmsh::model::geo::synchronize();
122 | 
123 | 	// We can then generate a 2D mesh...
124 | 	gmsh::model::mesh::generate(2);
125 | 	//gmsh::model::mesh::refine();
126 | 
127 | 	// get mesh information
128 | 	// Nodes
129 | 	vector<double> nodes,y;
130 | 	vector<std::size_t> nodeTags;
131 | 	gmsh::model::mesh::getNodes(nodeTags, nodes, y, -2, -1, false, true);
132 | 	cout<<"The number of nodes: "<<nodes.size()/3<<endl;
133 | #if 0
134 | 	cout<<endl;
135 | 	for(int i = 0; i < nodes.size(); i++)
136 | 	{
137 | 		cout<<nodes[i]<<endl;
138 | 	}
139 | #endif
140 | 	writeNodes(nodes, "nodes");
141 | 
142 | 	//Elements
143 | 	vector<int> elementTypes;
144 | 	vector<std::vector<std::size_t> > elementTags, nodeTags2;
145 | 	gmsh::model::mesh::getElements(elementTypes, elementTags, nodeTags2, -2, -1);
146 | 	cout<<"The number of elements: "<<nodeTags2[1].size()/3<<endl;
147 | #if 0
148 | 	cout<<endl;
149 | 	for(int j = 0; j < nodeTags2[1].size(); j++)
150 | 	{
151 | 		cout<<nodeTags2[1][j]<<endl;
152 | 		
153 | 	}
154 | #endif
155 | 	writeElements(nodeTags2, "elements");
156 | 
157 | 	// ... and save it to disk
158 | 	gmsh::write("demo.msh");
159 | 
160 | 	// Remember that by default, if physical groups are defined, Gmsh will export
161 | 	// in the output mesh file only those elements that belong to at least one
162 | 	// physical group. To force Gmsh to save all elements, you can use
163 | 	//
164 | 	// gmsh::option::setNumber("Mesh.SaveAll", 1);
165 | 
166 | 	// This should be called at the end:
167 | 	gmsh::finalize();
168 | 	return 0;
169 | }
170 | 


--------------------------------------------------------------------------------
/demo/oneDExample.cpp:
--------------------------------------------------------------------------------
  1 | #include <gmsh.h>
  2 | #include <iostream>
  3 | #include <vector>
  4 | #include <fstream>
  5 | #include <iomanip>
  6 | #include <string>
  7 | #include <map>
  8 | #include <set>
  9 | 
 10 | using namespace std;
 11 | 
 12 | class Point
 13 | {
 14 | public:
 15 |     Point(const double x = 0, const double y = 0, const double z = 0):_x(x), _y(y), _z(z) { }
 16 |     ~Point() { }
 17 |     const double getX() const {return _x;}
 18 |     const double getY() const {return _y;}
 19 |     const double getZ() const {return _z;}
 20 |     void setX(const double x) {_x = x;}
 21 |     void setY(const double y) {_y = y;}
 22 |     void setZ(const double z) {_z = z;}
 23 |     bool operator < (const Point& p) const;
 24 |     
 25 | protected:
 26 |     double _x, _y, _z;
 27 | };
 28 | 
 29 | bool Point::operator < (const Point& p) const
 30 | {
 31 |     double px = p.getX();
 32 |     double py = p.getY();
 33 |     if(px == _x) {
 34 |         return py < _y;
 35 |     } else {
 36 |         return px < _x;
 37 |     }
 38 | }
 39 | 
 40 | 
 41 | class Polygon
 42 | {
 43 | public:
 44 |     Polygon() : _vertex(), _name() {}
 45 |     Polygon(const vector<Point>& vertex, const string& name = "poly") :_vertex(vertex), _name(name) {}
 46 |     ~Polygon() {}
 47 |     const vector<Point>& getVertex() const {return _vertex;}
 48 |     void clear() {_vertex.clear();}
 49 |     void append(const Point& p) {_vertex.push_back(p);}
 50 |     const string& getName() const {return _name;}
 51 |     void setPoints(const vector<Point>& vertex) {_vertex = vertex;}
 52 |     void setName(const string& name) {_name = name;}
 53 |   
 54 | private:
 55 |     vector<Point> _vertex;
 56 |     string _name;
 57 | };
 58 | 
 59 | class Topology
 60 | {
 61 | public:
 62 |     Topology(const int firstId = -3, const int secondId = -3) : _firstId(firstId), _secondId(secondId) {}
 63 |     ~Topology() {}
 64 |     const int getFirstId() const {return _firstId;}
 65 |     const int getSecondId() const {return _secondId;}
 66 |     void setFirstId(const int firstId) {_firstId = firstId;}
 67 |     void setSecondId(const int secondId) {_secondId= secondId;}
 68 |     bool operator < (const Topology& topo) const;
 69 |     
 70 | private:
 71 |     int _firstId, _secondId;
 72 | };
 73 | 
 74 | bool Topology::operator < (const Topology& topo) const
 75 | {
 76 |     int tof = topo.getFirstId();
 77 |     int tos = topo.getSecondId();
 78 |     if(tof == _firstId) {
 79 |         return tos < _secondId;
 80 |     } else {
 81 |         return tof < _firstId;
 82 |     }
 83 | }
 84 | 
 85 | class Segment
 86 | {
 87 | public:
 88 |     Segment() : _startPt(), _endPt() {}
 89 |     Segment(const Point& startPt, const Point& endPt) : _startPt(startPt), _endPt(endPt) {}
 90 |     ~Segment() {}
 91 |     const Point& getStartPt() const {return _startPt;}
 92 |     const Point& getEndPt() const {return _endPt;}
 93 |     void setStartPt(const Point& startPt) {_startPt = startPt;}
 94 |     void setEndPt(const Point& endPt) {_endPt = endPt;}
 95 |   
 96 | protected:
 97 |     Point _startPt, _endPt;
 98 | };
 99 | 
100 | class Element : public Segment
101 | {
102 | public:
103 |     Element() : Segment(), _topo() {}
104 |     Element(const Point& startPt, const Point& endPt, const Topology& topo) : Segment(startPt, endPt), _topo(topo) {}
105 |     ~Element() {}
106 |     
107 |     const Topology& getTopology() const {return _topo;}
108 |     void setTopology(const Topology& topo) {_topo = topo;}
109 |     void setTopology(const int firstId, const int secondId) {_topo = Topology(firstId, secondId);}
110 |     
111 | private:
112 |     Topology _topo;
113 | };
114 | 
115 | void generatePolygons(vector<Polygon>& region)
116 | {
117 |     Polygon poly;
118 |     double left, right, bot, top;
119 |     double x, y, z = 0;
120 |     Point p;
121 |     
122 |     poly.clear();
123 |     p = Point(0, 0, z);
124 |     poly.append(p);
125 |     p = Point(8, 0, z);
126 |     poly.append(p);
127 |     p = Point(8, 1, z);
128 |     poly.append(p);
129 |     p = Point(6, 1, z);
130 |     poly.append(p);
131 |     p = Point(5, 1, z);
132 |     poly.append(p);
133 |     p = Point(3, 1, z);
134 |     poly.append(p);
135 |     p = Point(2, 1, z);
136 |     poly.append(p);
137 |     p = Point(0, 1, z);
138 |     poly.append(p);
139 |     poly.setName("poly_1");
140 |     region.push_back(poly);
141 |     
142 |     poly.clear();
143 |     p = Point(0, 1, z);
144 |     poly.append(p);
145 |     p = Point(2, 1, z);
146 |     poly.append(p);
147 |     p = Point(2, 2, z);
148 |     poly.append(p);
149 |     p = Point(0, 2, z);
150 |     poly.append(p);
151 |     poly.setName("poly_2");
152 |     region.push_back(poly);
153 |     
154 |     poly.clear();
155 |     p = Point(3, 1, z);
156 |     poly.append(p);
157 |     p = Point(5, 1, z);
158 |     poly.append(p);
159 |     p = Point(5, 2, z);
160 |     poly.append(p);
161 |     p = Point(3, 2, z);
162 |     poly.append(p);
163 |     poly.setName("poly_3");
164 |     region.push_back(poly);
165 | 
166 |     poly.clear();
167 |     p = Point(6, 1, z);
168 |     poly.append(p);
169 |     p = Point(8, 1, z);
170 |     poly.append(p);
171 |     p = Point(8, 2, z);
172 |     poly.append(p);
173 |     p = Point(6, 2, z);
174 |     poly.append(p);
175 |     poly.setName("poly_4");
176 |     region.push_back(poly);
177 | 
178 |     poly.clear();
179 |     p = Point(0, 2, z);
180 |     poly.append(p);
181 |     p = Point(2, 2, z);
182 |     poly.append(p);
183 |     p = Point(3, 2, z);
184 |     poly.append(p);
185 |     p = Point(5, 2, z);
186 |     poly.append(p);
187 |     p = Point(6, 2, z);
188 |     poly.append(p);
189 |     p = Point(8, 2, z);
190 |     poly.append(p);
191 |     p = Point(8, 3, z);
192 |     poly.append(p);
193 |     p = Point(0, 3, z);
194 |     poly.append(p);
195 |     poly.setName("poly_5");
196 |     region.push_back(poly);
197 | }
198 | 
199 | int main(int argc, char **argv)
200 | {
201 | 	clock_t start, finish;
202 | 	start = clock();
203 | 
204 | 	gmsh::initialize();
205 |     gmsh::option::setNumber("General.Terminal", 1);
206 |     gmsh::model::add("demo");
207 |     
208 |     double lc = 0.8;
209 |     vector<Polygon> region;
210 |     generatePolygons(region);
211 |     
212 |     int dim = 1;
213 |     int tag = 1;
214 |     double x, y, z;
215 |     vector<int> loop, curveloop, planeSurface;
216 |     
217 |     int ct_point = 1;
218 |     int ct_line = 1;
219 |     
220 |     // record the points where the previous polygon was inserted
221 |     map<Point, int> recPoints;
222 |     // record the line where the previous polygon wsa inserted
223 |     map<Topology, int> recLines;
224 |     vector<vector<int> > domElelIds(region.size());
225 |     for(int i = 0; i < region.size(); i++)
226 |     {
227 |         loop.clear();
228 |         Polygon poly = region[i];
229 |         const vector<Point>& vertex = poly.getVertex();
230 |         for(int j = 0; j < vertex.size(); j++)
231 |         {
232 |             Point p = vertex[j];
233 |             auto iter = recPoints.find(p);
234 |             if(iter == recPoints.end())
235 |             {
236 |                 x = p.getX();
237 |                 y = p.getY();
238 |                 z = p.getZ();
239 |                 gmsh::model::geo::addPoint(x, y, z, lc, ct_point);
240 |                 recPoints[p] = ct_point;
241 |                 loop.push_back(ct_point);
242 |                 ct_point++;
243 |             } else {
244 |                 loop.push_back(iter->second);
245 |             }
246 |         }
247 |         curveloop.clear();
248 |         for(int j = 0; j < loop.size(); j++)
249 |         {
250 |             int l1 = loop[j];
251 |             int l2 = j == loop.size()-1 ? loop[0] : loop[j+1];
252 |             Topology topo(l1, l2);
253 |             Topology topo2(l2, l1);
254 |             auto iterl = recLines.find(topo);
255 |             auto iterl2 = recLines.find(topo2);
256 |             bool flag = (iterl == recLines.end()) && (iterl2 == recLines.end());
257 |             if(flag)
258 |             {
259 |                 gmsh::model::geo::addLine(l1, l2, ct_line);
260 |                 recLines[topo] = ct_line;
261 |                 recLines[topo2] = -ct_line;
262 |                 curveloop.push_back(ct_line);
263 |                 ct_line++;
264 |             } else {
265 |                 if(iterl != recLines.end())
266 |                     curveloop.push_back(iterl->second);
267 |                 else
268 |                     curveloop.push_back(iterl2->second);
269 |             }
270 |         }
271 |         planeSurface.push_back(gmsh::model::geo::addCurveLoop(curveloop));
272 |         tag = gmsh::model::addPhysicalGroup(dim, curveloop);
273 |         string name = "Polygon_" + to_string(i);
274 |         gmsh::model::setPhysicalName(dim, tag, name);
275 |         domElelIds[i] = curveloop;
276 |     }
277 |     
278 |     gmsh::model::geo::addPlaneSurface(planeSurface);
279 |     
280 |     dim = 1;
281 |     gmsh::model::addPhysicalGroup(dim, planeSurface);
282 |     
283 |     gmsh::model::geo::synchronize();
284 |     
285 |     gmsh::model::mesh::generate(1);
286 | 
287 | #if 1
288 |     clock_t info_s, info_e;
289 |     info_s = clock();
290 |     // get Mesh information
291 |     
292 |     // get Nodes
293 |     dim = -1;
294 |     tag = -1;
295 |     bool includeBoundary = false;
296 |     vector<double> coord, parametricCoord;
297 |     vector<size_t> nodeTags;
298 |     gmsh::model::mesh::getNodes(nodeTags, coord, parametricCoord, dim, tag, includeBoundary, true);
299 |     cout<<"The number of nodes: "<<coord.size()/3<<endl;
300 |     vector<Point> nodes;
301 |     nodes.reserve(1000);
302 |     Point p;
303 |     for(int i = 0; i < coord.size(); i+=3)
304 |     {
305 |         x = coord[i];
306 |         y = coord[i+1];
307 |         z = coord[i+2];
308 |         p = Point(x, y, z);
309 |         nodes.push_back(p);
310 |     }
311 |     
312 |     // get Elemetnts
313 |     dim = 1;
314 |     tag = -1;
315 |     vector<int> elementTypes;
316 |     vector<vector<size_t> > elemetTags, nodeTagss;
317 |     gmsh::model::mesh::getElements(elementTypes, elemetTags, nodeTagss, dim, tag);
318 |     cout<<"The number of elemnets: "<<nodeTagss[0].size()/2<<endl;
319 |     Topology topo;
320 |     Element ele;
321 |     Point p1, p2;
322 |     unsigned int id1, id2;
323 |     for(int i = 0; i < nodeTagss[0].size(); i+=2)
324 |     {
325 |         id1 = nodeTagss[0][i] - 1;
326 |         id2 = nodeTagss[0][i+1] - 1;
327 |         p1 = nodes[id1];
328 |         p2 = nodes[id2];
329 |         ele = Element(p1, p2, topo);
330 |     }
331 |     
332 |     // get polygon nodes id and it's coordinate information
333 |     // coord is coordinate information (x1, y1, z1, .......)
334 |     dim = 1;
335 |     vector<size_t> bouNodesIds;
336 |     // boundaryTag = tag = gmsh::model::addPhysicalGroup(dim, curveloop)
337 |     // then we can get one polygon points id on edges
338 |     int boundaryTag = 1;
339 |     gmsh::model::mesh::getNodesForPhysicalGroup(dim, boundaryTag, bouNodesIds, coord);
340 |     vector<unsigned int> boundaryNodesIds(bouNodesIds.size());
341 |     // just because index from star 0
342 |     for(int i = 0; i < bouNodesIds.size(); i++)
343 |     {
344 |         boundaryNodesIds[i] = bouNodesIds[i]-1;
345 |         //cout<<boundaryNodesIds[i]<<endl;
346 |     }
347 |     
348 |     // get polygon arround elemnt id
349 |     // bcurveLoop is line's tag(id) of one polygon
350 |     // then we can get elemnts id of one polygon
351 |     vector<int> bcurveLoop = domElelIds[0];
352 |     vector<int> boundaryElementsIds;
353 |     dim = 1;
354 |     int dim2 = 1;
355 |     bool strict = true;
356 |     double xx, yy, zz;
357 |     vector<size_t> elemetTag;
358 |     for(int i = 0; i < bcurveLoop.size(); i++)
359 |     {
360 |         tag = abs(bcurveLoop[i]);
361 |         gmsh::model::mesh::getElements(elementTypes, elemetTags, nodeTagss, dim, tag);
362 |         for(int k = 0; k < nodeTagss[0].size(); k+=2)
363 |         {
364 |             id1 = nodeTagss[0][k]-1;
365 |             id2 = nodeTagss[0][k+1]-1;
366 |             xx = (nodes[id1].getX() + nodes[id2].getX())/2.0;
367 |             yy = (nodes[id1].getY() + nodes[id2].getY())/2.0;
368 |             zz = (nodes[id1].getZ() + nodes[id2].getZ())/2.0;
369 |             gmsh::model::mesh::getElementsByCoordinates(xx, yy, zz, elemetTag, dim2, strict);
370 |             for(int kk= 0; kk < elemetTag.size(); kk++)
371 |             {
372 |                 boundaryElementsIds.push_back(elemetTag[kk]-1);
373 |             }
374 |         }
375 |         if(bcurveLoop[i] < 0)
376 |         {
377 |             reverse(boundaryElementsIds.begin(), boundaryElementsIds.end());
378 |         }
379 |     }
380 |     info_e = clock();
381 | #endif
382 |     gmsh::write("demo.msh");
383 |     gmsh::finalize();
384 |     finish = clock();
385 |     cout<<"The run time of get information: "<<(double)(info_e - info_s)/CLOCKS_PER_SEC<<"s"<<endl;
386 |     cout<<"The run time is: "<<(double)(finish - start)/CLOCKS_PER_SEC<<"s"<<endl;
387 |     return 0;
388 | }
389 | 


--------------------------------------------------------------------------------
/demo/threeDDemo.cpp:
--------------------------------------------------------------------------------
  1 | #include <iostream>
  2 | #include <gmsh.h>
  3 | #include <vector> 
  4 | #include <fstream>
  5 | #include <iomanip>
  6 | using namespace std;
  7 | 
  8 | 
  9 | void writeNodes(vector<double>& nodes, string s)
 10 | {
 11 | 	string filename = s + ".txt";
 12 | 	ofstream outfileb(filename.c_str(), ios::out);
 13 | 	for(int i = 0; i < nodes.size(); i++)
 14 | 	{
 15 | 		outfileb<<(nodes[i])<<endl;
 16 | 	}
 17 | 	outfileb.close();
 18 | }
 19 | 
 20 | 
 21 | void writeElements(vector<std::vector<std::size_t> >& elements, string s)
 22 | {
 23 | 	string filename = s + ".txt";
 24 | 	ofstream outfileb(filename.c_str(), ios::out);
 25 | 	for(int i = 0; i < elements[1].size(); i++)
 26 | 	{
 27 | 		outfileb<<(elements[1][i])<<endl;
 28 | 	}
 29 | 	outfileb.close();
 30 | }
 31 | 
 32 | 
 33 | int main(int argc, char **argv)
 34 | {
 35 | 	// Before using any functions in the C++ API, Gmsh must be initialized.
 36 | 	gmsh::initialize();
 37 | 	// By default Gmsh will not print out any messages: in order to output
 38 | 	// messages on the terminal, just set the standard Gmsh option
 39 | 	// "General.Terminal" (same format and meaning as in .geo files) using
 40 | 	// gmsh::option::setNumber():
 41 | 	gmsh::option::setNumber("General.Terminal", 1);
 42 | 
 43 | 	// This adds a new model, named "t1". If gmsh::model::add() is not called, a
 44 | 	// new default (unnamed) model will be created on the fly, if necessary.
 45 | 	gmsh::model::add("Test1");
 46 | 	// The C++ API provides direct access to the internal CAD kernels. The
 47 |     // built-in CAD kernel was used in t1.geo: the corresponding API functions
 48 |     // live in the "gmsh::model::geo" namespace. To create geometrical points with
 49 |     // the built-in CAD kernel, one thus uses gmsh::model::geo::addPoint():
 50 |     //
 51 |     // - the first 3 arguments are the point coordinates (x, y, z)
 52 |     //
 53 |     // - the next (optional) argument is the target mesh size close to the point
 54 |     //
 55 |     // - the last (optional) argument is the point tag
 56 | 
 57 | 
 58 | 
 59 | 
 60 | 
 61 | 
 62 | 
 63 |   	double lc = 1;//2e-1;
 64 |     int pIndex = 1;
 65 |   	gmsh::model::geo::addPoint(-4, -4,  -4, lc, pIndex++);
 66 |   	gmsh::model::geo::addPoint(4,  -4,  -4, lc, pIndex++);
 67 |   	gmsh::model::geo::addPoint(4,   4,  -4, lc, pIndex++);
 68 |   	gmsh::model::geo::addPoint(-4,  4,  -4, lc, pIndex++);
 69 |   	gmsh::model::geo::addPoint(-4, -4,   4, lc,  pIndex++);
 70 |   	gmsh::model::geo::addPoint(4,  -4,   4, lc,  pIndex++);
 71 |   	gmsh::model::geo::addPoint(4,   4,   4, lc,  pIndex++);
 72 |   	gmsh::model::geo::addPoint(-4,  4,   4, lc,  pIndex++);
 73 | 
 74 |   	// The API to create lines with the built-in kernel follows the same
 75 |   	// conventions: the first 2 arguments are point tags, the last (optional one)
 76 |   	// is the line tag.
 77 |     int lIndex = 1;
 78 |     int p, q;
 79 |     p = 1;
 80 |     for(int i = 0; i < 3; i++)
 81 |     {
 82 |       gmsh::model::geo::addLine(p, p+1, lIndex++);
 83 |       p++;
 84 |     }
 85 |     gmsh::model::geo::addLine(p, 1, lIndex++);
 86 | 
 87 |     p = 5;
 88 |     for(int i = 0; i < 3; i++)
 89 |     {
 90 |       gmsh::model::geo::addLine(p, p+1, lIndex++);
 91 |       p++;
 92 |     }
 93 |     gmsh::model::geo::addLine(p, 5, lIndex++);
 94 | 
 95 |     p = 1;
 96 |     q = 5;
 97 |     for(int i = 0; i < 4; i++)
 98 |     {
 99 |        gmsh::model::geo::addLine(p, q, lIndex++);
100 |        p++;
101 |        q++;
102 |     }
103 | 
104 | 
105 |   	// The philosophy to construct curve loops and surfaces is similar: the first
106 |   	// argument is now a vector of integers.
107 |   	gmsh::model::geo::addCurveLoop({1, 2, 3, 4},     1);
108 |   	gmsh::model::geo::addCurveLoop({5, 6, 7, 8},     2);
109 |   	gmsh::model::geo::addCurveLoop({1, 10, -5, -9},  3);
110 |   	gmsh::model::geo::addCurveLoop({2, 11, -6, -10}, 4);
111 |   	gmsh::model::geo::addCurveLoop({3, 12, -7, -11}, 5);
112 |   	gmsh::model::geo::addCurveLoop({-4,12, 8,   -9}, 6);
113 | 
114 | 
115 | 
116 |   	gmsh::model::geo::addPlaneSurface({1}, 1);
117 |   	gmsh::model::geo::addPlaneSurface({2}, 2);
118 |   	gmsh::model::geo::addPlaneSurface({3}, 3);
119 |   	gmsh::model::geo::addPlaneSurface({4}, 4);
120 |   	gmsh::model::geo::addPlaneSurface({5}, 5);
121 |   	gmsh::model::geo::addPlaneSurface({6}, 6);
122 | 
123 |     gmsh::model::geo::addSurfaceLoop({1, 2, 3, 4, 5, 6}, 1);
124 |     gmsh::model::geo::addVolume({1}, 1);
125 |   	// Physical groups are defined by providing the dimension of the group (0 for
126 |   	// physical points, 1 for physical curves, 2 for physical surfaces and 3 for
127 |   	// phsyical volumes) followed by a vector of entity tags. The last (optional)
128 |   	// argument is the tag of the new group to create.
129 | 
130 | 	  gmsh::model::setPhysicalName(3, 1, "My volume");
131 | 
132 |     // Before it can be meshed, the internal CAD representation must be
133 |     // synchronized with the Gmsh model, which will create the relevant Gmsh data
134 |     // structures. This is achieved by the gmsh::model::geo::synchronize() API
135 |     // call for the built-in CAD kernel. Synchronizations can be called at any
136 |     // time, but they involve a non trivial amount of processing; so while you
137 |     // could synchronize the internal CAD data after every CAD command, it is
138 |     // usually better to minimize the number of synchronization points.
139 |     gmsh::model::geo::synchronize();
140 | 
141 |     // We can then generate a 2D mesh...
142 |     gmsh::model::mesh::generate();
143 |     //gmsh::model::mesh::refine();
144 | 
145 |     
146 |     //gmsh::model::mesh::recombine();
147 | #if 0
148 |     // get mesh information
149 |     // Nodes
150 |     vector<double> nodes,y;
151 |     vector<std::size_t> nodeTags;
152 |     gmsh::model::mesh::getNodes(nodeTags, nodes, y, -2, -1, false, true);
153 |     cout<<"The number of nodes: "<<nodes.size()/3<<endl;
154 |   #if 0
155 |     cout<<endl;
156 |     for(int i = 0; i < nodes.size(); i++)
157 |     {
158 |       cout<<nodes[i]<<endl;
159 |     }
160 |   #endif
161 |     writeNodes(nodes, "gmshNodes");
162 | 
163 |     //Elements
164 |     vector<int> elementTypes;
165 |     vector<std::vector<std::size_t> > elementTags, nodeTags2;
166 |     gmsh::model::mesh::getElements(elementTypes, elementTags, nodeTags2, -2, -1);
167 |     cout<<"The number of elements: "<<nodeTags2[1].size()/3<<endl;
168 |   #if 0
169 |     cout<<endl;
170 |     for(int j = 0; j < nodeTags2[1].size(); j++)
171 |     {
172 |       cout<<nodeTags2[1][j]<<endl;
173 | 
174 |     }
175 |   #endif
176 |     writeElements(nodeTags2, "gmshElements");
177 | 
178 | #endif
179 | 
180 |     // ... and save it to disk
181 |     gmsh::write("demo.msh");
182 | 
183 |     // Remember that by default, if physical groups are defined, Gmsh will export
184 |     // in the output mesh file only those elements that belong to at least one
185 |     // physical group. To force Gmsh to save all elements, you can use
186 |     //
187 |     // gmsh::option::setNumber("Mesh.SaveAll", 1);
188 | 
189 |     // This should be called at the end:
190 |     gmsh::finalize();
191 |     return 0;
192 | }
193 | 


--------------------------------------------------------------------------------
/demo/twoDExample.cpp:
--------------------------------------------------------------------------------
  1 | #include <gmsh.h>
  2 | #include <iostream>
  3 | #include <vector> 
  4 | #include <fstream>
  5 | #include <iomanip>
  6 | #include <string>
  7 | using namespace std;
  8 | 
  9 | 
 10 | class Point
 11 | {
 12 |   public:
 13 |     Point(const double x = 0, const double y = 0, const double z = 0):_x(x), _y(y), _z(z) {}
 14 |     ~Point() {}
 15 |     const double getX() const {return _x;}
 16 |     const double getY() const {return _y;}
 17 |     const double getZ() const {return _z;}
 18 |     void setX(const double x) {_x = x;}
 19 |     void setY(const double y) {_y = y;}
 20 |     void setZ(const double z) {_z = z;}
 21 | 
 22 |   private:
 23 |     double _x, _y, _z;
 24 | };
 25 | 
 26 | 
 27 | 
 28 | class Triangle
 29 | {
 30 |   /**
 31 |    * Triangle vertex index
 32 |    * _id[0] ~ _id[2] : three vertices index of triangle
 33 |    * _id[3] is the midpoint between _id[0] and _id[1]
 34 |    * _id[4] is the midpoint between _id[1] and _id[2]
 35 |    * _id[5] is the midpoint between _id[2] and _id[0]
 36 |    */
 37 |   public:
 38 |     Triangle() {}
 39 |     ~Triangle() {}
 40 |     const vector<unsigned int>& getIds() const {return _ids;}
 41 |     void clear() {_ids.clear();}
 42 |     void append(unsigned int id) {_ids.push_back(id);}
 43 |     void resize(const unsigned int size) {_ids.resize(size);}
 44 | 
 45 |   private:
 46 |     vector<unsigned int> _ids;
 47 | };
 48 | 
 49 | 
 50 | 
 51 | void generateBoudaryAndHoles(vector<Point>& boundary, vector<vector<Point> >& holes)
 52 | {
 53 |   double left, right, bot, top;
 54 |   double x, y, z;
 55 |   left = 0;
 56 |   right = 5;
 57 |   bot = 0;
 58 |   top = 4;
 59 |   z = 0;
 60 |   Point p;
 61 |   p = Point(left, bot, z);
 62 |   boundary.push_back(p);
 63 |   p = Point(right, bot, z);
 64 |   boundary.push_back(p);
 65 |   p = Point(right, top, z);
 66 |   boundary.push_back(p);
 67 |   p = Point(left, top, z);
 68 |   boundary.push_back(p);
 69 |   p = Point(left, (bot+top)/2.0, z);
 70 |   boundary.push_back(p);
 71 | 
 72 | 
 73 |   vector<Point> hole;
 74 |   left = 1;
 75 |   right = 2;
 76 |   bot = 1;
 77 |   top = 2;
 78 |   p = Point(left, bot, z);
 79 |   hole.push_back(p);
 80 |   p = Point(right, bot, z);
 81 |   hole.push_back(p);
 82 |   p = Point(right, top, z);
 83 |   hole.push_back(p);
 84 |   p = Point(left, top, z);
 85 |   hole.push_back(p);
 86 |   p = Point(left, (bot+top)/2.0, z);
 87 |   hole.push_back(p);
 88 |   holes.push_back(hole);
 89 | 
 90 | 
 91 |   hole.clear();
 92 |   left = 3;
 93 |   right = 4;
 94 |   bot = 1;
 95 |   top = 2;
 96 |   p = Point(left, bot, z);
 97 |   hole.push_back(p);
 98 |   p = Point(right, bot, z);
 99 |   hole.push_back(p);
100 |   p = Point(right, top, z);
101 |   hole.push_back(p);
102 |   p = Point(left, top, z);
103 |   hole.push_back(p);
104 |   holes.push_back(hole);
105 | }
106 | 
107 | 
108 | int main(int argc, char **argv)
109 | {
110 |   gmsh::initialize();
111 |   gmsh::option::setNumber("General.Terminal", 1);
112 |     
113 |   gmsh::model::add("demo");
114 | 
115 |   double boundartLc = 1, holeLc = 0.4;
116 |   vector<Point> boundary;
117 |   vector<vector<Point> > holes;
118 |   generateBoudaryAndHoles(boundary, holes);
119 | 
120 |   int dim = 1;
121 |   int tag = 1;
122 |   double x, y, z;
123 |   vector<int> loop, curveLoop, planeSurface;
124 | 
125 |   int boundaryTag;
126 |   for(int i = 0; i < boundary.size(); i++)
127 |   {
128 |     x = boundary[i].getX();
129 |     y = boundary[i].getY();
130 |     z = boundary[i].getZ();
131 |     loop.push_back(gmsh::model::geo::addPoint(x, y, z, boundartLc));
132 |   }
133 | 
134 |   for(int i = 0; i < loop.size(); i++)
135 |   {
136 |     curveLoop.push_back(gmsh::model::geo::addLine(loop[i], i == loop.size()-1 ? loop[0] : loop[i+1]));
137 |   }
138 |   planeSurface.push_back(gmsh::model::geo::addCurveLoop(curveLoop));
139 |   boundaryTag = gmsh::model::addPhysicalGroup(dim, curveLoop);
140 |   gmsh::model::setPhysicalName(dim, boundaryTag, "Boundary");
141 | 
142 |   vector<int> bcurveLoop = curveLoop;
143 | 
144 |   vector<vector<int> > recCurveLoop;
145 |   vector<int> holesTags;
146 |   for(int i = 0; i < holes.size(); i++)
147 |   {
148 |     loop.clear();
149 |     curveLoop.clear();
150 |     for(int j = 0; j < holes[i].size(); j++)
151 |     {
152 |       x = holes[i][j].getX();
153 |       y = holes[i][j].getY();
154 |       z = holes[i][j].getZ();
155 |       loop.push_back(gmsh::model::geo::addPoint(x, y, z, holeLc));
156 |     }
157 | 
158 |     for(int j = 0; j < loop.size(); j++)
159 |     {
160 |       curveLoop.push_back(gmsh::model::geo::addLine(loop[j], j == loop.size()-1 ? loop[0] : loop[j+1]));
161 |     }
162 |     planeSurface.push_back(-gmsh::model::geo::addCurveLoop(curveLoop));
163 |     
164 |     tag = gmsh::model::addPhysicalGroup(dim, curveLoop);
165 |     holesTags.push_back(tag);
166 |     string name = "hole" + to_string(i);
167 |     gmsh::model::setPhysicalName(dim, tag, name);
168 |     recCurveLoop.push_back(curveLoop);
169 |   }
170 | 
171 |   gmsh::model::geo::addPlaneSurface(planeSurface);
172 | 
173 |   dim = 2;
174 |   gmsh::model::addPhysicalGroup(dim, planeSurface);
175 | 
176 |   gmsh::model::geo::synchronize();
177 | 
178 |   gmsh::model::mesh::generate(2);
179 | 
180 | 
181 |   // get mesh information
182 |   // Nodes
183 |   dim = -1;
184 |   tag = -1;
185 |   bool includeBoundary = false;
186 |   vector<double> coord, parametricCoord;
187 |   vector<size_t> nodeTags;
188 |   gmsh::model::mesh::getNodes(nodeTags, coord, parametricCoord, dim, tag, includeBoundary, true);
189 |   cout<<"The number of nodes: "<<coord.size()/3<<endl;
190 |   vector<Point> nodes;
191 |   nodes.reserve(1000);
192 |   Point p;
193 |   for(int i = 0; i < coord.size(); i+=3)
194 |   {
195 |     x = coord[i];
196 |     y = coord[i+1];
197 |     z = coord[i+2];
198 |     p = Point(x, y, z);
199 |     nodes.push_back(p);
200 |   }
201 | 
202 | 
203 | 
204 |   //Elements
205 |   dim = -1;
206 |   tag = -1;
207 |   vector<int> elementTypes;
208 |   vector<vector<size_t> > elementTags, nodeTagss;
209 |   gmsh::model::mesh::getElements(elementTypes, elementTags, nodeTagss, dim, tag);
210 |   cout<<"The number of elements: "<<nodeTagss[1].size()/3<<endl;
211 |   // boundaryNodesId include holes boundary nodes id
212 |   // vector<size_t> boundaryNodesId = elementTags[0];
213 |   vector<Triangle> triangles;
214 |   triangles.reserve(1000);
215 |   unsigned int id1, id2, id3;
216 |   for(int i = 0; i < nodeTagss[1].size(); i += 3)
217 |   {
218 |     Triangle tri;
219 |     id1 = nodeTagss[1][i]-1;
220 |     id2 = nodeTagss[1][i+1]-1;
221 |     id3 = nodeTagss[1][i+2]-1;
222 |     tri.append(id1);
223 |     tri.append(id2);
224 |     tri.append(id3);
225 |     triangles.push_back(tri);
226 |   }
227 |   
228 | 
229 | 
230 |   // get boundary nodes id and it's coordinate information
231 |   // coord is coordinate information (x1, y1, z1, .....)
232 |   dim = 1;
233 |   vector<size_t> bouNodesIds;
234 |   gmsh::model::mesh::getNodesForPhysicalGroup(dim, boundaryTag, bouNodesIds, coord);
235 |   vector<unsigned int> boundaryNodesIds(bouNodesIds.size());
236 |   // just because index from start 0
237 |   for(int i = 0; i < bouNodesIds.size(); i++)
238 |   {
239 |     boundaryNodesIds[i] = bouNodesIds[i]-1;
240 |   }
241 | 
242 | 
243 |   // get holes nodes id and it's coordinate information
244 |   vector<vector<size_t> > hosNodesIds(holesTags.size());
245 |   vector<vector<unsigned int> > holesNodesIds(holesTags.size());
246 |   for(int i = 0; i < holesTags.size(); i++)
247 |   {
248 |     gmsh::model::mesh::getNodesForPhysicalGroup(dim, holesTags[i], hosNodesIds[i], coord);
249 |   }
250 |   for(int i = 0; i < hosNodesIds.size(); i++)
251 |   {
252 |     holesNodesIds[i].resize(hosNodesIds[i].size());
253 |     for(int j = 0; j < hosNodesIds[i].size(); j++)
254 |     {
255 |       holesNodesIds[i][j] = hosNodesIds[i][j]-1;
256 |     }
257 |   }
258 | 
259 | #if 0
260 |   for(int i = 0; i < holesNodesIds.size(); i++)
261 |   {
262 |     for(int j = 0; j < holesNodesIds[i].size(); j++)
263 |     {
264 |       cout<<holesNodesIds[i][j]<<" ";
265 |     }
266 |     cout<<endl;
267 |   }
268 | #endif
269 | 
270 | 
271 | 
272 |   // get holes arround triangle id
273 |   vector<vector<int> > holesElementsIds(recCurveLoop.size());
274 |   dim = 1;
275 |   int dim2 = 2;
276 |   vector<size_t> elementTag;
277 |   bool strict = true;
278 |   double xx, yy, zz;
279 |   vector<int> holeElementsIds;
280 |   for(int i = 0; i < recCurveLoop.size(); i++)
281 |   {
282 |     holeElementsIds.clear();
283 |     for(int j = 0; j < recCurveLoop[i].size(); j++)
284 |     {
285 |       //elementTypes.clear();
286 |       //elementTags.clear();
287 |       //nodeTagss.clear();
288 |       tag = recCurveLoop[i][j];
289 |       gmsh::model::mesh::getElements(elementTypes, elementTags, nodeTagss, dim, tag);
290 |       for(int k = 0; k < nodeTagss[0].size(); k += 2)
291 |       {
292 |         id1 = nodeTagss[0][k]-1;
293 |         id2 = nodeTagss[0][k+1]-1;
294 |         xx = (nodes[id1].getX() + nodes[id2].getX())/2.0;
295 |         yy = (nodes[id1].getY() + nodes[id2].getY())/2.0;
296 |         zz = (nodes[id1].getZ() + nodes[id2].getZ())/2.0;
297 |         gmsh::model::mesh::getElementsByCoordinates(xx, yy, zz, elementTag, dim2, strict);
298 |         for(int kk = 0; kk < elementTag.size(); kk++)
299 |         {
300 |           holeElementsIds.push_back(elementTag[kk]-1);
301 |         }
302 |       }
303 |     }
304 |     holesElementsIds.push_back(holeElementsIds);
305 |   }
306 | #if 0
307 |   for(int i = 0; i < holesElementsIds.size(); i++)
308 |   {
309 |     for(int j = 0; j < holesElementsIds[i].size(); j++)
310 |     {
311 |       cout<<holesElementsIds[i][j]<<" ";
312 |     }
313 |     cout<<endl;
314 |   }
315 | #endif
316 | 
317 |   // get boundary arround triangle id
318 |   vector<int> boundaryElementsIds;
319 |   dim = 1;
320 |   dim2 = 2;
321 |   strict = true;
322 |   for(int i = 0; i < bcurveLoop.size(); i++)
323 |   {
324 |     tag = bcurveLoop [i];
325 |     gmsh::model::mesh::getElements(elementTypes, elementTags, nodeTagss, dim, tag);
326 |     for(int k = 0; k < nodeTagss[0].size(); k += 2)
327 |     {
328 |       id1 = nodeTagss[0][k]-1;
329 |       id2 = nodeTagss[0][k+1]-1;
330 |       xx = (nodes[id1].getX() + nodes[id2].getX())/2.0;
331 |       yy = (nodes[id1].getY() + nodes[id2].getY())/2.0;
332 |       zz = (nodes[id1].getZ() + nodes[id2].getZ())/2.0;
333 |       gmsh::model::mesh::getElementsByCoordinates(xx, yy, zz, elementTag, dim2, strict);
334 |       for(int kk = 0; kk < elementTag.size(); kk++)
335 |       {
336 |         boundaryElementsIds.push_back(elementTag[kk]-1);
337 |       }
338 |     }
339 |   }
340 | 
341 |   gmsh::write("demo.msh");
342 |     
343 |   gmsh::finalize();
344 |   return 0;
345 | }
346 | 


--------------------------------------------------------------------------------