32 #define LAMBDA(E, NU) (E * NU / ((1. + NU) * (1. - 2. * NU)))
33 #define MU(E, NU) (0.5 * E / (1. + NU))
34 #define DELTA(NU_P, NU_PZ, E_P, E_Z) \
35 (((1 + NU_P) * (1 - NU_P - 2 * NU_PZ * (NU_PZ * E_Z / E_P))) / \
38 #define N_MBTET0(x, y, z) (1. - x - y - z)
39 #define N_MBTET1(x, y, z) (x)
40 #define N_MBTET2(x, y, z) (y)
41 #define N_MBTET3(x, y, z) (z)
42 #define diffN_MBTET0x (-1.)
43 #define diffN_MBTET0y (-1.)
44 #define diffN_MBTET0z (-1.)
45 #define diffN_MBTET1x (1.)
46 #define diffN_MBTET1y (0.)
47 #define diffN_MBTET1z (0.)
48 #define diffN_MBTET2x (0.)
49 #define diffN_MBTET2y (1.)
50 #define diffN_MBTET2z (0.)
51 #define diffN_MBTET3x (0.)
52 #define diffN_MBTET3y (0.)
53 #define diffN_MBTET3z (1.)
56 #define N_MBTRI0(x, y) (1. - x - y)
57 #define N_MBTRI1(x, y) (x)
58 #define N_MBTRI2(x, y) (y)
59 #define diffN_MBTRI0x (-1.)
60 #define diffN_MBTRI0y (-1.)
61 #define diffN_MBTRI1x (1.)
62 #define diffN_MBTRI1y (0.)
63 #define diffN_MBTRI2x (0.)
64 #define diffN_MBTRI2y (1.)
67 #define N_MBQUAD0(x, y) ((1. - x) * (1. - y))
68 #define N_MBQUAD1(x, y) ((x) * (1. - y))
69 #define N_MBQUAD2(x, y) ((x) * (y))
70 #define N_MBQUAD3(x, y) ((1. - x) * (y))
71 #define diffN_MBQUAD0x(y) (-(1. - y))
72 #define diffN_MBQUAD0y(x) (-(1. - x))
73 #define diffN_MBQUAD1x(y) ((1. - y))
74 #define diffN_MBQUAD1y(x) (-x)
75 #define diffN_MBQUAD2x(y) (y)
76 #define diffN_MBQUAD2y(x) (x)
77 #define diffN_MBQUAD3x(y) (-y)
78 #define diffN_MBQUAD3y(x) ((1. - x))
81 #define N_MBHEX0(x, y, z) (N_MBQUAD0(x, y) * (1 - z))
82 #define N_MBHEX1(x, y, z) (N_MBQUAD1(x, y) * (1 - z))
83 #define N_MBHEX2(x, y, z) (N_MBQUAD2(x, y) * (1 - z))
84 #define N_MBHEX3(x, y, z) (N_MBQUAD3(x, y) * (1 - z))
85 #define N_MBHEX4(x, y, z) (N_MBQUAD0(x, y) * (z))
86 #define N_MBHEX5(x, y, z) (N_MBQUAD1(x, y) * (z))
87 #define N_MBHEX6(x, y, z) (N_MBQUAD2(x, y) * (z))
88 #define N_MBHEX7(x, y, z) (N_MBQUAD3(x, y) * (z))
89 #define diffN_MBHEX0x(y, z) (diffN_MBQUAD0x(y) * (1 - z))
90 #define diffN_MBHEX1x(y, z) (diffN_MBQUAD1x(y) * (1 - z))
91 #define diffN_MBHEX2x(y, z) (diffN_MBQUAD2x(y) * (1 - z))
92 #define diffN_MBHEX3x(y, z) (diffN_MBQUAD3x(y) * (1 - z))
93 #define diffN_MBHEX4x(y, z) (diffN_MBQUAD0x(y) * (z))
94 #define diffN_MBHEX5x(y, z) (diffN_MBQUAD1x(y) * (z))
95 #define diffN_MBHEX6x(y, z) (diffN_MBQUAD2x(y) * (z))
96 #define diffN_MBHEX7x(y, z) (diffN_MBQUAD3x(y) * (z))
97 #define diffN_MBHEX0y(x, z) (diffN_MBQUAD0y(x) * (1 - z))
98 #define diffN_MBHEX1y(x, z) (diffN_MBQUAD1y(x) * (1 - z))
99 #define diffN_MBHEX2y(x, z) (diffN_MBQUAD2y(x) * (1 - z))
100 #define diffN_MBHEX3y(x, z) (diffN_MBQUAD3y(x) * (1 - z))
101 #define diffN_MBHEX4y(x, z) (diffN_MBQUAD0y(x) * (z))
102 #define diffN_MBHEX5y(x, z) (diffN_MBQUAD1y(x) * (z))
103 #define diffN_MBHEX6y(x, z) (diffN_MBQUAD2y(x) * (z))
104 #define diffN_MBHEX7y(x, z) (diffN_MBQUAD3y(x) * (z))
105 #define diffN_MBHEX0z(x, y) (-N_MBQUAD0(x, y))
106 #define diffN_MBHEX1z(x, y) (-N_MBQUAD1(x, y))
107 #define diffN_MBHEX2z(x, y) (-N_MBQUAD2(x, y))
108 #define diffN_MBHEX3z(x, y) (-N_MBQUAD3(x, y))
109 #define diffN_MBHEX4z(x, y) (N_MBQUAD0(x, y))
110 #define diffN_MBHEX5z(x, y) (N_MBQUAD1(x, y))
111 #define diffN_MBHEX6z(x, y) (N_MBQUAD2(x, y))
112 #define diffN_MBHEX7z(x, y) (N_MBQUAD3(x, y))
115 #define N_MBEDGE0(x) (1. - (x))
116 #define N_MBEDGE1(x) (x)
117 #define diffN_MBEDGE0 (-1.)
118 #define diffN_MBEDGE1 (1.)
121 #define N_MBTRIQ0(x, y) ((1. - x - y) * (2 * (1. - x - y) - 1.))
122 #define N_MBTRIQ1(x, y) (x * (2. * x - 1.))
123 #define N_MBTRIQ2(x, y) (y * (2. * y - 1.))
124 #define N_MBTRIQ3(x, y) (4. * (1. - x - y) * x)
125 #define N_MBTRIQ4(x, y) (4. * x * y)
126 #define N_MBTRIQ5(x, y) (4. * (1. - x - y) * y)
127 #define diffN_MBTRIQ0x(x, y) (x + y - 3. * (1. - x - y))
128 #define diffN_MBTRIQ0y(x, y) (x + y - 3. * (1. - x - y))
129 #define diffN_MBTRIQ1x(x, y) (-1. + 4. * x)
130 #define diffN_MBTRIQ1y(x, y) (0.)
131 #define diffN_MBTRIQ2x(x, y) (0.)
132 #define diffN_MBTRIQ2y(x, y) (-1. + 4. * y)
133 #define diffN_MBTRIQ3x(x, y) (4. * ((1. - x - y) - x))
134 #define diffN_MBTRIQ3y(x, y) (-4. * x)
135 #define diffN_MBTRIQ4x(x, y) (4. * y)
136 #define diffN_MBTRIQ4y(x, y) (4. * x)
137 #define diffN_MBTRIQ5x(x, y) (-4. * y)
138 #define diffN_MBTRIQ5y(x, y) (4. * ((1. - x - y) - y))
156 PetscErrorCode
ShapeMBTRI(
double *
N,
const double *X,
const double *Y,
168 double *
normal,
double *s1,
double *s2);
172 double *diff_normal);
178 PetscErrorCode
ShapeMBTET(
double *
N,
const double *G_X,
const double *G_Y,
179 const double *G_Z,
int DIM);
185 PetscErrorCode
ShapeJacMBTET(
double *diffN,
const double *coords,
double *jac);
192 double *diffNinvJac);
197 PetscErrorCode
Spin(
double *spinOmega,
double *vecOmega);
205 int *order_edge,
double *diffN,
double *diffN_face,
206 double *diffN_edge[],
double *dofs,
double *dofs_edge[],
207 double *dofs_face,
double *idofs,
double *idofs_edge[],
219 const double *elem_coords,
220 const double *glob_coords,
250 const double *elem_coords,
251 const double *glob_coords,
258 PetscErrorCode
ShapeMBEDGE(
double *
N,
const double *G_X,
int DIM);
262 PetscErrorCode
ShapeMBTRIQ(
double *
N,
const double *X,
const double *Y,
264 PetscErrorCode
ShapeDiffMBTRIQ(
double *diffN,
const double *X,
const double *Y,
266 PetscErrorCode
ShapeMBTETQ(
double *
N,
const double x,
const double y,
268 PetscErrorCode
ShapeDiffMBTETQ(
double *diffN,
const double x,
const double y,
271 const double *Z,
const int G_DIM);
273 const double *Y,
const double *Z,
275 PetscErrorCode
ShapeJacMBTETQ(
const double *diffN,
const double *coords,
279 const double *diffN,
const double *coords,
284 const double *elem_coords,
285 const double *glob_coords,
286 double *loc_coords,
const double eps);
291 CBLAS_TRANSPOSE Trans);
322 static const double G_TRI_X1[] = {3.3333333333333331e-01};
323 static const double G_TRI_Y1[] = {3.3333333333333331e-01};
328 3.3333333333333331e-01, 3.3333333333333331e-01, 3.3333333333333331e-01};
330 7.503111022260811058e-02, 1.785587282636164064e-01,
331 2.800199154990741235e-01, 6.663902460147014262e-01};
333 2.800199154990741235e-01, 6.663902460147014262e-01,
334 7.503111022260811058e-02, 1.785587282636164064e-01};
336 1.8195861825602258066e-01, 3.1804138174397683647e-01,
337 1.8195861825602258066e-01, 3.1804138174397683647e-01};
339 0.333333333333333, 0.736712498968435, 0.736712498968435, 0.237932366472434,
340 0.237932366472434, 0.025355134551932, 0.025355134551932};
342 0.333333333333333, 0.237932366472434, 0.025355134551932, 0.736712498968435,
343 0.025355134551932, 0.736712498968435, 0.237932366472434};
345 0.375000000000000, 0.104166666666667, 0.104166666666667, 0.104166666666667,
346 0.104166666666667, 0.104166666666667, 0.104166666666667};
348 0.333333333333333, 0.479308067841923, 0.260345966079038, 0.260345966079038,
349 0.869739794195568, 0.065130102902216, 0.065130102902216, 0.638444188569809,
350 0.638444188569809, 0.312865496004875, 0.312865496004875, 0.048690315425316,
353 0.333333333333333, 0.260345966079038, 0.479308067841923, 0.260345966079038,
354 0.065130102902216, 0.869739794195568, 0.065130102902216, 0.312865496004875,
355 0.048690315425316, 0.638444188569809, 0.048690315425316, 0.638444188569809,
358 -0.149570044467670, 0.175615257433204, 0.175615257433204, 0.175615257433204,
359 0.053347235608839, 0.053347235608839, 0.053347235608839, 0.077113760890257,
360 0.077113760890257, 0.077113760890257, 0.077113760890257, 0.077113760890257,
364 0.333333333333333, 0.797426985353087, 0.101286507323456, 0.101286507323456,
365 0.059715871789770, 0.470142064105115, 0.470142064105115, 0.535795346449899,
366 0.232102326775050, 0.232102326775050, 0.941038278231121, 0.029480860884440,
367 0.029480860884440, 0.738416812340510, 0.738416812340510, 0.232102326775050,
368 0.232102326775050, 0.029480860884440, 0.029480860884440};
370 0.333333333333333, 0.101286507323456, 0.797426985353087, 0.101286507323456,
371 0.470142064105115, 0.059715871789770, 0.470142064105115, 0.232102326775050,
372 0.535795346449899, 0.232102326775050, 0.029480860884440, 0.941038278231121,
373 0.029480860884440, 0.232102326775050, 0.029480860884440, 0.738416812340510,
374 0.029480860884440, 0.738416812340510, 0.232102326775050};
376 9.71357962827961025E-002, 3.13347002271398278E-002,
377 3.13347002271398278E-002, 3.13347002271398278E-002,
378 7.78275410047754301E-002, 7.78275410047754301E-002,
379 7.78275410047754301E-002, 7.96477389272090969E-002,
380 7.96477389272090969E-002, 7.96477389272090969E-002,
381 2.55776756586981006E-002, 2.55776756586981006E-002,
382 2.55776756586981006E-002, 4.32835393772893970E-002,
383 4.32835393772893970E-002, 4.32835393772893970E-002,
384 4.32835393772893970E-002, 4.32835393772893970E-002,
385 4.32835393772893970E-002};
388 0.333333333333333, 0.948021718143423, 0.025989140928288,
389 0.025989140928288, 0.811424994704155, 0.094287502647923,
390 0.094287502647923, 0.010726449965571, 0.494636775017215,
391 0.494636775017215, 0.585313234770972, 0.207343382614514,
392 0.207343382614514, 0.122184388599019, 0.438907805700491,
393 0.438907805700491, 0.677937654882590, 0.677937654882590,
394 0.044841677589131, 0.044841677589131, 0.277220667528279,
395 0.277220667528279, 0.858870281282636, 0.858870281282636,
396 0.0000000000000000, 0.0000000000000000, 0.141129718717364,
399 0.333333333333333, 0.025989140928288, 0.948021718143423, 0.025989140928288,
400 0.094287502647923, 0.811424994704155, 0.094287502647923, 0.494636775017215,
401 0.010726449965571, 0.494636775017215, 0.207343382614514, 0.585313234770972,
402 0.207343382614514, 0.438907805700491, 0.122184388599019, 0.438907805700491,
403 0.044841677589131, 0.277220667528279, 0.677937654882590, 0.277220667528279,
404 0.677937654882590, 0.044841677589131, 0.000000000000000, 0.141129718717364,
405 0.858870281282636, 0.141129718717364, 0.858870281282636, 0.000000000000000};
407 0.08797730116222190, 0.008744311553736190, 0.008744311553736190,
408 0.008744311553736190, 0.03808157199393533, 0.03808157199393533,
409 0.03808157199393533, 0.01885544805613125, 0.01885544805613125,
410 0.01885544805613125, 0.07215969754474100, 0.07215969754474100,
411 0.07215969754474100, 0.06932913870553720, 0.06932913870553720,
412 0.06932913870553720, 0.04105631542928860, 0.04105631542928860,
413 0.04105631542928860, 0.04105631542928860, 0.04105631542928860,
414 0.04105631542928860, 0.007362383783300573, 0.007362383783300573,
415 0.007362383783300573, 0.007362383783300573, 0.007362383783300573,
416 0.007362383783300573};
419 0.333333333333333, 0.950275662924106, 0.024862168537947, 0.024862168537947,
420 0.171614914923835, 0.414192542538082, 0.414192542538082, 0.539412243677190,
421 0.230293878161405, 0.230293878161405, 0.772160036676533, 0.113919981661734,
422 0.113919981661734, 0.009085399949835, 0.495457300025082, 0.495457300025082,
423 0.062277290305887, 0.468861354847056, 0.468861354847056, 0.022076289653624,
424 0.022076289653624, 0.851306504174348, 0.851306504174348, 0.126617206172027,
425 0.126617206172027, 0.018620522802521, 0.018620522802521, 0.689441970728591,
426 0.689441970728591, 0.291937506468888, 0.291937506468888, 0.096506481292159,
427 0.096506481292159, 0.635867859433873, 0.635867859433873, 0.267625659273968,
430 0.333333333333333, 0.024862168537947, 0.950275662924106, 0.024862168537947,
431 0.414192542538082, 0.171614914923835, 0.414192542538082, 0.230293878161405,
432 0.539412243677190, 0.230293878161405, 0.113919981661734, 0.772160036676533,
433 0.113919981661734, 0.495457300025082, 0.009085399949835, 0.495457300025082,
434 0.468861354847056, 0.062277290305887, 0.468861354847056, 0.851306504174348,
435 0.126617206172027, 0.022076289653624, 0.126617206172027, 0.022076289653624,
436 0.851306504174348, 0.689441970728591, 0.291937506468888, 0.018620522802521,
437 0.291937506468888, 0.018620522802521, 0.689441970728591, 0.635867859433873,
438 0.267625659273968, 0.096506481292159, 0.267625659273968, 0.096506481292159,
441 0.051739766065744, 0.008007799555565, 0.008007799555565, 0.008007799555565,
442 0.046868898981822, 0.046868898981822, 0.046868898981822, 0.046590940183976,
443 0.046590940183976, 0.046590940183976, 0.031016943313796, 0.031016943313796,
444 0.031016943313796, 0.010791612736631, 0.010791612736631, 0.010791612736631,
445 0.032195534242432, 0.032195534242432, 0.032195534242432, 0.015445834210702,
446 0.015445834210702, 0.015445834210702, 0.015445834210702, 0.015445834210702,
447 0.015445834210702, 0.017822989923179, 0.017822989923179, 0.017822989923179,
448 0.017822989923179, 0.017822989923179, 0.017822989923179, 0.037038683681385,
449 0.037038683681385, 0.037038683681385, 0.037038683681385, 0.037038683681385,
1022 0.3333333333333333};
1024 2.912193380035668, 2.912193380035668, 2.912193380035668,
1025 2.912193380035668, 2.912193380035668, 2.912193380035668,
1026 2.912193380035668, 2.912193380035668, 2.912193380035668,
1027 2.912193380035668, 2.912193380035668, 2.912193380035668,
1028 2.912193380035668, 2.912193380035668, 2.912193380035668,
1029 2.912193380035668, 2.912193380035668, 2.912193380035668,
1030 2.912193380035668, 2.912193380035668, 2.912193380035668,
1031 2.912193380035668, 2.912193380035668, 2.912193380035668,
1032 2.912193380035668, 2.912193380035668, 2.912193380035668,
1033 2.912193380035668, 2.912193380035668, 2.912193380035668,
1034 2.912193380035668, 2.912193380035668, 2.912193380035668,
1035 2.912193380035668, 2.912193380035668, 2.912193380035668,
1036 2.912193380035668, 2.912193380035668, 2.912193380035668,
1037 2.912193380035668, 2.912193380035668, 2.912193380035668,
1038 2.912193380035668, 2.912193380035668, 2.912193380035668,
1039 2.912193380035668, 2.912193380035668, 2.912193380035668,
1040 2.912193380035668, 2.912193380035668, 2.912193380035668,
1041 2.912193380035668, 2.912193380035668, 2.912193380035668,
1042 2.912193380035668, 2.912193380035668, 2.912193380035668,
1043 2.912193380035668, 2.912193380035668, 2.912193380035668,
1044 2.912193380035668, 2.912193380035668, 2.912193380035668,
1045 2.912193380035668, 2.912193380035668, 2.912193380035668,
1046 -9.914451197589852, -9.914451197589852, -9.914451197589852,
1047 -9.914451197589852, -9.914451197589852, -9.914451197589852,
1048 -9.914451197589852, -9.914451197589852, -9.914451197589852,
1049 -9.914451197589852, -9.914451197589852, -9.914451197589852,
1050 -9.914451197589852, -9.914451197589852, -9.914451197589852,
1051 -9.914451197589852, -9.914451197589852, -9.914451197589852,
1052 -9.914451197589852, -9.914451197589852, -9.914451197589852,
1053 -9.914451197589852, -9.914451197589852, -9.914451197589852,
1054 -9.914451197589852, -9.914451197589852, -9.914451197589852,
1055 -9.914451197589852, -9.914451197589852, -9.914451197589852,
1056 -9.914451197589852, -9.914451197589852, -9.914451197589852,
1057 -9.914451197589852, -9.914451197589852, -9.914451197589852,
1058 -9.914451197589852, -9.914451197589852, -9.914451197589852,
1059 -9.914451197589852, -9.914451197589852, -9.914451197589852,
1060 -9.914451197589852, -9.914451197589852, -9.914451197589852,
1061 -9.914451197589852, -9.914451197589852, -9.914451197589852,
1062 -9.914451197589852, -9.914451197589852, -9.914451197589852,
1063 -9.914451197589852, -9.914451197589852, -9.914451197589852,
1064 -9.914451197589852, 13.33158527957992, 13.33158527957992,
1065 13.33158527957992, 13.33158527957992, 13.33158527957992,
1066 13.33158527957992, 13.33158527957992, 13.33158527957992,
1067 13.33158527957992, 13.33158527957992, 13.33158527957992,
1068 13.33158527957992, 13.33158527957992, 13.33158527957992,
1069 13.33158527957992, 13.33158527957992, 13.33158527957992,
1070 13.33158527957992, 13.33158527957992, 13.33158527957992,
1071 13.33158527957992, 13.33158527957992, 13.33158527957992,
1072 13.33158527957992, 13.33158527957992, 13.33158527957992,
1073 13.33158527957992, 13.33158527957992, 13.33158527957992,
1074 13.33158527957992, 13.33158527957992, 13.33158527957992,
1075 13.33158527957992, 13.33158527957992, 13.33158527957992,
1076 13.33158527957992, 13.33158527957992, 13.33158527957992,
1077 13.33158527957992, 13.33158527957992, 13.33158527957992,
1078 13.33158527957992, 13.33158527957992, 13.33158527957992,
1079 13.33158527957992, -9.027792408986382, -9.027792408986382,
1080 -9.027792408986382, -9.027792408986382, -9.027792408986382,
1081 -9.027792408986382, -9.027792408986382, -9.027792408986382,
1082 -9.027792408986382, -9.027792408986382, -9.027792408986382,
1083 -9.027792408986382, -9.027792408986382, -9.027792408986382,
1084 -9.027792408986382, -9.027792408986382, -9.027792408986382,
1085 -9.027792408986382, -9.027792408986382, -9.027792408986382,
1086 -9.027792408986382, -9.027792408986382, -9.027792408986382,
1087 -9.027792408986382, -9.027792408986382, -9.027792408986382,
1088 -9.027792408986382, -9.027792408986382, -9.027792408986382,
1089 -9.027792408986382, -9.027792408986382, -9.027792408986382,
1090 -9.027792408986382, -9.027792408986382, -9.027792408986382,
1091 -9.027792408986382, 3.258672079964582, 3.258672079964582,
1092 3.258672079964582, 3.258672079964582, 3.258672079964582,
1093 3.258672079964582, 3.258672079964582, 3.258672079964582,
1094 3.258672079964582, 3.258672079964582, 3.258672079964582,
1095 3.258672079964582, 3.258672079964582, 3.258672079964582,
1096 3.258672079964582, 3.258672079964582, 3.258672079964582,
1097 3.258672079964582, 3.258672079964582, 3.258672079964582,
1098 3.258672079964582, 3.258672079964582, 3.258672079964582,
1099 3.258672079964582, 3.258672079964582, 3.258672079964582,
1100 3.258672079964582, 3.258672079964582, -0.6133639040302452,
1101 -0.6133639040302452, -0.6133639040302452, -0.6133639040302452,
1102 -0.6133639040302452, -0.6133639040302452, -0.6133639040302452,
1103 -0.6133639040302452, -0.6133639040302452, -0.6133639040302452,
1104 -0.6133639040302452, -0.6133639040302452, -0.6133639040302452,
1105 -0.6133639040302452, -0.6133639040302452, -0.6133639040302452,
1106 -0.6133639040302452, -0.6133639040302452, -0.6133639040302452,
1107 -0.6133639040302452, -0.6133639040302452, 0.05511571669513555,
1108 0.05511571669513555, 0.05511571669513555, 0.05511571669513555,
1109 0.05511571669513555, 0.05511571669513555, 0.05511571669513555,
1110 0.05511571669513555, 0.05511571669513555, 0.05511571669513555,
1111 0.05511571669513555, 0.05511571669513555, 0.05511571669513555,
1112 0.05511571669513555, 0.05511571669513555, -0.001979122382447095,
1113 -0.001979122382447095, -0.001979122382447095, -0.001979122382447095,
1114 -0.001979122382447095, -0.001979122382447095, -0.001979122382447095,
1115 -0.001979122382447095, -0.001979122382447095, -0.001979122382447095,
1116 2.02054621415273e-05, 2.02054621415273e-05, 2.02054621415273e-05,
1117 2.02054621415273e-05, 2.02054621415273e-05, 2.02054621415273e-05,
1118 -2.874940020535803e-08, -2.874940020535803e-08, -2.874940020535803e-08,
1119 8.829438425435718e-13};
1126 static const double G_TET_X4[] = {0.1757281246520584, 0.2445310270213291,
1127 0.5556470949048655, 0.0240937534217468};
1128 static const double G_TET_Y4[] = {0.5656137776620919, 0.0501800797762026,
1129 0.1487681308666864, 0.2354380116950194};
1130 static const double G_TET_Z4[] = {0.2180665126782654, 0.5635595064952189,
1131 0.0350112499848832, 0.1833627308416330};
1132 static const double G_TET_W4[] = {0.25, 0.25, 0.25, 0.25};
1133 static const double G_TET_X5[] = {0.25000000000000000, 0.50000000000000000,
1134 0.16666666666666667, 0.16666666666666667,
1135 0.16666666666666667};
1136 static const double G_TET_Y5[] = {0.25000000000000000, 0.16666666666666667,
1137 0.50000000000000000, 0.16666666666666667,
1138 0.16666666666666667};
1139 static const double G_TET_Z5[] = {0.25000000000000000, 0.16666666666666667,
1140 0.16666666666666667, 0.50000000000000000,
1141 0.16666666666666667};
1142 static const double G_TET_W5[] = {-0.80000000000000000, 0.45000000000000000,
1143 0.45000000000000000, 0.45000000000000000,
1144 0.45000000000000000};
1145 static const double G_TET_X10[] = {0.5684305841968444, 0.1438564719343852,
1146 0.1438564719343852, 0.1438564719343852,
1147 0.0000000000000000, 0.5000000000000000,
1148 0.5000000000000000, 0.5000000000000000,
1149 0.0000000000000000, 0.0000000000000000};
1150 static const double G_TET_Y10[] = {0.1438564719343852, 0.1438564719343852,
1151 0.1438564719343852, 0.5684305841968444,
1152 0.5000000000000000, 0.0000000000000000,
1153 0.5000000000000000, 0.0000000000000000,
1154 0.5000000000000000, 0.0000000000000000};
1155 static const double G_TET_Z10[] = {0.1438564719343852, 0.1438564719343852,
1156 0.5684305841968444, 0.1438564719343852,
1157 0.5000000000000000, 0.5000000000000000,
1158 0.0000000000000000, 0.0000000000000000,
1159 0.0000000000000000, 0.5000000000000000};
1160 static const double G_TET_W10[] = {0.2177650698804054, 0.2177650698804054,
1161 0.2177650698804054, 0.2177650698804054,
1162 0.0214899534130631, 0.0214899534130631,
1163 0.0214899534130631, 0.0214899534130631,
1164 0.0214899534130631, 0.0214899534130631};
1167 0.2500000000000000, 0.6175871903000830, 0.1274709365666390,
1168 0.1274709365666390, 0.1274709365666390, 0.9037635088221031,
1169 0.0320788303926323, 0.0320788303926323, 0.0320788303926323,
1170 0.4502229043567190, 0.0497770956432810, 0.0497770956432810,
1171 0.0497770956432810, 0.4502229043567190, 0.4502229043567190,
1172 0.3162695526014501, 0.1837304473985499, 0.1837304473985499,
1173 0.1837304473985499, 0.3162695526014501, 0.3162695526014501,
1174 0.0229177878448171, 0.2319010893971509, 0.2319010893971509,
1175 0.5132800333608811, 0.2319010893971509, 0.2319010893971509,
1176 0.2319010893971509, 0.0229177878448171, 0.5132800333608811,
1177 0.2319010893971509, 0.0229177878448171, 0.5132800333608811,
1178 0.7303134278075384, 0.0379700484718286, 0.0379700484718286,
1179 0.1937464752488044, 0.0379700484718286, 0.0379700484718286,
1180 0.0379700484718286, 0.7303134278075384, 0.1937464752488044,
1181 0.0379700484718286, 0.7303134278075384, 0.1937464752488044};
1183 0.2500000000000000, 0.1274709365666390, 0.1274709365666390,
1184 0.1274709365666390, 0.6175871903000830, 0.0320788303926323,
1185 0.0320788303926323, 0.0320788303926323, 0.9037635088221031,
1186 0.0497770956432810, 0.4502229043567190, 0.0497770956432810,
1187 0.4502229043567190, 0.0497770956432810, 0.4502229043567190,
1188 0.1837304473985499, 0.3162695526014501, 0.1837304473985499,
1189 0.3162695526014501, 0.1837304473985499, 0.3162695526014501,
1190 0.2319010893971509, 0.0229177878448171, 0.2319010893971509,
1191 0.2319010893971509, 0.5132800333608811, 0.2319010893971509,
1192 0.0229177878448171, 0.5132800333608811, 0.2319010893971509,
1193 0.5132800333608811, 0.2319010893971509, 0.0229177878448171,
1194 0.0379700484718286, 0.7303134278075384, 0.0379700484718286,
1195 0.0379700484718286, 0.1937464752488044, 0.0379700484718286,
1196 0.7303134278075384, 0.1937464752488044, 0.0379700484718286,
1197 0.1937464752488044, 0.0379700484718286, 0.7303134278075384};
1199 0.2500000000000000, 0.1274709365666390, 0.1274709365666390,
1200 0.6175871903000830, 0.1274709365666390, 0.0320788303926323,
1201 0.0320788303926323, 0.9037635088221031, 0.0320788303926323,
1202 0.0497770956432810, 0.0497770956432810, 0.4502229043567190,
1203 0.4502229043567190, 0.4502229043567190, 0.0497770956432810,
1204 0.1837304473985499, 0.1837304473985499, 0.3162695526014501,
1205 0.3162695526014501, 0.3162695526014501, 0.1837304473985499,
1206 0.2319010893971509, 0.2319010893971509, 0.0229177878448171,
1207 0.2319010893971509, 0.2319010893971509, 0.5132800333608811,
1208 0.5132800333608811, 0.2319010893971509, 0.0229177878448171,
1209 0.0229177878448171, 0.5132800333608811, 0.2319010893971509,
1210 0.0379700484718286, 0.0379700484718286, 0.7303134278075384,
1211 0.0379700484718286, 0.0379700484718286, 0.1937464752488044,
1212 0.1937464752488044, 0.0379700484718286, 0.7303134278075384,
1213 0.7303134278075384, 0.1937464752488044, 0.0379700484718286};
1215 -0.2359620398477557, 0.0244878963560562, 0.0244878963560562,
1216 0.0244878963560562, 0.0244878963560562, 0.0039485206398261,
1217 0.0039485206398261, 0.0039485206398261, 0.0039485206398261,
1218 0.0263055529507371, 0.0263055529507371, 0.0263055529507371,
1219 0.0263055529507371, 0.0263055529507371, 0.0263055529507371,
1220 0.0829803830550589, 0.0829803830550589, 0.0829803830550589,
1221 0.0829803830550589, 0.0829803830550589, 0.0829803830550589,
1222 0.0254426245481023, 0.0254426245481023, 0.0254426245481023,
1223 0.0254426245481023, 0.0254426245481023, 0.0254426245481023,
1224 0.0254426245481023, 0.0254426245481023, 0.0254426245481023,
1225 0.0254426245481023, 0.0254426245481023, 0.0254426245481023,
1226 0.0134324384376852, 0.0134324384376852, 0.0134324384376852,
1227 0.0134324384376852, 0.0134324384376852, 0.0134324384376852,
1228 0.0134324384376852, 0.0134324384376852, 0.0134324384376852,
1229 0.0134324384376852, 0.0134324384376852, 0.0134324384376852};
1231 0.1000000000000000, 0.1000000000000000, 0.1000000000000000,
1232 0.7000000000000000, 0.1000000000000000, 0.1000000000000000,
1233 0.1000000000000000, 0.1000000000000000, 0.1000000000000000,
1234 0.1000000000000000, 0.6000000000000000, 0.6000000000000000,
1235 0.6000000000000000, 0.2000000000000000, 0.2000000000000000,
1236 0.2000000000000000, 0.1000000000000000, 0.1000000000000000,
1237 0.1000000000000000, 0.1000000000000000, 0.1000000000000000,
1238 0.1000000000000000, 0.5000000000000000, 0.5000000000000000,
1239 0.5000000000000000, 0.3000000000000000, 0.3000000000000000,
1240 0.3000000000000000, 0.2000000000000000, 0.2000000000000000,
1241 0.2000000000000000, 0.2000000000000000, 0.2000000000000000,
1242 0.2000000000000000, 0.1000000000000000, 0.1000000000000000,
1243 0.1000000000000000, 0.5000000000000000, 0.5000000000000000,
1244 0.5000000000000000, 0.4000000000000000, 0.4000000000000000,
1245 0.4000000000000000, 0.1000000000000000, 0.1000000000000000,
1246 0.1000000000000000, 0.4000000000000000, 0.4000000000000000,
1247 0.4000000000000000, 0.4000000000000000, 0.4000000000000000,
1248 0.4000000000000000, 0.3000000000000000, 0.3000000000000000,
1249 0.3000000000000000, 0.3000000000000000, 0.3000000000000000,
1250 0.3000000000000000, 0.2000000000000000, 0.2000000000000000,
1251 0.2000000000000000, 0.2000000000000000, 0.2000000000000000,
1252 0.2000000000000000, 0.1000000000000000, 0.1000000000000000,
1253 0.1000000000000000, 0.1000000000000000, 0.1000000000000000,
1254 0.1000000000000000, 0.2000000000000000, 0.2000000000000000,
1255 0.2000000000000000, 0.4000000000000000, 0.3000000000000000,
1256 0.3000000000000000, 0.3000000000000000, 0.1000000000000000,
1257 0.3000000000000000, 0.3000000000000000, 0.3000000000000000,
1258 0.2000000000000000, 0.2000000000000000, 0.2000000000000000};
1260 0.1000000000000000, 0.1000000000000000, 0.7000000000000000,
1261 0.1000000000000000, 0.1000000000000000, 0.1000000000000000,
1262 0.6000000000000000, 0.6000000000000000, 0.2000000000000000,
1263 0.2000000000000000, 0.1000000000000000, 0.1000000000000000,
1264 0.2000000000000000, 0.1000000000000000, 0.1000000000000000,
1265 0.6000000000000000, 0.1000000000000000, 0.1000000000000000,
1266 0.5000000000000000, 0.5000000000000000, 0.3000000000000000,
1267 0.3000000000000000, 0.1000000000000000, 0.1000000000000000,
1268 0.3000000000000000, 0.1000000000000000, 0.1000000000000000,
1269 0.5000000000000000, 0.2000000000000000, 0.2000000000000000,
1270 0.1000000000000000, 0.1000000000000000, 0.5000000000000000,
1271 0.5000000000000000, 0.2000000000000000, 0.2000000000000000,
1272 0.5000000000000000, 0.2000000000000000, 0.2000000000000000,
1273 0.1000000000000000, 0.4000000000000000, 0.1000000000000000,
1274 0.1000000000000000, 0.4000000000000000, 0.4000000000000000,
1275 0.1000000000000000, 0.3000000000000000, 0.3000000000000000,
1276 0.2000000000000000, 0.2000000000000000, 0.1000000000000000,
1277 0.1000000000000000, 0.4000000000000000, 0.4000000000000000,
1278 0.2000000000000000, 0.2000000000000000, 0.1000000000000000,
1279 0.1000000000000000, 0.4000000000000000, 0.4000000000000000,
1280 0.3000000000000000, 0.3000000000000000, 0.1000000000000000,
1281 0.1000000000000000, 0.4000000000000000, 0.4000000000000000,
1282 0.3000000000000000, 0.3000000000000000, 0.2000000000000000,
1283 0.2000000000000000, 0.2000000000000000, 0.2000000000000000,
1284 0.4000000000000000, 0.2000000000000000, 0.3000000000000000,
1285 0.3000000000000000, 0.1000000000000000, 0.3000000000000000,
1286 0.3000000000000000, 0.2000000000000000, 0.2000000000000000,
1287 0.3000000000000000, 0.3000000000000000, 0.2000000000000000};
1289 0.1000000000000000, 0.7000000000000000, 0.1000000000000000,
1290 0.1000000000000000, 0.6000000000000000, 0.2000000000000000,
1291 0.1000000000000000, 0.2000000000000000, 0.1000000000000000,
1292 0.6000000000000000, 0.1000000000000000, 0.2000000000000000,
1293 0.1000000000000000, 0.1000000000000000, 0.6000000000000000,
1294 0.1000000000000000, 0.5000000000000000, 0.3000000000000000,
1295 0.1000000000000000, 0.3000000000000000, 0.1000000000000000,
1296 0.5000000000000000, 0.1000000000000000, 0.3000000000000000,
1297 0.1000000000000000, 0.1000000000000000, 0.5000000000000000,
1298 0.1000000000000000, 0.1000000000000000, 0.5000000000000000,
1299 0.2000000000000000, 0.5000000000000000, 0.2000000000000000,
1300 0.1000000000000000, 0.2000000000000000, 0.5000000000000000,
1301 0.2000000000000000, 0.2000000000000000, 0.1000000000000000,
1302 0.2000000000000000, 0.1000000000000000, 0.4000000000000000,
1303 0.1000000000000000, 0.4000000000000000, 0.1000000000000000,
1304 0.4000000000000000, 0.2000000000000000, 0.1000000000000000,
1305 0.3000000000000000, 0.1000000000000000, 0.3000000000000000,
1306 0.2000000000000000, 0.2000000000000000, 0.1000000000000000,
1307 0.4000000000000000, 0.1000000000000000, 0.4000000000000000,
1308 0.2000000000000000, 0.3000000000000000, 0.1000000000000000,
1309 0.4000000000000000, 0.1000000000000000, 0.4000000000000000,
1310 0.3000000000000000, 0.3000000000000000, 0.2000000000000000,
1311 0.4000000000000000, 0.2000000000000000, 0.4000000000000000,
1312 0.3000000000000000, 0.2000000000000000, 0.4000000000000000,
1313 0.2000000000000000, 0.2000000000000000, 0.3000000000000000,
1314 0.1000000000000000, 0.3000000000000000, 0.3000000000000000,
1315 0.2000000000000000, 0.3000000000000000, 0.2000000000000000,
1316 0.3000000000000000, 0.2000000000000000, 0.3000000000000000};
1318 0.2843915343915344, 0.2843915343915344, 0.2843915343915344,
1319 0.2843915343915344, -0.3882275132275133, -0.3882275132275133,
1320 -0.3882275132275133, -0.3882275132275133, -0.3882275132275133,
1321 -0.3882275132275133, -0.3882275132275133, -0.3882275132275133,
1322 -0.3882275132275133, -0.3882275132275133, -0.3882275132275133,
1323 -0.3882275132275133, 0.8776455026455027, 0.8776455026455027,
1324 0.8776455026455027, 0.8776455026455027, 0.8776455026455027,
1325 0.8776455026455027, 0.8776455026455027, 0.8776455026455027,
1326 0.8776455026455027, 0.8776455026455027, 0.8776455026455027,
1327 0.8776455026455027, 0.1236772486772487, 0.1236772486772487,
1328 0.1236772486772487, 0.1236772486772487, 0.1236772486772487,
1329 0.1236772486772487, 0.1236772486772487, 0.1236772486772487,
1330 0.1236772486772487, 0.1236772486772487, 0.1236772486772487,
1331 0.1236772486772487, -0.8584656084656085, -0.8584656084656085,
1332 -0.8584656084656085, -0.8584656084656085, -0.8584656084656085,
1333 -0.8584656084656085, -0.2632275132275133, -0.2632275132275133,
1334 -0.2632275132275133, -0.2632275132275133, -0.2632275132275133,
1335 -0.2632275132275133, -0.2632275132275133, -0.2632275132275133,
1336 -0.2632275132275133, -0.2632275132275133, -0.2632275132275133,
1337 -0.2632275132275133, -0.2632275132275133, -0.2632275132275133,
1338 -0.2632275132275133, -0.2632275132275133, -0.2632275132275133,
1339 -0.2632275132275133, -0.2632275132275133, -0.2632275132275133,
1340 -0.2632275132275133, -0.2632275132275133, -0.2632275132275133,
1341 -0.2632275132275133, 0.0145502645502645, 0.0145502645502645,
1342 0.0145502645502645, 0.0145502645502645, 1.0165343915343916,
1343 1.0165343915343916, 1.0165343915343916, 1.0165343915343916,
1344 -0.0251322751322751, -0.0251322751322751, -0.0251322751322751,
1345 -0.0251322751322751, -0.0251322751322751, -0.0251322751322751};
FTensor::Index< 'm', 3 > m