ParMooN
 All Classes Functions Variables Friends Pages
BF_C_Q_UL5SE_2D.h
1 // ***********************************************************************
2 // Q5 element with bubbles, conforming, 2D
3 // ***********************************************************************
4 
5 // base function values
6 static void C_Q_UL5SE_2D_Funct(double xi, double eta, double *values)
7 {
8  double t1 = eta*eta;
9  double t2 = xi*t1;
10  double t3 = 359.0/6144.0*t2;
11  double t4 = xi*xi;
12  double t5 = t4*t1;
13  double t6 = 31.0/384.0*t5;
14  double t7 = t4*eta;
15  double t8 = 359.0/6144.0*t7;
16  double t9 = t4*xi;
17  double t10 = t9*eta;
18  double t11 = t10/96.0;
19  double t12 = t9*t1;
20  double t13 = 453.0/2048.0*t12;
21  double t14 = t4*t4;
22  double t15 = t14*t1;
23  double t16 = 371.0/1536.0*t15;
24  double t17 = t1*eta;
25  double t18 = xi*t17;
26  double t19 = t18/96.0;
27  double t20 = t1/512.0;
28  double t21 = xi*eta;
29  double t22 = 191.0/768.0*t21;
30  double t23 = 923.0/6144.0*t9;
31  double t24 = 887.0/6144.0*xi;
32  double t25 = 887.0/6144.0*eta;
33  double t26 = 923.0/6144.0*t17;
34  double t27 = -t3+t6-t8+t11+t13-t16+t19-t20-t22-t23+t24+t25-t26;
35  double t28 = t4/512.0;
36  double t29 = t14/128.0;
37  double t31 = 1.0-t4;
38  double t35 = 5.0/2.0*t9-3.0/2.0*xi;
39  double t36 = (1.0-eta)*t31*t35;
40  double t37 = 125.0/768.0*t36;
41  double t38 = t4*t17;
42  double t39 = 453.0/2048.0*t38;
43  double t40 = t9*t17;
44  double t41 = 367.0/768.0*t40;
45  double t42 = t14*t17;
46  double t43 = 859.0/1536.0*t42;
47  double t44 = t1*t1;
48  double t45 = xi*t44;
49  double t46 = 39.0/256.0*t45;
50  double t47 = t4*t44;
51  double t48 = 371.0/1536.0*t47;
52  double t49 = t9*t44;
53  double t50 = 859.0/1536.0*t49;
54  double t51 = t14*t44;
55  double t52 = 41.0/64.0*t51;
56  double t53 = t44/128.0;
57  double t54 = t14*eta;
58  double t55 = 39.0/256.0*t54;
59  double t57 = 1.0-t1;
60  double t61 = 5.0/2.0*t17-3.0/2.0*eta;
61  double t62 = (1.0-xi)*t57*t61;
62  double t63 = 125.0/768.0*t62;
63  double t64 = -t28+t29+t37+t39+t41-t43+t46-t48-t50+t52+t53+t55+t63;
64  double t66 = 4375.0/6144.0*xi;
65  double t67 = 25.0/1536.0*eta;
66  double t68 = 875.0/1536.0*t21;
67  double t69 = 325.0/768.0*t7;
68  double t70 = 875.0/1536.0*t10;
69  double t71 = 625.0/1536.0*t54;
70  double t72 = 875.0/6144.0*t2;
71  double t73 = 325.0/768.0*t5;
72  double t74 = 875.0/6144.0*t12;
73  double t75 = 625.0/1536.0*t15;
74  double t76 = 875.0/1536.0*t18;
75  double t77 = -t66-t67+t68+t69-t70-t71+t72-t73-t74+t75+t76;
76  double t78 = 325.0/192.0*t38;
77  double t79 = 875.0/1536.0*t40;
78  double t80 = 625.0/384.0*t42;
79  double t81 = 875.0/1536.0*t45;
80  double t82 = 325.0/192.0*t47;
81  double t83 = 875.0/1536.0*t49;
82  double t84 = 625.0/384.0*t51;
83  double t85 = 4375.0/6144.0*t9;
84  double t86 = 25.0/1536.0*t1;
85  double t87 = 25.0/384.0*t17;
86  double t88 = 25.0/384.0*t44;
87  double t89 = 625.0/768.0*t36;
88  double t90 = -t78-t79+t80-t81+t82+t83-t84+t85+t86+t87-t88-t89;
89  double t92 = 4375.0/3072.0*xi;
90  double t93 = 25.0/256.0*eta;
91  double t94 = 3625.0/2304.0*t21;
92  double t95 = 425.0/1152.0*t7;
93  double t96 = 3625.0/2304.0*t10;
94  double t97 = 625.0/2304.0*t54;
95  double t98 = 1375.0/9216.0*t2;
96  double t99 = 425.0/1152.0*t5;
97  double t100 = 1375.0/9216.0*t12;
98  double t101 = 625.0/2304.0*t15;
99  double t102 = 1375.0/2304.0*t18;
100  double t103 = t92+t93-t94-t95+t96+t97+t98+t99-t100-t101+t102;
101  double t104 = 425.0/288.0*t38;
102  double t105 = 1375.0/2304.0*t40;
103  double t106 = 625.0/576.0*t42;
104  double t107 = 1375.0/2304.0*t45;
105  double t108 = 425.0/288.0*t47;
106  double t109 = 1375.0/2304.0*t49;
107  double t110 = 625.0/576.0*t51;
108  double t111 = 4375.0/3072.0*t9;
109  double t112 = 25.0/256.0*t1;
110  double t113 = 25.0/64.0*t17;
111  double t114 = 25.0/64.0*t44;
112  double t115 = 625.0/384.0*t36;
113  double t116 = t104-t105-t106-t107-t108+t109+t110-t111-t112-t113+t114+t115;
114  double t118 = -t92+t93+t94-t95-t96+t97-t98+t99+t100-t101-t102;
115  double t119 = t104+t105-t106+t107-t108-t109+t110+t111-t112-t113+t114-t115;
116  double t121 = t66-t67-t68+t69+t70-t71-t72-t73+t74+t75-t76;
117  double t122 = -t78+t79+t80+t81+t82-t83-t84-t85+t86+t87-t88+t89;
118  double t124 = t3+t6-t8-t11-t13-t16-t19-t20+t22+t23-t24+t25-t26;
119  double t127 = (1.0+xi)*t57*t61;
120  double t128 = 125.0/768.0*t127;
121  double t129 = -t28+t29-t37+t39-t41-t43-t46-t48+t50+t52+t53+t55+t128;
122  double t131 = 25.0/1536.0*xi;
123  double t132 = 4375.0/6144.0*eta;
124  double t133 = 875.0/6144.0*t7;
125  double t134 = 875.0/1536.0*t54;
126  double t135 = 325.0/768.0*t2;
127  double t136 = 325.0/192.0*t12;
128  double t137 = 325.0/192.0*t15;
129  double t138 = t131-t132-t68+t133-t70-t134-t135-t73+t136+t137+t76;
130  double t139 = 875.0/6144.0*t38;
131  double t140 = 875.0/1536.0*t42;
132  double t141 = 625.0/1536.0*t45;
133  double t142 = 625.0/1536.0*t47;
134  double t143 = 625.0/384.0*t49;
135  double t144 = 25.0/1536.0*t4;
136  double t145 = 25.0/384.0*t9;
137  double t146 = 25.0/384.0*t14;
138  double t147 = 4375.0/6144.0*t17;
139  double t148 = 625.0/768.0*t127;
140  double t149 = -t139+t79+t140+t141+t142-t143-t84+t144-t145-t146+t147-t148;
141  double t151 = 25.0/256.0*xi;
142  double t152 = 4375.0/3072.0*eta;
143  double t153 = 1375.0/9216.0*t7;
144  double t154 = 1375.0/2304.0*t10;
145  double t155 = 1375.0/2304.0*t54;
146  double t156 = 425.0/1152.0*t2;
147  double t157 = 425.0/288.0*t12;
148  double t158 = 425.0/288.0*t15;
149  double t159 = 3625.0/2304.0*t18;
150  double t160 = -t151+t152+t94+t153-t154-t155+t156+t99-t157-t158-t159;
151  double t161 = 1375.0/9216.0*t38;
152  double t162 = 1375.0/2304.0*t42;
153  double t163 = 625.0/2304.0*t45;
154  double t164 = 625.0/2304.0*t47;
155  double t165 = 625.0/576.0*t49;
156  double t166 = 25.0/256.0*t4;
157  double t167 = 25.0/64.0*t9;
158  double t168 = 25.0/64.0*t14;
159  double t169 = 4375.0/3072.0*t17;
160  double t170 = 625.0/384.0*t127;
161  double t171 = -t161+t105+t162-t163-t164+t165+t110-t166+t167+t168-t169+t170;
162  double t173 = -t151-t152-t94-t153+t154+t155+t156+t99-t157-t158+t159;
163  double t174 = t161-t105-t162-t163-t164+t165+t110-t166+t167+t168+t169-t170;
164  double t176 = t131+t132+t68-t133+t70+t134-t135-t73+t136+t137-t76;
165  double t177 = t139-t79-t140+t141+t142-t143-t84+t144-t145-t146-t147+t148;
166  double t179 = t3+t6+t8+t11-t13-t16+t19-t20-t22+t23-t24-t25+t26;
167  double t182 = (1.0+eta)*t31*t35;
168  double t183 = 125.0/768.0*t182;
169  double t184 = -t28+t29-t39+t41+t43-t46-t48+t50+t52+t53-t55-t128-t183;
170  double t186 = t66+t67+t68-t69-t70+t71-t72-t73+t74+t75+t76;
171  double t187 = 625.0/768.0*t182;
172  double t188 = t78-t79-t80+t81+t82-t83-t84-t85+t86-t87-t88+t187;
173  double t190 = -t92-t93-t94+t95+t96-t97-t98+t99+t100-t101+t102;
174  double t191 = 625.0/384.0*t182;
175  double t192 = -t104-t105+t106+t107-t108-t109+t110+t111-t112+t113+t114-t191;
176  double t194 = t92-t93+t94+t95-t96-t97+t98+t99-t100-t101-t102;
177  double t195 = -t104+t105+t106-t107-t108+t109+t110-t111-t112+t113+t114+t191;
178  double t197 = -t66+t67-t68-t69+t70+t71+t72-t73-t74+t75-t76;
179  double t198 = t78+t79-t80-t81+t82+t83-t84+t85+t86-t87-t88-t187;
180  double t200 = -t3+t6+t8-t11+t13-t16-t19-t20+t22-t23+t24-t25+t26;
181  double t201 = -t28+t29-t39-t41+t43+t46-t48-t50+t52+t53-t55+t183-t63;
182  double t203 = -t131+t132-t68-t133-t70+t134+t135-t73-t136+t137+t76;
183  double t204 = 625.0/768.0*t62;
184  double t205 = t139+t79-t140-t141+t142+t143-t84+t144+t145-t146-t147+t204;
185  double t207 = t151-t152+t94-t153-t154+t155-t156+t99+t157-t158-t159;
186  double t208 = 625.0/384.0*t62;
187  double t209 = t161+t105-t162+t163-t164-t165+t110-t166-t167+t168+t169-t208;
188  double t211 = t151+t152-t94+t153+t154-t155-t156+t99+t157-t158+t159;
189  double t212 = -t161-t105+t162+t163-t164-t165+t110-t166-t167+t168-t169+t208;
190  double t214 = -t131-t132+t68+t133+t70-t134+t135-t73-t136+t137-t76;
191  double t215 = -t139-t79+t140-t141+t142+t143-t84+t144+t145-t146+t147-t204;
192  double t217 = 16.0/9.0*t21;
193  double t218 = 32.0/9.0*t7;
194  double t219 = 16.0/9.0*t10;
195  double t220 = 32.0/9.0*t54;
196  double t221 = 32.0/9.0*t2;
197  double t222 = 64.0/9.0*t5;
198  double t223 = 32.0/9.0*t12;
199  double t224 = 64.0/9.0*t15;
200  double t225 = 16.0/9.0*t18;
201  double t226 = 32.0/9.0*t38;
202  double t227 = 16.0/9.0*t40;
203  double t228 = 32.0/9.0*t42;
204  double t229 = 32.0/9.0*t45;
205  double t230 = 64.0/9.0*t47;
206  double t231 = 32.0/9.0*t49;
207  double t232 = 64.0/9.0*t51;
208  double t233 = t217-t218-t219+t220-t221+t222+t223-t224-t225+t226+t227-t228+t229-t230-t231+t232;
209  double t234 = 4.0/3.0*xi;
210  double t235 = 20.0/3.0*t2;
211  double t236 = 40.0/3.0*t5;
212  double t237 = 20.0/3.0*t12;
213  double t238 = 40.0/3.0*t15;
214  double t239 = 16.0/3.0*t45;
215  double t240 = 32.0/3.0*t47;
216  double t241 = 16.0/3.0*t49;
217  double t242 = 32.0/3.0*t51;
218  double t243 = 8.0/3.0*t4;
219  double t244 = 4.0/3.0*t9;
220  double t245 = 8.0/3.0*t14;
221  double t246 = -t234+t235-t236-t237+t238-t239+t240+t241-t242+t243+t244-t245;
222  double t247 = -t217+t218+t219-t220-t221+t222+t223-t224+t225-t226-t227+t228+t229-t230-t231+t232;
223  double t248 = 4.0/3.0*eta;
224  double t249 = 20.0/3.0*t7;
225  double t250 = 16.0/3.0*t54;
226  double t251 = 32.0/3.0*t15;
227  double t252 = 20.0/3.0*t38;
228  double t253 = 16.0/3.0*t42;
229  double t254 = 40.0/3.0*t47;
230  double t255 = 8.0/3.0*t1;
231  double t256 = 4.0/3.0*t17;
232  double t257 = 8.0/3.0*t44;
233  double t258 = -t248+t249-t250-t236+t251-t252+t253+t254-t242+t255+t256-t257;
234  double t268 = t248-t249+t250-t236+t251+t252-t253+t254-t242+t255-t256-t257;
235  double t269 = -t217-t218+t219+t220+t221+t222-t223-t224+t225+t226-t227-t228-t229-t230+t231+t232;
236  double t270 = t234-t235-t236+t237+t238+t239+t240-t241-t242+t243-t244-t245;
237  double t271 = t217+t218-t219-t220+t221+t222-t223-t224-t225-t226+t227+t228-t229-t230+t231+t232;
238  values[0] = t27+t64;
239  values[1] = t77+t90;
240  values[2] = t103+t116;
241  values[3] = t118+t119;
242  values[4] = t121+t122;
243  values[5] = t124+t129;
244  values[6] = t138+t149;
245  values[7] = t160+t171;
246  values[8] = t173+t174;
247  values[9] = t176+t177;
248  values[10] = t179+t184;
249  values[11] = t186+t188;
250  values[12] = t190+t192;
251  values[13] = t194+t195;
252  values[14] = t197+t198;
253  values[15] = t200+t201;
254  values[16] = t203+t205;
255  values[17] = t207+t209;
256  values[18] = t211+t212;
257  values[19] = t214+t215;
258  values[20] = t233;
259  values[21] = t246;
260  values[22] = t247;
261  values[23] = t258;
262  values[24] = 1.0+25.0*t5-20.0*t15-20.0*t47+16.0*t51-5.0*t4+4.0*t14-5.0*t1+4.0*t44;
263  values[25] = t268;
264  values[26] = t269;
265  values[27] = t270;
266  values[28] = t271;
267 }
268 
269 // values of the derivatives in xi direction
270 static void C_Q_UL5SE_2D_DeriveXi(double xi, double eta, double *values)
271 {
272  double t1 = eta*eta;
273  double t3 = t1*eta;
274  double t7 = (1.0-t1)*(5.0/2.0*t3-3.0/2.0*eta);
275  double t8 = 125.0/768.0*t7;
276  double t9 = 1.0-eta;
277  double t11 = xi*xi;
278  double t12 = t11*xi;
279  double t15 = 5.0/2.0*t12-3.0/2.0*xi;
280  double t16 = t9*xi*t15;
281  double t17 = 125.0/384.0*t16;
282  double t18 = 1.0-t11;
283  double t21 = 15.0/2.0*t11-3.0/2.0;
284  double t22 = t9*t18*t21;
285  double t23 = 125.0/768.0*t22;
286  double t24 = 191.0/768.0*eta;
287  double t25 = xi*eta;
288  double t26 = 359.0/3072.0*t25;
289  double t27 = t11*eta;
290  double t28 = t27/32.0;
291  double t29 = t12*eta;
292  double t30 = 39.0/64.0*t29;
293  double t31 = 359.0/6144.0*t1;
294  double t32 = xi*t1;
295  double t33 = 31.0/192.0*t32;
296  double t34 = t11*t1;
297  double t35 = 1359.0/2048.0*t34;
298  double t36 = t12*t1;
299  double t37 = 371.0/384.0*t36;
300  double t38 = -t8-t17+t23-t24-t26+t28+t30-t31+t33+t35-t37;
301  double t39 = t3/96.0;
302  double t40 = xi*t3;
303  double t41 = 453.0/1024.0*t40;
304  double t42 = t11*t3;
305  double t43 = 367.0/256.0*t42;
306  double t44 = t12*t3;
307  double t45 = 859.0/384.0*t44;
308  double t46 = t1*t1;
309  double t47 = 39.0/256.0*t46;
310  double t48 = xi*t46;
311  double t49 = 371.0/768.0*t48;
312  double t50 = t11*t46;
313  double t51 = 859.0/512.0*t50;
314  double t52 = t12*t46;
315  double t53 = 41.0/16.0*t52;
316  double t54 = xi/256.0;
317  double t55 = 923.0/2048.0*t11;
318  double t56 = t12/32.0;
319  double t57 = 887.0/6144.0+t39+t41+t43-t45+t47-t49-t51+t53-t54-t55+t56;
320  double t59 = 875.0/1536.0*eta;
321  double t60 = 325.0/384.0*t25;
322  double t61 = 875.0/512.0*t27;
323  double t62 = 625.0/384.0*t29;
324  double t63 = 875.0/6144.0*t1;
325  double t64 = 325.0/384.0*t32;
326  double t65 = 875.0/2048.0*t34;
327  double t66 = 625.0/384.0*t36;
328  double t67 = 875.0/1536.0*t3;
329  double t68 = 325.0/96.0*t40;
330  double t69 = 875.0/512.0*t42;
331  double t70 = 625.0/96.0*t44;
332  double t71 = 875.0/1536.0*t46;
333  double t72 = 325.0/96.0*t48;
334  double t73 = 875.0/512.0*t50;
335  double t74 = 625.0/96.0*t52;
336  double t75 = 4375.0/2048.0*t11;
337  double t76 = 625.0/384.0*t16;
338  double t77 = 625.0/768.0*t22;
339  double t78 = -4375.0/6144.0+t59+t60-t61-t62+t63-t64-t65+t66+t67-t68-t69+t70-t71+t72+t73-t74+t75+t76-t77;
340  double t79 = 3625.0/2304.0*eta;
341  double t80 = 425.0/576.0*t25;
342  double t81 = 3625.0/768.0*t27;
343  double t82 = 625.0/576.0*t29;
344  double t83 = 1375.0/9216.0*t1;
345  double t84 = 425.0/576.0*t32;
346  double t85 = 1375.0/3072.0*t34;
347  double t86 = 625.0/576.0*t36;
348  double t87 = 1375.0/2304.0*t3;
349  double t88 = 425.0/144.0*t40;
350  double t89 = 1375.0/768.0*t42;
351  double t90 = 625.0/144.0*t44;
352  double t91 = 1375.0/2304.0*t46;
353  double t92 = 425.0/144.0*t48;
354  double t93 = 1375.0/768.0*t50;
355  double t94 = 625.0/144.0*t52;
356  double t95 = 4375.0/1024.0*t11;
357  double t96 = 625.0/192.0*t16;
358  double t97 = 625.0/384.0*t22;
359  double t98 = 4375.0/3072.0-t79-t80+t81+t82+t83+t84-t85-t86+t87+t88-t89-t90-t91-t92+t93+t94-t95-t96+t97;
360  double t99 = -4375.0/3072.0+t79-t80-t81+t82-t83+t84+t85-t86-t87+t88+t89-t90+t91-t92-t93+t94+t95+t96-t97;
361  double t100 = 4375.0/6144.0-t59+t60+t61-t62-t63-t64+t65+t66-t67-t68+t69+t70+t71+t72-t73-t74-t75-t76+t77;
362  double t101 = t17-t23+t8+t24-t26-t28+t30+t31+t33-t35-t37;
363  double t102 = -887.0/6144.0-t39+t41-t43-t45-t47-t49+t51+t53-t54+t55+t56;
364  double t104 = 875.0/3072.0*t25;
365  double t105 = 875.0/384.0*t29;
366  double t106 = 325.0/768.0*t1;
367  double t107 = 325.0/64.0*t34;
368  double t108 = 325.0/48.0*t36;
369  double t110 = 875.0/3072.0*t40;
370  double t111 = 875.0/384.0*t44;
371  double t112 = 625.0/1536.0*t46;
372  double t113 = 625.0/768.0*t48;
373  double t114 = 625.0/128.0*t50;
374  double t115 = 25.0/768.0*xi;
375  double t116 = 25.0/128.0*t11;
376  double t117 = 25.0/96.0*t12;
377  double t118 = 625.0/768.0*t7;
378  double t119 = -t110+t69+t111+t112+t113-t114-t74+t115-t116-t117-t118;
379  double t121 = 1375.0/4608.0*t25;
380  double t122 = 1375.0/768.0*t27;
381  double t123 = 1375.0/576.0*t29;
382  double t124 = 425.0/1152.0*t1;
383  double t125 = 425.0/96.0*t34;
384  double t126 = 425.0/72.0*t36;
385  double t127 = 3625.0/2304.0*t3;
386  double t129 = 1375.0/4608.0*t40;
387  double t130 = 1375.0/576.0*t44;
388  double t131 = 625.0/2304.0*t46;
389  double t132 = 625.0/1152.0*t48;
390  double t133 = 625.0/192.0*t50;
391  double t134 = 25.0/128.0*xi;
392  double t135 = 75.0/64.0*t11;
393  double t136 = 25.0/16.0*t12;
394  double t137 = 625.0/384.0*t7;
395  double t138 = -t129+t89+t130-t131-t132+t133+t94-t134+t135+t136+t137;
396  double t141 = t129-t89-t130-t131-t132+t133+t94-t134+t135+t136-t137;
397  double t144 = t110-t69-t111+t112+t113-t114-t74+t115-t116-t117+t118;
398  double t146 = 1.0+eta;
399  double t148 = t146*xi*t15;
400  double t149 = 125.0/384.0*t148;
401  double t151 = t146*t18*t21;
402  double t152 = 125.0/768.0*t151;
403  double t153 = t149-t152-t8-t24+t26+t28-t30+t31+t33-t35-t37;
404  double t154 = -887.0/6144.0+t39-t41+t43+t45-t47-t49+t51+t53-t54+t55+t56;
405  double t156 = 625.0/384.0*t148;
406  double t157 = 625.0/768.0*t151;
407  double t158 = 4375.0/6144.0+t59-t60-t61+t62-t63-t64+t65+t66+t67+t68-t69-t70+t71+t72-t73-t74-t75-t156+t157;
408  double t159 = 625.0/192.0*t148;
409  double t160 = 625.0/384.0*t151;
410  double t161 = -4375.0/3072.0-t79+t80+t81-t82-t83+t84+t85-t86+t87-t88-t89+t90+t91-t92-t93+t94+t95+t159-t160;
411  double t162 = 4375.0/3072.0+t79+t80-t81-t82+t83+t84-t85-t86-t87-t88+t89+t90-t91-t92+t93+t94-t95-t159+t160;
412  double t163 = -4375.0/6144.0-t59-t60+t61+t62+t63-t64-t65+t66-t67+t68+t69-t70-t71+t72+t73-t74+t75+t156-t157;
413  double t164 = t8-t149+t152+t24+t26-t28-t30-t31+t33+t35-t37;
414  double t165 = 887.0/6144.0-t39-t41-t43+t45+t47-t49-t51+t53-t54-t55+t56;
415  double t168 = t110+t69-t111-t112+t113+t114-t74+t115+t116-t117-t118;
416  double t171 = t129+t89-t130+t131-t132-t133+t94-t134-t135+t136+t137;
417  double t174 = -t129-t89+t130+t131-t132-t133+t94-t134-t135+t136-t137;
418  double t177 = -t110-t69+t111-t112+t113+t114-t74+t115+t116-t117+t118;
419  double t179 = 16.0/9.0*eta;
420  double t180 = 64.0/9.0*t25;
421  double t181 = 16.0/3.0*t27;
422  double t182 = 128.0/9.0*t29;
423  double t183 = 32.0/9.0*t1;
424  double t184 = 128.0/9.0*t32;
425  double t185 = 32.0/3.0*t34;
426  double t186 = 256.0/9.0*t36;
427  double t187 = 16.0/9.0*t3;
428  double t188 = 64.0/9.0*t40;
429  double t189 = 16.0/3.0*t42;
430  double t190 = 128.0/9.0*t44;
431  double t191 = 32.0/9.0*t46;
432  double t192 = 128.0/9.0*t48;
433  double t193 = 32.0/3.0*t50;
434  double t194 = 256.0/9.0*t52;
435  double t195 = t179-t180-t181+t182-t183+t184+t185-t186-t187+t188+t189-t190+t191-t192-t193+t194;
436  double t196 = 20.0/3.0*t1;
437  double t197 = 80.0/3.0*t32;
438  double t198 = 20.0*t34;
439  double t199 = 160.0/3.0*t36;
440  double t200 = 16.0/3.0*t46;
441  double t201 = 64.0/3.0*t48;
442  double t202 = 16.0*t50;
443  double t203 = 128.0/3.0*t52;
444  double t204 = 16.0/3.0*xi;
445  double t205 = 4.0*t11;
446  double t206 = 32.0/3.0*t12;
447  double t207 = -4.0/3.0+t196-t197-t198+t199-t200+t201+t202-t203+t204+t205-t206;
448  double t208 = -t179+t180+t181-t182-t183+t184+t185-t186+t187-t188-t189+t190+t191-t192-t193+t194;
449  double t209 = 40.0/3.0*t25;
450  double t210 = 64.0/3.0*t29;
451  double t211 = 128.0/3.0*t36;
452  double t212 = 40.0/3.0*t40;
453  double t213 = 64.0/3.0*t44;
454  double t214 = 80.0/3.0*t48;
455  double t224 = -t179-t180+t181+t182+t183+t184-t185-t186+t187+t188-t189-t190-t191-t192+t193+t194;
456  double t225 = 4.0/3.0-t196-t197+t198+t199+t200+t201-t202-t203+t204-t205-t206;
457  double t226 = t179+t180-t181-t182+t183+t184-t185-t186-t187-t188+t189+t190-t191-t192+t193+t194;
458  values[0] = t38+t57;
459  values[1] = t78;
460  values[2] = t98;
461  values[3] = t99;
462  values[4] = t100;
463  values[5] = t101+t102;
464  values[6] = 25.0/1536.0-t59+t104-t61-t105-t106-t64+t107+t108+t67+t119;
465  values[7] = -25.0/256.0+t79+t121-t122-t123+t124+t84-t125-t126-t127+t138;
466  values[8] = -25.0/256.0-t79-t121+t122+t123+t124+t84-t125-t126+t127+t141;
467  values[9] = 25.0/1536.0+t59-t104+t61+t105-t106-t64+t107+t108-t67+t144;
468  values[10] = t153+t154;
469  values[11] = t158;
470  values[12] = t161;
471  values[13] = t162;
472  values[14] = t163;
473  values[15] = t164+t165;
474  values[16] = -25.0/1536.0-t59-t104-t61+t105+t106-t64-t107+t108+t67+t168;
475  values[17] = 25.0/256.0+t79-t121-t122+t123-t124+t84+t125-t126-t127+t171;
476  values[18] = 25.0/256.0-t79+t121+t122-t123-t124+t84+t125-t126+t127+t174;
477  values[19] = -25.0/1536.0+t59+t104+t61-t105+t106-t64-t107+t108-t67+t177;
478  values[20] = t195;
479  values[21] = t207;
480  values[22] = t208;
481  values[23] = t209-t210-t197+t211-t212+t213+t214-t203;
482  values[24] = 50.0*t32-80.0*t36-40.0*t48+64.0*t52-10.0*xi+16.0*t12;
483  values[25] = -t209+t210-t197+t211+t212-t213+t214-t203;
484  values[26] = t224;
485  values[27] = t225;
486  values[28] = t226;
487 }
488 
489 // values of the derivatives in eta direction
490 static void C_Q_UL5SE_2D_DeriveEta(double xi, double eta, double *values)
491 {
492  double t1 = 1.0-xi;
493  double t3 = eta*eta;
494  double t4 = t3*eta;
495  double t7 = 5.0/2.0*t4-3.0/2.0*eta;
496  double t8 = t1*eta*t7;
497  double t9 = 125.0/384.0*t8;
498  double t10 = 1.0-t3;
499  double t13 = 15.0/2.0*t3-3.0/2.0;
500  double t14 = t1*t10*t13;
501  double t15 = 125.0/768.0*t14;
502  double t16 = xi*xi;
503  double t18 = t16*xi;
504  double t22 = (1.0-t16)*(5.0/2.0*t18-3.0/2.0*xi);
505  double t23 = 125.0/768.0*t22;
506  double t24 = 191.0/768.0*xi;
507  double t25 = 359.0/6144.0*t16;
508  double t26 = t18/96.0;
509  double t27 = t16*t16;
510  double t28 = 39.0/256.0*t27;
511  double t29 = xi*eta;
512  double t30 = 359.0/3072.0*t29;
513  double t31 = t16*eta;
514  double t32 = 31.0/192.0*t31;
515  double t33 = t18*eta;
516  double t34 = 453.0/1024.0*t33;
517  double t35 = t27*eta;
518  double t36 = 371.0/768.0*t35;
519  double t37 = -t9+t15-t23-t24-t25+t26+t28-t30+t32+t34-t36;
520  double t38 = xi*t3;
521  double t39 = t38/32.0;
522  double t40 = t16*t3;
523  double t41 = 1359.0/2048.0*t40;
524  double t42 = t18*t3;
525  double t43 = 367.0/256.0*t42;
526  double t44 = t27*t3;
527  double t45 = 859.0/512.0*t44;
528  double t46 = xi*t4;
529  double t47 = 39.0/64.0*t46;
530  double t48 = t16*t4;
531  double t49 = 371.0/384.0*t48;
532  double t50 = t18*t4;
533  double t51 = 859.0/384.0*t50;
534  double t52 = t27*t4;
535  double t53 = 41.0/16.0*t52;
536  double t54 = eta/256.0;
537  double t55 = 923.0/2048.0*t3;
538  double t56 = t4/32.0;
539  double t57 = 887.0/6144.0+t39+t41+t43-t45+t47-t49-t51+t53-t54-t55+t56;
540  double t59 = 875.0/1536.0*xi;
541  double t60 = 325.0/768.0*t16;
542  double t61 = 875.0/1536.0*t18;
543  double t62 = 625.0/1536.0*t27;
544  double t63 = 875.0/3072.0*t29;
545  double t64 = 325.0/384.0*t31;
546  double t65 = 875.0/3072.0*t33;
547  double t66 = 625.0/768.0*t35;
548  double t67 = 875.0/512.0*t38;
549  double t69 = 325.0/64.0*t40;
550  double t70 = 875.0/512.0*t42;
551  double t71 = 625.0/128.0*t44;
552  double t72 = 875.0/384.0*t46;
553  double t73 = 325.0/48.0*t48;
554  double t74 = 875.0/384.0*t50;
555  double t75 = 625.0/96.0*t52;
556  double t76 = 25.0/768.0*eta;
557  double t77 = 25.0/128.0*t3;
558  double t78 = 25.0/96.0*t4;
559  double t79 = 625.0/768.0*t22;
560  double t80 = -t69-t70+t71-t72+t73+t74-t75+t76+t77-t78+t79;
561  double t82 = 3625.0/2304.0*xi;
562  double t83 = 425.0/1152.0*t16;
563  double t84 = 3625.0/2304.0*t18;
564  double t85 = 625.0/2304.0*t27;
565  double t86 = 1375.0/4608.0*t29;
566  double t87 = 425.0/576.0*t31;
567  double t88 = 1375.0/4608.0*t33;
568  double t89 = 625.0/1152.0*t35;
569  double t90 = 1375.0/768.0*t38;
570  double t92 = 425.0/96.0*t40;
571  double t93 = 1375.0/768.0*t42;
572  double t94 = 625.0/192.0*t44;
573  double t95 = 1375.0/576.0*t46;
574  double t96 = 425.0/72.0*t48;
575  double t97 = 1375.0/576.0*t50;
576  double t98 = 625.0/144.0*t52;
577  double t99 = 25.0/128.0*eta;
578  double t100 = 75.0/64.0*t3;
579  double t101 = 25.0/16.0*t4;
580  double t102 = 625.0/384.0*t22;
581  double t103 = t92-t93-t94-t95-t96+t97+t98-t99-t100+t101-t102;
582  double t106 = t92+t93-t94+t95-t96-t97+t98-t99-t100+t101+t102;
583  double t109 = -t69+t70+t71+t72+t73-t74-t75+t76+t77-t78-t79;
584  double t111 = 1.0+xi;
585  double t113 = t111*eta*t7;
586  double t114 = 125.0/384.0*t113;
587  double t116 = t111*t10*t13;
588  double t117 = 125.0/768.0*t116;
589  double t118 = t23-t114+t117+t24-t25-t26+t28+t30+t32-t34-t36;
590  double t119 = 887.0/6144.0-t39+t41-t43-t45-t47-t49+t51+t53-t54-t55+t56;
591  double t121 = 875.0/6144.0*t16;
592  double t122 = 875.0/1536.0*t27;
593  double t123 = 325.0/384.0*t29;
594  double t124 = 325.0/96.0*t33;
595  double t125 = 325.0/96.0*t35;
596  double t126 = 875.0/2048.0*t40;
597  double t127 = 875.0/512.0*t44;
598  double t128 = 625.0/384.0*t46;
599  double t129 = 625.0/384.0*t48;
600  double t130 = 625.0/96.0*t50;
601  double t131 = 4375.0/2048.0*t3;
602  double t132 = 625.0/384.0*t113;
603  double t133 = 625.0/768.0*t116;
604  double t134 = -4375.0/6144.0-t59+t121-t61-t122-t123-t64+t124+t125+t67-t126+t70+t127+t128+t129-t130-t75+t131+t132-t133;
605  double t135 = 1375.0/9216.0*t16;
606  double t136 = 1375.0/2304.0*t18;
607  double t137 = 1375.0/2304.0*t27;
608  double t138 = 425.0/576.0*t29;
609  double t139 = 425.0/144.0*t33;
610  double t140 = 425.0/144.0*t35;
611  double t141 = 3625.0/768.0*t38;
612  double t142 = 1375.0/3072.0*t40;
613  double t143 = 1375.0/768.0*t44;
614  double t144 = 625.0/576.0*t46;
615  double t145 = 625.0/576.0*t48;
616  double t146 = 625.0/144.0*t50;
617  double t147 = 4375.0/1024.0*t3;
618  double t148 = 625.0/192.0*t113;
619  double t149 = 625.0/384.0*t116;
620  double t150 = 4375.0/3072.0+t82+t135-t136-t137+t138+t87-t139-t140-t141-t142+t93+t143-t144-t145+t146+t98-t147-t148+t149;
621  double t151 = -4375.0/3072.0-t82-t135+t136+t137+t138+t87-t139-t140+t141+t142-t93-t143-t144-t145+t146+t98+t147+t148-t149;
622  double t152 = 4375.0/6144.0+t59-t121+t61+t122-t123-t64+t124+t125-t67+t126-t70-t127+t128+t129-t130-t75-t131-t132+t133;
623  double t153 = -t23+t114-t117-t24+t25+t26-t28+t30+t32-t34-t36;
624  double t154 = -887.0/6144.0+t39-t41+t43+t45-t47-t49+t51+t53-t54+t55+t56;
625  double t157 = t69-t70-t71+t72+t73-t74-t75+t76-t77-t78+t79;
626  double t160 = -t92-t93+t94+t95-t96-t97+t98-t99+t100+t101-t102;
627  double t163 = -t92+t93+t94-t95-t96+t97+t98-t99+t100+t101+t102;
628  double t166 = t69+t70-t71-t72+t73+t74-t75+t76-t77-t78-t79;
629  double t168 = t9-t15+t23+t24+t25-t26-t28-t30+t32+t34-t36;
630  double t169 = -887.0/6144.0-t39-t41-t43+t45+t47-t49-t51+t53-t54+t55+t56;
631  double t171 = 625.0/384.0*t8;
632  double t172 = 625.0/768.0*t14;
633  double t173 = 4375.0/6144.0-t59-t121-t61+t122+t123-t64-t124+t125+t67+t126+t70-t127-t128+t129+t130-t75-t131-t171+t172;
634  double t174 = 625.0/192.0*t8;
635  double t175 = 625.0/384.0*t14;
636  double t176 = -4375.0/3072.0+t82-t135-t136+t137-t138+t87+t139-t140-t141+t142+t93-t143+t144-t145-t146+t98+t147+t174-t175;
637  double t177 = 4375.0/3072.0-t82+t135+t136-t137-t138+t87+t139-t140+t141-t142-t93+t143+t144-t145-t146+t98-t147-t174+t175;
638  double t178 = -4375.0/6144.0+t59+t121+t61-t122+t123-t64-t124+t125-t67-t126-t70+t127-t128+t129+t130-t75+t131+t171-t172;
639  double t179 = 16.0/9.0*xi;
640  double t180 = 32.0/9.0*t16;
641  double t181 = 16.0/9.0*t18;
642  double t182 = 32.0/9.0*t27;
643  double t183 = 64.0/9.0*t29;
644  double t184 = 128.0/9.0*t31;
645  double t185 = 64.0/9.0*t33;
646  double t186 = 128.0/9.0*t35;
647  double t187 = 16.0/3.0*t38;
648  double t188 = 32.0/3.0*t40;
649  double t189 = 16.0/3.0*t42;
650  double t190 = 32.0/3.0*t44;
651  double t191 = 128.0/9.0*t46;
652  double t192 = 256.0/9.0*t48;
653  double t193 = 128.0/9.0*t50;
654  double t194 = 256.0/9.0*t52;
655  double t195 = t179-t180-t181+t182-t183+t184+t185-t186-t187+t188+t189-t190+t191-t192-t193+t194;
656  double t196 = 40.0/3.0*t29;
657  double t197 = 80.0/3.0*t31;
658  double t198 = 40.0/3.0*t33;
659  double t199 = 80.0/3.0*t35;
660  double t200 = 64.0/3.0*t46;
661  double t201 = 128.0/3.0*t48;
662  double t202 = 64.0/3.0*t50;
663  double t203 = 128.0/3.0*t52;
664  double t205 = -t179+t180+t181-t182-t183+t184+t185-t186+t187-t188-t189+t190+t191-t192-t193+t194;
665  double t206 = 20.0/3.0*t16;
666  double t207 = 16.0/3.0*t27;
667  double t208 = 64.0/3.0*t35;
668  double t209 = 20.0*t40;
669  double t210 = 16.0*t44;
670  double t211 = 160.0/3.0*t48;
671  double t212 = 16.0/3.0*eta;
672  double t213 = 4.0*t3;
673  double t214 = 32.0/3.0*t4;
674  double t215 = -4.0/3.0+t206-t207-t197+t208-t209+t210+t211-t203+t212+t213-t214;
675  double t223 = 4.0/3.0-t206+t207-t197+t208+t209-t210+t211-t203+t212-t213-t214;
676  double t224 = -t179-t180+t181+t182+t183+t184-t185-t186+t187+t188-t189-t190-t191-t192+t193+t194;
677  double t226 = t179+t180-t181-t182+t183+t184-t185-t186-t187-t188+t189+t190-t191-t192+t193+t194;
678  values[0] = t37+t57;
679  values[1] = -25.0/1536.0+t59+t60-t61-t62+t63-t64-t65+t66+t67+t80;
680  values[2] = 25.0/256.0-t82-t83+t84+t85+t86+t87-t88-t89+t90+t103;
681  values[3] = 25.0/256.0+t82-t83-t84+t85-t86+t87+t88-t89-t90+t106;
682  values[4] = -25.0/1536.0-t59+t60+t61-t62-t63-t64+t65+t66-t67+t109;
683  values[5] = t118+t119;
684  values[6] = t134;
685  values[7] = t150;
686  values[8] = t151;
687  values[9] = t152;
688  values[10] = t153+t154;
689  values[11] = 25.0/1536.0+t59-t60-t61+t62-t63-t64+t65+t66+t67+t157;
690  values[12] = -25.0/256.0-t82+t83+t84-t85-t86+t87+t88-t89+t90+t160;
691  values[13] = -25.0/256.0+t82+t83-t84-t85+t86+t87-t88-t89-t90+t163;
692  values[14] = 25.0/1536.0-t59-t60+t61+t62+t63-t64-t65+t66-t67+t166;
693  values[15] = t168+t169;
694  values[16] = t173;
695  values[17] = t176;
696  values[18] = t177;
697  values[19] = t178;
698  values[20] = t195;
699  values[21] = t196-t197-t198+t199-t200+t201+t202-t203;
700  values[22] = t205;
701  values[23] = t215;
702  values[24] = 50.0*t31-40.0*t35-80.0*t48+64.0*t52-10.0*eta+16.0*t4;
703  values[25] = t223;
704  values[26] = t224;
705  values[27] = -t196-t197+t198+t199+t200+t201-t202-t203;
706  values[28] = t226;
707 }
708 
709 // values of the derivatives in xi-xi direction
710 static void C_Q_UL5SE_2D_DeriveXiXi(double xi, double eta, double *values)
711 {
712  double t1 = 1.0-eta;
713  double t2 = xi*xi;
714  double t6 = 5.0/2.0*t2*xi-3.0/2.0*xi;
715  double t7 = t1*t6;
716  double t8 = 125.0/384.0*t7;
717  double t11 = 15.0/2.0*t2-3.0/2.0;
718  double t12 = t1*xi*t11;
719  double t13 = 125.0/192.0*t12;
720  double t14 = 1.0-t2;
721  double t16 = t1*t14*xi;
722  double t17 = 625.0/256.0*t16;
723  double t18 = 359.0/3072.0*eta;
724  double t19 = xi*eta;
725  double t20 = t19/16.0;
726  double t21 = t2*eta;
727  double t22 = 117.0/64.0*t21;
728  double t23 = eta*eta;
729  double t24 = 31.0/192.0*t23;
730  double t25 = xi*t23;
731  double t26 = 1359.0/1024.0*t25;
732  double t27 = t2*t23;
733  double t28 = 371.0/128.0*t27;
734  double t29 = t23*eta;
735  double t30 = 453.0/1024.0*t29;
736  double t31 = xi*t29;
737  double t32 = 367.0/128.0*t31;
738  double t33 = t2*t29;
739  double t34 = 859.0/128.0*t33;
740  double t35 = t23*t23;
741  double t36 = 371.0/768.0*t35;
742  double t37 = xi*t35;
743  double t38 = 859.0/256.0*t37;
744  double t39 = t2*t35;
745  double t40 = 123.0/16.0*t39;
746  double t41 = 923.0/1024.0*xi;
747  double t42 = 3.0/32.0*t2;
748  double t43 = -t8-t13+t17-t18+t20+t22+t24+t26-t28+t30+t32-t34-t36-t38+t40-1.0/256.0-t41+t42;
749  double t44 = 325.0/384.0*eta;
750  double t45 = 875.0/256.0*t19;
751  double t46 = 625.0/128.0*t21;
752  double t47 = 325.0/384.0*t23;
753  double t48 = 875.0/1024.0*t25;
754  double t49 = 625.0/128.0*t27;
755  double t50 = 325.0/96.0*t29;
756  double t51 = 875.0/256.0*t31;
757  double t52 = 625.0/32.0*t33;
758  double t53 = 325.0/96.0*t35;
759  double t54 = 875.0/256.0*t37;
760  double t55 = 625.0/32.0*t39;
761  double t56 = 4375.0/1024.0*xi;
762  double t57 = 625.0/384.0*t7;
763  double t58 = 625.0/192.0*t12;
764  double t59 = 3125.0/256.0*t16;
765  double t60 = t44-t45-t46-t47-t48+t49-t50-t51+t52+t53+t54-t55+t56+t57+t58-t59;
766  double t61 = 425.0/576.0*eta;
767  double t62 = 3625.0/384.0*t19;
768  double t63 = 625.0/192.0*t21;
769  double t64 = 425.0/576.0*t23;
770  double t65 = 1375.0/1536.0*t25;
771  double t66 = 625.0/192.0*t27;
772  double t67 = 425.0/144.0*t29;
773  double t68 = 1375.0/384.0*t31;
774  double t69 = 625.0/48.0*t33;
775  double t70 = 425.0/144.0*t35;
776  double t71 = 1375.0/384.0*t37;
777  double t72 = 625.0/48.0*t39;
778  double t73 = 4375.0/512.0*xi;
779  double t74 = 625.0/192.0*t7;
780  double t75 = 625.0/96.0*t12;
781  double t76 = 3125.0/128.0*t16;
782  double t77 = -t61+t62+t63+t64-t65-t66+t67-t68-t69-t70+t71+t72-t73-t74-t75+t76;
783  double t78 = -t61-t62+t63+t64+t65-t66+t67+t68-t69-t70-t71+t72+t73+t74+t75-t76;
784  double t79 = t44+t45-t46-t47+t48+t49-t50+t51+t52+t53-t54-t55-t56-t57-t58+t59;
785  double t80 = t8+t13-t17-t18-t20+t22+t24-t26-t28+t30-t32-t34-t36+t38+t40-1.0/256.0+t41+t42;
786  double t81 = 875.0/3072.0*eta;
787  double t82 = 875.0/128.0*t21;
788  double t83 = 325.0/32.0*t25;
789  double t84 = 325.0/16.0*t27;
790  double t85 = 875.0/3072.0*t29;
791  double t86 = 875.0/128.0*t33;
792  double t87 = 625.0/768.0*t35;
793  double t88 = 625.0/64.0*t37;
794  double t89 = 25.0/64.0*xi;
795  double t90 = 25.0/32.0*t2;
796  double t91 = t81-t45-t82-t47+t83+t84-t85+t51+t86+t87-t88-t55+25.0/768.0-t89-t90;
797  double t92 = 1375.0/4608.0*eta;
798  double t93 = 1375.0/384.0*t19;
799  double t94 = 1375.0/192.0*t21;
800  double t95 = 425.0/48.0*t25;
801  double t96 = 425.0/24.0*t27;
802  double t97 = 1375.0/4608.0*t29;
803  double t98 = 1375.0/192.0*t33;
804  double t99 = 625.0/1152.0*t35;
805  double t100 = 625.0/96.0*t37;
806  double t101 = 75.0/32.0*xi;
807  double t102 = 75.0/16.0*t2;
808  double t103 = t92-t93-t94+t64-t95-t96-t97+t68+t98-t99+t100+t72-25.0/128.0+t101+t102;
809  double t104 = -t92+t93+t94+t64-t95-t96+t97-t68-t98-t99+t100+t72-25.0/128.0+t101+t102;
810  double t105 = -t81+t45+t82-t47+t83+t84+t85-t51-t86+t87-t88-t55+25.0/768.0-t89-t90;
811  double t106 = 1.0+eta;
812  double t107 = t106*t6;
813  double t108 = 125.0/384.0*t107;
814  double t110 = t106*xi*t11;
815  double t111 = 125.0/192.0*t110;
816  double t113 = t106*t14*xi;
817  double t114 = 625.0/256.0*t113;
818  double t115 = t108+t111-t114+t18+t20-t22+t24-t26-t28-t30+t32+t34-t36+t38+t40-1.0/256.0+t41+t42;
819  double t116 = 625.0/384.0*t107;
820  double t117 = 625.0/192.0*t110;
821  double t118 = 3125.0/256.0*t113;
822  double t119 = -t44-t45+t46-t47+t48+t49+t50-t51-t52+t53-t54-t55-t56-t116-t117+t118;
823  double t120 = 625.0/192.0*t107;
824  double t121 = 625.0/96.0*t110;
825  double t122 = 3125.0/128.0*t113;
826  double t123 = t61+t62-t63+t64+t65-t66-t67-t68+t69-t70-t71+t72+t73+t120+t121-t122;
827  double t124 = t61-t62-t63+t64-t65-t66-t67+t68+t69-t70+t71+t72-t73-t120-t121+t122;
828  double t125 = -t44+t45+t46-t47-t48+t49+t50+t51-t52+t53+t54-t55+t56+t116+t117-t118;
829  double t126 = -t108-t111+t114+t18-t20-t22+t24+t26-t28-t30-t32+t34-t36-t38+t40-1.0/256.0-t41+t42;
830  double t127 = -t81-t45+t82-t47-t83+t84+t85+t51-t86+t87+t88-t55+25.0/768.0+t89-t90;
831  double t128 = -t92-t93+t94+t64+t95-t96+t97+t68-t98-t99-t100+t72-25.0/128.0-t101+t102;
832  double t129 = t92+t93-t94+t64+t95-t96-t97-t68+t98-t99-t100+t72-25.0/128.0-t101+t102;
833  double t130 = t81+t45-t82-t47-t83+t84-t85-t51+t86+t87+t88-t55+25.0/768.0+t89-t90;
834  double t131 = 64.0/9.0*eta;
835  double t132 = 32.0/3.0*t19;
836  double t133 = 128.0/3.0*t21;
837  double t134 = 128.0/9.0*t23;
838  double t135 = 64.0/3.0*t25;
839  double t136 = 256.0/3.0*t27;
840  double t137 = 64.0/9.0*t29;
841  double t138 = 32.0/3.0*t31;
842  double t139 = 128.0/3.0*t33;
843  double t140 = 128.0/9.0*t35;
844  double t141 = 64.0/3.0*t37;
845  double t142 = 256.0/3.0*t39;
846  double t143 = -t131-t132+t133+t134+t135-t136+t137+t138-t139-t140-t141+t142;
847  double t144 = 80.0/3.0*t23;
848  double t145 = 40.0*t25;
849  double t146 = 160.0*t27;
850  double t147 = 64.0/3.0*t35;
851  double t148 = 32.0*t37;
852  double t149 = 128.0*t39;
853  double t150 = 8.0*xi;
854  double t151 = 32.0*t2;
855  double t153 = t131+t132-t133+t134+t135-t136-t137-t138+t139-t140-t141+t142;
856  double t154 = 40.0/3.0*eta;
857  double t155 = 64.0*t21;
858  double t156 = 128.0*t27;
859  double t157 = 40.0/3.0*t29;
860  double t158 = 64.0*t33;
861  double t159 = 80.0/3.0*t35;
862  double t168 = -t131+t132+t133+t134-t135-t136+t137-t138-t139-t140+t141+t142;
863  double t170 = t131-t132-t133+t134-t135-t136-t137+t138+t139-t140+t141+t142;
864  values[0] = t43;
865  values[1] = t60;
866  values[2] = t77;
867  values[3] = t78;
868  values[4] = t79;
869  values[5] = t80;
870  values[6] = t91;
871  values[7] = t103;
872  values[8] = t104;
873  values[9] = t105;
874  values[10] = t115;
875  values[11] = t119;
876  values[12] = t123;
877  values[13] = t124;
878  values[14] = t125;
879  values[15] = t126;
880  values[16] = t127;
881  values[17] = t128;
882  values[18] = t129;
883  values[19] = t130;
884  values[20] = t143;
885  values[21] = -t144-t145+t146+t147+t148-t149+16.0/3.0+t150-t151;
886  values[22] = t153;
887  values[23] = t154-t155-t144+t156-t157+t158+t159-t149;
888  values[24] = 50.0*t23-240.0*t27-40.0*t35+192.0*t39-10.0+48.0*t2;
889  values[25] = -t154+t155-t144+t156+t157-t158+t159-t149;
890  values[26] = t168;
891  values[27] = -t144+t145+t146+t147-t148-t149+16.0/3.0-t150-t151;
892  values[28] = t170;
893 }
894 
895 // values of the derivatives in xi-eta direction
896 static void C_Q_UL5SE_2D_DeriveXiEta(double xi, double eta, double *values)
897 {
898  double t1 = eta*eta;
899  double t2 = t1*eta;
900  double t6 = eta*(5.0/2.0*t2-3.0/2.0*eta);
901  double t7 = 125.0/384.0*t6;
902  double t11 = (1.0-t1)*(15.0/2.0*t1-3.0/2.0);
903  double t12 = 125.0/768.0*t11;
904  double t13 = xi*xi;
905  double t14 = t13*xi;
906  double t18 = xi*(5.0/2.0*t14-3.0/2.0*xi);
907  double t19 = 125.0/384.0*t18;
908  double t23 = (1.0-t13)*(15.0/2.0*t13-3.0/2.0);
909  double t24 = 125.0/768.0*t23;
910  double t25 = 359.0/3072.0*xi;
911  double t26 = t13/32.0;
912  double t27 = 39.0/64.0*t14;
913  double t28 = 359.0/3072.0*eta;
914  double t29 = xi*eta;
915  double t30 = 31.0/96.0*t29;
916  double t31 = t13*eta;
917  double t32 = 1359.0/1024.0*t31;
918  double t33 = t14*eta;
919  double t34 = 371.0/192.0*t33;
920  double t35 = t1/32.0;
921  double t36 = xi*t1;
922  double t37 = 1359.0/1024.0*t36;
923  double t38 = t13*t1;
924  double t39 = 1101.0/256.0*t38;
925  double t40 = t14*t1;
926  double t41 = 859.0/128.0*t40;
927  double t42 = 39.0/64.0*t2;
928  double t43 = t2*xi;
929  double t44 = 371.0/192.0*t43;
930  double t45 = t13*t2;
931  double t46 = 859.0/128.0*t45;
932  double t47 = t14*t2;
933  double t48 = 41.0/4.0*t47;
934  double t49 = t7-t12+t19-t24-191.0/768.0-t25+t26+t27-t28+t30+t32-t34+t35+t37+t39-t41+t42-t44-t46+t48;
935  double t50 = 325.0/384.0*xi;
936  double t51 = 875.0/512.0*t13;
937  double t52 = 625.0/384.0*t14;
938  double t53 = 875.0/3072.0*eta;
939  double t54 = 325.0/192.0*t29;
940  double t55 = 875.0/1024.0*t31;
941  double t56 = 625.0/192.0*t33;
942  double t57 = 875.0/512.0*t1;
943  double t58 = 325.0/32.0*t36;
944  double t59 = 2625.0/512.0*t38;
945  double t60 = 625.0/32.0*t40;
946  double t61 = 875.0/384.0*t2;
947  double t62 = 325.0/24.0*t43;
948  double t63 = 875.0/128.0*t45;
949  double t64 = 625.0/24.0*t47;
950  double t65 = 625.0/384.0*t18;
951  double t66 = 625.0/768.0*t23;
952  double t67 = 875.0/1536.0+t50-t51-t52+t53-t54-t55+t56+t57-t58-t59+t60-t61+t62+t63-t64-t65+t66;
953  double t68 = 425.0/576.0*xi;
954  double t69 = 3625.0/768.0*t13;
955  double t70 = 625.0/576.0*t14;
956  double t71 = 1375.0/4608.0*eta;
957  double t72 = 425.0/288.0*t29;
958  double t73 = 1375.0/1536.0*t31;
959  double t74 = 625.0/288.0*t33;
960  double t75 = 1375.0/768.0*t1;
961  double t76 = 425.0/48.0*t36;
962  double t77 = 1375.0/256.0*t38;
963  double t78 = 625.0/48.0*t40;
964  double t79 = 1375.0/576.0*t2;
965  double t80 = 425.0/36.0*t43;
966  double t81 = 1375.0/192.0*t45;
967  double t82 = 625.0/36.0*t47;
968  double t83 = 625.0/192.0*t18;
969  double t84 = 625.0/384.0*t23;
970  double t85 = -3625.0/2304.0-t68+t69+t70+t71+t72-t73-t74+t75+t76-t77-t78-t79-t80+t81+t82+t83-t84;
971  double t86 = 3625.0/2304.0-t68-t69+t70-t71+t72+t73-t74-t75+t76+t77-t78+t79-t80-t81+t82-t83+t84;
972  double t87 = -875.0/1536.0+t50+t51-t52-t53-t54+t55+t56-t57-t58+t59+t60+t61+t62-t63-t64+t65-t66;
973  double t88 = -t19+t24-t7+t12+191.0/768.0-t25-t26+t27+t28+t30-t32-t34-t35+t37-t39-t41-t42-t44+t46+t48;
974  double t89 = 875.0/3072.0*xi;
975  double t90 = 875.0/384.0*t14;
976  double t91 = 325.0/384.0*eta;
977  double t92 = 325.0/32.0*t31;
978  double t93 = 325.0/24.0*t33;
979  double t94 = 875.0/1024.0*t36;
980  double t95 = 875.0/128.0*t40;
981  double t96 = 625.0/384.0*t2;
982  double t97 = 625.0/192.0*t43;
983  double t98 = 625.0/32.0*t45;
984  double t99 = 625.0/384.0*t6;
985  double t100 = 625.0/768.0*t11;
986  double t101 = -875.0/1536.0+t89-t51-t90-t91-t54+t92+t93+t57-t94+t59+t95+t96+t97-t98-t64+t99-t100;
987  double t102 = 1375.0/4608.0*xi;
988  double t103 = 1375.0/768.0*t13;
989  double t104 = 1375.0/576.0*t14;
990  double t105 = 425.0/576.0*eta;
991  double t106 = 425.0/48.0*t31;
992  double t107 = 425.0/36.0*t33;
993  double t108 = 3625.0/768.0*t1;
994  double t109 = 1375.0/1536.0*t36;
995  double t110 = 1375.0/192.0*t40;
996  double t111 = 625.0/576.0*t2;
997  double t112 = 625.0/288.0*t43;
998  double t113 = 625.0/48.0*t45;
999  double t114 = 625.0/192.0*t6;
1000  double t115 = 625.0/384.0*t11;
1001  double t116 = 3625.0/2304.0+t102-t103-t104+t105+t72-t106-t107-t108-t109+t77+t110-t111-t112+t113+t82-t114+t115;
1002  double t117 = -3625.0/2304.0-t102+t103+t104+t105+t72-t106-t107+t108+t109-t77-t110-t111-t112+t113+t82+t114-t115;
1003  double t118 = 875.0/1536.0-t89+t51+t90-t91-t54+t92+t93-t57+t94-t59-t95+t96+t97-t98-t64-t99+t100;
1004  double t119 = t19-t24+t7-t12-191.0/768.0+t25+t26-t27+t28+t30-t32-t34+t35-t37+t39+t41-t42-t44+t46+t48;
1005  double t120 = 875.0/1536.0-t50-t51+t52-t53-t54+t55+t56+t57+t58-t59-t60+t61+t62-t63-t64-t65+t66;
1006  double t121 = -3625.0/2304.0+t68+t69-t70-t71+t72+t73-t74+t75-t76-t77+t78+t79-t80-t81+t82+t83-t84;
1007  double t122 = 3625.0/2304.0+t68-t69-t70+t71+t72-t73-t74-t75-t76+t77+t78-t79-t80+t81+t82-t83+t84;
1008  double t123 = -875.0/1536.0-t50+t51+t52+t53-t54-t55+t56-t57+t58+t59-t60-t61+t62+t63-t64+t65-t66;
1009  double t124 = -t7+t12-t19+t24+191.0/768.0+t25-t26-t27-t28+t30+t32-t34-t35-t37-t39+t41+t42-t44-t46+t48;
1010  double t125 = -875.0/1536.0-t89-t51+t90+t91-t54-t92+t93+t57+t94+t59-t95-t96+t97+t98-t64+t99-t100;
1011  double t126 = 3625.0/2304.0-t102-t103+t104-t105+t72+t106-t107-t108+t109+t77-t110+t111-t112-t113+t82-t114+t115;
1012  double t127 = -3625.0/2304.0+t102+t103-t104-t105+t72+t106-t107+t108-t109-t77+t110+t111-t112-t113+t82+t114-t115;
1013  double t128 = 875.0/1536.0+t89+t51-t90+t91-t54-t92+t93-t57-t94-t59+t95-t96+t97+t98-t64-t99+t100;
1014  double t129 = 64.0/9.0*xi;
1015  double t130 = 16.0/3.0*t13;
1016  double t131 = 128.0/9.0*t14;
1017  double t132 = 64.0/9.0*eta;
1018  double t133 = 256.0/9.0*t29;
1019  double t134 = 64.0/3.0*t31;
1020  double t135 = 512.0/9.0*t33;
1021  double t136 = 16.0/3.0*t1;
1022  double t137 = 64.0/3.0*t36;
1023  double t138 = 16.0*t38;
1024  double t139 = 128.0/3.0*t40;
1025  double t140 = 128.0/9.0*t2;
1026  double t141 = 512.0/9.0*t43;
1027  double t142 = 128.0/3.0*t45;
1028  double t143 = 1024.0/9.0*t47;
1029  double t144 = 16.0/9.0-t129-t130+t131-t132+t133+t134-t135-t136+t137+t138-t139+t140-t141-t142+t143;
1030  double t145 = 40.0/3.0*eta;
1031  double t146 = 160.0/3.0*t29;
1032  double t147 = 40.0*t31;
1033  double t148 = 320.0/3.0*t33;
1034  double t149 = 64.0/3.0*t2;
1035  double t150 = 256.0/3.0*t43;
1036  double t151 = 64.0*t45;
1037  double t152 = 512.0/3.0*t47;
1038  double t154 = -16.0/9.0+t129+t130-t131-t132+t133+t134-t135+t136-t137-t138+t139+t140-t141-t142+t143;
1039  double t155 = 40.0/3.0*xi;
1040  double t156 = 64.0/3.0*t14;
1041  double t157 = 256.0/3.0*t33;
1042  double t158 = 40.0*t36;
1043  double t159 = 64.0*t40;
1044  double t160 = 320.0/3.0*t43;
1045  double t168 = -16.0/9.0-t129+t130+t131+t132+t133-t134-t135+t136+t137-t138-t139-t140-t141+t142+t143;
1046  double t170 = 16.0/9.0+t129-t130-t131+t132+t133-t134-t135-t136-t137+t138+t139-t140-t141+t142+t143;
1047  values[0] = t49;
1048  values[1] = t67;
1049  values[2] = t85;
1050  values[3] = t86;
1051  values[4] = t87;
1052  values[5] = t88;
1053  values[6] = t101;
1054  values[7] = t116;
1055  values[8] = t117;
1056  values[9] = t118;
1057  values[10] = t119;
1058  values[11] = t120;
1059  values[12] = t121;
1060  values[13] = t122;
1061  values[14] = t123;
1062  values[15] = t124;
1063  values[16] = t125;
1064  values[17] = t126;
1065  values[18] = t127;
1066  values[19] = t128;
1067  values[20] = t144;
1068  values[21] = t145-t146-t147+t148-t149+t150+t151-t152;
1069  values[22] = t154;
1070  values[23] = t155-t156-t146+t157-t158+t159+t160-t152;
1071  values[24] = 100.0*t29-160.0*t33-160.0*t43+256.0*t47;
1072  values[25] = -t155+t156-t146+t157+t158-t159+t160-t152;
1073  values[26] = t168;
1074  values[27] = -t145-t146+t147+t148+t149+t150-t151-t152;
1075  values[28] = t170;
1076 }
1077 
1078 // values of the derivatives in eta-eta direction
1079 static void C_Q_UL5SE_2D_DeriveEtaEta(double xi, double eta, double *values)
1080 {
1081  double t1 = 1.0-xi;
1082  double t2 = eta*eta;
1083  double t6 = 5.0/2.0*t2*eta-3.0/2.0*eta;
1084  double t7 = t1*t6;
1085  double t8 = 125.0/384.0*t7;
1086  double t11 = 15.0/2.0*t2-3.0/2.0;
1087  double t12 = t1*eta*t11;
1088  double t13 = 125.0/192.0*t12;
1089  double t14 = 1.0-t2;
1090  double t16 = t14*t1*eta;
1091  double t17 = 625.0/256.0*t16;
1092  double t18 = 359.0/3072.0*xi;
1093  double t19 = xi*xi;
1094  double t20 = 31.0/192.0*t19;
1095  double t21 = t19*xi;
1096  double t22 = 453.0/1024.0*t21;
1097  double t23 = t19*t19;
1098  double t24 = 371.0/768.0*t23;
1099  double t25 = xi*eta;
1100  double t26 = t25/16.0;
1101  double t27 = t19*eta;
1102  double t28 = 1359.0/1024.0*t27;
1103  double t29 = t21*eta;
1104  double t30 = 367.0/128.0*t29;
1105  double t31 = t23*eta;
1106  double t32 = 859.0/256.0*t31;
1107  double t33 = t2*xi;
1108  double t34 = 117.0/64.0*t33;
1109  double t35 = t19*t2;
1110  double t36 = 371.0/128.0*t35;
1111  double t37 = t21*t2;
1112  double t38 = 859.0/128.0*t37;
1113  double t39 = t2*t23;
1114  double t40 = 123.0/16.0*t39;
1115  double t41 = 923.0/1024.0*eta;
1116  double t42 = 3.0/32.0*t2;
1117  double t43 = -t8-t13+t17-t18+t20+t22-t24+t26+t28+t30-t32+t34-t36-t38+t40-1.0/256.0-t41+t42;
1118  double t44 = 875.0/3072.0*xi;
1119  double t45 = 325.0/384.0*t19;
1120  double t46 = 875.0/3072.0*t21;
1121  double t47 = 625.0/768.0*t23;
1122  double t48 = 875.0/256.0*t25;
1123  double t49 = 325.0/32.0*t27;
1124  double t50 = 875.0/256.0*t29;
1125  double t51 = 625.0/64.0*t31;
1126  double t52 = 875.0/128.0*t33;
1127  double t53 = 325.0/16.0*t35;
1128  double t54 = 875.0/128.0*t37;
1129  double t55 = 625.0/32.0*t39;
1130  double t56 = 25.0/64.0*eta;
1131  double t57 = 25.0/32.0*t2;
1132  double t58 = t44-t45-t46+t47+t48-t49-t50+t51-t52+t53+t54-t55+25.0/768.0+t56-t57;
1133  double t59 = 1375.0/4608.0*xi;
1134  double t60 = 425.0/576.0*t19;
1135  double t61 = 1375.0/4608.0*t21;
1136  double t62 = 625.0/1152.0*t23;
1137  double t63 = 1375.0/384.0*t25;
1138  double t64 = 425.0/48.0*t27;
1139  double t65 = 1375.0/384.0*t29;
1140  double t66 = 625.0/96.0*t31;
1141  double t67 = 1375.0/192.0*t33;
1142  double t68 = 425.0/24.0*t35;
1143  double t69 = 1375.0/192.0*t37;
1144  double t70 = 625.0/48.0*t39;
1145  double t71 = 75.0/32.0*eta;
1146  double t72 = 75.0/16.0*t2;
1147  double t73 = t59+t60-t61-t62+t63+t64-t65-t66-t67-t68+t69+t70-25.0/128.0-t71+t72;
1148  double t74 = -t59+t60+t61-t62-t63+t64+t65-t66+t67-t68-t69+t70-25.0/128.0-t71+t72;
1149  double t75 = -t44-t45+t46+t47-t48-t49+t50+t51+t52+t53-t54-t55+25.0/768.0+t56-t57;
1150  double t76 = 1.0+xi;
1151  double t77 = t76*t6;
1152  double t78 = 125.0/384.0*t77;
1153  double t80 = t76*eta*t11;
1154  double t81 = 125.0/192.0*t80;
1155  double t83 = t76*t14*eta;
1156  double t84 = 625.0/256.0*t83;
1157  double t85 = -t78-t81+t84+t18+t20-t22-t24-t26+t28-t30-t32-t34-t36+t38+t40-1.0/256.0-t41+t42;
1158  double t86 = 325.0/384.0*xi;
1159  double t87 = 325.0/96.0*t21;
1160  double t88 = 325.0/96.0*t23;
1161  double t89 = 875.0/1024.0*t27;
1162  double t90 = 875.0/256.0*t31;
1163  double t91 = 625.0/128.0*t33;
1164  double t92 = 625.0/128.0*t35;
1165  double t93 = 625.0/32.0*t37;
1166  double t94 = 4375.0/1024.0*eta;
1167  double t95 = 625.0/384.0*t77;
1168  double t96 = 625.0/192.0*t80;
1169  double t97 = 3125.0/256.0*t83;
1170  double t98 = -t86-t45+t87+t88+t48-t89+t50+t90+t91+t92-t93-t55+t94+t95+t96-t97;
1171  double t99 = 425.0/576.0*xi;
1172  double t100 = 425.0/144.0*t21;
1173  double t101 = 425.0/144.0*t23;
1174  double t102 = 3625.0/384.0*t25;
1175  double t103 = 1375.0/1536.0*t27;
1176  double t104 = 1375.0/384.0*t31;
1177  double t105 = 625.0/192.0*t33;
1178  double t106 = 625.0/192.0*t35;
1179  double t107 = 625.0/48.0*t37;
1180  double t108 = 4375.0/512.0*eta;
1181  double t109 = 625.0/192.0*t77;
1182  double t110 = 625.0/96.0*t80;
1183  double t111 = 3125.0/128.0*t83;
1184  double t112 = t99+t60-t100-t101-t102-t103+t65+t104-t105-t106+t107+t70-t108-t109-t110+t111;
1185  double t113 = t99+t60-t100-t101+t102+t103-t65-t104-t105-t106+t107+t70+t108+t109+t110-t111;
1186  double t114 = -t86-t45+t87+t88-t48+t89-t50-t90+t91+t92-t93-t55-t94-t95-t96+t97;
1187  double t115 = t78+t81-t84+t18+t20-t22-t24+t26-t28+t30+t32-t34-t36+t38+t40-1.0/256.0+t41+t42;
1188  double t116 = -t44-t45+t46+t47+t48+t49-t50-t51+t52+t53-t54-t55+25.0/768.0-t56-t57;
1189  double t117 = -t59+t60+t61-t62+t63-t64-t65+t66+t67-t68-t69+t70-25.0/128.0+t71+t72;
1190  double t118 = t59+t60-t61-t62-t63-t64+t65+t66-t67-t68+t69+t70-25.0/128.0+t71+t72;
1191  double t119 = t44-t45-t46+t47-t48+t49+t50-t51-t52+t53+t54-t55+25.0/768.0-t56-t57;
1192  double t120 = t8+t13-t17-t18+t20+t22-t24-t26-t28-t30+t32+t34-t36-t38+t40-1.0/256.0+t41+t42;
1193  double t121 = 625.0/384.0*t7;
1194  double t122 = 625.0/192.0*t12;
1195  double t123 = 3125.0/256.0*t16;
1196  double t124 = t86-t45-t87+t88+t48+t89+t50-t90-t91+t92+t93-t55-t94-t121-t122+t123;
1197  double t125 = 625.0/192.0*t7;
1198  double t126 = 625.0/96.0*t12;
1199  double t127 = 3125.0/128.0*t16;
1200  double t128 = -t99+t60+t100-t101-t102+t103+t65-t104+t105-t106-t107+t70+t108+t125+t126-t127;
1201  double t129 = -t99+t60+t100-t101+t102-t103-t65+t104+t105-t106-t107+t70-t108-t125-t126+t127;
1202  double t130 = t86-t45-t87+t88-t48-t89-t50+t90-t91+t92+t93-t55+t94+t121+t122-t123;
1203  double t131 = 64.0/9.0*xi;
1204  double t132 = 128.0/9.0*t19;
1205  double t133 = 64.0/9.0*t21;
1206  double t134 = 128.0/9.0*t23;
1207  double t135 = 32.0/3.0*t25;
1208  double t136 = 64.0/3.0*t27;
1209  double t137 = 32.0/3.0*t29;
1210  double t138 = 64.0/3.0*t31;
1211  double t139 = 128.0/3.0*t33;
1212  double t140 = 256.0/3.0*t35;
1213  double t141 = 128.0/3.0*t37;
1214  double t142 = 256.0/3.0*t39;
1215  double t143 = -t131+t132+t133-t134-t135+t136+t137-t138+t139-t140-t141+t142;
1216  double t144 = 40.0/3.0*xi;
1217  double t145 = 80.0/3.0*t19;
1218  double t146 = 40.0/3.0*t21;
1219  double t147 = 80.0/3.0*t23;
1220  double t148 = 64.0*t33;
1221  double t149 = 128.0*t35;
1222  double t150 = 64.0*t37;
1223  double t151 = 128.0*t39;
1224  double t153 = -t131+t132+t133-t134+t135-t136-t137+t138+t139-t140-t141+t142;
1225  double t154 = 64.0/3.0*t23;
1226  double t155 = 40.0*t27;
1227  double t156 = 32.0*t31;
1228  double t157 = 160.0*t35;
1229  double t158 = 8.0*eta;
1230  double t159 = 32.0*t2;
1231  double t168 = t131+t132-t133-t134+t135+t136-t137-t138-t139-t140+t141+t142;
1232  double t170 = t131+t132-t133-t134-t135-t136+t137+t138-t139-t140+t141+t142;
1233  values[0] = t43;
1234  values[1] = t58;
1235  values[2] = t73;
1236  values[3] = t74;
1237  values[4] = t75;
1238  values[5] = t85;
1239  values[6] = t98;
1240  values[7] = t112;
1241  values[8] = t113;
1242  values[9] = t114;
1243  values[10] = t115;
1244  values[11] = t116;
1245  values[12] = t117;
1246  values[13] = t118;
1247  values[14] = t119;
1248  values[15] = t120;
1249  values[16] = t124;
1250  values[17] = t128;
1251  values[18] = t129;
1252  values[19] = t130;
1253  values[20] = t143;
1254  values[21] = t144-t145-t146+t147-t148+t149+t150-t151;
1255  values[22] = t153;
1256  values[23] = -t145+t154-t155+t156+t157-t151+16.0/3.0+t158-t159;
1257  values[24] = 50.0*t19-40.0*t23-240.0*t35+192.0*t39-10.0+48.0*t2;
1258  values[25] = -t145+t154+t155-t156+t157-t151+16.0/3.0-t158-t159;
1259  values[26] = t168;
1260  values[27] = -t144-t145+t146+t147+t148+t149-t150-t151;
1261  values[28] = t170;
1262 }
1263 
1264 // ***********************************************************************
1265 
1266 TBaseFunct2D *BF_C_Q_UL5SE_2D_Obj = new TBaseFunct2D
1267  (29, BF_C_Q_UL5SE_2D, BFUnitSquare,
1268  C_Q_UL5SE_2D_Funct, C_Q_UL5SE_2D_DeriveXi,
1269  C_Q_UL5SE_2D_DeriveEta, C_Q_UL5SE_2D_DeriveXiXi,
1270  C_Q_UL5SE_2D_DeriveXiEta, C_Q_UL5SE_2D_DeriveEtaEta, 5, 5,
1271  0, NULL);
TBaseFunct2D
Definition: BaseFunct2D.h:27
Generated on Tue May 10 2016 09:07:58 for ParMooN by   doxygen 1.8.5