ParMooN
 All Classes Functions Variables Friends Pages
BF_N_H_Q3_3D.h
1 // ***********************************************************************
2 // Q3Rot element, nonconforming, 3D
3 // ***********************************************************************
4 
5 static void N_H_Q3_3D_Funct(double xi, double eta, double zeta,
6  double *values)
7 {
8  double t1, t2, t3, t4, t5, t6, t7, t8, t9, t10, t11, t12, t13, t16, t18;
9  double t19, t20, t21, t24, t26, t28, t29, t30, t31, t32, t35, t37, t39;
10  double t40, t41, t43, t44, t46, t47, t48, t50, t51, t52, t53, t58, t60;
11  double t62, t63, t65, t67, t68, t69, t70, t71, t72, t74, t75, t76, t78;
12  double t79, t80, t81, t82, t83, t84, t86, t88, t89, t91, t93, t95, t96;
13  double t97, t98, t101, t102, t103, t104, t105, t106, t107, t108, t109;
14  double t110, t111, t113, t114, t115, t116, t119, t120, t121, t133, t134;
15  double t135, t136, t137, t138, t139, t141, t143, t144, t145, t147, t148;
16 
17  t1 = zeta*zeta;
18  t2 = 3.0/4.0*t1;
19  t3 = eta*eta;
20  t4 = 3.0/2.0*t3;
21  t5 = -1.0/2.0+t4;
22  t6 = 3.0/2.0*t1;
23  t7 = -1.0/2.0+t6;
24  t8 = t5*t7;
25  t9 = t8/4.0;
26  t10 = t1*zeta;
27  t11 = 5.0/4.0*t10;
28  t12 = xi*xi;
29  t13 = t12*t12;
30  t16 = 3.0/8.0+35.0/8.0*t13-15.0/4.0*t12;
31  t18 = t16*t7/4.0;
32  t19 = 3.0/2.0*t12;
33  t20 = -1.0/2.0+t19;
34  t21 = t1*t1;
35  t24 = 3.0/8.0+35.0/8.0*t21-15.0/4.0*t1;
36  t26 = t20*t24/4.0;
37  t28 = t5*t24/4.0;
38  t29 = t20*t7;
39  t30 = t29/4.0;
40  t31 = 3.0/4.0*zeta;
41  t32 = t3*t3;
42  t35 = 3.0/8.0+35.0/8.0*t32-15.0/4.0*t3;
43  t37 = t35*t7/4.0;
44  t39 = t3*eta;
45  t40 = 5.0/4.0*t39;
46  t41 = 3.0/4.0*eta;
47  t43 = t16*t5/4.0;
48  t44 = 3.0/4.0*t3;
49  t46 = t20*t35/4.0;
50  t47 = t20*t5;
51  t48 = t47/4.0;
52  t50 = 3.0/4.0*xi;
53  t51 = t12*xi;
54  t52 = 5.0/4.0*t51;
55  t53 = 3.0/4.0*t12;
56  t58 = 5.0/2.0*t51;
57  t60 = t58-3.0/2.0*xi;
58  t62 = t60*zeta/4.0;
59  t63 = 5.0/2.0*t10;
60  t65 = t63-3.0/2.0*zeta;
61  t67 = xi*t65/4.0;
62  t68 = t60*t7;
63  t69 = t68/2.0;
64  t70 = xi*t5;
65  t71 = t70*zeta;
66  t72 = t71/4.0;
67  t74 = xi*zeta/4.0;
68  t75 = xi*t7;
69  t76 = t75/2.0;
70  t78 = t20*eta;
71  t79 = t78*zeta;
72  t80 = t79/4.0;
73  t81 = t5*t65;
74  t82 = t81/2.0;
75  t83 = t5*zeta;
76  t84 = t83/2.0;
77  t86 = eta*zeta/4.0;
78  t88 = eta*t65/4.0;
79  t89 = 5.0/2.0*t39;
80  t91 = t89-3.0/2.0*eta;
81  t93 = t91*zeta/4.0;
82  t95 = t20*t65;
83  t96 = t95/2.0;
84  t97 = t20*zeta;
85  t98 = t97/2.0;
86  t101 = t60*eta;
87  t102 = t101/4.0;
88  t103 = xi*t91;
89  t104 = t103/4.0;
90  t105 = xi*eta;
91  t106 = t105/4.0;
92  t107 = t105*t7;
93  t108 = t107/4.0;
94  t109 = t78/2.0;
95  t110 = t20*t91;
96  t111 = t110/2.0;
97  t113 = eta*t7;
98  t114 = t113/2.0;
99  t115 = t91*t7;
100  t116 = t115/2.0;
101  t119 = t60*t5;
102  t120 = t119/2.0;
103  t121 = t70/2.0;
104  t133 = t105*t65;
105  t134 = t133/3.0;
106  t135 = t107/2.0;
107  t136 = t101*zeta;
108  t137 = t136/6.0;
109  t138 = t103*zeta;
110  t139 = t138/6.0;
111  t141 = t105*zeta/6.0;
112  t143 = t133/6.0;
113  t144 = t71/2.0;
114  t145 = t138/3.0;
115  t147 = t79/2.0;
116  t148 = t136/3.0;
117 
118  values[0] = t2-t9-t11-t18+t26+t28-t30-1.0/4.0+t31-t37;
119  values[1] = -t40-t9+t41-t28-t43+t44+t46-t48-1.0/4.0+t37;
120  values[2] = -t50+t52+t18-t26+t43-t46-t30-t48-1.0/4.0+t53;
121  values[3] = t40-t9-t41-t28-t43+t44+t46-t48-1.0/4.0+t37;
122  values[4] = t50-t52+t18-t26+t43-t46-t30-t48-1.0/4.0+t53;
123  values[5] = t2-t9+t11-t18+t26+t28-t30-1.0/4.0-t31-t37;
124  values[6] = t62-t67-t69+t72-t74+t76;
125  values[7] = t80-t82+t84-t86+t88-t93;
126  values[8] = t62-t67-t96-t72+t74+t98;
127  values[9] = -t80-t82+t84+t86-t88+t93;
128  values[10] = -t102+t104-t106+t108+t109-t111;
129  values[11] = t114-t80-t116+t86+t88-t93;
130  values[12] = t114+t80-t116-t86-t88+t93;
131  values[13] = t102-t104-t106-t120+t108+t121;
132  values[14] = t102-t104+t106-t108+t109-t111;
133  values[15] = t102-t104-t106+t120+t108-t121;
134  values[16] = -t62+t67-t96+t72-t74+t98;
135  values[17] = -t62+t67-t69-t72+t74+t76;
136  values[18] = -t96-t18+t26+t30;
137  values[19] = t9-t28-t116+t37;
138  values[20] = t69+t18-t26+t30;
139  values[21] = t9-t28+t116+t37;
140  values[22] = -t120+t43-t46+t48;
141  values[23] = t9+t28+t82-t37;
142  values[24] = -t134+t135+t137+t139-t141;
143  values[25] = t143+t144+t137-t145-t141;
144  values[26] = -t143+t147+t148-t139+t141;
145  values[27] = t143-t144+t137-t145-t141;
146  values[28] = t143+t147-t148+t139-t141;
147  values[29] = t134+t135-t137-t139+t141;
148  values[30] = t9+t28-t82-t37;
149  values[31] = -t43+t46-t111+t48;
150  values[32] = t120+t43-t46+t48;
151  values[33] = -t43+t46+t111+t48;
152  values[34] = -t69+t18-t26+t30;
153  values[35] = t96-t18+t26+t30;
154  values[36] = -t6+t8+5.0/2.0-t4+t29+t47-t19;
155  values[37] = t119+5.0/2.0*xi-t58+t68-t75-t70;
156  values[38] = -t89-t113+5.0/2.0*eta+t115-t78+t110;
157  values[39] = t95-t63+t81-t97+5.0/2.0*zeta-t83;
158 }
159 
160 static void N_H_Q3_3D_DeriveXi(double xi, double eta, double zeta,
161  double *values)
162 {
163  double t1, t5, t6, t7, t8, t10, t11, t16, t17, t18, t19, t20, t21, t22;
164  double t24, t25, t30, t31, t32, t33, t34, t35, t38, t39, t41, t42, t43;
165  double t44, t45, t46, t47, t48, t49, t51, t52, t53, t54, t57, t58, t59;
166  double t60, t61, t63, t64, t65, t66, t67, t68, t69, t70, t73, t74, t75;
167  double t77, t78, t79, t80, t89, t90, t91, t92, t93, t94, t95, t97, t99;
168  double t100, t101, t103, t104;
169 
170  t1 = xi*xi;
171  t5 = 35.0/2.0*t1*xi-15.0/2.0*xi;
172  t6 = zeta*zeta;
173  t7 = 3.0/2.0*t6;
174  t8 = -1.0/2.0+t7;
175  t10 = t5*t8/4.0;
176  t11 = t6*t6;
177  t16 = 3.0/4.0*xi*(3.0/8.0+35.0/8.0*t11-15.0/4.0*t6);
178  t17 = xi*t8;
179  t18 = 3.0/4.0*t17;
180  t19 = -t10+t16-t18;
181  t20 = eta*eta;
182  t21 = 3.0/2.0*t20;
183  t22 = -1.0/2.0+t21;
184  t24 = t5*t22/4.0;
185  t25 = t20*t20;
186  t30 = 3.0/4.0*xi*(3.0/8.0+35.0/8.0*t25-15.0/4.0*t20);
187  t31 = xi*t22;
188  t32 = 3.0/4.0*t31;
189  t33 = -t24+t30-t32;
190  t34 = 15.0/4.0*t1;
191  t35 = 3.0/2.0*xi;
192  t38 = 15.0/2.0*t1;
193  t39 = t38-3.0/2.0;
194  t41 = t39*zeta/4.0;
195  t42 = t6*zeta;
196  t43 = 5.0/8.0*t42;
197  t44 = zeta/8.0;
198  t45 = t39*t8;
199  t46 = t45/2.0;
200  t47 = t22*zeta;
201  t48 = t47/4.0;
202  t49 = 3.0/4.0*t6;
203  t51 = xi*eta;
204  t52 = t51*zeta;
205  t53 = 3.0/4.0*t52;
206  t54 = 5.0/8.0*zeta;
207  t57 = 5.0/2.0*t42-3.0/2.0*zeta;
208  t58 = xi*t57;
209  t59 = 3.0/2.0*t58;
210  t60 = xi*zeta;
211  t61 = 3.0/2.0*t60;
212  t63 = t39*eta;
213  t64 = t63/4.0;
214  t65 = t20*eta;
215  t66 = 5.0/8.0*t65;
216  t67 = 5.0/8.0*eta;
217  t68 = eta*t8;
218  t69 = t68/4.0;
219  t70 = 3.0/2.0*t51;
220  t73 = 5.0/2.0*t65-3.0/2.0*eta;
221  t74 = xi*t73;
222  t75 = 3.0/2.0*t74;
223  t77 = eta/8.0;
224  t78 = t39*t22;
225  t79 = t78/2.0;
226  t80 = 3.0/4.0*t20;
227  t89 = eta*t57;
228  t90 = t89/3.0;
229  t91 = t68/2.0;
230  t92 = t63*zeta;
231  t93 = t92/6.0;
232  t94 = t73*zeta;
233  t95 = t94/6.0;
234  t97 = eta*zeta/6.0;
235  t99 = t89/6.0;
236  t100 = t47/2.0;
237  t101 = t94/3.0;
238  t103 = 3.0/2.0*t52;
239  t104 = t92/3.0;
240 
241  values[0] = t19;
242  values[1] = t33;
243  values[2] = -3.0/4.0+t34+t10-t16+t24-t30-t18-t32+t35;
244  values[3] = t33;
245  values[4] = 3.0/4.0-t34+t10-t16+t24-t30-t18-t32+t35;
246  values[5] = t19;
247  values[6] = t41-t43+t44-t46+t48-1.0/4.0+t49;
248  values[7] = t53;
249  values[8] = t41-t43+t54-t59-t48+t61;
250  values[9] = -t53;
251  values[10] = -t64+t66-t67+t69+t70-t75;
252  values[11] = -t53;
253  values[12] = t53;
254  values[13] = t64-t66+t77-t79+t69-1.0/4.0+t80;
255  values[14] = t64-t66+t67-t69+t70-t75;
256  values[15] = t64-t66+t77+t79+t69+1.0/4.0-t80;
257  values[16] = -t41+t43-t54-t59+t48+t61;
258  values[17] = -t41+t43-t44-t46-t48-1.0/4.0+t49;
259  values[18] = -t59-t10+t16+t18;
260  values[19] = 0.0;
261  values[20] = t46+t10-t16+t18;
262  values[21] = 0.0;
263  values[22] = -t79+t24-t30+t32;
264  values[23] = 0.0;
265  values[24] = -t90+t91+t93+t95-t97;
266  values[25] = t99+t100+t93-t101-t97;
267  values[26] = -t99+t103+t104-t95+t97;
268  values[27] = t99-t100+t93-t101-t97;
269  values[28] = t99+t103-t104+t95-t97;
270  values[29] = t90+t91-t93-t95+t97;
271  values[30] = 0.0;
272  values[31] = -t24+t30-t75+t32;
273  values[32] = t79+t24-t30+t32;
274  values[33] = -t24+t30+t75+t32;
275  values[34] = -t46+t10-t16+t18;
276  values[35] = t59-t10+t16+t18;
277  values[36] = 3.0*t17+3.0*t31-3.0*xi;
278  values[37] = t78+7.0/2.0-t38+t45-t7-t21;
279  values[38] = -3.0*t51+3.0*t74;
280  values[39] = 3.0*t58-3.0*t60;
281 }
282 
283 static void N_H_Q3_3D_DeriveEta(double xi, double eta, double zeta,
284  double *values)
285 {
286  double t1, t2, t3, t4, t5, t6, t11, t12, t16, t18, t19, t20, t21, t22;
287  double t27, t28, t29, t30, t32, t33, t34, t36, t38, t39, t40, t41, t42;
288  double t43, t46, t47, t48, t49, t50, t51, t52, t53, t54, t56, t59, t60;
289  double t61, t62, t63, t64, t65, t66, t67, t68, t70, t71, t72, t73, t76;
290  double t79, t80, t81, t82, t90, t91, t92, t93, t94, t95, t96, t98, t100;
291  double t101, t102, t104, t105;
292 
293  t1 = zeta*zeta;
294  t2 = 3.0/2.0*t1;
295  t3 = -1.0/2.0+t2;
296  t4 = t3*eta;
297  t5 = 3.0/4.0*t4;
298  t6 = t1*t1;
299  t11 = 3.0/4.0*eta*(3.0/8.0+35.0/8.0*t6-15.0/4.0*t1);
300  t12 = eta*eta;
301  t16 = 35.0/2.0*t12*eta-15.0/2.0*eta;
302  t18 = t16*t3/4.0;
303  t19 = -t5+t11-t18;
304  t20 = 15.0/4.0*t12;
305  t21 = xi*xi;
306  t22 = t21*t21;
307  t27 = 3.0/4.0*(3.0/8.0+35.0/8.0*t22-15.0/4.0*t21)*eta;
308  t28 = 3.0/2.0*eta;
309  t29 = 3.0/2.0*t21;
310  t30 = -1.0/2.0+t29;
311  t32 = t30*t16/4.0;
312  t33 = t30*eta;
313  t34 = 3.0/4.0*t33;
314  t36 = t27-t32-t34;
315  t38 = xi*eta;
316  t39 = t38*zeta;
317  t40 = 3.0/4.0*t39;
318  t41 = t30*zeta;
319  t42 = t41/4.0;
320  t43 = t1*zeta;
321  t46 = 5.0/2.0*t43-3.0/2.0*zeta;
322  t47 = eta*t46;
323  t48 = 3.0/2.0*t47;
324  t49 = eta*zeta;
325  t50 = 3.0/2.0*t49;
326  t51 = 5.0/8.0*zeta;
327  t52 = 5.0/8.0*t43;
328  t53 = 15.0/2.0*t12;
329  t54 = t53-3.0/2.0;
330  t56 = t54*zeta/4.0;
331  t59 = t21*xi;
332  t60 = 5.0/8.0*t59;
333  t61 = xi/8.0;
334  t62 = xi*t54;
335  t63 = t62/4.0;
336  t64 = xi*t3;
337  t65 = t64/4.0;
338  t66 = 3.0/4.0*t21;
339  t67 = t30*t54;
340  t68 = t67/2.0;
341  t70 = 3.0/4.0*t1;
342  t71 = t54*t3;
343  t72 = t71/2.0;
344  t73 = zeta/8.0;
345  t76 = 5.0/8.0*xi;
346  t79 = 5.0/2.0*t59-3.0/2.0*xi;
347  t80 = t79*eta;
348  t81 = 3.0/2.0*t80;
349  t82 = 3.0/2.0*t38;
350  t90 = xi*t46;
351  t91 = t90/3.0;
352  t92 = t64/2.0;
353  t93 = t79*zeta;
354  t94 = t93/6.0;
355  t95 = t62*zeta;
356  t96 = t95/6.0;
357  t98 = xi*zeta/6.0;
358  t100 = t90/6.0;
359  t101 = 3.0/2.0*t39;
360  t102 = t95/3.0;
361  t104 = t41/2.0;
362  t105 = t93/3.0;
363 
364  values[0] = t19;
365  values[1] = -t20-t5+3.0/4.0-t11-t27+t28+t32-t34+t18;
366  values[2] = t36;
367  values[3] = t20-t5-3.0/4.0-t11-t27+t28+t32-t34+t18;
368  values[4] = t36;
369  values[5] = t19;
370  values[6] = t40;
371  values[7] = t42-t48+t50-t51+t52-t56;
372  values[8] = -t40;
373  values[9] = -t42-t48+t50+t51-t52+t56;
374  values[10] = -t60+t61+t63+t65-1.0/4.0+t66-t68;
375  values[11] = -1.0/4.0+t70-t42-t72-t73+t52-t56;
376  values[12] = -1.0/4.0+t70+t42-t72+t73-t52+t56;
377  values[13] = t60-t76-t63-t81+t65+t82;
378  values[14] = t60-t61-t63-t65-1.0/4.0+t66-t68;
379  values[15] = t60-t76-t63+t81+t65-t82;
380  values[16] = t40;
381  values[17] = -t40;
382  values[18] = 0.0;
383  values[19] = t5-t11-t72+t18;
384  values[20] = 0.0;
385  values[21] = t5-t11+t72+t18;
386  values[22] = -t81+t27-t32+t34;
387  values[23] = t5+t11+t48-t18;
388  values[24] = -t91+t92+t94+t96-t98;
389  values[25] = t100+t101+t94-t102-t98;
390  values[26] = -t100+t104+t105-t96+t98;
391  values[27] = t100-t101+t94-t102-t98;
392  values[28] = t100+t104-t105+t96-t98;
393  values[29] = t91+t92-t94-t96+t98;
394  values[30] = t5+t11-t48-t18;
395  values[31] = -t27+t32-t68+t34;
396  values[32] = t81+t27-t32+t34;
397  values[33] = -t27+t32+t68+t34;
398  values[34] = 0.0;
399  values[35] = 0.0;
400  values[36] = 3.0*t4-3.0*eta+3.0*t33;
401  values[37] = 3.0*t80-3.0*t38;
402  values[38] = -t53+7.0/2.0-t2+t71-t29+t67;
403  values[39] = 3.0*t47-3.0*t49;
404 }
405 
406 static void N_H_Q3_3D_DeriveZeta(double xi, double eta, double zeta,
407  double *values)
408 {
409  double t1, t2, t3, t4, t5, t6, t7, t8, t9, t10, t15, t16, t17, t21, t23;
410  double t25, t26, t27, t28, t33, t35, t36, t38, t39, t40, t41, t42, t44;
411  double t47, t48, t49, t50, t51, t52, t53, t55, t56, t57, t58, t59, t60;
412  double t62, t63, t64, t66, t67, t68, t69, t72, t73, t74, t75, t76, t79;
413  double t80, t81, t82, t92, t93, t94, t95, t96, t97, t98, t99, t101, t102;
414  double t103, t105, t106;
415 
416  t1 = 3.0/2.0*zeta;
417  t2 = eta*eta;
418  t3 = 3.0/2.0*t2;
419  t4 = -1.0/2.0+t3;
420  t5 = t4*zeta;
421  t6 = 3.0/4.0*t5;
422  t7 = zeta*zeta;
423  t8 = 15.0/4.0*t7;
424  t9 = xi*xi;
425  t10 = t9*t9;
426  t15 = 3.0/4.0*(3.0/8.0+35.0/8.0*t10-15.0/4.0*t9)*zeta;
427  t16 = 3.0/2.0*t9;
428  t17 = -1.0/2.0+t16;
429  t21 = 35.0/2.0*t7*zeta-15.0/2.0*zeta;
430  t23 = t17*t21/4.0;
431  t25 = t4*t21/4.0;
432  t26 = t17*zeta;
433  t27 = 3.0/4.0*t26;
434  t28 = t2*t2;
435  t33 = 3.0/4.0*(3.0/8.0+35.0/8.0*t28-15.0/4.0*t2)*zeta;
436  t35 = -t6-t25+t33;
437  t36 = t15-t23-t27;
438  t38 = t9*xi;
439  t39 = 5.0/8.0*t38;
440  t40 = 5.0/8.0*xi;
441  t41 = 15.0/2.0*t7;
442  t42 = t41-3.0/2.0;
443  t44 = xi*t42/4.0;
444  t47 = 5.0/2.0*t38-3.0/2.0*xi;
445  t48 = t47*zeta;
446  t49 = 3.0/2.0*t48;
447  t50 = xi*t4;
448  t51 = t50/4.0;
449  t52 = xi*zeta;
450  t53 = 3.0/2.0*t52;
451  t55 = t17*eta;
452  t56 = t55/4.0;
453  t57 = t4*t42;
454  t58 = t57/2.0;
455  t59 = 3.0/4.0*t2;
456  t60 = eta/8.0;
457  t62 = eta*t42/4.0;
458  t63 = t2*eta;
459  t64 = 5.0/8.0*t63;
460  t66 = xi/8.0;
461  t67 = t17*t42;
462  t68 = t67/2.0;
463  t69 = 3.0/4.0*t9;
464  t72 = xi*eta;
465  t73 = t72*zeta;
466  t74 = 3.0/4.0*t73;
467  t75 = eta*zeta;
468  t76 = 3.0/2.0*t75;
469  t79 = 5.0/2.0*t63-3.0/2.0*eta;
470  t80 = t79*zeta;
471  t81 = 3.0/2.0*t80;
472  t82 = 5.0/8.0*eta;
473  t92 = t72*t42;
474  t93 = t92/3.0;
475  t94 = 3.0/2.0*t73;
476  t95 = t47*eta;
477  t96 = t95/6.0;
478  t97 = xi*t79;
479  t98 = t97/6.0;
480  t99 = t72/6.0;
481  t101 = t92/6.0;
482  t102 = t50/2.0;
483  t103 = t97/3.0;
484  t105 = t55/2.0;
485  t106 = t95/3.0;
486 
487  values[0] = t1-t6-t8-t15+t23+t25-t27+3.0/4.0-t33;
488  values[1] = t35;
489  values[2] = t36;
490  values[3] = t35;
491  values[4] = t36;
492  values[5] = t1-t6+t8-t15+t23+t25-t27-3.0/4.0-t33;
493  values[6] = t39-t40-t44-t49+t51+t53;
494  values[7] = t56-t58-1.0/4.0+t59+t60+t62-t64;
495  values[8] = t39-t66-t44-t68-t51-1.0/4.0+t69;
496  values[9] = -t56-t58-1.0/4.0+t59-t60-t62+t64;
497  values[10] = t74;
498  values[11] = t76-t56-t81+t82+t62-t64;
499  values[12] = t76+t56-t81-t82-t62+t64;
500  values[13] = t74;
501  values[14] = -t74;
502  values[15] = t74;
503  values[16] = -t39+t66+t44-t68+t51-1.0/4.0+t69;
504  values[17] = -t39+t40+t44-t49-t51+t53;
505  values[18] = -t68-t15+t23+t27;
506  values[19] = t6-t25-t81+t33;
507  values[20] = t49+t15-t23+t27;
508  values[21] = t6-t25+t81+t33;
509  values[22] = 0.0;
510  values[23] = t6+t25+t58-t33;
511  values[24] = -t93+t94+t96+t98-t99;
512  values[25] = t101+t102+t96-t103-t99;
513  values[26] = -t101+t105+t106-t98+t99;
514  values[27] = t101-t102+t96-t103-t99;
515  values[28] = t101+t105-t106+t98-t99;
516  values[29] = t93+t94-t96-t98+t99;
517  values[30] = t6+t25-t58-t33;
518  values[31] = 0.0;
519  values[32] = 0.0;
520  values[33] = 0.0;
521  values[34] = -t49+t15-t23+t27;
522  values[35] = t68-t15+t23+t27;
523  values[36] = -3.0*zeta+3.0*t5+3.0*t26;
524  values[37] = 3.0*t48-3.0*t52;
525  values[38] = -3.0*t75+3.0*t80;
526  values[39] = t67-t41+t57+7.0/2.0-t16-t3;
527 }
528 
529 static void N_H_Q3_3D_DeriveXiXi(double xi, double eta, double zeta,
530  double *values)
531 {
532  double t1, t3, t4, t6, t8, t9, t10, t11, t12, t13, t15, t17, t18, t19;
533  double t20, t21, t22, t23, t24, t27, t28, t29, t30, t32, t33, t34, t36;
534  double t37, t39, t40, t41, t47, t48, t52, t53, t54, t55, t58, t59;
535 
536  t1 = xi*xi;
537  t3 = 105.0/2.0*t1-15.0/2.0;
538  t4 = zeta*zeta;
539  t6 = -1.0/2.0+3.0/2.0*t4;
540  t8 = t3*t6/4.0;
541  t9 = t4*t4;
542  t10 = 105.0/32.0*t9;
543  t11 = 63.0/16.0*t4;
544  t12 = -t8+21.0/32.0+t10-t11;
545  t13 = eta*eta;
546  t15 = -1.0/2.0+3.0/2.0*t13;
547  t17 = t3*t15/4.0;
548  t18 = t13*t13;
549  t19 = 105.0/32.0*t18;
550  t20 = 63.0/16.0*t13;
551  t21 = -t17+21.0/32.0+t19-t20;
552  t22 = 15.0/2.0*xi;
553  t23 = 27.0/16.0*t4;
554  t24 = 27.0/16.0*t13;
555  t27 = xi*zeta;
556  t28 = 15.0/4.0*t27;
557  t29 = xi*t6;
558  t30 = 15.0/2.0*t29;
559  t32 = eta*zeta;
560  t33 = 3.0/4.0*t32;
561  t34 = t4*zeta;
562  t36 = xi*eta;
563  t37 = t13*eta;
564  t39 = 15.0/4.0*t36;
565  t40 = xi*t15;
566  t41 = 15.0/2.0*t40;
567  t47 = 15.0/4.0*t34;
568  t48 = 9.0/4.0*zeta;
569  t52 = t36*zeta;
570  t53 = 5.0/2.0*t52;
571  t54 = 3.0/2.0*t32;
572  t55 = 5.0*t52;
573  t58 = 15.0/4.0*t37;
574  t59 = 9.0/4.0*eta;
575 
576  values[0] = t12;
577  values[1] = t21;
578  values[2] = t22+t8+27.0/16.0-t10+t23+t17-t19+t24;
579  values[3] = t21;
580  values[4] = -t22+t8+27.0/16.0-t10+t23+t17-t19+t24;
581  values[5] = t12;
582  values[6] = t28-t30;
583  values[7] = t33;
584  values[8] = 15.0/4.0*t27-15.0/4.0*t34+15.0/4.0*zeta;
585  values[9] = -t33;
586  values[10] = -15.0/4.0*t36+15.0/4.0*eta-15.0/4.0*t37;
587  values[11] = -t33;
588  values[12] = t33;
589  values[13] = t39-t41;
590  values[14] = 15.0/4.0*t36+15.0/4.0*eta-15.0/4.0*t37;
591  values[15] = t39+t41;
592  values[16] = -15.0/4.0*t27-15.0/4.0*t34+15.0/4.0*zeta;
593  values[17] = -t28-t30;
594  values[18] = -t47+t48-t8-3.0/32.0+t10-t23;
595  values[19] = 0.0;
596  values[20] = t30+t8-21.0/32.0-t10+t11;
597  values[21] = 0.0;
598  values[22] = -t41+t17-21.0/32.0-t19+t20;
599  values[23] = 0.0;
600  values[24] = t53;
601  values[25] = t53;
602  values[26] = t54+t55;
603  values[27] = t53;
604  values[28] = t54-t55;
605  values[29] = -t53;
606  values[30] = 0.0;
607  values[31] = -t17-3.0/32.0+t19-t24-t58+t59;
608  values[32] = t41+t17-21.0/32.0-t19+t20;
609  values[33] = -t17-3.0/32.0+t19-t24+t58-t59;
610  values[34] = -t30+t8-21.0/32.0-t10+t11;
611  values[35] = t47-t48-t8-3.0/32.0+t10-t23;
612  values[36] = -6.0+9.0/2.0*t4+9.0/2.0*t13;
613  values[37] = 15.0*t40-15.0*xi+15.0*t29;
614  values[38] = -15.0/2.0*eta+15.0/2.0*t37;
615  values[39] = 15.0/2.0*t34-15.0/2.0*zeta;
616 }
617 
618 static void N_H_Q3_3D_DeriveXiEta(double xi, double eta, double zeta,
619  double *values)
620 {
621  double t1, t7, t8, t14, t15, t16, t17, t18, t19, t20, t21, t22, t23, t24;
622  double t25, t26, t27, t29, t30, t31, t34, t35, t36, t37, t42, t43, t44;
623  double t45, t46, t47, t48, t49, t51, t52, t53, t54, t56, t57;
624 
625  t1 = xi*xi;
626  t7 = 3.0/4.0*(35.0/2.0*t1*xi-15.0/2.0*xi)*eta;
627  t8 = eta*eta;
628  t14 = 3.0/4.0*xi*(35.0/2.0*t8*eta-15.0/2.0*eta);
629  t15 = xi*eta;
630  t16 = 9.0/4.0*t15;
631  t17 = -t7+t14-t16;
632  t18 = t7-t14-t16;
633  t19 = eta*zeta;
634  t20 = 3.0/4.0*t19;
635  t21 = xi*zeta;
636  t22 = 3.0/4.0*t21;
637  t23 = 15.0/8.0*t1;
638  t24 = 15.0/8.0*t8;
639  t25 = zeta*zeta;
640  t26 = 3.0/8.0*t25;
641  t27 = 3.0/2.0*xi;
642  t29 = 15.0/2.0*t8-3.0/2.0;
643  t30 = xi*t29;
644  t31 = 3.0/2.0*t30;
645  t34 = 15.0/2.0*t1-3.0/2.0;
646  t35 = t34*eta;
647  t36 = 3.0/2.0*t35;
648  t37 = 3.0/2.0*eta;
649  t42 = t25*zeta;
650  t43 = 5.0/6.0*t42;
651  t44 = zeta/3.0;
652  t45 = 3.0/4.0*t25;
653  t46 = t34*zeta;
654  t47 = t46/6.0;
655  t48 = t29*zeta;
656  t49 = t48/6.0;
657  t51 = 5.0/12.0*t42;
658  t52 = 5.0/12.0*zeta;
659  t53 = 3.0/2.0*t19;
660  t54 = t48/3.0;
661  t56 = 3.0/2.0*t21;
662  t57 = t46/3.0;
663 
664  values[0] = 0.0;
665  values[1] = t17;
666  values[2] = t18;
667  values[3] = t17;
668  values[4] = t18;
669  values[5] = 0.0;
670  values[6] = t20;
671  values[7] = t22;
672  values[8] = -t20;
673  values[9] = -t22;
674  values[10] = -t23-3.0/8.0+t24+t26+t27-t31;
675  values[11] = -t22;
676  values[12] = t22;
677  values[13] = t23-3.0/8.0-t24-t36+t26+t37;
678  values[14] = t23+3.0/8.0-t24-t26+t27-t31;
679  values[15] = t23-3.0/8.0-t24+t36+t26-t37;
680  values[16] = t20;
681  values[17] = -t20;
682  values[18] = 0.0;
683  values[19] = 0.0;
684  values[20] = 0.0;
685  values[21] = 0.0;
686  values[22] = -t36+t7-t14+t16;
687  values[23] = 0.0;
688  values[24] = -t43+t44-1.0/4.0+t45+t47+t49;
689  values[25] = t51-t52+t53+t47-t54;
690  values[26] = -t51+t52+t56+t57-t49;
691  values[27] = t51-t52-t53+t47-t54;
692  values[28] = t51-t52+t56-t57+t49;
693  values[29] = t43-t44-1.0/4.0+t45-t47-t49;
694  values[30] = 0.0;
695  values[31] = -t7+t14-t31+t16;
696  values[32] = t36+t7-t14+t16;
697  values[33] = -t7+t14+t31+t16;
698  values[34] = 0.0;
699  values[35] = 0.0;
700  values[36] = 9.0*t15;
701  values[37] = 3.0*t35-3.0*eta;
702  values[38] = -3.0*xi+3.0*t30;
703  values[39] = 0.0;
704 }
705 
706 static void N_H_Q3_3D_DeriveXiZeta(double xi, double eta, double zeta,
707  double *values)
708 {
709  double t1, t7, t8, t14, t15, t16, t17, t18, t19, t20, t22, t23, t24, t25;
710  double t26, t27, t29, t30, t32, t33, t34, t35, t37, t38, t43, t44, t45;
711  double t46, t47, t48, t49, t50, t52, t53, t54, t55, t57, t58;
712 
713  t1 = xi*xi;
714  t7 = 3.0/4.0*(35.0/2.0*t1*xi-15.0/2.0*xi)*zeta;
715  t8 = zeta*zeta;
716  t14 = 3.0/4.0*xi*(35.0/2.0*t8*zeta-15.0/2.0*zeta);
717  t15 = xi*zeta;
718  t16 = 9.0/4.0*t15;
719  t17 = -t7+t14-t16;
720  t18 = t7-t14-t16;
721  t19 = 15.0/8.0*t1;
722  t20 = 15.0/8.0*t8;
723  t22 = 15.0/2.0*t1-3.0/2.0;
724  t23 = t22*zeta;
725  t24 = 3.0/2.0*t23;
726  t25 = eta*eta;
727  t26 = 3.0/8.0*t25;
728  t27 = 3.0/2.0*zeta;
729  t29 = xi*eta;
730  t30 = 3.0/4.0*t29;
731  t32 = 15.0/2.0*t8-3.0/2.0;
732  t33 = xi*t32;
733  t34 = 3.0/2.0*t33;
734  t35 = 3.0/2.0*xi;
735  t37 = eta*zeta;
736  t38 = 3.0/4.0*t37;
737  t43 = eta*t32;
738  t44 = t43/3.0;
739  t45 = 3.0/2.0*t37;
740  t46 = t22*eta;
741  t47 = t46/6.0;
742  t48 = t25*eta;
743  t49 = 5.0/12.0*t48;
744  t50 = 5.0/12.0*eta;
745  t52 = t43/6.0;
746  t53 = 3.0/4.0*t25;
747  t54 = 5.0/6.0*t48;
748  t55 = eta/3.0;
749  t57 = 3.0/2.0*t29;
750  t58 = t46/3.0;
751 
752  values[0] = t17;
753  values[1] = 0.0;
754  values[2] = t18;
755  values[3] = 0.0;
756  values[4] = t18;
757  values[5] = t17;
758  values[6] = t19-3.0/8.0-t20-t24+t26+t27;
759  values[7] = t30;
760  values[8] = t19+3.0/8.0-t20-t34-t26+t35;
761  values[9] = -t30;
762  values[10] = t38;
763  values[11] = -t30;
764  values[12] = t30;
765  values[13] = t38;
766  values[14] = -t38;
767  values[15] = t38;
768  values[16] = -t19-3.0/8.0+t20-t34+t26+t35;
769  values[17] = -t19+3.0/8.0+t20-t24-t26+t27;
770  values[18] = -t34-t7+t14+t16;
771  values[19] = 0.0;
772  values[20] = t24+t7-t14+t16;
773  values[21] = 0.0;
774  values[22] = 0.0;
775  values[23] = 0.0;
776  values[24] = -t44+t45+t47+t49-t50;
777  values[25] = t52-1.0/4.0+t53+t47-t54+t55;
778  values[26] = -t52+t57+t58-t49+t50;
779  values[27] = t52+1.0/4.0-t53+t47-t54+t55;
780  values[28] = t52+t57-t58+t49-t50;
781  values[29] = t44+t45-t47-t49+t50;
782  values[30] = 0.0;
783  values[31] = 0.0;
784  values[32] = 0.0;
785  values[33] = 0.0;
786  values[34] = -t24+t7-t14+t16;
787  values[35] = t34-t7+t14+t16;
788  values[36] = 9.0*t15;
789  values[37] = 3.0*t23-3.0*zeta;
790  values[38] = 0.0;
791  values[39] = 3.0*t33-3.0*xi;
792 }
793 
794 static void N_H_Q3_3D_DeriveEtaEta(double xi, double eta, double zeta,
795  double *values)
796 {
797  double t1, t2, t3, t4, t5, t7, t9, t11, t12, t13, t14, t15, t16, t17;
798  double t18, t20, t22, t24, t25, t27, t28, t29, t30, t33, t34, t35, t36;
799  double t38, t39, t40, t43, t49, t50, t52, t53, t55, t56, t57, t58;
800 
801  t1 = zeta*zeta;
802  t2 = 63.0/16.0*t1;
803  t3 = t1*t1;
804  t4 = 105.0/32.0*t3;
805  t5 = eta*eta;
806  t7 = 105.0/2.0*t5-15.0/2.0;
807  t9 = -1.0/2.0+3.0/2.0*t1;
808  t11 = t7*t9/4.0;
809  t12 = 21.0/32.0-t2+t4-t11;
810  t13 = 15.0/2.0*eta;
811  t14 = 27.0/16.0*t1;
812  t15 = xi*xi;
813  t16 = t15*t15;
814  t17 = 105.0/32.0*t16;
815  t18 = 27.0/16.0*t15;
816  t20 = -1.0/2.0+3.0/2.0*t15;
817  t22 = t20*t7/4.0;
818  t24 = 63.0/16.0*t15;
819  t25 = 21.0/32.0+t17-t24-t22;
820  t27 = xi*zeta;
821  t28 = 3.0/4.0*t27;
822  t29 = t1*zeta;
823  t30 = eta*zeta;
824  t33 = xi*eta;
825  t34 = 15.0/4.0*t33;
826  t35 = t20*eta;
827  t36 = 15.0/2.0*t35;
828  t38 = eta*t9;
829  t39 = 15.0/2.0*t38;
830  t40 = 15.0/4.0*t30;
831  t43 = t15*xi;
832  t49 = 15.0/4.0*t43;
833  t50 = 9.0/4.0*xi;
834  t52 = 15.0/4.0*t29;
835  t53 = 9.0/4.0*zeta;
836  t55 = t33*zeta;
837  t56 = 5.0/2.0*t55;
838  t57 = 3.0/2.0*t27;
839  t58 = 5.0*t55;
840 
841  values[0] = t12;
842  values[1] = -t13+27.0/16.0+t14-t4-t17+t18+t22+t11;
843  values[2] = t25;
844  values[3] = t13+27.0/16.0+t14-t4-t17+t18+t22+t11;
845  values[4] = t25;
846  values[5] = t12;
847  values[6] = t28;
848  values[7] = -15.0/4.0*t29+15.0/4.0*zeta-15.0/4.0*t30;
849  values[8] = -t28;
850  values[9] = -15.0/4.0*t29+15.0/4.0*zeta+15.0/4.0*t30;
851  values[10] = t34-t36;
852  values[11] = -t39-t40;
853  values[12] = -t39+t40;
854  values[13] = -15.0/4.0*t33-15.0/4.0*t43+15.0/4.0*xi;
855  values[14] = -t34-t36;
856  values[15] = -15.0/4.0*t33+15.0/4.0*t43-15.0/4.0*xi;
857  values[16] = t28;
858  values[17] = -t28;
859  values[18] = 0.0;
860  values[19] = -21.0/32.0+t2-t4-t39+t11;
861  values[20] = 0.0;
862  values[21] = -21.0/32.0+t2-t4+t39+t11;
863  values[22] = -t49+t50-3.0/32.0+t17-t18-t22;
864  values[23] = -3.0/32.0-t14+t4+t52-t53-t11;
865  values[24] = t56;
866  values[25] = t57-t58;
867  values[26] = -t56;
868  values[27] = -t57-t58;
869  values[28] = t56;
870  values[29] = -t56;
871  values[30] = -3.0/32.0-t14+t4-t52+t53-t11;
872  values[31] = -21.0/32.0-t17+t24+t22-t36;
873  values[32] = t49-t50-3.0/32.0+t17-t18-t22;
874  values[33] = -21.0/32.0-t17+t24+t22+t36;
875  values[34] = 0.0;
876  values[35] = 0.0;
877  values[36] = -6.0+9.0/2.0*t1+9.0/2.0*t15;
878  values[37] = 15.0/2.0*t43-15.0/2.0*xi;
879  values[38] = -15.0*eta+15.0*t38+15.0*t35;
880  values[39] = 15.0/2.0*t29-15.0/2.0*zeta;
881 }
882 
883 static void N_H_Q3_3D_DeriveEtaZeta(double xi, double eta, double zeta,
884  double *values)
885 {
886  double t1, t2, t3, t9, t10, t16, t17, t18, t19, t20, t21, t22, t24, t25;
887  double t26, t27, t28, t29, t32, t33, t34, t36, t37, t38, t44, t45, t46;
888  double t47, t48, t49, t50, t51, t53, t54, t55, t57, t58, t59;
889 
890  t1 = eta*zeta;
891  t2 = 9.0/4.0*t1;
892  t3 = zeta*zeta;
893  t9 = 3.0/4.0*eta*(35.0/2.0*t3*zeta-15.0/2.0*zeta);
894  t10 = eta*eta;
895  t16 = 3.0/4.0*(35.0/2.0*t10*eta-15.0/2.0*eta)*zeta;
896  t17 = -t2+t9-t16;
897  t18 = -t2-t9+t16;
898  t19 = xi*eta;
899  t20 = 3.0/4.0*t19;
900  t21 = xi*xi;
901  t22 = 3.0/8.0*t21;
902  t24 = 15.0/2.0*t3-3.0/2.0;
903  t25 = eta*t24;
904  t26 = 3.0/2.0*t25;
905  t27 = 3.0/2.0*eta;
906  t28 = 15.0/8.0*t3;
907  t29 = 15.0/8.0*t10;
908  t32 = xi*zeta;
909  t33 = 3.0/4.0*t32;
910  t34 = 3.0/2.0*zeta;
911  t36 = 15.0/2.0*t10-3.0/2.0;
912  t37 = t36*zeta;
913  t38 = 3.0/2.0*t37;
914  t44 = xi*t24;
915  t45 = t44/3.0;
916  t46 = 3.0/2.0*t32;
917  t47 = t21*xi;
918  t48 = 5.0/12.0*t47;
919  t49 = 5.0/12.0*xi;
920  t50 = xi*t36;
921  t51 = t50/6.0;
922  t53 = t44/6.0;
923  t54 = 3.0/2.0*t19;
924  t55 = t50/3.0;
925  t57 = 3.0/4.0*t21;
926  t58 = 5.0/6.0*t47;
927  t59 = xi/3.0;
928 
929  values[0] = t17;
930  values[1] = t18;
931  values[2] = 0.0;
932  values[3] = t18;
933  values[4] = 0.0;
934  values[5] = t17;
935  values[6] = t20;
936  values[7] = -3.0/8.0+t22-t26+t27+t28-t29;
937  values[8] = -t20;
938  values[9] = 3.0/8.0-t22-t26+t27-t28+t29;
939  values[10] = t33;
940  values[11] = t34+3.0/8.0-t22-t38+t28-t29;
941  values[12] = t34-3.0/8.0+t22-t38-t28+t29;
942  values[13] = t33;
943  values[14] = -t33;
944  values[15] = t33;
945  values[16] = t20;
946  values[17] = -t20;
947  values[18] = 0.0;
948  values[19] = t2-t9-t38+t16;
949  values[20] = 0.0;
950  values[21] = t2-t9+t38+t16;
951  values[22] = 0.0;
952  values[23] = t2+t9+t26-t16;
953  values[24] = -t45+t46+t48-t49+t51;
954  values[25] = t53+t54+t48-t49-t55;
955  values[26] = -t53-1.0/4.0+t57+t58-t59-t51;
956  values[27] = t53-t54+t48-t49-t55;
957  values[28] = t53-1.0/4.0+t57-t58+t59+t51;
958  values[29] = t45+t46-t48+t49-t51;
959  values[30] = t2+t9-t26-t16;
960  values[31] = 0.0;
961  values[32] = 0.0;
962  values[33] = 0.0;
963  values[34] = 0.0;
964  values[35] = 0.0;
965  values[36] = 9.0*t1;
966  values[37] = 0.0;
967  values[38] = -3.0*zeta+3.0*t37;
968  values[39] = 3.0*t25-3.0*eta;
969 }
970 
971 static void N_H_Q3_3D_DeriveZetaZeta(double xi, double eta, double zeta,
972  double *values)
973 {
974  double t1, t2, t3, t4, t5, t6, t7, t9, t10, t12, t14, t16, t18, t19, t20;
975  double t22, t23, t24, t25, t27, t28, t30, t31, t32, t33, t35, t36, t37;
976  double t40, t41, t42, t48, t49, t51, t52, t56, t57, t58, t60;
977 
978  t1 = eta*eta;
979  t2 = 27.0/16.0*t1;
980  t3 = 15.0/2.0*zeta;
981  t4 = xi*xi;
982  t5 = t4*t4;
983  t6 = 105.0/32.0*t5;
984  t7 = 27.0/16.0*t4;
985  t9 = -1.0/2.0+3.0/2.0*t4;
986  t10 = zeta*zeta;
987  t12 = 105.0/2.0*t10-15.0/2.0;
988  t14 = t9*t12/4.0;
989  t16 = -1.0/2.0+3.0/2.0*t1;
990  t18 = t16*t12/4.0;
991  t19 = t1*t1;
992  t20 = 105.0/32.0*t19;
993  t22 = 63.0/16.0*t1;
994  t23 = 21.0/32.0-t22-t18+t20;
995  t24 = 63.0/16.0*t4;
996  t25 = 21.0/32.0+t6-t24-t14;
997  t27 = xi*zeta;
998  t28 = t4*xi;
999  t30 = t16*zeta;
1000  t31 = 15.0/2.0*t30;
1001  t32 = eta*zeta;
1002  t33 = 15.0/4.0*t32;
1003  t35 = 15.0/4.0*t27;
1004  t36 = t9*zeta;
1005  t37 = 15.0/2.0*t36;
1006  t40 = xi*eta;
1007  t41 = 3.0/4.0*t40;
1008  t42 = t1*eta;
1009  t48 = 15.0/4.0*t42;
1010  t49 = 9.0/4.0*eta;
1011  t51 = 15.0/4.0*t28;
1012  t52 = 9.0/4.0*xi;
1013  t56 = t40*zeta;
1014  t57 = 5.0*t56;
1015  t58 = 3.0/2.0*t40;
1016  t60 = 5.0/2.0*t56;
1017 
1018  values[0] = 27.0/16.0+t2-t3-t6+t7+t14+t18-t20;
1019  values[1] = t23;
1020  values[2] = t25;
1021  values[3] = t23;
1022  values[4] = t25;
1023  values[5] = 27.0/16.0+t2+t3-t6+t7+t14+t18-t20;
1024  values[6] = -15.0/4.0*t27-15.0/4.0*t28+15.0/4.0*xi;
1025  values[7] = -t31+t33;
1026  values[8] = -t35-t37;
1027  values[9] = -t31-t33;
1028  values[10] = t41;
1029  values[11] = 15.0/4.0*eta-15.0/4.0*t42+15.0/4.0*t32;
1030  values[12] = 15.0/4.0*eta-15.0/4.0*t42-15.0/4.0*t32;
1031  values[13] = t41;
1032  values[14] = -t41;
1033  values[15] = t41;
1034  values[16] = t35-t37;
1035  values[17] = 15.0/4.0*t27-15.0/4.0*t28+15.0/4.0*xi;
1036  values[18] = -t37-21.0/32.0-t6+t24+t14;
1037  values[19] = -3.0/32.0-t2-t18-t48+t49+t20;
1038  values[20] = t51-t52-3.0/32.0+t6-t7-t14;
1039  values[21] = -3.0/32.0-t2-t18+t48-t49+t20;
1040  values[22] = 0.0;
1041  values[23] = -21.0/32.0+t22+t18+t31-t20;
1042  values[24] = -t57+t58;
1043  values[25] = t60;
1044  values[26] = -t60;
1045  values[27] = t60;
1046  values[28] = t60;
1047  values[29] = t57+t58;
1048  values[30] = -21.0/32.0+t22+t18-t31-t20;
1049  values[31] = 0.0;
1050  values[32] = 0.0;
1051  values[33] = 0.0;
1052  values[34] = -t51+t52-3.0/32.0+t6-t7-t14;
1053  values[35] = t37-21.0/32.0-t6+t24+t14;
1054  values[36] = -6.0+9.0/2.0*t1+9.0/2.0*t4;
1055  values[37] = 15.0/2.0*t28-15.0/2.0*xi;
1056  values[38] = -15.0/2.0*eta+15.0/2.0*t42;
1057  values[39] = 15.0*t36-15.0*zeta+15.0*t30;
1058 }
1059 
1060 static int N_H_Q3_3D_ChangeU0[2] = { 6, 24 };
1061 static int N_H_Q3_3D_ChangeU1[2] = { 7, 25 };
1062 static int N_H_Q3_3D_ChangeU2[2] = { 8, 26 };
1063 static int N_H_Q3_3D_ChangeU3[2] = { 9, 27 };
1064 static int N_H_Q3_3D_ChangeU4[2] = { 10, 28 };
1065 static int N_H_Q3_3D_ChangeU5[2] = { 11, 29 };
1066 
1067 static int N_H_Q3_3D_ChangeV0[2] = { 12, 24 };
1068 static int N_H_Q3_3D_ChangeV1[2] = { 13, 25 };
1069 static int N_H_Q3_3D_ChangeV2[2] = { 14, 26 };
1070 static int N_H_Q3_3D_ChangeV3[2] = { 15, 27 };
1071 static int N_H_Q3_3D_ChangeV4[2] = { 16, 28 };
1072 static int N_H_Q3_3D_ChangeV5[2] = { 17, 29 };
1073 
1074 static int *N_H_Q3_3D_ChangeU[6] = { N_H_Q3_3D_ChangeU0, N_H_Q3_3D_ChangeU1,
1075  N_H_Q3_3D_ChangeU2, N_H_Q3_3D_ChangeU3,
1076  N_H_Q3_3D_ChangeU4, N_H_Q3_3D_ChangeU5 };
1077 
1078 static int *N_H_Q3_3D_ChangeV[6] = { N_H_Q3_3D_ChangeV0, N_H_Q3_3D_ChangeV1,
1079  N_H_Q3_3D_ChangeV2, N_H_Q3_3D_ChangeV3,
1080  N_H_Q3_3D_ChangeV4, N_H_Q3_3D_ChangeV5 };
1081 
1082 static int **N_H_Q3_3D_Change[2] = { N_H_Q3_3D_ChangeU, N_H_Q3_3D_ChangeV };
1083 
1084 TBaseFunct3D *BF_N_H_Q3_3D_Obj =
1085 new TBaseFunct3D(40, BF_N_H_Q3_3D, BFUnitHexahedron,
1086  N_H_Q3_3D_Funct, N_H_Q3_3D_DeriveXi,
1087  N_H_Q3_3D_DeriveEta, N_H_Q3_3D_DeriveZeta,
1088  N_H_Q3_3D_DeriveXiXi, N_H_Q3_3D_DeriveXiEta,
1089  N_H_Q3_3D_DeriveXiZeta, N_H_Q3_3D_DeriveEtaEta,
1090  N_H_Q3_3D_DeriveEtaZeta, N_H_Q3_3D_DeriveZetaZeta,
1091  4, 3,
1092  2, N_H_Q3_3D_Change);
TBaseFunct3D
Definition: BaseFunct3D.h:27
Generated on Tue May 10 2016 09:08:00 for ParMooN by   doxygen 1.8.5