12#ifndef MFEM_NCMESH_TABLES
13#define MFEM_NCMESH_TABLES
18static constexpr int ref_type_num_children[8] = { 0, 2, 2, 4, 2, 4, 4, 8 };
26static constexpr int quad_deref_table[3][4 + 4] =
28 { 0, 1, 1, 0, 1, 1, 0, 0 },
29 { 0, 0, 1, 1, 0, 0, 1, 1 },
30 { 0, 1, 2, 3, 1, 1, 3, 3 }
33static constexpr int hex_deref_table[7][8 + 6] =
35 { 0, 1, 1, 0, 0, 1, 1, 0, 1, 1, 1, 0, 0, 0 },
36 { 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 0, 1, 1, 1 },
37 { 0, 1, 2, 3, 0, 1, 2, 3, 1, 1, 1, 3, 3, 3 },
38 { 0, 0, 0, 0, 1, 1, 1, 1, 0, 0, 0, 1, 1, 1 },
39 { 0, 1, 1, 0, 3, 2, 2, 3, 1, 1, 1, 3, 3, 3 },
40 { 0, 0, 1, 1, 2, 2, 3, 3, 0, 0, 0, 3, 3, 3 },
41 { 0, 1, 2, 3, 4, 5, 6, 7, 1, 1, 1, 7, 7, 7 }
44static constexpr int prism_deref_table[7][6 + 5] =
46 {-1,-1,-1,-1,-1,-1, -1,-1,-1,-1,-1 },
47 {-1,-1,-1,-1,-1,-1, -1,-1,-1,-1,-1 },
48 { 0, 1, 2, 0, 1, 2, 0, 0, 0, 1, 0 },
49 { 0, 0, 0, 1, 1, 1, 0, 1, 0, 0, 0 },
50 {-1,-1,-1,-1,-1,-1, -1,-1,-1,-1,-1 },
51 {-1,-1,-1,-1,-1,-1, -1,-1,-1,-1,-1 },
52 { 0, 1, 2, 4, 5, 6, 0, 5, 0, 5, 0 }
55static constexpr int pyramid_deref_table[7][5 + 5] =
57 {-1,-1,-1,-1,-1, -1,-1,-1,-1,-1 },
58 {-1,-1,-1,-1,-1, -1,-1,-1,-1,-1 },
59 {-1,-1,-1,-1,-1, -1,-1,-1,-1,-1 },
60 {-1,-1,-1,-1,-1, -1,-1,-1,-1,-1 },
61 {-1,-1,-1,-1,-1, -1,-1,-1,-1,-1 },
62 {-1,-1,-1,-1,-1, -1,-1,-1,-1,-1 },
63 { 0, 1, 2, 3, 5, 0, 5, 5, 5, 5 }
68static constexpr char quad_hilbert_child_order[8][4] =
70 {0,1,2,3}, {0,3,2,1}, {1,2,3,0}, {1,0,3,2},
71 {2,3,0,1}, {2,1,0,3}, {3,0,1,2}, {3,2,1,0}
74static constexpr char quad_hilbert_child_state[8][4] =
76 {1,0,0,5}, {0,1,1,4}, {3,2,2,7}, {2,3,3,6},
77 {5,4,4,1}, {4,5,5,0}, {7,6,6,3}, {6,7,7,2}
80static constexpr char hex_hilbert_child_order[24][8] =
82 {0,1,2,3,7,6,5,4}, {0,3,7,4,5,6,2,1}, {0,4,5,1,2,6,7,3},
83 {1,0,3,2,6,7,4,5}, {1,2,6,5,4,7,3,0}, {1,5,4,0,3,7,6,2},
84 {2,1,5,6,7,4,0,3}, {2,3,0,1,5,4,7,6}, {2,6,7,3,0,4,5,1},
85 {3,0,4,7,6,5,1,2}, {3,2,1,0,4,5,6,7}, {3,7,6,2,1,5,4,0},
86 {4,0,1,5,6,2,3,7}, {4,5,6,7,3,2,1,0}, {4,7,3,0,1,2,6,5},
87 {5,1,0,4,7,3,2,6}, {5,4,7,6,2,3,0,1}, {5,6,2,1,0,3,7,4},
88 {6,2,3,7,4,0,1,5}, {6,5,1,2,3,0,4,7}, {6,7,4,5,1,0,3,2},
89 {7,3,2,6,5,1,0,4}, {7,4,0,3,2,1,5,6}, {7,6,5,4,0,1,2,3}
92static constexpr char hex_hilbert_child_state[24][8] =
94 {1,2,2,7,7,21,21,17}, {2,0,0,22,22,16,16,8}, {0,1,1,15,15,6,6,23},
95 {4,5,5,10,10,18,18,14}, {5,3,3,19,19,13,13,11}, {3,4,4,12,12,9,9,20},
96 {8,7,7,17,17,23,23,2}, {6,8,8,0,0,15,15,22}, {7,6,6,21,21,1,1,16},
97 {11,10,10,14,14,20,20,5}, {9,11,11,3,3,12,12,19}, {10,9,9,18,18,4,4,13},
98 {13,14,14,5,5,19,19,10}, {14,12,12,20,20,11,11,4}, {12,13,13,9,9,3,3,18},
99 {16,17,17,2,2,22,22,7}, {17,15,15,23,23,8,8,1}, {15,16,16,6,6,0,0,21},
100 {20,19,19,11,11,14,14,3}, {18,20,20,4,4,10,10,12}, {19,18,18,13,13,5,5,9},
101 {23,22,22,8,8,17,17,0}, {21,23,23,1,1,7,7,15}, {22,21,21,16,16,2,2,6}
109static constexpr RefCoord T_HALF = (1ll << 59);
110static constexpr RefCoord T_ONE = (1ll << 60);
111static constexpr RefCoord T_TWO = (1ll << 61);
114static constexpr RefCoord S_HALF = 1;
118static constexpr RefCoord tri_corners[3][3] =
125static constexpr RefCoord quad_corners[4][3] =
133static constexpr RefCoord hex_corners[8][3] =
141 {T_ONE, T_ONE, T_ONE},
145static constexpr RefCoord prism_corners[6][3] =
155static constexpr RefCoord pyramid_corners[5][3] =
165static const RefPoint* geom_corners[8] =
184 for (
int i = 0; i < 3; i++)
186 dst[i] = (src[i]*
s[i] >> 1) +
t[i];
191static constexpr RefTrf quad_parent_rt1[2] =
193 { {S_HALF, S_ONE, 0}, { 0, 0, 0} },
194 { {S_HALF, S_ONE, 0}, {T_HALF, 0, 0} }
197static constexpr RefTrf quad_child_rt1[2] =
199 { {S_TWO, S_ONE, 0}, { 0, 0, 0} },
200 { {S_TWO, S_ONE, 0}, {-T_ONE, 0, 0} }
203static constexpr RefTrf quad_parent_rt2[2] =
205 { {S_ONE, S_HALF, 0}, {0, 0, 0} },
206 { {S_ONE, S_HALF, 0}, {0, T_HALF, 0} }
209static constexpr RefTrf quad_child_rt2[2] =
211 { {S_ONE, S_TWO, 0}, {0, 0, 0} },
212 { {S_ONE, S_TWO, 0}, {0, -T_ONE, 0} }
215static constexpr RefTrf quad_parent_rt3[4] =
217 { {S_HALF, S_HALF, 0}, { 0, 0, 0} },
218 { {S_HALF, S_HALF, 0}, {T_HALF, 0, 0} },
219 { {S_HALF, S_HALF, 0}, {T_HALF, T_HALF, 0} },
220 { {S_HALF, S_HALF, 0}, { 0, T_HALF, 0} }
223static constexpr RefTrf quad_child_rt3[4] =
225 { {S_TWO, S_TWO, 0}, { 0, 0, 0} },
226 { {S_TWO, S_TWO, 0}, {-T_ONE, 0, 0} },
227 { {S_TWO, S_TWO, 0}, {-T_ONE, -T_ONE, 0} },
228 { {S_TWO, S_TWO, 0}, { 0, -T_ONE, 0} }
231static const RefTrf* quad_parent[4] =
239static const RefTrf* quad_child[4] =
247static constexpr RefTrf hex_parent_rt1[2] =
249 { {S_HALF, S_ONE, S_ONE}, { 0, 0, 0} },
250 { {S_HALF, S_ONE, S_ONE}, {T_HALF, 0, 0} }
253static constexpr RefTrf hex_child_rt1[2] =
255 { {S_TWO, S_ONE, S_ONE}, { 0, 0, 0} },
256 { {S_TWO, S_ONE, S_ONE}, {-T_ONE, 0, 0} }
259static constexpr RefTrf hex_parent_rt2[2] =
261 { {S_ONE, S_HALF, S_ONE}, {0, 0, 0} },
262 { {S_ONE, S_HALF, S_ONE}, {0, T_HALF, 0} }
265static constexpr RefTrf hex_child_rt2[2] =
267 { {S_ONE, S_TWO, S_ONE}, {0, 0, 0} },
268 { {S_ONE, S_TWO, S_ONE}, {0, -T_ONE, 0} }
271static constexpr RefTrf hex_parent_rt3[4] =
273 { {S_HALF, S_HALF, S_ONE}, { 0, 0, 0} },
274 { {S_HALF, S_HALF, S_ONE}, {T_HALF, 0, 0} },
275 { {S_HALF, S_HALF, S_ONE}, {T_HALF, T_HALF, 0} },
276 { {S_HALF, S_HALF, S_ONE}, { 0, T_HALF, 0} }
279static constexpr RefTrf hex_child_rt3[4] =
281 { {S_TWO, S_TWO, S_ONE}, { 0, 0, 0} },
282 { {S_TWO, S_TWO, S_ONE}, {-T_ONE, 0, 0} },
283 { {S_TWO, S_TWO, S_ONE}, {-T_ONE, -T_ONE, 0} },
284 { {S_TWO, S_TWO, S_ONE}, { 0, -T_ONE, 0} }
287static constexpr RefTrf hex_parent_rt4[2] =
289 { {S_ONE, S_ONE, S_HALF}, {0, 0, 0} },
290 { {S_ONE, S_ONE, S_HALF}, {0, 0, T_HALF} }
293static constexpr RefTrf hex_child_rt4[2] =
295 { {S_ONE, S_ONE, S_TWO}, {0, 0, 0} },
296 { {S_ONE, S_ONE, S_TWO}, {0, 0, -T_ONE} }
299static constexpr RefTrf hex_parent_rt5[4] =
301 { {S_HALF, S_ONE, S_HALF}, { 0, 0, 0} },
302 { {S_HALF, S_ONE, S_HALF}, {T_HALF, 0, 0} },
303 { {S_HALF, S_ONE, S_HALF}, {T_HALF, 0, T_HALF} },
304 { {S_HALF, S_ONE, S_HALF}, { 0, 0, T_HALF} }
307static constexpr RefTrf hex_child_rt5[4] =
309 { {S_TWO, S_ONE, S_TWO}, { 0, 0, 0} },
310 { {S_TWO, S_ONE, S_TWO}, {-T_ONE, 0, 0} },
311 { {S_TWO, S_ONE, S_TWO}, {-T_ONE, 0, -T_ONE} },
312 { {S_TWO, S_ONE, S_TWO}, { 0, 0, -T_ONE} }
315static constexpr RefTrf hex_parent_rt6[4] =
317 { {S_ONE, S_HALF, S_HALF}, {0, 0, 0} },
318 { {S_ONE, S_HALF, S_HALF}, {0, T_HALF, 0} },
319 { {S_ONE, S_HALF, S_HALF}, {0, 0, T_HALF} },
320 { {S_ONE, S_HALF, S_HALF}, {0, T_HALF, T_HALF} }
323static constexpr RefTrf hex_child_rt6[4] =
325 { {S_ONE, S_TWO, S_TWO}, {0, 0, 0} },
326 { {S_ONE, S_TWO, S_TWO}, {0, -T_ONE, 0} },
327 { {S_ONE, S_TWO, S_TWO}, {0, 0, -T_ONE} },
328 { {S_ONE, S_TWO, S_TWO}, {0, -T_ONE, -T_ONE} }
331static constexpr RefTrf hex_parent_rt7[8] =
333 { {S_HALF, S_HALF, S_HALF}, { 0, 0, 0} },
334 { {S_HALF, S_HALF, S_HALF}, {T_HALF, 0, 0} },
335 { {S_HALF, S_HALF, S_HALF}, {T_HALF, T_HALF, 0} },
336 { {S_HALF, S_HALF, S_HALF}, { 0, T_HALF, 0} },
337 { {S_HALF, S_HALF, S_HALF}, { 0, 0, T_HALF} },
338 { {S_HALF, S_HALF, S_HALF}, {T_HALF, 0, T_HALF} },
339 { {S_HALF, S_HALF, S_HALF}, {T_HALF, T_HALF, T_HALF} },
340 { {S_HALF, S_HALF, S_HALF}, { 0, T_HALF, T_HALF} }
343static constexpr RefTrf hex_child_rt7[8] =
345 { {S_TWO, S_TWO, S_TWO}, { 0, 0, 0} },
346 { {S_TWO, S_TWO, S_TWO}, {-T_ONE, 0, 0} },
347 { {S_TWO, S_TWO, S_TWO}, {-T_ONE, -T_ONE, 0} },
348 { {S_TWO, S_TWO, S_TWO}, { 0, -T_ONE, 0} },
349 { {S_TWO, S_TWO, S_TWO}, { 0, 0, -T_ONE} },
350 { {S_TWO, S_TWO, S_TWO}, {-T_ONE, 0, -T_ONE} },
351 { {S_TWO, S_TWO, S_TWO}, {-T_ONE, -T_ONE, -T_ONE} },
352 { {S_TWO, S_TWO, S_TWO}, { 0, -T_ONE, -T_ONE} }
355static const RefTrf* hex_parent[8] =
367static const RefTrf* hex_child[8] =
379static constexpr RefTrf tri_parent_rt3[4] =
381 { { S_HALF, S_HALF, 0}, { 0, 0, 0} },
382 { { S_HALF, S_HALF, 0}, {T_HALF, 0, 0} },
383 { { S_HALF, S_HALF, 0}, { 0, T_HALF, 0} },
384 { {-S_HALF, -S_HALF, 0}, {T_HALF, T_HALF, 0} }
387static constexpr RefTrf tri_child_rt3[4] =
389 { { S_TWO, S_TWO, 0}, { 0, 0, 0} },
390 { { S_TWO, S_TWO, 0}, {-T_ONE, 0, 0} },
391 { { S_TWO, S_TWO, 0}, { 0, -T_ONE, 0} },
392 { {-S_TWO, -S_TWO, 0}, { T_ONE, T_ONE, 0} }
395static const RefTrf* tri_parent[4] =
401static const RefTrf* tri_child[4] =
407static constexpr RefTrf prism_parent_rt3[4] =
409 { { S_HALF, S_HALF, S_ONE}, { 0, 0, 0} },
410 { { S_HALF, S_HALF, S_ONE}, {T_HALF, 0, 0} },
411 { { S_HALF, S_HALF, S_ONE}, { 0, T_HALF, 0} },
412 { {-S_HALF, -S_HALF, S_ONE}, {T_HALF, T_HALF, 0} }
415static constexpr RefTrf prism_child_rt3[4] =
417 { { S_TWO, S_TWO, S_ONE}, { 0, 0, 0} },
418 { { S_TWO, S_TWO, S_ONE}, {-T_ONE, 0, 0} },
419 { { S_TWO, S_TWO, S_ONE}, { 0, -T_ONE, 0} },
420 { {-S_TWO, -S_TWO, S_ONE}, { T_ONE, T_ONE, 0} }
423static constexpr RefTrf prism_parent_rt4[2] =
425 { {S_ONE, S_ONE, S_HALF}, {0, 0, 0} },
426 { {S_ONE, S_ONE, S_HALF}, {0, 0, T_HALF} }
429static constexpr RefTrf prism_child_rt4[2] =
431 { {S_ONE, S_ONE, S_TWO}, {0, 0, 0} },
432 { {S_ONE, S_ONE, S_TWO}, {0, 0, -T_ONE} }
435static constexpr RefTrf prism_parent_rt7[8] =
437 { { S_HALF, S_HALF, S_HALF}, { 0, 0, 0} },
438 { { S_HALF, S_HALF, S_HALF}, {T_HALF, 0, 0} },
439 { { S_HALF, S_HALF, S_HALF}, { 0, T_HALF, 0} },
440 { {-S_HALF, -S_HALF, S_HALF}, {T_HALF, T_HALF, 0} },
441 { { S_HALF, S_HALF, S_HALF}, { 0, 0, T_HALF} },
442 { { S_HALF, S_HALF, S_HALF}, {T_HALF, 0, T_HALF} },
443 { { S_HALF, S_HALF, S_HALF}, { 0, T_HALF, T_HALF} },
444 { {-S_HALF, -S_HALF, S_HALF}, {T_HALF, T_HALF, T_HALF} }
447static constexpr RefTrf prism_child_rt7[8] =
449 { { S_TWO, S_TWO, S_TWO}, { 0, 0, 0} },
450 { { S_TWO, S_TWO, S_TWO}, {-T_ONE, 0, 0} },
451 { { S_TWO, S_TWO, S_TWO}, { 0, -T_ONE, 0} },
452 { {-S_TWO, -S_TWO, S_TWO}, { T_ONE, T_ONE, 0} },
453 { { S_TWO, S_TWO, S_TWO}, { 0, 0, -T_ONE} },
454 { { S_TWO, S_TWO, S_TWO}, {-T_ONE, 0, -T_ONE} },
455 { { S_TWO, S_TWO, S_TWO}, { 0, -T_ONE, -T_ONE} },
456 { {-S_TWO, -S_TWO, S_TWO}, { T_ONE, T_ONE, -T_ONE} }
459static const RefTrf* prism_parent[8] =
468static const RefTrf* prism_child[8] =
477static const RefTrf** geom_parent[7] =
488static const RefTrf** geom_child[7] =
NCMesh::RefCoord RefCoord
void Apply(const RefCoord src[3], RefCoord dst[3]) const