22#define LAMBDA(E, NU) (E * NU / ((1. + NU) * (1. - 2. * NU)))
23#define MU(E, NU) (0.5 * E / (1. + NU))
24#define DELTA(NU_P, NU_PZ, E_P, E_Z) \
25 (((1 + NU_P) * (1 - NU_P - 2 * NU_PZ * (NU_PZ * E_Z / E_P))) / \
28#define N_MBTET0(x, y, z) (1. - x - y - z)
29#define N_MBTET1(x, y, z) (x)
30#define N_MBTET2(x, y, z) (y)
31#define N_MBTET3(x, y, z) (z)
32#define diffN_MBTET0x (-1.)
33#define diffN_MBTET0y (-1.)
34#define diffN_MBTET0z (-1.)
35#define diffN_MBTET1x (1.)
36#define diffN_MBTET1y (0.)
37#define diffN_MBTET1z (0.)
38#define diffN_MBTET2x (0.)
39#define diffN_MBTET2y (1.)
40#define diffN_MBTET2z (0.)
41#define diffN_MBTET3x (0.)
42#define diffN_MBTET3y (0.)
43#define diffN_MBTET3z (1.)
46#define N_MBTRI0(x, y) (1. - x - y)
47#define N_MBTRI1(x, y) (x)
48#define N_MBTRI2(x, y) (y)
49#define diffN_MBTRI0x (-1.)
50#define diffN_MBTRI0y (-1.)
51#define diffN_MBTRI1x (1.)
52#define diffN_MBTRI1y (0.)
53#define diffN_MBTRI2x (0.)
54#define diffN_MBTRI2y (1.)
57#define N_MBQUAD0(x, y) ((1. - x) * (1. - y))
58#define N_MBQUAD1(x, y) ((x) * (1. - y))
59#define N_MBQUAD2(x, y) ((x) * (y))
60#define N_MBQUAD3(x, y) ((1. - x) * (y))
61#define diffN_MBQUAD0x(y) (-(1. - y))
62#define diffN_MBQUAD0y(x) (-(1. - x))
63#define diffN_MBQUAD1x(y) ((1. - y))
64#define diffN_MBQUAD1y(x) (-x)
65#define diffN_MBQUAD2x(y) (y)
66#define diffN_MBQUAD2y(x) (x)
67#define diffN_MBQUAD3x(y) (-y)
68#define diffN_MBQUAD3y(x) ((1. - x))
71#define N_MBHEX0(x, y, z) (N_MBQUAD0(x, y) * (1 - z))
72#define N_MBHEX1(x, y, z) (N_MBQUAD1(x, y) * (1 - z))
73#define N_MBHEX2(x, y, z) (N_MBQUAD2(x, y) * (1 - z))
74#define N_MBHEX3(x, y, z) (N_MBQUAD3(x, y) * (1 - z))
75#define N_MBHEX4(x, y, z) (N_MBQUAD0(x, y) * (z))
76#define N_MBHEX5(x, y, z) (N_MBQUAD1(x, y) * (z))
77#define N_MBHEX6(x, y, z) (N_MBQUAD2(x, y) * (z))
78#define N_MBHEX7(x, y, z) (N_MBQUAD3(x, y) * (z))
79#define diffN_MBHEX0x(y, z) (diffN_MBQUAD0x(y) * (1 - z))
80#define diffN_MBHEX1x(y, z) (diffN_MBQUAD1x(y) * (1 - z))
81#define diffN_MBHEX2x(y, z) (diffN_MBQUAD2x(y) * (1 - z))
82#define diffN_MBHEX3x(y, z) (diffN_MBQUAD3x(y) * (1 - z))
83#define diffN_MBHEX4x(y, z) (diffN_MBQUAD0x(y) * (z))
84#define diffN_MBHEX5x(y, z) (diffN_MBQUAD1x(y) * (z))
85#define diffN_MBHEX6x(y, z) (diffN_MBQUAD2x(y) * (z))
86#define diffN_MBHEX7x(y, z) (diffN_MBQUAD3x(y) * (z))
87#define diffN_MBHEX0y(x, z) (diffN_MBQUAD0y(x) * (1 - z))
88#define diffN_MBHEX1y(x, z) (diffN_MBQUAD1y(x) * (1 - z))
89#define diffN_MBHEX2y(x, z) (diffN_MBQUAD2y(x) * (1 - z))
90#define diffN_MBHEX3y(x, z) (diffN_MBQUAD3y(x) * (1 - z))
91#define diffN_MBHEX4y(x, z) (diffN_MBQUAD0y(x) * (z))
92#define diffN_MBHEX5y(x, z) (diffN_MBQUAD1y(x) * (z))
93#define diffN_MBHEX6y(x, z) (diffN_MBQUAD2y(x) * (z))
94#define diffN_MBHEX7y(x, z) (diffN_MBQUAD3y(x) * (z))
95#define diffN_MBHEX0z(x, y) (-N_MBQUAD0(x, y))
96#define diffN_MBHEX1z(x, y) (-N_MBQUAD1(x, y))
97#define diffN_MBHEX2z(x, y) (-N_MBQUAD2(x, y))
98#define diffN_MBHEX3z(x, y) (-N_MBQUAD3(x, y))
99#define diffN_MBHEX4z(x, y) (N_MBQUAD0(x, y))
100#define diffN_MBHEX5z(x, y) (N_MBQUAD1(x, y))
101#define diffN_MBHEX6z(x, y) (N_MBQUAD2(x, y))
102#define diffN_MBHEX7z(x, y) (N_MBQUAD3(x, y))
105#define N_MBEDGE0(x) (1. - (x))
106#define N_MBEDGE1(x) (x)
107#define diffN_MBEDGE0 (-1.)
108#define diffN_MBEDGE1 (1.)
111#define N_MBTRIQ0(x, y) ((1. - x - y) * (2 * (1. - x - y) - 1.))
112#define N_MBTRIQ1(x, y) (x * (2. * x - 1.))
113#define N_MBTRIQ2(x, y) (y * (2. * y - 1.))
114#define N_MBTRIQ3(x, y) (4. * (1. - x - y) * x)
115#define N_MBTRIQ4(x, y) (4. * x * y)
116#define N_MBTRIQ5(x, y) (4. * (1. - x - y) * y)
117#define diffN_MBTRIQ0x(x, y) (x + y - 3. * (1. - x - y))
118#define diffN_MBTRIQ0y(x, y) (x + y - 3. * (1. - x - y))
119#define diffN_MBTRIQ1x(x, y) (-1. + 4. * x)
120#define diffN_MBTRIQ1y(x, y) (0.)
121#define diffN_MBTRIQ2x(x, y) (0.)
122#define diffN_MBTRIQ2y(x, y) (-1. + 4. * y)
123#define diffN_MBTRIQ3x(x, y) (4. * ((1. - x - y) - x))
124#define diffN_MBTRIQ3y(x, y) (-4. * x)
125#define diffN_MBTRIQ4x(x, y) (4. * y)
126#define diffN_MBTRIQ4y(x, y) (4. * x)
127#define diffN_MBTRIQ5x(x, y) (-4. * y)
128#define diffN_MBTRIQ5y(x, y) (4. * ((1. - x - y) - y))
146PetscErrorCode
ShapeMBTRI(
double *
N,
const double *X,
const double *Y,
158 double *normal,
double *s1,
double *s2);
162 double *diff_normal);
168PetscErrorCode
ShapeMBTET(
double *
N,
const double *G_X,
const double *G_Y,
169 const double *G_Z,
int DIM);
175PetscErrorCode
ShapeJacMBTET(
double *diffN,
const double *coords,
double *jac);
182 double *diffNinvJac);
187PetscErrorCode
Spin(
double *spinOmega,
double *vecOmega);
195 int *order_edge,
double *diffN,
double *diffN_face,
196 double *diffN_edge[],
double *dofs,
double *dofs_edge[],
197 double *dofs_face,
double *idofs,
double *idofs_edge[],
209 const double *elem_coords,
210 const double *glob_coords,
240 const double *elem_coords,
241 const double *glob_coords,
248PetscErrorCode
ShapeMBEDGE(
double *
N,
const double *G_X,
int DIM);
252PetscErrorCode
ShapeMBTRIQ(
double *
N,
const double *X,
const double *Y,
254PetscErrorCode
ShapeDiffMBTRIQ(
double *diffN,
const double *X,
const double *Y,
256PetscErrorCode
ShapeMBTETQ(
double *
N,
const double x,
const double y,
258PetscErrorCode
ShapeDiffMBTETQ(
double *diffN,
const double x,
const double y,
261 const double *Z,
const int G_DIM);
263 const double *Y,
const double *Z,
265PetscErrorCode
ShapeJacMBTETQ(
const double *diffN,
const double *coords,
269 const double *diffN,
const double *coords,
274 const double *elem_coords,
275 const double *glob_coords,
276 double *loc_coords,
const double eps);
281 enum CBLAS_TRANSPOSE Trans);
312static const double G_TRI_X1[] = {3.3333333333333331e-01};
313static const double G_TRI_Y1[] = {3.3333333333333331e-01};
318 3.3333333333333331e-01, 3.3333333333333331e-01, 3.3333333333333331e-01};
320 7.503111022260811058e-02, 1.785587282636164064e-01,
321 2.800199154990741235e-01, 6.663902460147014262e-01};
323 2.800199154990741235e-01, 6.663902460147014262e-01,
324 7.503111022260811058e-02, 1.785587282636164064e-01};
326 1.8195861825602258066e-01, 3.1804138174397683647e-01,
327 1.8195861825602258066e-01, 3.1804138174397683647e-01};
329 0.333333333333333, 0.736712498968435, 0.736712498968435, 0.237932366472434,
330 0.237932366472434, 0.025355134551932, 0.025355134551932};
332 0.333333333333333, 0.237932366472434, 0.025355134551932, 0.736712498968435,
333 0.025355134551932, 0.736712498968435, 0.237932366472434};
335 0.375000000000000, 0.104166666666667, 0.104166666666667, 0.104166666666667,
336 0.104166666666667, 0.104166666666667, 0.104166666666667};
338 0.333333333333333, 0.479308067841923, 0.260345966079038, 0.260345966079038,
339 0.869739794195568, 0.065130102902216, 0.065130102902216, 0.638444188569809,
340 0.638444188569809, 0.312865496004875, 0.312865496004875, 0.048690315425316,
343 0.333333333333333, 0.260345966079038, 0.479308067841923, 0.260345966079038,
344 0.065130102902216, 0.869739794195568, 0.065130102902216, 0.312865496004875,
345 0.048690315425316, 0.638444188569809, 0.048690315425316, 0.638444188569809,
348 -0.149570044467670, 0.175615257433204, 0.175615257433204, 0.175615257433204,
349 0.053347235608839, 0.053347235608839, 0.053347235608839, 0.077113760890257,
350 0.077113760890257, 0.077113760890257, 0.077113760890257, 0.077113760890257,
354 0.333333333333333, 0.797426985353087, 0.101286507323456, 0.101286507323456,
355 0.059715871789770, 0.470142064105115, 0.470142064105115, 0.535795346449899,
356 0.232102326775050, 0.232102326775050, 0.941038278231121, 0.029480860884440,
357 0.029480860884440, 0.738416812340510, 0.738416812340510, 0.232102326775050,
358 0.232102326775050, 0.029480860884440, 0.029480860884440};
360 0.333333333333333, 0.101286507323456, 0.797426985353087, 0.101286507323456,
361 0.470142064105115, 0.059715871789770, 0.470142064105115, 0.232102326775050,
362 0.535795346449899, 0.232102326775050, 0.029480860884440, 0.941038278231121,
363 0.029480860884440, 0.232102326775050, 0.029480860884440, 0.738416812340510,
364 0.029480860884440, 0.738416812340510, 0.232102326775050};
366 9.71357962827961025E-002, 3.13347002271398278E-002,
367 3.13347002271398278E-002, 3.13347002271398278E-002,
368 7.78275410047754301E-002, 7.78275410047754301E-002,
369 7.78275410047754301E-002, 7.96477389272090969E-002,
370 7.96477389272090969E-002, 7.96477389272090969E-002,
371 2.55776756586981006E-002, 2.55776756586981006E-002,
372 2.55776756586981006E-002, 4.32835393772893970E-002,
373 4.32835393772893970E-002, 4.32835393772893970E-002,
374 4.32835393772893970E-002, 4.32835393772893970E-002,
375 4.32835393772893970E-002};
378 0.333333333333333, 0.948021718143423, 0.025989140928288,
379 0.025989140928288, 0.811424994704155, 0.094287502647923,
380 0.094287502647923, 0.010726449965571, 0.494636775017215,
381 0.494636775017215, 0.585313234770972, 0.207343382614514,
382 0.207343382614514, 0.122184388599019, 0.438907805700491,
383 0.438907805700491, 0.677937654882590, 0.677937654882590,
384 0.044841677589131, 0.044841677589131, 0.277220667528279,
385 0.277220667528279, 0.858870281282636, 0.858870281282636,
386 0.0000000000000000, 0.0000000000000000, 0.141129718717364,
389 0.333333333333333, 0.025989140928288, 0.948021718143423, 0.025989140928288,
390 0.094287502647923, 0.811424994704155, 0.094287502647923, 0.494636775017215,
391 0.010726449965571, 0.494636775017215, 0.207343382614514, 0.585313234770972,
392 0.207343382614514, 0.438907805700491, 0.122184388599019, 0.438907805700491,
393 0.044841677589131, 0.277220667528279, 0.677937654882590, 0.277220667528279,
394 0.677937654882590, 0.044841677589131, 0.000000000000000, 0.141129718717364,
395 0.858870281282636, 0.141129718717364, 0.858870281282636, 0.000000000000000};
397 0.08797730116222190, 0.008744311553736190, 0.008744311553736190,
398 0.008744311553736190, 0.03808157199393533, 0.03808157199393533,
399 0.03808157199393533, 0.01885544805613125, 0.01885544805613125,
400 0.01885544805613125, 0.07215969754474100, 0.07215969754474100,
401 0.07215969754474100, 0.06932913870553720, 0.06932913870553720,
402 0.06932913870553720, 0.04105631542928860, 0.04105631542928860,
403 0.04105631542928860, 0.04105631542928860, 0.04105631542928860,
404 0.04105631542928860, 0.007362383783300573, 0.007362383783300573,
405 0.007362383783300573, 0.007362383783300573, 0.007362383783300573,
406 0.007362383783300573};
409 0.333333333333333, 0.950275662924106, 0.024862168537947, 0.024862168537947,
410 0.171614914923835, 0.414192542538082, 0.414192542538082, 0.539412243677190,
411 0.230293878161405, 0.230293878161405, 0.772160036676533, 0.113919981661734,
412 0.113919981661734, 0.009085399949835, 0.495457300025082, 0.495457300025082,
413 0.062277290305887, 0.468861354847056, 0.468861354847056, 0.022076289653624,
414 0.022076289653624, 0.851306504174348, 0.851306504174348, 0.126617206172027,
415 0.126617206172027, 0.018620522802521, 0.018620522802521, 0.689441970728591,
416 0.689441970728591, 0.291937506468888, 0.291937506468888, 0.096506481292159,
417 0.096506481292159, 0.635867859433873, 0.635867859433873, 0.267625659273968,
420 0.333333333333333, 0.024862168537947, 0.950275662924106, 0.024862168537947,
421 0.414192542538082, 0.171614914923835, 0.414192542538082, 0.230293878161405,
422 0.539412243677190, 0.230293878161405, 0.113919981661734, 0.772160036676533,
423 0.113919981661734, 0.495457300025082, 0.009085399949835, 0.495457300025082,
424 0.468861354847056, 0.062277290305887, 0.468861354847056, 0.851306504174348,
425 0.126617206172027, 0.022076289653624, 0.126617206172027, 0.022076289653624,
426 0.851306504174348, 0.689441970728591, 0.291937506468888, 0.018620522802521,
427 0.291937506468888, 0.018620522802521, 0.689441970728591, 0.635867859433873,
428 0.267625659273968, 0.096506481292159, 0.267625659273968, 0.096506481292159,
431 0.051739766065744, 0.008007799555565, 0.008007799555565, 0.008007799555565,
432 0.046868898981822, 0.046868898981822, 0.046868898981822, 0.046590940183976,
433 0.046590940183976, 0.046590940183976, 0.031016943313796, 0.031016943313796,
434 0.031016943313796, 0.010791612736631, 0.010791612736631, 0.010791612736631,
435 0.032195534242432, 0.032195534242432, 0.032195534242432, 0.015445834210702,
436 0.015445834210702, 0.015445834210702, 0.015445834210702, 0.015445834210702,
437 0.015445834210702, 0.017822989923179, 0.017822989923179, 0.017822989923179,
438 0.017822989923179, 0.017822989923179, 0.017822989923179, 0.037038683681385,
439 0.037038683681385, 0.037038683681385, 0.037038683681385, 0.037038683681385,
1012 0.3333333333333333};
1014 2.912193380035668, 2.912193380035668, 2.912193380035668,
1015 2.912193380035668, 2.912193380035668, 2.912193380035668,
1016 2.912193380035668, 2.912193380035668, 2.912193380035668,
1017 2.912193380035668, 2.912193380035668, 2.912193380035668,
1018 2.912193380035668, 2.912193380035668, 2.912193380035668,
1019 2.912193380035668, 2.912193380035668, 2.912193380035668,
1020 2.912193380035668, 2.912193380035668, 2.912193380035668,
1021 2.912193380035668, 2.912193380035668, 2.912193380035668,
1022 2.912193380035668, 2.912193380035668, 2.912193380035668,
1023 2.912193380035668, 2.912193380035668, 2.912193380035668,
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 -9.914451197589852, -9.914451197589852, -9.914451197589852,
1037 -9.914451197589852, -9.914451197589852, -9.914451197589852,
1038 -9.914451197589852, -9.914451197589852, -9.914451197589852,
1039 -9.914451197589852, -9.914451197589852, -9.914451197589852,
1040 -9.914451197589852, -9.914451197589852, -9.914451197589852,
1041 -9.914451197589852, -9.914451197589852, -9.914451197589852,
1042 -9.914451197589852, -9.914451197589852, -9.914451197589852,
1043 -9.914451197589852, -9.914451197589852, -9.914451197589852,
1044 -9.914451197589852, -9.914451197589852, -9.914451197589852,
1045 -9.914451197589852, -9.914451197589852, -9.914451197589852,
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, 13.33158527957992, 13.33158527957992,
1055 13.33158527957992, 13.33158527957992, 13.33158527957992,
1056 13.33158527957992, 13.33158527957992, 13.33158527957992,
1057 13.33158527957992, 13.33158527957992, 13.33158527957992,
1058 13.33158527957992, 13.33158527957992, 13.33158527957992,
1059 13.33158527957992, 13.33158527957992, 13.33158527957992,
1060 13.33158527957992, 13.33158527957992, 13.33158527957992,
1061 13.33158527957992, 13.33158527957992, 13.33158527957992,
1062 13.33158527957992, 13.33158527957992, 13.33158527957992,
1063 13.33158527957992, 13.33158527957992, 13.33158527957992,
1064 13.33158527957992, 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, -9.027792408986382, -9.027792408986382,
1070 -9.027792408986382, -9.027792408986382, -9.027792408986382,
1071 -9.027792408986382, -9.027792408986382, -9.027792408986382,
1072 -9.027792408986382, -9.027792408986382, -9.027792408986382,
1073 -9.027792408986382, -9.027792408986382, -9.027792408986382,
1074 -9.027792408986382, -9.027792408986382, -9.027792408986382,
1075 -9.027792408986382, -9.027792408986382, -9.027792408986382,
1076 -9.027792408986382, -9.027792408986382, -9.027792408986382,
1077 -9.027792408986382, -9.027792408986382, -9.027792408986382,
1078 -9.027792408986382, -9.027792408986382, -9.027792408986382,
1079 -9.027792408986382, -9.027792408986382, -9.027792408986382,
1080 -9.027792408986382, -9.027792408986382, -9.027792408986382,
1081 -9.027792408986382, 3.258672079964582, 3.258672079964582,
1082 3.258672079964582, 3.258672079964582, 3.258672079964582,
1083 3.258672079964582, 3.258672079964582, 3.258672079964582,
1084 3.258672079964582, 3.258672079964582, 3.258672079964582,
1085 3.258672079964582, 3.258672079964582, 3.258672079964582,
1086 3.258672079964582, 3.258672079964582, 3.258672079964582,
1087 3.258672079964582, 3.258672079964582, 3.258672079964582,
1088 3.258672079964582, 3.258672079964582, 3.258672079964582,
1089 3.258672079964582, 3.258672079964582, 3.258672079964582,
1090 3.258672079964582, 3.258672079964582, -0.6133639040302452,
1091 -0.6133639040302452, -0.6133639040302452, -0.6133639040302452,
1092 -0.6133639040302452, -0.6133639040302452, -0.6133639040302452,
1093 -0.6133639040302452, -0.6133639040302452, -0.6133639040302452,
1094 -0.6133639040302452, -0.6133639040302452, -0.6133639040302452,
1095 -0.6133639040302452, -0.6133639040302452, -0.6133639040302452,
1096 -0.6133639040302452, -0.6133639040302452, -0.6133639040302452,
1097 -0.6133639040302452, -0.6133639040302452, 0.05511571669513555,
1098 0.05511571669513555, 0.05511571669513555, 0.05511571669513555,
1099 0.05511571669513555, 0.05511571669513555, 0.05511571669513555,
1100 0.05511571669513555, 0.05511571669513555, 0.05511571669513555,
1101 0.05511571669513555, 0.05511571669513555, 0.05511571669513555,
1102 0.05511571669513555, 0.05511571669513555, -0.001979122382447095,
1103 -0.001979122382447095, -0.001979122382447095, -0.001979122382447095,
1104 -0.001979122382447095, -0.001979122382447095, -0.001979122382447095,
1105 -0.001979122382447095, -0.001979122382447095, -0.001979122382447095,
1106 2.02054621415273e-05, 2.02054621415273e-05, 2.02054621415273e-05,
1107 2.02054621415273e-05, 2.02054621415273e-05, 2.02054621415273e-05,
1108 -2.874940020535803e-08, -2.874940020535803e-08, -2.874940020535803e-08,
1109 8.829438425435718e-13};
1116static const double G_TET_X4[] = {0.1757281246520584, 0.2445310270213291,
1117 0.5556470949048655, 0.0240937534217468};
1118static const double G_TET_Y4[] = {0.5656137776620919, 0.0501800797762026,
1119 0.1487681308666864, 0.2354380116950194};
1120static const double G_TET_Z4[] = {0.2180665126782654, 0.5635595064952189,
1121 0.0350112499848832, 0.1833627308416330};
1122static const double G_TET_W4[] = {0.25, 0.25, 0.25, 0.25};
1123static const double G_TET_X5[] = {0.25000000000000000, 0.50000000000000000,
1124 0.16666666666666667, 0.16666666666666667,
1125 0.16666666666666667};
1126static const double G_TET_Y5[] = {0.25000000000000000, 0.16666666666666667,
1127 0.50000000000000000, 0.16666666666666667,
1128 0.16666666666666667};
1129static const double G_TET_Z5[] = {0.25000000000000000, 0.16666666666666667,
1130 0.16666666666666667, 0.50000000000000000,
1131 0.16666666666666667};
1132static const double G_TET_W5[] = {-0.80000000000000000, 0.45000000000000000,
1133 0.45000000000000000, 0.45000000000000000,
1134 0.45000000000000000};
1135static const double G_TET_X10[] = {0.5684305841968444, 0.1438564719343852,
1136 0.1438564719343852, 0.1438564719343852,
1137 0.0000000000000000, 0.5000000000000000,
1138 0.5000000000000000, 0.5000000000000000,
1139 0.0000000000000000, 0.0000000000000000};
1140static const double G_TET_Y10[] = {0.1438564719343852, 0.1438564719343852,
1141 0.1438564719343852, 0.5684305841968444,
1142 0.5000000000000000, 0.0000000000000000,
1143 0.5000000000000000, 0.0000000000000000,
1144 0.5000000000000000, 0.0000000000000000};
1145static const double G_TET_Z10[] = {0.1438564719343852, 0.1438564719343852,
1146 0.5684305841968444, 0.1438564719343852,
1147 0.5000000000000000, 0.5000000000000000,
1148 0.0000000000000000, 0.0000000000000000,
1149 0.0000000000000000, 0.5000000000000000};
1150static const double G_TET_W10[] = {0.2177650698804054, 0.2177650698804054,
1151 0.2177650698804054, 0.2177650698804054,
1152 0.0214899534130631, 0.0214899534130631,
1153 0.0214899534130631, 0.0214899534130631,
1154 0.0214899534130631, 0.0214899534130631};
1157 0.2500000000000000, 0.6175871903000830, 0.1274709365666390,
1158 0.1274709365666390, 0.1274709365666390, 0.9037635088221031,
1159 0.0320788303926323, 0.0320788303926323, 0.0320788303926323,
1160 0.4502229043567190, 0.0497770956432810, 0.0497770956432810,
1161 0.0497770956432810, 0.4502229043567190, 0.4502229043567190,
1162 0.3162695526014501, 0.1837304473985499, 0.1837304473985499,
1163 0.1837304473985499, 0.3162695526014501, 0.3162695526014501,
1164 0.0229177878448171, 0.2319010893971509, 0.2319010893971509,
1165 0.5132800333608811, 0.2319010893971509, 0.2319010893971509,
1166 0.2319010893971509, 0.0229177878448171, 0.5132800333608811,
1167 0.2319010893971509, 0.0229177878448171, 0.5132800333608811,
1168 0.7303134278075384, 0.0379700484718286, 0.0379700484718286,
1169 0.1937464752488044, 0.0379700484718286, 0.0379700484718286,
1170 0.0379700484718286, 0.7303134278075384, 0.1937464752488044,
1171 0.0379700484718286, 0.7303134278075384, 0.1937464752488044};
1173 0.2500000000000000, 0.1274709365666390, 0.1274709365666390,
1174 0.1274709365666390, 0.6175871903000830, 0.0320788303926323,
1175 0.0320788303926323, 0.0320788303926323, 0.9037635088221031,
1176 0.0497770956432810, 0.4502229043567190, 0.0497770956432810,
1177 0.4502229043567190, 0.0497770956432810, 0.4502229043567190,
1178 0.1837304473985499, 0.3162695526014501, 0.1837304473985499,
1179 0.3162695526014501, 0.1837304473985499, 0.3162695526014501,
1180 0.2319010893971509, 0.0229177878448171, 0.2319010893971509,
1181 0.2319010893971509, 0.5132800333608811, 0.2319010893971509,
1182 0.0229177878448171, 0.5132800333608811, 0.2319010893971509,
1183 0.5132800333608811, 0.2319010893971509, 0.0229177878448171,
1184 0.0379700484718286, 0.7303134278075384, 0.0379700484718286,
1185 0.0379700484718286, 0.1937464752488044, 0.0379700484718286,
1186 0.7303134278075384, 0.1937464752488044, 0.0379700484718286,
1187 0.1937464752488044, 0.0379700484718286, 0.7303134278075384};
1189 0.2500000000000000, 0.1274709365666390, 0.1274709365666390,
1190 0.6175871903000830, 0.1274709365666390, 0.0320788303926323,
1191 0.0320788303926323, 0.9037635088221031, 0.0320788303926323,
1192 0.0497770956432810, 0.0497770956432810, 0.4502229043567190,
1193 0.4502229043567190, 0.4502229043567190, 0.0497770956432810,
1194 0.1837304473985499, 0.1837304473985499, 0.3162695526014501,
1195 0.3162695526014501, 0.3162695526014501, 0.1837304473985499,
1196 0.2319010893971509, 0.2319010893971509, 0.0229177878448171,
1197 0.2319010893971509, 0.2319010893971509, 0.5132800333608811,
1198 0.5132800333608811, 0.2319010893971509, 0.0229177878448171,
1199 0.0229177878448171, 0.5132800333608811, 0.2319010893971509,
1200 0.0379700484718286, 0.0379700484718286, 0.7303134278075384,
1201 0.0379700484718286, 0.0379700484718286, 0.1937464752488044,
1202 0.1937464752488044, 0.0379700484718286, 0.7303134278075384,
1203 0.7303134278075384, 0.1937464752488044, 0.0379700484718286};
1205 -0.2359620398477557, 0.0244878963560562, 0.0244878963560562,
1206 0.0244878963560562, 0.0244878963560562, 0.0039485206398261,
1207 0.0039485206398261, 0.0039485206398261, 0.0039485206398261,
1208 0.0263055529507371, 0.0263055529507371, 0.0263055529507371,
1209 0.0263055529507371, 0.0263055529507371, 0.0263055529507371,
1210 0.0829803830550589, 0.0829803830550589, 0.0829803830550589,
1211 0.0829803830550589, 0.0829803830550589, 0.0829803830550589,
1212 0.0254426245481023, 0.0254426245481023, 0.0254426245481023,
1213 0.0254426245481023, 0.0254426245481023, 0.0254426245481023,
1214 0.0254426245481023, 0.0254426245481023, 0.0254426245481023,
1215 0.0254426245481023, 0.0254426245481023, 0.0254426245481023,
1216 0.0134324384376852, 0.0134324384376852, 0.0134324384376852,
1217 0.0134324384376852, 0.0134324384376852, 0.0134324384376852,
1218 0.0134324384376852, 0.0134324384376852, 0.0134324384376852,
1219 0.0134324384376852, 0.0134324384376852, 0.0134324384376852};
1221 0.1000000000000000, 0.1000000000000000, 0.1000000000000000,
1222 0.7000000000000000, 0.1000000000000000, 0.1000000000000000,
1223 0.1000000000000000, 0.1000000000000000, 0.1000000000000000,
1224 0.1000000000000000, 0.6000000000000000, 0.6000000000000000,
1225 0.6000000000000000, 0.2000000000000000, 0.2000000000000000,
1226 0.2000000000000000, 0.1000000000000000, 0.1000000000000000,
1227 0.1000000000000000, 0.1000000000000000, 0.1000000000000000,
1228 0.1000000000000000, 0.5000000000000000, 0.5000000000000000,
1229 0.5000000000000000, 0.3000000000000000, 0.3000000000000000,
1230 0.3000000000000000, 0.2000000000000000, 0.2000000000000000,
1231 0.2000000000000000, 0.2000000000000000, 0.2000000000000000,
1232 0.2000000000000000, 0.1000000000000000, 0.1000000000000000,
1233 0.1000000000000000, 0.5000000000000000, 0.5000000000000000,
1234 0.5000000000000000, 0.4000000000000000, 0.4000000000000000,
1235 0.4000000000000000, 0.1000000000000000, 0.1000000000000000,
1236 0.1000000000000000, 0.4000000000000000, 0.4000000000000000,
1237 0.4000000000000000, 0.4000000000000000, 0.4000000000000000,
1238 0.4000000000000000, 0.3000000000000000, 0.3000000000000000,
1239 0.3000000000000000, 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.1000000000000000, 0.1000000000000000,
1244 0.1000000000000000, 0.2000000000000000, 0.2000000000000000,
1245 0.2000000000000000, 0.4000000000000000, 0.3000000000000000,
1246 0.3000000000000000, 0.3000000000000000, 0.1000000000000000,
1247 0.3000000000000000, 0.3000000000000000, 0.3000000000000000,
1248 0.2000000000000000, 0.2000000000000000, 0.2000000000000000};
1250 0.1000000000000000, 0.1000000000000000, 0.7000000000000000,
1251 0.1000000000000000, 0.1000000000000000, 0.1000000000000000,
1252 0.6000000000000000, 0.6000000000000000, 0.2000000000000000,
1253 0.2000000000000000, 0.1000000000000000, 0.1000000000000000,
1254 0.2000000000000000, 0.1000000000000000, 0.1000000000000000,
1255 0.6000000000000000, 0.1000000000000000, 0.1000000000000000,
1256 0.5000000000000000, 0.5000000000000000, 0.3000000000000000,
1257 0.3000000000000000, 0.1000000000000000, 0.1000000000000000,
1258 0.3000000000000000, 0.1000000000000000, 0.1000000000000000,
1259 0.5000000000000000, 0.2000000000000000, 0.2000000000000000,
1260 0.1000000000000000, 0.1000000000000000, 0.5000000000000000,
1261 0.5000000000000000, 0.2000000000000000, 0.2000000000000000,
1262 0.5000000000000000, 0.2000000000000000, 0.2000000000000000,
1263 0.1000000000000000, 0.4000000000000000, 0.1000000000000000,
1264 0.1000000000000000, 0.4000000000000000, 0.4000000000000000,
1265 0.1000000000000000, 0.3000000000000000, 0.3000000000000000,
1266 0.2000000000000000, 0.2000000000000000, 0.1000000000000000,
1267 0.1000000000000000, 0.4000000000000000, 0.4000000000000000,
1268 0.2000000000000000, 0.2000000000000000, 0.1000000000000000,
1269 0.1000000000000000, 0.4000000000000000, 0.4000000000000000,
1270 0.3000000000000000, 0.3000000000000000, 0.1000000000000000,
1271 0.1000000000000000, 0.4000000000000000, 0.4000000000000000,
1272 0.3000000000000000, 0.3000000000000000, 0.2000000000000000,
1273 0.2000000000000000, 0.2000000000000000, 0.2000000000000000,
1274 0.4000000000000000, 0.2000000000000000, 0.3000000000000000,
1275 0.3000000000000000, 0.1000000000000000, 0.3000000000000000,
1276 0.3000000000000000, 0.2000000000000000, 0.2000000000000000,
1277 0.3000000000000000, 0.3000000000000000, 0.2000000000000000};
1279 0.1000000000000000, 0.7000000000000000, 0.1000000000000000,
1280 0.1000000000000000, 0.6000000000000000, 0.2000000000000000,
1281 0.1000000000000000, 0.2000000000000000, 0.1000000000000000,
1282 0.6000000000000000, 0.1000000000000000, 0.2000000000000000,
1283 0.1000000000000000, 0.1000000000000000, 0.6000000000000000,
1284 0.1000000000000000, 0.5000000000000000, 0.3000000000000000,
1285 0.1000000000000000, 0.3000000000000000, 0.1000000000000000,
1286 0.5000000000000000, 0.1000000000000000, 0.3000000000000000,
1287 0.1000000000000000, 0.1000000000000000, 0.5000000000000000,
1288 0.1000000000000000, 0.1000000000000000, 0.5000000000000000,
1289 0.2000000000000000, 0.5000000000000000, 0.2000000000000000,
1290 0.1000000000000000, 0.2000000000000000, 0.5000000000000000,
1291 0.2000000000000000, 0.2000000000000000, 0.1000000000000000,
1292 0.2000000000000000, 0.1000000000000000, 0.4000000000000000,
1293 0.1000000000000000, 0.4000000000000000, 0.1000000000000000,
1294 0.4000000000000000, 0.2000000000000000, 0.1000000000000000,
1295 0.3000000000000000, 0.1000000000000000, 0.3000000000000000,
1296 0.2000000000000000, 0.2000000000000000, 0.1000000000000000,
1297 0.4000000000000000, 0.1000000000000000, 0.4000000000000000,
1298 0.2000000000000000, 0.3000000000000000, 0.1000000000000000,
1299 0.4000000000000000, 0.1000000000000000, 0.4000000000000000,
1300 0.3000000000000000, 0.3000000000000000, 0.2000000000000000,
1301 0.4000000000000000, 0.2000000000000000, 0.4000000000000000,
1302 0.3000000000000000, 0.2000000000000000, 0.4000000000000000,
1303 0.2000000000000000, 0.2000000000000000, 0.3000000000000000,
1304 0.1000000000000000, 0.3000000000000000, 0.3000000000000000,
1305 0.2000000000000000, 0.3000000000000000, 0.2000000000000000,
1306 0.3000000000000000, 0.2000000000000000, 0.3000000000000000};
1308 0.2843915343915344, 0.2843915343915344, 0.2843915343915344,
1309 0.2843915343915344, -0.3882275132275133, -0.3882275132275133,
1310 -0.3882275132275133, -0.3882275132275133, -0.3882275132275133,
1311 -0.3882275132275133, -0.3882275132275133, -0.3882275132275133,
1312 -0.3882275132275133, -0.3882275132275133, -0.3882275132275133,
1313 -0.3882275132275133, 0.8776455026455027, 0.8776455026455027,
1314 0.8776455026455027, 0.8776455026455027, 0.8776455026455027,
1315 0.8776455026455027, 0.8776455026455027, 0.8776455026455027,
1316 0.8776455026455027, 0.8776455026455027, 0.8776455026455027,
1317 0.8776455026455027, 0.1236772486772487, 0.1236772486772487,
1318 0.1236772486772487, 0.1236772486772487, 0.1236772486772487,
1319 0.1236772486772487, 0.1236772486772487, 0.1236772486772487,
1320 0.1236772486772487, 0.1236772486772487, 0.1236772486772487,
1321 0.1236772486772487, -0.8584656084656085, -0.8584656084656085,
1322 -0.8584656084656085, -0.8584656084656085, -0.8584656084656085,
1323 -0.8584656084656085, -0.2632275132275133, -0.2632275132275133,
1324 -0.2632275132275133, -0.2632275132275133, -0.2632275132275133,
1325 -0.2632275132275133, -0.2632275132275133, -0.2632275132275133,
1326 -0.2632275132275133, -0.2632275132275133, -0.2632275132275133,
1327 -0.2632275132275133, -0.2632275132275133, -0.2632275132275133,
1328 -0.2632275132275133, -0.2632275132275133, -0.2632275132275133,
1329 -0.2632275132275133, -0.2632275132275133, -0.2632275132275133,
1330 -0.2632275132275133, -0.2632275132275133, -0.2632275132275133,
1331 -0.2632275132275133, 0.0145502645502645, 0.0145502645502645,
1332 0.0145502645502645, 0.0145502645502645, 1.0165343915343916,
1333 1.0165343915343916, 1.0165343915343916, 1.0165343915343916,
1334 -0.0251322751322751, -0.0251322751322751, -0.0251322751322751,
1335 -0.0251322751322751, -0.0251322751322751, -0.0251322751322751};
FTensor::Index< 'n', SPACE_DIM > n
FTensor::Index< 'm', SPACE_DIM > m