ParMooN
 All Classes Functions Variables Friends Pages
BF_C_Q_UL5_2D.h
1 // ***********************************************************************
2 // Q5 element with bubbles, conforming, 2D
3 // ***********************************************************************
4 
5 // base function values
6 static void C_Q_UL5_2D_Funct(double xi, double eta, double *values)
7 {
8  double t1, t2, t3, t4, t5, t6, t7, t8, t9, t10, t11, t12, t13, t14, t15;
9  double t16, t17, t18, t19, t20, t21, t22, t23, t24, t25, t26, t27, t28;
10  double t29, t30, t31, t32, t33, t34, t35, t36, t37, t38, t39, t40, t41;
11  double t42, t43, t44, t45, t46, t47, t48, t49, t50, t51, t52, t53, t54;
12  double t55, t56, t57, t58, t59, t60, t61, t62, t63, t64, t65, t66, t67;
13  double t68, t69, t70, t71, t72, t73, t74, t75, t76, t77, t78, t80, t81;
14  double t82, t83, t84, t85, t86, t87, t88, t89, t90, t91, t92, t93, t94;
15  double t95, t96, t97, t98, t99, t100, t101, t102, t103, t104, t105, t106;
16  double t107, t108, t109, t110, t111, t112, t113, t114, t115, t116, t117;
17  double t118, t120, t121, t122, t123, t124, t125, t126, t127, t128, t129;
18  double t130, t131, t132, t133, t134, t135, t136, t137, t138, t139, t140;
19  double t141, t142, t143, t144, t145, t146, t147, t148, t149, t150, t151;
20  double t152, t153, t154, t155, t156, t157, t158, t160, t161, t163, t164;
21  double t166, t167, t169, t170, t171, t172, t173, t174, t175, t176, t177;
22  double t178, t179, t180, t181, t182, t183, t184, t185, t186, t187, t188;
23  double t189, t190, t191, t192, t193, t194, t195, t196, t197, t198, t199;
24  double t200, t201, t202, t204, t205, t206, t207, t208, t209, t210, t211;
25  double t212, t213, t214, t215, t216, t217, t218, t219, t220, t221, t222;
26  double t223, t224, t225, t226, t227, t228, t229, t230, t231, t232, t233;
27  double t234, t235, t237, t238, t240, t241, t243, t244, t246, t247, t249;
28  double t250, t252, t253, t255, t256, t258, t259, t261, t262, t264, t265;
29  double t267, t268, t270, t271, t285, t296, t297, t306, t349, t352, t354;
30  double t365, t436;
31 
32  t1 = 333.0/2048.0*xi;
33  t2 = 333.0/2048.0*eta;
34  t3 = xi*xi;
35  t4 = t3*t3;
36  t5 = t4*xi;
37  t6 = t5*eta;
38  t7 = 5903.0/2048.0*t6;
39  t8 = t3*xi;
40  t9 = t8*eta;
41  t10 = 3835.0/1024.0*t9;
42  t11 = eta*eta;
43  t12 = t11*eta;
44  t13 = xi*t12;
45  t14 = 3835.0/1024.0*t13;
46  t15 = xi*eta;
47  t16 = 1755.0/2048.0*t15;
48  t17 = t3*t11;
49  t18 = 3395.0/1024.0*t17;
50  t19 = xi*t11;
51  t20 = 1465.0/1024.0*t19;
52  t21 = 725.0/1024.0*t8;
53  t22 = 725.0/1024.0*t12;
54  t23 = 1105.0/2048.0*t5;
55  t24 = 2497.0/3072.0*t3;
56  t25 = 2497.0/3072.0*t11;
57  t26 = t3*eta;
58  t27 = 1465.0/1024.0*t26;
59  t28 = t3*t12;
60  t29 = 9475.0/1536.0*t28;
61  t30 = t11*t11;
62  t31 = t30*eta;
63  t32 = xi*t31;
64  t33 = 5903.0/2048.0*t32;
65  t34 = t4*t30;
66  t35 = 7525.0/6144.0*t34;
67  t36 = t8*t11;
68  t37 = 9475.0/1536.0*t36;
69  t38 = 1105.0/2048.0*t31;
70  t39 = -607.0/6144.0+t1+t2-t7+t10+t14-t16-t18-t20-t21-t22+t23+t24+t25-t27+t29-t33-t35+t37+t38;
71  t40 = 2185.0/1536.0*t4;
72  t41 = 2185.0/1536.0*t30;
73  t42 = t4*t3;
74  t43 = 1463.0/2048.0*t42;
75  t44 = t3*t31;
76  t45 = 4685.0/1024.0*t44;
77  t46 = t8*t30;
78  t47 = 5425.0/1024.0*t46;
79  t48 = t5*t11;
80  t49 = 4685.0/1024.0*t48;
81  t50 = t5*t12;
82  t51 = 38045.0/3072.0*t50;
83  t52 = t8*t31;
84  t53 = 38045.0/3072.0*t52;
85  t54 = t42*t11;
86  t55 = 1463.0/2048.0*t54;
87  t56 = t5*t30;
88  t57 = 22295.0/6144.0*t56;
89  t58 = t5*t31;
90  t59 = 18627.0/2048.0*t58;
91  t60 = t30*t11;
92  t61 = t3*t60;
93  t62 = 1463.0/2048.0*t61;
94  t63 = t4*t31;
95  t64 = 22295.0/6144.0*t63;
96  t65 = t3*t30;
97  t66 = 6255.0/2048.0*t65;
98  t67 = t8*t12;
99  t68 = 25025.0/1536.0*t67;
100  t69 = t4*t11;
101  t70 = 6255.0/2048.0*t69;
102  t71 = xi*t30;
103  t72 = 2585.0/2048.0*t71;
104  t73 = t4*eta;
105  t74 = 2585.0/2048.0*t73;
106  t75 = t4*t12;
107  t76 = 5425.0/1024.0*t75;
108  t77 = 1463.0/2048.0*t60;
109  t78 = -t40-t41+t43-t45-t47-t49+t51+t53-t55+t57-t59-t62+t64+t66-t68+t70+t72+t74-t76+t77;
110  t80 = 125.0/4096.0*xi;
111  t81 = 375.0/4096.0*eta;
112  t82 = 15625.0/4096.0*t6;
113  t83 = 8125.0/2048.0*t9;
114  t84 = 4375.0/6144.0*t13;
115  t85 = 625.0/4096.0*t15;
116  t86 = 13875.0/2048.0*t17;
117  t87 = 625.0/2048.0*t19;
118  t88 = 1625.0/2048.0*t8;
119  t89 = 875.0/2048.0*t12;
120  t90 = 3125.0/4096.0*t5;
121  t91 = 625.0/1024.0*t3;
122  t92 = 1125.0/512.0*t11;
123  t93 = 4875.0/2048.0*t26;
124  t94 = 11375.0/1024.0*t28;
125  t95 = 2625.0/4096.0*t32;
126  t96 = 21875.0/4096.0*t34;
127  t97 = 8125.0/1024.0*t36;
128  t98 = 1575.0/4096.0*t31;
129  t99 = t80+t81-t82+t83+t84-t85-t86-t87-t88-t89+t90+t91+t92-t93+t94-t95-t96+t97+t98;
130  t100 = 1875.0/4096.0*t4;
131  t101 = 20125.0/4096.0*t30;
132  t102 = 20475.0/2048.0*t44;
133  t103 = 56875.0/6144.0*t46;
134  t104 = 15625.0/2048.0*t48;
135  t105 = 109375.0/6144.0*t50;
136  t106 = 34125.0/2048.0*t52;
137  t107 = 109375.0/12288.0*t56;
138  t108 = 65625.0/4096.0*t58;
139  t109 = 5775.0/2048.0*t61;
140  t110 = 39375.0/4096.0*t63;
141  t111 = 2625.0/256.0*t65;
142  t112 = 56875.0/3072.0*t67;
143  t113 = 9375.0/2048.0*t69;
144  t114 = 4375.0/12288.0*t71;
145  t115 = 9375.0/4096.0*t73;
146  t116 = 21875.0/2048.0*t75;
147  t117 = 5775.0/2048.0*t60;
148  t118 = -625.0/4096.0-t100-t101-t102-t103-t104+t105+t106+t107-t108-t109+t110+t111-t112+t113+t114+t115-t116+t117;
149  t120 = 1125.0/2048.0*xi;
150  t121 = 1125.0/2048.0*eta;
151  t122 = 15625.0/2048.0*t6;
152  t123 = 10625.0/1024.0*t9;
153  t124 = 13125.0/1024.0*t13;
154  t125 = 5625.0/2048.0*t15;
155  t126 = 2875.0/1024.0*t17;
156  t127 = 5625.0/1024.0*t19;
157  t128 = 2125.0/1024.0*t8;
158  t129 = 2625.0/1024.0*t12;
159  t130 = 3125.0/2048.0*t5;
160  t131 = 125.0/384.0*t3;
161  t132 = 125.0/512.0*t11;
162  t133 = 2125.0/1024.0*t26;
163  t134 = 14875.0/1536.0*t28;
164  t135 = 23625.0/2048.0*t32;
165  t136 = 21875.0/6144.0*t34;
166  t137 = 10625.0/512.0*t36;
167  t138 = 125.0/6144.0-t120-t121+t122-t123-t124+t125+t126+t127+t128+t129-t130-t131+t132+t133-t134+t135+t136-t137;
168  t139 = 4725.0/2048.0*t31;
169  t140 = 625.0/2048.0*t4;
170  t141 = 11375.0/6144.0*t30;
171  t142 = 8925.0/1024.0*t44;
172  t143 = 74375.0/3072.0*t46;
173  t144 = 15625.0/1024.0*t48;
174  t145 = 109375.0/3072.0*t50;
175  t146 = 44625.0/1024.0*t52;
176  t147 = 109375.0/6144.0*t56;
177  t148 = 65625.0/2048.0*t58;
178  t149 = 1925.0/1024.0*t61;
179  t150 = 13125.0/2048.0*t63;
180  t151 = 875.0/512.0*t65;
181  t152 = 74375.0/1536.0*t67;
182  t153 = 3125.0/1024.0*t69;
183  t154 = 13125.0/2048.0*t71;
184  t155 = 3125.0/2048.0*t73;
185  t156 = 21875.0/3072.0*t75;
186  t157 = 1925.0/1024.0*t60;
187  t158 = -t139+t140-t141+t142+t143+t144-t145-t146-t147+t148-t149-t150-t151+t152-t153-t154-t155+t156+t157;
188  t160 = 125.0/6144.0+t120-t121-t122+t123+t124-t125+t126-t127-t128+t129+t130-t131+t132+t133-t134-t135+t136+t137;
189  t161 = -t139+t140-t141+t142-t143-t144+t145+t146+t147-t148-t149-t150-t151-t152-t153+t154-t155+t156+t157;
190  t163 = -t80+t81+t82-t83-t84+t85-t86+t87+t88-t89-t90+t91+t92-t93+t94+t95-t96-t97+t98;
191  t164 = -625.0/4096.0-t100-t101-t102+t103+t104-t105-t106-t107+t108-t109+t110+t111+t112+t113-t114+t115-t116+t117;
192  t166 = -607.0/6144.0-t1+t2+t7-t10-t14+t16-t18+t20+t21-t22-t23+t24+t25-t27+t29+t33-t35-t37+t38;
193  t167 = -t40-t41+t43-t45+t47+t49-t51-t53-t55-t57+t59-t62+t64+t66+t68+t70-t72+t74-t76+t77;
194  t169 = 375.0/4096.0*xi;
195  t170 = 125.0/4096.0*eta;
196  t171 = 2625.0/4096.0*t6;
197  t172 = 4375.0/6144.0*t9;
198  t173 = 8125.0/2048.0*t13;
199  t174 = 4875.0/2048.0*t19;
200  t175 = 875.0/2048.0*t8;
201  t176 = 1625.0/2048.0*t12;
202  t177 = 1575.0/4096.0*t5;
203  t178 = 1125.0/512.0*t3;
204  t179 = 625.0/1024.0*t11;
205  t180 = 625.0/2048.0*t26;
206  t181 = 8125.0/1024.0*t28;
207  t182 = 15625.0/4096.0*t32;
208  t183 = 11375.0/1024.0*t36;
209  t184 = 3125.0/4096.0*t31;
210  t185 = -t169+t170+t171-t172-t173+t85-t86+t174+t175-t176-t177+t178+t179-t180+t181+t182-t96-t183+t184;
211  t186 = 20125.0/4096.0*t4;
212  t187 = 1875.0/4096.0*t30;
213  t188 = 5775.0/2048.0*t42;
214  t189 = 15625.0/2048.0*t44;
215  t190 = 21875.0/2048.0*t46;
216  t191 = 20475.0/2048.0*t48;
217  t192 = 34125.0/2048.0*t50;
218  t193 = 109375.0/6144.0*t52;
219  t194 = 5775.0/2048.0*t54;
220  t195 = 39375.0/4096.0*t56;
221  t196 = 109375.0/12288.0*t63;
222  t197 = 9375.0/2048.0*t65;
223  t198 = 2625.0/256.0*t69;
224  t199 = 9375.0/4096.0*t71;
225  t200 = 4375.0/12288.0*t73;
226  t201 = 56875.0/6144.0*t75;
227  t202 = -625.0/4096.0-t186-t187+t188-t189+t190+t191-t192-t193-t194-t195+t108+t196+t197+t112+t198-t199+t200-t201;
228  t204 = 23625.0/2048.0*t6;
229  t205 = 13125.0/1024.0*t9;
230  t206 = 10625.0/1024.0*t13;
231  t207 = 2125.0/1024.0*t19;
232  t208 = 2625.0/1024.0*t8;
233  t209 = 2125.0/1024.0*t12;
234  t210 = 4725.0/2048.0*t5;
235  t211 = 125.0/512.0*t3;
236  t212 = 125.0/384.0*t11;
237  t213 = 5625.0/1024.0*t26;
238  t214 = 10625.0/512.0*t28;
239  t215 = 15625.0/2048.0*t32;
240  t216 = 14875.0/1536.0*t36;
241  t217 = 125.0/6144.0+t120-t121-t204+t205+t206-t125+t126-t207-t208+t209+t210+t211-t212+t213-t214-t215+t136+t216;
242  t218 = 3125.0/2048.0*t31;
243  t219 = 11375.0/6144.0*t4;
244  t220 = 625.0/2048.0*t30;
245  t221 = 1925.0/1024.0*t42;
246  t222 = 15625.0/1024.0*t44;
247  t223 = 21875.0/3072.0*t46;
248  t224 = 8925.0/1024.0*t48;
249  t225 = 44625.0/1024.0*t50;
250  t226 = 109375.0/3072.0*t52;
251  t227 = 1925.0/1024.0*t54;
252  t228 = 13125.0/2048.0*t56;
253  t229 = 109375.0/6144.0*t63;
254  t230 = 3125.0/1024.0*t65;
255  t231 = 875.0/512.0*t69;
256  t232 = 3125.0/2048.0*t71;
257  t233 = 13125.0/2048.0*t73;
258  t234 = 74375.0/3072.0*t75;
259  t235 = -t218-t219+t220+t221+t222-t223-t224+t225+t226-t227+t228-t148-t229-t230-t152-t231+t232-t233+t234;
260  t237 = 125.0/6144.0+t120+t121+t204-t205-t206+t125+t126-t207-t208-t209+t210+t211-t212-t213+t214+t215+t136+t216;
261  t238 = t218-t219+t220+t221-t222-t223-t224-t225-t226-t227+t228+t148+t229-t230+t152-t231+t232+t233-t234;
262  t240 = -t169-t170-t171+t172+t173-t85-t86+t174+t175+t176-t177+t178+t179+t180-t181-t182-t96-t183-t184;
263  t241 = -625.0/4096.0-t186-t187+t188+t189+t190+t191+t192+t193-t194-t195-t108-t196+t197-t112+t198-t199-t200+t201;
264  t243 = -607.0/6144.0-t1-t2-t7+t10+t14-t16-t18+t20+t21+t22-t23+t24+t25+t27-t29-t33-t35-t37-t38;
265  t244 = -t40-t41+t43+t45+t47+t49+t51+t53-t55-t57-t59-t62-t64+t66-t68+t70-t72-t74+t76+t77;
266  t246 = -t80-t81-t82+t83+t84-t85-t86+t87+t88+t89-t90+t91+t92+t93-t94-t95-t96-t97-t98;
267  t247 = -625.0/4096.0-t100-t101+t102+t103+t104+t105+t106-t107-t108-t109-t110+t111-t112+t113-t114-t115+t116+t117;
268  t249 = 125.0/6144.0+t120+t121+t122-t123-t124+t125+t126-t127-t128-t129+t130-t131+t132-t133+t134+t135+t136+t137;
269  t250 = t139+t140-t141-t142-t143-t144-t145-t146+t147+t148-t149+t150-t151+t152-t153+t154+t155-t156+t157;
270  t252 = 125.0/6144.0-t120+t121-t122+t123+t124-t125+t126+t127+t128-t129-t130-t131+t132-t133+t134-t135+t136-t137;
271  t253 = t139+t140-t141-t142+t143+t144+t145+t146-t147-t148-t149+t150-t151-t152-t153-t154+t155-t156+t157;
272  t255 = t80-t81+t82-t83-t84+t85-t86-t87-t88+t89+t90+t91+t92+t93-t94+t95-t96+t97-t98;
273  t256 = -625.0/4096.0-t100-t101+t102-t103-t104-t105-t106+t107+t108-t109-t110+t111+t112+t113+t114-t115+t116+t117;
274  t258 = -607.0/6144.0+t1-t2+t7-t10-t14+t16-t18-t20-t21+t22+t23+t24+t25+t27-t29+t33-t35+t37-t38;
275  t259 = -t40-t41+t43+t45-t47-t49-t51-t53-t55+t57+t59-t62-t64+t66+t68+t70+t72-t74+t76+t77;
276  t261 = t169-t170+t171-t172-t173+t85-t86-t174-t175+t176+t177+t178+t179+t180-t181+t182-t96+t183-t184;
277  t262 = -625.0/4096.0-t186-t187+t188+t189-t190-t191-t192-t193-t194+t195+t108-t196+t197+t112+t198+t199-t200+t201;
278  t264 = 125.0/6144.0-t120+t121-t204+t205+t206-t125+t126+t207+t208-t209-t210+t211-t212-t213+t214-t215+t136-t216;
279  t265 = t218-t219+t220+t221-t222+t223+t224+t225+t226-t227-t228-t148+t229-t230-t152-t231-t232+t233-t234;
280  t267 = 125.0/6144.0-t120-t121+t204-t205-t206+t125+t126+t207+t208+t209-t210+t211-t212+t213-t214+t215+t136-t216;
281  t268 = -t218-t219+t220+t221+t222+t223+t224-t225-t226-t227-t228+t148-t229-t230+t152-t231-t232-t233+t234;
282  t270 = t169+t170-t171+t172+t173-t85-t86-t174-t175-t176+t177+t178+t179-t180+t181-t182-t96+t183+t184;
283  t271 = -625.0/4096.0-t186-t187+t188-t189-t190-t191+t192+t193-t194+t195-t108+t196+t197-t112+t198+t199+t200-t201;
284  t285 = 45.0/4.0*t17+1225.0/256.0*t34+2135.0/256.0*t4-453.0/128.0*t3-453.0/128.0*t11-105.0/8.0*t65+693.0/128.0*t54+2135.0/256.0*t30-693.0/128.0*t42-693.0/128.0*t60+157.0/256.0+693.0/128.0*t61-105.0/8.0*t69;
285  t296 = 5625.0/128.0*t17;
286  t297 = 6125.0/256.0*t34;
287  t306 = t296+t297+1225.0/256.0*t4-345.0/64.0*t3-975.0/64.0*t11-525.0/8.0*t65+10675.0/256.0*t30+155.0/256.0-3465.0/128.0*t60+3465.0/128.0*t61-3675.0/128.0*t69;
288  t349 = 3969.0/256.0*t15;
289  t352 = 33075.0/64.0*t67;
290  t354 = 83349.0/256.0*t58;
291  t365 = t296+t297+10675.0/256.0*t4-975.0/64.0*t3-345.0/64.0*t11-3675.0/128.0*t65+3465.0/128.0*t54+1225.0/256.0*t30-3465.0/128.0*t42+155.0/256.0-525.0/8.0*t69;
292  t436 = 4455.0/128.0*t17;
293 
294  values[0] = t39+t78;
295  values[1] = t99+t118;
296  values[2] = t138+t158;
297  values[3] = t160+t161;
298  values[4] = t163+t164;
299  values[5] = t166+t167;
300  values[6] = t185+t202;
301  values[7] = t217+t235;
302  values[8] = t237+t238;
303  values[9] = t240+t241;
304  values[10] = t243+t244;
305  values[11] = t246+t247;
306  values[12] = t249+t250;
307  values[13] = t252+t253;
308  values[14] = t255+t256;
309  values[15] = t258+t259;
310  values[16] = t261+t262;
311  values[17] = t264+t265;
312  values[18] = t267+t268;
313  values[19] = t270+t271;
314  values[20] = t285;
315  values[21] = -105.0/256.0*eta-945.0/256.0*t31+525.0/128.0*t12+735.0/256.0*t73-315.0/128.0*t26+6615.0/256.0*t63-3675.0/128.0*t75-2835.0/128.0*t44+1575.0/64.0*t28;
316  values[22] = t306;
317  values[23] = -945.0/256.0*eta-2205.0/256.0*t31+1575.0/128.0*t12-11025.0/128.0*t75+6615.0/256.0*t73+4725.0/64.0*t28-2835.0/128.0*t26+15435.0/256.0*t63-6615.0/128.0*t44;
318  values[24] = -105.0/256.0*xi-945.0/256.0*t5+735.0/256.0*t71-315.0/128.0*t19+525.0/128.0*t8+6615.0/256.0*t56-2835.0/128.0*t48-3675.0/128.0*t46+1575.0/64.0*t36;
319  values[25] = 441.0/256.0*t15+3969.0/256.0*t32-2205.0/128.0*t13+3969.0/256.0*t6-2205.0/128.0*t9+35721.0/256.0*t58-19845.0/128.0*t50-19845.0/128.0*t52+11025.0/64.0*t67;
320  values[26] = 735.0/256.0*xi-2205.0/128.0*t19+3675.0/256.0*t71+6615.0/256.0*t5-19845.0/128.0*t48-3675.0/128.0*t8+11025.0/64.0*t36+33075.0/256.0*t56-18375.0/128.0*t46;
321  values[27] = 9261.0/256.0*t32-6615.0/128.0*t13+t349-59535.0/128.0*t50+35721.0/256.0*t6+t352-19845.0/128.0*t9+t354-46305.0/128.0*t52;
322  values[28] = t365;
323  values[29] = 735.0/256.0*eta+3675.0/256.0*t73-2205.0/128.0*t26+6615.0/256.0*t31-3675.0/128.0*t12-19845.0/128.0*t44+11025.0/64.0*t28+33075.0/256.0*t63-18375.0/128.0*t75;
324  values[30] = 11025.0/64.0*t17+30625.0/256.0*t34+1225.0/256.0+6125.0/256.0*t4-3675.0/128.0*t3-3675.0/128.0*t11-18375.0/128.0*t65+6125.0/256.0*t30-18375.0/128.0*t69;
325  values[31] = 6615.0/256.0*eta-46305.0/128.0*t44-55125.0/128.0*t75+33075.0/256.0*t73+15435.0/256.0*t31-11025.0/128.0*t12-19845.0/128.0*t26+33075.0/64.0*t28+77175.0/256.0*t63;
326  values[32] = -945.0/256.0*xi-2205.0/256.0*t5+1575.0/128.0*t8-11025.0/128.0*t46+4725.0/64.0*t36+6615.0/256.0*t71-2835.0/128.0*t19+15435.0/256.0*t56-6615.0/128.0*t48;
327  values[33] = -59535.0/128.0*t52+t352-6615.0/128.0*t9+35721.0/256.0*t32-19845.0/128.0*t13+t349+9261.0/256.0*t6+t354-46305.0/128.0*t50;
328  values[34] = 6615.0/256.0*xi+33075.0/256.0*t71-55125.0/128.0*t46-46305.0/128.0*t48+15435.0/256.0*t5-11025.0/128.0*t8+33075.0/64.0*t36-19845.0/128.0*t19+77175.0/256.0*t56;
329  values[35] = 83349.0/256.0*t32+35721.0/256.0*t15-138915.0/128.0*t52+194481.0/256.0*t58-59535.0/128.0*t13+83349.0/256.0*t6-138915.0/128.0*t50-59535.0/128.0*t9+99225.0/64.0*t67;
330  values[36] = 297.0/128.0+10395.0/128.0*t30-4455.0/128.0*t11-297.0/128.0*t3-10395.0/128.0*t65+t436-6237.0/128.0*t60+6237.0/128.0*t61;
331  values[37] = 10395.0/128.0*t4-4455.0/128.0*t3+297.0/128.0-297.0/128.0*t11-10395.0/128.0*t69+t436-6237.0/128.0*t42+6237.0/128.0*t54;
332 }
333 
334 
335 // values of the derivatives in xi direction
336 static void C_Q_UL5_2D_DeriveXi(double xi, double eta, double *values)
337 {
338  double t1, t2, t3, t4, t5, t6, t7, t8, t9, t10, t11, t12, t13, t14, t15;
339  double t16, t17, t18, t19, t20, t21, t22, t23, t24, t25, t26, t27, t28;
340  double t29, t30, t31, t32, t33, t34, t35, t36, t37, t38, t39, t40, t41;
341  double t43, t44, t45, t46, t47, t48, t49, t50, t51, t52, t53, t54, t55;
342  double t56, t57, t58, t59, t60, t61, t62, t63, t64, t65, t67, t68, t69;
343  double t70, t71, t72, t73, t74, t75, t76, t77, t78, t79, t80, t81, t82;
344  double t83, t84, t85, t86, t87, t88, t89, t90, t91, t92, t93, t94, t95;
345  double t96, t97, t98, t100, t101, t102, t103, t104, t105, t106, t107;
346  double t108, t109, t110, t111, t112, t113, t114, t115, t116, t117, t118;
347  double t119, t120, t121, t122, t123, t124, t125, t126, t127, t128, t129;
348  double t130, t131, t133, t134, t136, t137, t139, t140, t142, t143, t144;
349  double t145, t146, t147, t148, t149, t150, t151, t152, t153, t154, t155;
350  double t156, t157, t158, t159, t160, t161, t162, t163, t164, t165, t166;
351  double t167, t168, t169, t171, t172, t173, t174, t175, t176, t177, t178;
352  double t179, t180, t181, t182, t183, t184, t185, t186, t187, t188, t189;
353  double t190, t191, t192, t193, t194, t195, t196, t197, t198, t200, t201;
354  double t203, t204, t206, t207, t209, t210, t212, t213, t215, t216, t218;
355  double t219, t221, t222, t224, t225, t227, t228, t230, t231, t233, t234;
356  double t253, t254, t298, t301, t303, t371;
357 
358  t1 = 2497.0/1536.0*xi;
359  t2 = xi*xi;
360  t3 = t2*xi;
361  t4 = 2185.0/384.0*t3;
362  t5 = t3*eta;
363  t6 = 2585.0/512.0*t5;
364  t7 = eta*eta;
365  t8 = t7*eta;
366  t9 = xi*t8;
367  t10 = 9475.0/768.0*t9;
368  t11 = t2*t7;
369  t12 = 9475.0/512.0*t11;
370  t13 = xi*t7;
371  t14 = 3395.0/512.0*t13;
372  t15 = t2*eta;
373  t16 = 11505.0/1024.0*t15;
374  t17 = t3*t7;
375  t18 = 6255.0/512.0*t17;
376  t19 = t2*t8;
377  t20 = 25025.0/512.0*t19;
378  t21 = t7*t7;
379  t22 = t21*eta;
380  t23 = xi*t22;
381  t24 = 4685.0/512.0*t23;
382  t25 = t2*t2;
383  t26 = t25*t21;
384  t27 = 111475.0/6144.0*t26;
385  t28 = t2*t22;
386  t29 = 38045.0/1024.0*t28;
387  t30 = t25*t8;
388  t31 = 190225.0/3072.0*t30;
389  t32 = t3*t21;
390  t33 = 7525.0/1536.0*t32;
391  t34 = t25*xi;
392  t35 = t34*t7;
393  t36 = 4389.0/1024.0*t35;
394  t37 = 333.0/2048.0+t1-t4+t6+t10+t12-t14+t16+t18-t20-t24+t27+t29+t31-t33-t36;
395  t38 = t3*t22;
396  t39 = 22295.0/1536.0*t38;
397  t40 = t25*t22;
398  t41 = 93135.0/2048.0*t40;
399  t43 = xi*t21*t7;
400  t44 = 1463.0/1024.0*t43;
401  t45 = t2*t21;
402  t46 = 16275.0/1024.0*t45;
403  t47 = t3*t8;
404  t48 = 5425.0/256.0*t47;
405  t49 = t25*t7;
406  t50 = 23425.0/1024.0*t49;
407  t51 = t25*eta;
408  t52 = 29515.0/2048.0*t51;
409  t53 = 5903.0/2048.0*t22;
410  t54 = 3835.0/1024.0*t8;
411  t55 = xi*t21;
412  t56 = 6255.0/1024.0*t55;
413  t57 = 1755.0/2048.0*eta;
414  t58 = 2175.0/1024.0*t2;
415  t59 = 1465.0/1024.0*t7;
416  t60 = xi*eta;
417  t61 = 1465.0/512.0*t60;
418  t62 = 5525.0/2048.0*t25;
419  t63 = 2585.0/2048.0*t21;
420  t64 = 4389.0/1024.0*t34;
421  t65 = t39-t41-t44-t46-t48-t50-t52-t53+t54+t56-t57-t58-t59-t61+t62+t63+t64;
422  t67 = 625.0/512.0*xi;
423  t68 = 1875.0/1024.0*t3;
424  t69 = 9375.0/1024.0*t5;
425  t70 = 11375.0/512.0*t9;
426  t71 = 24375.0/1024.0*t11;
427  t72 = 13875.0/1024.0*t13;
428  t73 = 24375.0/2048.0*t15;
429  t74 = 9375.0/512.0*t17;
430  t75 = 56875.0/1024.0*t19;
431  t76 = 20475.0/1024.0*t23;
432  t77 = 546875.0/12288.0*t26;
433  t78 = 102375.0/2048.0*t28;
434  t79 = 546875.0/6144.0*t30;
435  t80 = 21875.0/1024.0*t32;
436  t81 = 39375.0/1024.0*t38;
437  t82 = t67-t68+t69+t70+t71-t72+t73+t74-t75-t76+t77+t78+t79-t80+t81;
438  t83 = 328125.0/4096.0*t40;
439  t84 = 5775.0/1024.0*t43;
440  t85 = 56875.0/2048.0*t45;
441  t86 = 21875.0/512.0*t47;
442  t87 = 78125.0/2048.0*t49;
443  t88 = 78125.0/4096.0*t51;
444  t89 = 2625.0/4096.0*t22;
445  t90 = 4375.0/6144.0*t8;
446  t91 = 2625.0/128.0*t55;
447  t92 = 625.0/4096.0*eta;
448  t93 = 4875.0/2048.0*t2;
449  t94 = 625.0/2048.0*t7;
450  t95 = 4875.0/1024.0*t60;
451  t96 = 15625.0/4096.0*t25;
452  t97 = 4375.0/12288.0*t21;
453  t98 = 125.0/4096.0-t83-t84-t85-t86-t87-t88-t89+t90+t91-t92-t93-t94-t95+t96+t97;
454  t100 = 125.0/192.0*xi;
455  t101 = 625.0/512.0*t3;
456  t102 = 3125.0/512.0*t5;
457  t103 = 14875.0/768.0*t9;
458  t104 = 31875.0/512.0*t11;
459  t105 = 2875.0/512.0*t13;
460  t106 = 31875.0/1024.0*t15;
461  t107 = 3125.0/256.0*t17;
462  t108 = 74375.0/512.0*t19;
463  t109 = 8925.0/512.0*t23;
464  t110 = 546875.0/6144.0*t26;
465  t111 = 133875.0/1024.0*t28;
466  t112 = 546875.0/3072.0*t30;
467  t113 = 21875.0/1536.0*t32;
468  t114 = -1125.0/2048.0-t100+t101-t102-t103-t104+t105-t106-t107+t108+t109-t110-t111-t112+t113;
469  t115 = 13125.0/512.0*t38;
470  t116 = 328125.0/2048.0*t40;
471  t117 = 1925.0/512.0*t43;
472  t118 = 74375.0/1024.0*t45;
473  t119 = 21875.0/768.0*t47;
474  t120 = 78125.0/1024.0*t49;
475  t121 = 78125.0/2048.0*t51;
476  t122 = 23625.0/2048.0*t22;
477  t123 = 13125.0/1024.0*t8;
478  t124 = 875.0/256.0*t55;
479  t125 = 5625.0/2048.0*eta;
480  t126 = 6375.0/1024.0*t2;
481  t127 = 5625.0/1024.0*t7;
482  t128 = 2125.0/512.0*t60;
483  t129 = 15625.0/2048.0*t25;
484  t130 = 13125.0/2048.0*t21;
485  t131 = -t115+t116-t117+t118+t119+t120+t121+t122-t123-t124+t125+t126+t127+t128-t129-t130;
486  t133 = -t100+t101-t102-t103+t104+t105+t106-t107-t108+t109+t110+t111+t112+t113-t115;
487  t134 = 1125.0/2048.0-t116-t117-t118+t119-t120-t121-t122+t123-t124-t125-t126-t127+t128+t129+t130;
488  t136 = t67-t68+t69+t70-t71-t72-t73+t74+t75-t76-t77-t78-t79-t80+t81;
489  t137 = -125.0/4096.0+t83-t84+t85-t86+t87+t88+t89-t90+t91+t92+t93+t94-t95-t96-t97;
490  t139 = t1-t4+t6+t10-t12-t14-t16+t18+t20-t24-t27-t29-t31-t33-t36+t39;
491  t140 = -333.0/2048.0+t41-t44+t46-t48+t50+t52+t53-t54+t56+t57+t58+t59-t61-t62-t63+t64;
492  t142 = 1125.0/256.0*xi;
493  t143 = 20125.0/1024.0*t3;
494  t144 = 4375.0/3072.0*t5;
495  t145 = 8125.0/512.0*t9;
496  t146 = 34125.0/1024.0*t11;
497  t147 = 4375.0/2048.0*t15;
498  t148 = 2625.0/64.0*t17;
499  t149 = 15625.0/1024.0*t23;
500  t150 = 196875.0/4096.0*t26;
501  t151 = 109375.0/2048.0*t28;
502  t152 = 170625.0/2048.0*t30;
503  t153 = 17325.0/1024.0*t35;
504  t154 = -375.0/4096.0+t142-t143+t144+t145-t146-t72-t147+t148+t75-t149-t150-t151-t152-t80-t153;
505  t155 = 109375.0/3072.0*t38;
506  t156 = 65625.0/2048.0*t45;
507  t157 = 56875.0/1536.0*t47;
508  t158 = 102375.0/2048.0*t49;
509  t159 = 13125.0/4096.0*t51;
510  t160 = 15625.0/4096.0*t22;
511  t161 = 8125.0/2048.0*t8;
512  t162 = 9375.0/1024.0*t55;
513  t163 = 2625.0/2048.0*t2;
514  t164 = 4875.0/2048.0*t7;
515  t165 = 625.0/1024.0*t60;
516  t166 = 7875.0/4096.0*t25;
517  t167 = 9375.0/4096.0*t21;
518  t168 = 17325.0/1024.0*t34;
519  t169 = t155+t83+t156-t157+t158+t159+t160-t161+t162+t92+t163+t164-t165-t166-t167+t168;
520  t171 = 125.0/256.0*xi;
521  t172 = 11375.0/1536.0*t3;
522  t173 = 13125.0/512.0*t5;
523  t174 = 10625.0/256.0*t9;
524  t175 = 14875.0/512.0*t11;
525  t176 = 39375.0/1024.0*t15;
526  t177 = 875.0/128.0*t17;
527  t178 = 15625.0/512.0*t23;
528  t179 = 65625.0/2048.0*t26;
529  t180 = 109375.0/1024.0*t28;
530  t181 = 223125.0/1024.0*t30;
531  t182 = 5775.0/512.0*t35;
532  t183 = 109375.0/1536.0*t38;
533  t184 = t171-t172-t173-t174+t175+t105+t176-t177-t108+t178+t179+t180+t181+t113-t182-t183;
534  t185 = 21875.0/1024.0*t45;
535  t186 = 74375.0/768.0*t47;
536  t187 = 44625.0/1024.0*t49;
537  t188 = 118125.0/2048.0*t51;
538  t189 = 15625.0/2048.0*t22;
539  t190 = 10625.0/1024.0*t8;
540  t191 = 3125.0/512.0*t55;
541  t192 = 7875.0/1024.0*t2;
542  t193 = 2125.0/1024.0*t7;
543  t194 = 5625.0/512.0*t60;
544  t195 = 23625.0/2048.0*t25;
545  t196 = 3125.0/2048.0*t21;
546  t197 = 5775.0/512.0*t34;
547  t198 = 1125.0/2048.0-t116-t185+t186-t187-t188-t189+t190-t191-t125-t192-t193+t194+t195+t196+t197;
548  t200 = t171-t172+t173+t174+t175+t105-t176-t177+t108-t178+t179-t180-t181+t113-t182+t183;
549  t201 = 1125.0/2048.0+t116-t185-t186-t187+t188+t189-t190-t191+t125-t192-t193-t194+t195+t196+t197;
550  t203 = -375.0/4096.0+t142-t143-t144-t145-t146-t72+t147+t148-t75+t149-t150+t151+t152-t80-t153;
551  t204 = -t155-t83+t156+t157+t158-t159-t160+t161+t162-t92+t163+t164+t165-t166-t167+t168;
552  t206 = t1-t4-t6-t10-t12-t14+t16+t18-t20+t24-t27+t29+t31-t33-t36-t39;
553  t207 = -333.0/2048.0-t41-t44+t46+t48+t50-t52-t53+t54+t56-t57+t58+t59+t61-t62-t63+t64;
554  t209 = t67-t68-t69-t70-t71-t72+t73+t74-t75+t76-t77+t78+t79-t80-t81;
555  t210 = -125.0/4096.0-t83-t84+t85+t86+t87-t88-t89+t90+t91-t92+t93+t94+t95-t96-t97;
556  t212 = -t100+t101+t102+t103+t104+t105-t106-t107+t108-t109+t110-t111-t112+t113+t115;
557  t213 = 1125.0/2048.0+t116-t117-t118-t119-t120+t121+t122-t123-t124+t125-t126-t127-t128+t129+t130;
558  t215 = -1125.0/2048.0-t100+t101+t102+t103-t104+t105+t106-t107-t108-t109-t110+t111+t112+t113;
559  t216 = t115-t116-t117+t118-t119+t120-t121-t122+t123-t124-t125+t126+t127-t128-t129-t130;
560  t218 = t67-t68-t69-t70+t71-t72-t73+t74+t75+t76+t77-t78-t79-t80-t81;
561  t219 = 125.0/4096.0+t83-t84-t85+t86-t87+t88+t89-t90+t91+t92-t93-t94+t95+t96+t97;
562  t221 = 333.0/2048.0+t1-t4-t6-t10+t12-t14-t16+t18+t20+t24+t27-t29-t31-t33-t36;
563  t222 = -t39+t41-t44-t46+t48-t50+t52+t53-t54+t56+t57-t58-t59+t61+t62+t63+t64;
564  t224 = t142-t143-t144-t145+t146-t72-t147+t148+t75+t149+t150-t151-t152-t80-t153-t155;
565  t225 = 375.0/4096.0+t83-t156+t157-t158+t159+t160-t161+t162+t92-t163-t164+t165+t166+t167+t168;
566  t227 = -1125.0/2048.0+t171-t172+t173+t174-t175+t105+t176-t177-t108-t178-t179+t180+t181+t113-t182;
567  t228 = t183-t116+t185-t186+t187-t188-t189+t190-t191-t125+t192+t193-t194-t195-t196+t197;
568  t230 = -1125.0/2048.0+t171-t172-t173-t174-t175+t105-t176-t177+t108+t178-t179-t180-t181+t113-t182;
569  t231 = -t183+t116+t185+t186+t187+t188+t189-t190-t191+t125+t192+t193+t194-t195-t196+t197;
570  t233 = t142-t143+t144+t145+t146-t72+t147+t148-t75-t149+t150+t151+t152-t80-t153+t155;
571  t234 = 375.0/4096.0-t83-t156-t157-t158-t159-t160+t161+t162-t92-t163-t164-t165+t166+t167+t168;
572  t253 = 5625.0/64.0*t13;
573  t254 = 6125.0/64.0*t32;
574  t298 = 3969.0/256.0*eta;
575  t301 = 99225.0/64.0*t19;
576  t303 = 416745.0/256.0*t40;
577  t371 = 4455.0/64.0*t13;
578 
579  values[0] = t37+t65;
580  values[1] = t82+t98;
581  values[2] = t114+t131;
582  values[3] = t133+t134;
583  values[4] = t136+t137;
584  values[5] = t139+t140;
585  values[6] = t154+t169;
586  values[7] = t184+t198;
587  values[8] = t200+t201;
588  values[9] = t203+t204;
589  values[10] = t206+t207;
590  values[11] = t209+t210;
591  values[12] = t212+t213;
592  values[13] = t215+t216;
593  values[14] = t218+t219;
594  values[15] = t221+t222;
595  values[16] = t224+t225;
596  values[17] = t227+t228;
597  values[18] = t230+t231;
598  values[19] = t233+t234;
599  values[20] = 45.0/2.0*t13+1225.0/64.0*t32+2135.0/64.0*t3-453.0/64.0*xi-105.0/4.0*t55+2079.0/64.0*t35-2079.0/64.0*t34+693.0/64.0*t43-105.0/2.0*t17;
600  values[21] = 735.0/64.0*t5-315.0/64.0*t60+6615.0/64.0*t38-3675.0/32.0*t47-2835.0/64.0*t23+1575.0/32.0*t9;
601  values[22] = t253+t254+1225.0/64.0*t3-345.0/32.0*xi-525.0/4.0*t55+3465.0/64.0*t43-3675.0/32.0*t17;
602  values[23] = -11025.0/32.0*t47+6615.0/64.0*t5+4725.0/32.0*t9-2835.0/64.0*t60+15435.0/64.0*t38-6615.0/64.0*t23;
603  values[24] = -105.0/256.0-4725.0/256.0*t25+735.0/256.0*t21-315.0/128.0*t7+1575.0/128.0*t2+33075.0/256.0*t26-14175.0/128.0*t49-11025.0/128.0*t45+4725.0/64.0*t11;
604  values[25] = 441.0/256.0*eta+3969.0/256.0*t22-2205.0/128.0*t8+19845.0/256.0*t51-6615.0/128.0*t15+178605.0/256.0*t40-99225.0/128.0*t30-59535.0/128.0*t28+33075.0/64.0*t19;
605  values[26] = 735.0/256.0-2205.0/128.0*t7+3675.0/256.0*t21+33075.0/256.0*t25-99225.0/128.0*t49-11025.0/128.0*t2+33075.0/64.0*t11+165375.0/256.0*t26-55125.0/128.0*t45;
606  values[27] = 9261.0/256.0*t22-6615.0/128.0*t8+t298-297675.0/128.0*t30+178605.0/256.0*t51+t301-59535.0/128.0*t15+t303-138915.0/128.0*t28;
607  values[28] = t253+t254+10675.0/64.0*t3-975.0/32.0*xi-3675.0/64.0*t55+10395.0/64.0*t35-10395.0/64.0*t34-525.0/2.0*t17;
608  values[29] = 3675.0/64.0*t5-2205.0/64.0*t60-19845.0/64.0*t23+11025.0/32.0*t9+33075.0/64.0*t38-18375.0/32.0*t47;
609  values[30] = 11025.0/32.0*t13+30625.0/64.0*t32+6125.0/64.0*t3-3675.0/64.0*xi-18375.0/64.0*t55-18375.0/32.0*t17;
610  values[31] = -46305.0/64.0*t23-55125.0/32.0*t47+33075.0/64.0*t5-19845.0/64.0*t60+33075.0/32.0*t9+77175.0/64.0*t38;
611  values[32] = -945.0/256.0-11025.0/256.0*t25+4725.0/128.0*t2-33075.0/128.0*t45+14175.0/64.0*t11+6615.0/256.0*t21-2835.0/128.0*t7+77175.0/256.0*t26-33075.0/128.0*t49;
612  values[33] = -178605.0/128.0*t28+t301-19845.0/128.0*t15+35721.0/256.0*t22-19845.0/128.0*t8+t298+46305.0/256.0*t51+t303-231525.0/128.0*t30;
613  values[34] = 6615.0/256.0+33075.0/256.0*t21-165375.0/128.0*t45-231525.0/128.0*t49+77175.0/256.0*t25-33075.0/128.0*t2+99225.0/64.0*t11-19845.0/128.0*t7+385875.0/256.0*t26;
614  values[35] = 83349.0/256.0*t22+35721.0/256.0*eta-416745.0/128.0*t28+972405.0/256.0*t40-59535.0/128.0*t8+416745.0/256.0*t51-694575.0/128.0*t30-178605.0/128.0*t15+297675.0/64.0*t19;
615  values[36] = -297.0/64.0*xi-10395.0/64.0*t55+t371+6237.0/64.0*t43;
616  values[37] = 10395.0/32.0*t3-4455.0/64.0*xi-10395.0/32.0*t17+t371-18711.0/64.0*t34+18711.0/64.0*t35;
617 }
618 
619 // values of the derivatives in eta direction
620 static void C_Q_UL5_2D_DeriveEta(double xi, double eta, double *values)
621 {
622  double t1, t2, t3, t4, t5, t6, t7, t8, t9, t10, t11, t12, t13, t14, t15;
623  double t16, t17, t18, t19, t20, t21, t22, t23, t24, t25, t26, t27, t28;
624  double t29, t30, t31, t32, t33, t34, t35, t36, t37, t38, t39, t40, t41;
625  double t42, t43, t44, t45, t46, t47, t48, t49, t50, t51, t53, t54, t55;
626  double t56, t57, t58, t59, t60, t61, t62, t63, t64, t65, t67, t68, t69;
627  double t70, t71, t72, t73, t74, t75, t76, t77, t78, t79, t80, t81, t82;
628  double t83, t84, t85, t86, t87, t88, t89, t90, t91, t92, t93, t94, t95;
629  double t96, t97, t98, t99, t101, t102, t103, t104, t105, t106, t107;
630  double t108, t109, t110, t111, t112, t113, t114, t115, t116, t117, t118;
631  double t119, t120, t121, t122, t123, t124, t125, t126, t127, t128, t129;
632  double t130, t131, t132, t133, t135, t136, t138, t139, t141, t142, t144;
633  double t145, t146, t147, t148, t149, t150, t151, t152, t153, t154, t155;
634  double t156, t157, t158, t159, t160, t161, t162, t163, t164, t165, t166;
635  double t167, t168, t169, t170, t172, t173, t174, t175, t176, t177, t178;
636  double t179, t180, t181, t182, t183, t184, t185, t186, t187, t188, t189;
637  double t190, t191, t192, t193, t194, t195, t196, t197, t198, t200, t201;
638  double t203, t204, t206, t207, t209, t210, t212, t213, t215, t216, t218;
639  double t219, t221, t222, t224, t225, t227, t228, t230, t231, t233, t234;
640  double t255, t256, t299, t302, t304, t372;
641 
642  t1 = 2497.0/1536.0*eta;
643  t2 = eta*eta;
644  t3 = t2*eta;
645  t4 = 2185.0/384.0*t3;
646  t5 = 1755.0/2048.0*xi;
647  t6 = xi*xi;
648  t7 = t6*xi;
649  t8 = 3835.0/1024.0*t7;
650  t9 = t6*t6;
651  t10 = t9*xi;
652  t11 = 5903.0/2048.0*t10;
653  t12 = t2*t2;
654  t13 = t10*t12;
655  t14 = 93135.0/2048.0*t13;
656  t15 = t10*t3;
657  t16 = 22295.0/1536.0*t15;
658  t17 = 2175.0/1024.0*t2;
659  t18 = 1465.0/1024.0*t6;
660  t19 = t10*t2;
661  t20 = 38045.0/1024.0*t19;
662  t21 = t7*t12;
663  t22 = 190225.0/3072.0*t21;
664  t23 = t12*eta;
665  t24 = t6*t23;
666  t25 = 4389.0/1024.0*t24;
667  t26 = t9*t3;
668  t27 = 7525.0/1536.0*t26;
669  t28 = t9*t12;
670  t29 = 111475.0/6144.0*t28;
671  t30 = 2585.0/2048.0*t9;
672  t31 = 5525.0/2048.0*t12;
673  t32 = t1-t4-t5+t8-t11-t14+t16-t17-t18+t20+t22-t25-t27+t29+t30+t31;
674  t33 = t6*t3;
675  t34 = 6255.0/512.0*t33;
676  t35 = t2*t7;
677  t36 = 25025.0/512.0*t35;
678  t37 = 4389.0/1024.0*t23;
679  t38 = t6*eta;
680  t39 = 3395.0/512.0*t38;
681  t40 = xi*t12;
682  t41 = 29515.0/2048.0*t40;
683  t42 = t2*xi;
684  t43 = 11505.0/1024.0*t42;
685  t44 = t6*t2;
686  t45 = 9475.0/512.0*t44;
687  t46 = xi*eta;
688  t47 = 1465.0/512.0*t46;
689  t48 = xi*t3;
690  t49 = 2585.0/512.0*t48;
691  t50 = t7*eta;
692  t51 = 9475.0/768.0*t50;
693  t53 = t9*t6*eta;
694  t54 = 1463.0/1024.0*t53;
695  t55 = t3*t7;
696  t56 = 5425.0/256.0*t55;
697  t57 = t10*eta;
698  t58 = 4685.0/512.0*t57;
699  t59 = t6*t12;
700  t60 = 23425.0/1024.0*t59;
701  t61 = t9*eta;
702  t62 = 6255.0/1024.0*t61;
703  t63 = t9*t2;
704  t64 = 16275.0/1024.0*t63;
705  t65 = 333.0/2048.0+t34-t36+t37-t39-t41+t43+t45-t47+t49+t51-t54-t56-t58-t60+t62-t64;
706  t67 = 1125.0/256.0*eta;
707  t68 = 20125.0/1024.0*t3;
708  t69 = 625.0/4096.0*xi;
709  t70 = 8125.0/2048.0*t7;
710  t71 = 15625.0/4096.0*t10;
711  t72 = 328125.0/4096.0*t13;
712  t73 = 109375.0/3072.0*t15;
713  t74 = 2625.0/2048.0*t2;
714  t75 = 4875.0/2048.0*t6;
715  t76 = 109375.0/2048.0*t19;
716  t77 = 170625.0/2048.0*t21;
717  t78 = 17325.0/1024.0*t24;
718  t79 = 21875.0/1024.0*t26;
719  t80 = 196875.0/4096.0*t28;
720  t81 = 9375.0/4096.0*t9;
721  t82 = 7875.0/4096.0*t12;
722  t83 = t67-t68-t69+t70-t71-t72+t73-t74-t75+t76+t77-t78-t79+t80+t81+t82;
723  t84 = 2625.0/64.0*t33;
724  t85 = 56875.0/1024.0*t35;
725  t86 = 17325.0/1024.0*t23;
726  t87 = 13875.0/1024.0*t38;
727  t88 = 13125.0/4096.0*t40;
728  t89 = 4375.0/2048.0*t42;
729  t90 = 34125.0/1024.0*t44;
730  t91 = 625.0/1024.0*t46;
731  t92 = 4375.0/3072.0*t48;
732  t93 = 8125.0/512.0*t50;
733  t94 = 56875.0/1536.0*t55;
734  t95 = 15625.0/1024.0*t57;
735  t96 = 102375.0/2048.0*t59;
736  t97 = 9375.0/1024.0*t61;
737  t98 = 65625.0/2048.0*t63;
738  t99 = 375.0/4096.0+t84-t85+t86-t87-t88+t89+t90-t91+t92+t93-t94-t95-t96+t97-t98;
739  t101 = 125.0/256.0*eta;
740  t102 = 11375.0/1536.0*t3;
741  t103 = 5625.0/2048.0*xi;
742  t104 = 10625.0/1024.0*t7;
743  t105 = 15625.0/2048.0*t10;
744  t106 = 328125.0/2048.0*t13;
745  t107 = 109375.0/1536.0*t15;
746  t108 = 7875.0/1024.0*t2;
747  t109 = 2125.0/1024.0*t6;
748  t110 = 109375.0/1024.0*t19;
749  t111 = 223125.0/1024.0*t21;
750  t112 = 5775.0/512.0*t24;
751  t113 = 21875.0/1536.0*t26;
752  t114 = 65625.0/2048.0*t28;
753  t115 = 3125.0/2048.0*t9;
754  t116 = 23625.0/2048.0*t12;
755  t117 = t101-t102+t103-t104+t105+t106-t107+t108+t109-t110-t111-t112+t113-t114-t115-t116;
756  t118 = 875.0/128.0*t33;
757  t119 = 74375.0/512.0*t35;
758  t120 = 5775.0/512.0*t23;
759  t121 = 2875.0/512.0*t38;
760  t122 = 118125.0/2048.0*t40;
761  t123 = 39375.0/1024.0*t42;
762  t124 = 14875.0/512.0*t44;
763  t125 = 5625.0/512.0*t46;
764  t126 = 13125.0/512.0*t48;
765  t127 = 10625.0/256.0*t50;
766  t128 = 74375.0/768.0*t55;
767  t129 = 15625.0/512.0*t57;
768  t130 = 44625.0/1024.0*t59;
769  t131 = 3125.0/512.0*t61;
770  t132 = 21875.0/1024.0*t63;
771  t133 = -1125.0/2048.0-t118+t119+t120+t121+t122-t123-t124+t125-t126-t127+t128+t129+t130-t131+t132;
772  t135 = t101-t102-t103+t104-t105-t106+t107+t108+t109+t110+t111-t112+t113-t114-t115-t116;
773  t136 = -1125.0/2048.0-t118-t119+t120+t121-t122+t123-t124-t125+t126+t127-t128-t129+t130-t131+t132;
774  t138 = t67-t68+t69-t70+t71+t72-t73-t74-t75-t76-t77-t78-t79+t80+t81+t82;
775  t139 = 375.0/4096.0+t84+t85+t86-t87+t88-t89+t90+t91-t92-t93+t94+t95-t96+t97-t98;
776  t141 = t1-t4+t5-t8+t11+t14-t16-t17-t18-t20-t22-t25-t27+t29+t30+t31;
777  t142 = 333.0/2048.0+t34+t36+t37-t39+t41-t43+t45+t47-t49-t51-t54+t56+t58-t60+t62-t64;
778  t144 = 625.0/512.0*eta;
779  t145 = 1875.0/1024.0*t3;
780  t146 = 4375.0/6144.0*t7;
781  t147 = 2625.0/4096.0*t10;
782  t148 = 39375.0/1024.0*t15;
783  t149 = 4875.0/2048.0*t2;
784  t150 = 625.0/2048.0*t6;
785  t151 = 102375.0/2048.0*t19;
786  t152 = 546875.0/6144.0*t21;
787  t153 = 546875.0/12288.0*t28;
788  t154 = 4375.0/12288.0*t9;
789  t155 = 15625.0/4096.0*t12;
790  t156 = t144-t145+t69-t146+t147+t72-t148-t149-t150-t151-t152-t79+t153+t154+t155;
791  t157 = 9375.0/512.0*t33;
792  t158 = 78125.0/4096.0*t40;
793  t159 = 24375.0/2048.0*t42;
794  t160 = 24375.0/1024.0*t44;
795  t161 = 4875.0/1024.0*t46;
796  t162 = 9375.0/1024.0*t48;
797  t163 = 11375.0/512.0*t50;
798  t164 = 5775.0/1024.0*t53;
799  t165 = 21875.0/512.0*t55;
800  t166 = 20475.0/1024.0*t57;
801  t167 = 78125.0/2048.0*t59;
802  t168 = 2625.0/128.0*t61;
803  t169 = 56875.0/2048.0*t63;
804  t170 = 125.0/4096.0+t157+t85-t87+t158-t159+t160+t161-t162-t163-t164+t165+t166-t167+t168-t169;
805  t172 = 125.0/192.0*eta;
806  t173 = 625.0/512.0*t3;
807  t174 = 13125.0/1024.0*t7;
808  t175 = 23625.0/2048.0*t10;
809  t176 = 13125.0/512.0*t15;
810  t177 = 6375.0/1024.0*t2;
811  t178 = 5625.0/1024.0*t6;
812  t179 = 133875.0/1024.0*t19;
813  t180 = 546875.0/3072.0*t21;
814  t181 = 546875.0/6144.0*t28;
815  t182 = 13125.0/2048.0*t9;
816  t183 = 15625.0/2048.0*t12;
817  t184 = -t172+t173-t103+t174-t175-t106+t176+t177+t178+t179+t180+t113-t181-t182-t183;
818  t185 = 3125.0/256.0*t33;
819  t186 = 78125.0/2048.0*t40;
820  t187 = 31875.0/1024.0*t42;
821  t188 = 31875.0/512.0*t44;
822  t189 = 2125.0/512.0*t46;
823  t190 = 3125.0/512.0*t48;
824  t191 = 14875.0/768.0*t50;
825  t192 = 1925.0/512.0*t53;
826  t193 = 21875.0/768.0*t55;
827  t194 = 8925.0/512.0*t57;
828  t195 = 78125.0/1024.0*t59;
829  t196 = 875.0/256.0*t61;
830  t197 = 74375.0/1024.0*t63;
831  t198 = -1125.0/2048.0-t185-t119+t121-t186+t187-t188-t189+t190+t191-t192-t193-t194+t195-t196+t197;
832  t200 = -t172+t173+t103-t174+t175+t106+t176-t177-t178-t179-t180+t113+t181+t182+t183;
833  t201 = 1125.0/2048.0-t185+t119+t121+t186-t187+t188-t189+t190+t191-t192-t193-t194-t195-t196-t197;
834  t203 = t144-t145-t69+t146-t147-t72-t148+t149+t150+t151+t152-t79-t153-t154-t155;
835  t204 = -125.0/4096.0+t157-t85-t87-t158+t159-t160+t161-t162-t163-t164+t165+t166+t167+t168+t169;
836  t206 = t1-t4-t5+t8-t11-t14-t16+t17+t18+t20+t22-t25-t27-t29-t30-t31;
837  t207 = -333.0/2048.0+t34-t36+t37-t39-t41+t43-t45+t47-t49-t51-t54+t56+t58+t60+t62+t64;
838  t209 = -375.0/4096.0+t67-t68-t69+t70-t71-t72-t73+t74+t75+t76+t77-t78-t79-t80-t81;
839  t210 = -t82+t84-t85+t86-t87-t88+t89-t90+t91-t92-t93+t94+t95+t96+t97+t98;
840  t212 = t101-t102+t103-t104+t105+t106+t107-t108-t109-t110-t111-t112+t113+t114+t115+t116;
841  t213 = 1125.0/2048.0-t118+t119+t120+t121+t122-t123+t124-t125+t126+t127-t128-t129-t130-t131-t132;
842  t215 = t101-t102-t103+t104-t105-t106-t107-t108-t109+t110+t111-t112+t113+t114+t115+t116;
843  t216 = 1125.0/2048.0-t118-t119+t120+t121-t122+t123+t124+t125-t126-t127+t128+t129-t130-t131-t132;
844  t218 = -375.0/4096.0+t67-t68+t69-t70+t71+t72+t73+t74+t75-t76-t77-t78-t79-t80-t81;
845  t219 = -t82+t84+t85+t86-t87+t88-t89-t90-t91+t92+t93-t94-t95+t96+t97+t98;
846  t221 = t1-t4+t5-t8+t11+t14+t16+t17+t18-t20-t22-t25-t27-t29-t30-t31;
847  t222 = -333.0/2048.0+t34+t36+t37-t39+t41-t43-t45-t47+t49+t51-t54-t56-t58+t60+t62+t64;
848  t224 = t144-t145+t69-t146+t147+t72+t148+t149+t150-t151-t152-t79-t153-t154-t155;
849  t225 = -125.0/4096.0+t157+t85-t87+t158-t159-t160-t161+t162+t163-t164-t165-t166+t167+t168+t169;
850  t227 = -t172+t173-t103+t174-t175-t106-t176-t177-t178+t179+t180+t113+t181+t182+t183;
851  t228 = 1125.0/2048.0-t185-t119+t121-t186+t187+t188+t189-t190-t191-t192+t193+t194-t195-t196-t197;
852  t230 = -t172+t173+t103-t174+t175+t106-t176+t177+t178-t179-t180+t113-t181-t182-t183;
853  t231 = -1125.0/2048.0-t185+t119+t121+t186-t187-t188+t189-t190-t191-t192+t193+t194+t195-t196+t197;
854  t233 = t144-t145-t69+t146-t147-t72+t148-t149-t150+t151+t152-t79+t153+t154+t155;
855  t234 = 125.0/4096.0+t157-t85-t87-t158+t159+t160-t161+t162+t163-t164-t165-t166-t167+t168-t169;
856  t255 = 5625.0/64.0*t38;
857  t256 = 6125.0/64.0*t26;
858  t299 = 3969.0/256.0*xi;
859  t302 = 99225.0/64.0*t35;
860  t304 = 416745.0/256.0*t13;
861  t372 = 4455.0/64.0*t38;
862 
863  values[0] = t32+t65;
864  values[1] = t83+t99;
865  values[2] = t117+t133;
866  values[3] = t135+t136;
867  values[4] = t138+t139;
868  values[5] = t141+t142;
869  values[6] = t156+t170;
870  values[7] = t184+t198;
871  values[8] = t200+t201;
872  values[9] = t203+t204;
873  values[10] = t206+t207;
874  values[11] = t209+t210;
875  values[12] = t212+t213;
876  values[13] = t215+t216;
877  values[14] = t218+t219;
878  values[15] = t221+t222;
879  values[16] = t224+t225;
880  values[17] = t227+t228;
881  values[18] = t230+t231;
882  values[19] = t233+t234;
883  values[20] = 45.0/2.0*t38+1225.0/64.0*t26-453.0/64.0*eta-105.0/2.0*t33+693.0/64.0*t53+2135.0/64.0*t3-2079.0/64.0*t23+2079.0/64.0*t24-105.0/4.0*t61;
884  values[21] = -105.0/256.0-4725.0/256.0*t12+1575.0/128.0*t2+735.0/256.0*t9-315.0/128.0*t6+33075.0/256.0*t28-11025.0/128.0*t63-14175.0/128.0*t59+4725.0/64.0*t44;
885  values[22] = t255+t256-975.0/32.0*eta-525.0/2.0*t33+10675.0/64.0*t3-10395.0/64.0*t23+10395.0/64.0*t24-3675.0/64.0*t61;
886  values[23] = -945.0/256.0-11025.0/256.0*t12+4725.0/128.0*t2-33075.0/128.0*t63+6615.0/256.0*t9+14175.0/64.0*t44-2835.0/128.0*t6+77175.0/256.0*t28-33075.0/128.0*t59;
887  values[24] = 735.0/64.0*t48-315.0/64.0*t46+6615.0/64.0*t15-2835.0/64.0*t57-3675.0/32.0*t55+1575.0/32.0*t50;
888  values[25] = 441.0/256.0*xi+19845.0/256.0*t40-6615.0/128.0*t42+3969.0/256.0*t10-2205.0/128.0*t7+178605.0/256.0*t13-59535.0/128.0*t19-99225.0/128.0*t21+33075.0/64.0*t35;
889  values[26] = -2205.0/64.0*t46+3675.0/64.0*t48-19845.0/64.0*t57+11025.0/32.0*t50+33075.0/64.0*t15-18375.0/32.0*t55;
890  values[27] = 46305.0/256.0*t40-19845.0/128.0*t42+t299-178605.0/128.0*t19+35721.0/256.0*t10+t302-19845.0/128.0*t7+t304-231525.0/128.0*t21;
891  values[28] = t255+t256-345.0/32.0*eta-3675.0/32.0*t33+3465.0/64.0*t53+1225.0/64.0*t3-525.0/4.0*t61;
892  values[29] = 735.0/256.0+3675.0/256.0*t9-2205.0/128.0*t6+33075.0/256.0*t12-11025.0/128.0*t2-99225.0/128.0*t59+33075.0/64.0*t44+165375.0/256.0*t28-55125.0/128.0*t63;
893  values[30] = 11025.0/32.0*t38+30625.0/64.0*t26-3675.0/64.0*eta-18375.0/32.0*t33+6125.0/64.0*t3-18375.0/64.0*t61;
894  values[31] = 6615.0/256.0-231525.0/128.0*t59-165375.0/128.0*t63+33075.0/256.0*t9+77175.0/256.0*t12-33075.0/128.0*t2-19845.0/128.0*t6+99225.0/64.0*t44+385875.0/256.0*t28;
895  values[32] = -11025.0/32.0*t55+4725.0/32.0*t50+6615.0/64.0*t48-2835.0/64.0*t46+15435.0/64.0*t15-6615.0/64.0*t57;
896  values[33] = -297675.0/128.0*t21+t302-6615.0/128.0*t7+178605.0/256.0*t40-59535.0/128.0*t42+t299+9261.0/256.0*t10+t304-138915.0/128.0*t19;
897  values[34] = 33075.0/64.0*t48-55125.0/32.0*t55-46305.0/64.0*t57+33075.0/32.0*t50-19845.0/64.0*t46+77175.0/64.0*t15;
898  values[35] = 416745.0/256.0*t40+35721.0/256.0*xi-694575.0/128.0*t21+972405.0/256.0*t13-178605.0/128.0*t42+83349.0/256.0*t10-416745.0/128.0*t19-59535.0/128.0*t7+297675.0/64.0*t35;
899  values[36] = 10395.0/32.0*t3-4455.0/64.0*eta-10395.0/32.0*t33+t372-18711.0/64.0*t23+18711.0/64.0*t24;
900  values[37] = -297.0/64.0*eta-10395.0/64.0*t61+t372+6237.0/64.0*t53;
901 }
902 
903 // values of the derivatives in xi-xi direction
904 static void C_Q_UL5_2D_DeriveXiXi(double xi, double eta, double *values)
905 {
906  double t1, t2, t3, t4, t5, t6, t7, t8, t9, t10, t11, t12, t13, t14, t15;
907  double t16, t17, t18, t19, t20, t21, t22, t23, t24, t25, t26, t27, t28;
908  double t29, t30, t31, t32, t33, t34, t35, t36, t37, t38, t39, t40, t41;
909  double t42, t43, t44, t45, t46, t47, t48, t49, t50, t51, t52, t54, t55;
910  double t56, t57, t58, t59, t60, t61, t62, t63, t64, t65, t66, t67, t68;
911  double t69, t70, t71, t72, t73, t74, t75, t76, t77, t78, t79, t81, t82;
912  double t83, t84, t85, t86, t87, t88, t89, t90, t91, t92, t93, t94, t95;
913  double t96, t97, t98, t99, t100, t101, t102, t103, t104, t105, t106, t108;
914  double t109, t111, t112, t114, t115, t117, t118, t119, t120, t121, t122;
915  double t123, t124, t125, t126, t127, t128, t129, t130, t131, t132, t133;
916  double t134, t135, t136, t137, t138, t139, t141, t142, t143, t144, t145;
917  double t146, t147, t148, t149, t150, t151, t152, t153, t154, t155, t156;
918  double t157, t158, t159, t160, t161, t162, t163, t165, t166, t168, t169;
919  double t171, t172, t174, t175, t177, t178, t180, t181, t183, t184, t186;
920  double t187, t189, t190, t192, t193, t195, t196, t198, t199, t217, t218;
921  double t254, t256, t312;
922 
923  t1 = xi*xi;
924  t2 = t1*t1;
925  t3 = 21945.0/1024.0*t2;
926  t4 = eta*eta;
927  t5 = t1*t4;
928  t6 = 18765.0/512.0*t5;
929  t7 = xi*t4;
930  t8 = 9475.0/256.0*t7;
931  t9 = t1*eta;
932  t10 = 7755.0/512.0*t9;
933  t11 = t1*xi;
934  t12 = t11*t4;
935  t13 = 23425.0/256.0*t12;
936  t14 = t4*t4;
937  t15 = t14*eta;
938  t16 = xi*t15;
939  t17 = 38045.0/512.0*t16;
940  t18 = t1*t15;
941  t19 = 22295.0/512.0*t18;
942  t20 = t11*t15;
943  t21 = 93135.0/512.0*t20;
944  t22 = t1*t14;
945  t23 = 7525.0/512.0*t22;
946  t24 = t4*eta;
947  t25 = t11*t24;
948  t26 = 190225.0/768.0*t25;
949  t27 = t2*t4;
950  t28 = 21945.0/1024.0*t27;
951  t29 = xi*t14;
952  t30 = 16275.0/512.0*t29;
953  t31 = t11*eta;
954  t32 = 29515.0/512.0*t31;
955  t33 = t3+t6+t8+t10-t13+t17+t19-t21-t23+t26-t28-t30-t32;
956  t34 = xi*t24;
957  t35 = 25025.0/256.0*t34;
958  t36 = xi*eta;
959  t37 = 11505.0/512.0*t36;
960  t38 = 5525.0/512.0*t11;
961  t39 = 2175.0/512.0*xi;
962  t40 = 1465.0/512.0*eta;
963  t41 = 9475.0/768.0*t24;
964  t42 = 2185.0/128.0*t1;
965  t43 = 4685.0/512.0*t15;
966  t44 = t1*t24;
967  t45 = 16275.0/256.0*t44;
968  t46 = 6255.0/1024.0*t14;
969  t47 = t14*t4;
970  t48 = 1463.0/1024.0*t47;
971  t49 = t11*t14;
972  t50 = 111475.0/1536.0*t49;
973  t51 = 3395.0/512.0*t4;
974  t52 = 2497.0/1536.0-t35+t37+t38-t39-t40+t41-t42-t43-t45+t46-t48+t50-t51;
975  t54 = 28125.0/512.0*t5;
976  t55 = 24375.0/512.0*t7;
977  t56 = 28125.0/1024.0*t9;
978  t57 = 78125.0/512.0*t12;
979  t58 = 102375.0/1024.0*t16;
980  t59 = 118125.0/1024.0*t18;
981  t60 = 328125.0/1024.0*t20;
982  t61 = 65625.0/1024.0*t22;
983  t62 = 546875.0/1536.0*t25;
984  t63 = 56875.0/1024.0*t29;
985  t64 = 78125.0/1024.0*t31;
986  t65 = 56875.0/512.0*t34;
987  t66 = t54+t55+t56-t57+t58+t59-t60-t61+t62-t63-t64-t65;
988  t67 = 24375.0/1024.0*t36;
989  t68 = 15625.0/1024.0*t11;
990  t69 = 4875.0/1024.0*xi;
991  t70 = 4875.0/1024.0*eta;
992  t71 = 11375.0/512.0*t24;
993  t72 = 5625.0/1024.0*t1;
994  t73 = 20475.0/1024.0*t15;
995  t74 = 65625.0/512.0*t44;
996  t75 = 2625.0/128.0*t14;
997  t76 = 5775.0/1024.0*t47;
998  t77 = 546875.0/3072.0*t49;
999  t78 = 13875.0/1024.0*t4;
1000  t79 = 625.0/512.0+t67+t68-t69-t70+t71-t72-t73-t74+t75-t76+t77-t78;
1001  t81 = 9375.0/256.0*t5;
1002  t82 = 31875.0/256.0*t7;
1003  t83 = 9375.0/512.0*t9;
1004  t84 = 78125.0/256.0*t12;
1005  t85 = 133875.0/512.0*t16;
1006  t86 = 39375.0/512.0*t18;
1007  t87 = 328125.0/512.0*t20;
1008  t88 = 21875.0/512.0*t22;
1009  t89 = 546875.0/768.0*t25;
1010  t90 = 74375.0/512.0*t29;
1011  t91 = 78125.0/512.0*t31;
1012  t92 = 74375.0/256.0*t34;
1013  t93 = -t81-t82-t83+t84-t85-t86+t87+t88-t89+t90+t91+t92;
1014  t94 = 31875.0/512.0*t36;
1015  t95 = 15625.0/512.0*t11;
1016  t96 = 6375.0/512.0*xi;
1017  t97 = 2125.0/512.0*eta;
1018  t98 = 14875.0/768.0*t24;
1019  t99 = 1875.0/512.0*t1;
1020  t100 = 8925.0/512.0*t15;
1021  t101 = 21875.0/256.0*t44;
1022  t102 = 875.0/256.0*t14;
1023  t103 = 1925.0/512.0*t47;
1024  t104 = 546875.0/1536.0*t49;
1025  t105 = 2875.0/512.0*t4;
1026  t106 = -125.0/192.0-t94-t95+t96+t97-t98+t99+t100+t101-t102-t103-t104+t105;
1027  t108 = -t81+t82-t83-t84+t85-t86-t87+t88+t89-t90-t91-t92;
1028  t109 = -125.0/192.0+t94+t95-t96+t97-t98+t99+t100+t101-t102-t103+t104+t105;
1029  t111 = t54-t55+t56+t57-t58+t59+t60-t61-t62+t63+t64+t65;
1030  t112 = 625.0/512.0-t67-t68+t69-t70+t71-t72-t73-t74+t75-t76-t77-t78;
1031  t114 = t3+t6-t8+t10+t13-t17+t19+t21-t23-t26-t28+t30+t32;
1032  t115 = 2497.0/1536.0+t35-t37-t38+t39-t40+t41-t42-t43-t45+t46-t48-t50-t51;
1033  t117 = 86625.0/1024.0*t2;
1034  t118 = 7875.0/64.0*t5;
1035  t119 = 34125.0/512.0*t7;
1036  t120 = 4375.0/1024.0*t9;
1037  t121 = 102375.0/512.0*t12;
1038  t122 = 109375.0/1024.0*t16;
1039  t123 = 109375.0/1024.0*t18;
1040  t124 = 170625.0/512.0*t25;
1041  t125 = 86625.0/1024.0*t27;
1042  t126 = 65625.0/1024.0*t29;
1043  t127 = 13125.0/1024.0*t31;
1044  t128 = t117+t118-t119+t120+t121-t122+t123+t60-t61-t124-t125+t126+t127;
1045  t129 = 4375.0/1024.0*t36;
1046  t130 = 7875.0/1024.0*t11;
1047  t131 = 2625.0/1024.0*xi;
1048  t132 = 625.0/1024.0*eta;
1049  t133 = 8125.0/512.0*t24;
1050  t134 = 60375.0/1024.0*t1;
1051  t135 = 15625.0/1024.0*t15;
1052  t136 = 56875.0/512.0*t44;
1053  t137 = 9375.0/1024.0*t14;
1054  t138 = 196875.0/1024.0*t49;
1055  t139 = 1125.0/256.0+t65-t129-t130+t131-t132+t133-t134-t135-t136+t137-t138-t78;
1056  t141 = 28875.0/512.0*t2;
1057  t142 = 2625.0/128.0*t5;
1058  t143 = 14875.0/256.0*t7;
1059  t144 = 39375.0/512.0*t9;
1060  t145 = 44625.0/256.0*t12;
1061  t146 = 109375.0/512.0*t16;
1062  t147 = 109375.0/512.0*t18;
1063  t148 = 223125.0/256.0*t25;
1064  t149 = 28875.0/512.0*t27;
1065  t150 = 21875.0/512.0*t29;
1066  t151 = 125.0/256.0+t141-t142+t143-t144-t145+t146-t147-t87+t88+t148-t149-t150;
1067  t152 = 118125.0/512.0*t31;
1068  t153 = 39375.0/512.0*t36;
1069  t154 = 23625.0/512.0*t11;
1070  t155 = 7875.0/512.0*xi;
1071  t156 = 5625.0/512.0*eta;
1072  t157 = 10625.0/256.0*t24;
1073  t158 = 11375.0/512.0*t1;
1074  t159 = 15625.0/512.0*t15;
1075  t160 = 74375.0/256.0*t44;
1076  t161 = 3125.0/512.0*t14;
1077  t162 = 65625.0/512.0*t49;
1078  t163 = -t152-t92+t153+t154-t155+t156-t157-t158+t159+t160-t161+t162+t105;
1079  t165 = 125.0/256.0+t141-t142+t143+t144-t145-t146+t147+t87+t88-t148-t149-t150;
1080  t166 = t152+t92-t153+t154-t155-t156+t157-t158-t159-t160-t161+t162+t105;
1081  t168 = t117+t118-t119-t120+t121+t122-t123-t60-t61+t124-t125+t126-t127;
1082  t169 = 1125.0/256.0-t65+t129-t130+t131+t132-t133-t134+t135+t136+t137-t138-t78;
1083  t171 = t3+t6-t8-t10+t13+t17-t19-t21-t23+t26-t28+t30-t32;
1084  t172 = 2497.0/1536.0-t35+t37-t38+t39+t40-t41-t42+t43+t45+t46-t48-t50-t51;
1085  t174 = t54-t55-t56+t57+t58-t59-t60-t61+t62+t63-t64-t65;
1086  t175 = 625.0/512.0+t67-t68+t69+t70-t71-t72+t73+t74+t75-t76-t77-t78;
1087  t177 = -t81+t82+t83-t84-t85+t86+t87+t88-t89-t90+t91+t92;
1088  t178 = -125.0/192.0-t94+t95-t96-t97+t98+t99-t100-t101-t102-t103+t104+t105;
1089  t180 = -t81-t82+t83+t84+t85+t86-t87+t88+t89+t90-t91-t92;
1090  t181 = -125.0/192.0+t94-t95+t96-t97+t98+t99-t100-t101-t102-t103-t104+t105;
1091  t183 = t54+t55-t56-t57-t58-t59+t60-t61-t62-t63+t64+t65;
1092  t184 = 625.0/512.0-t67+t68-t69+t70-t71-t72+t73+t74+t75-t76+t77-t78;
1093  t186 = t3+t6+t8-t10-t13-t17-t19+t21-t23-t26-t28-t30+t32;
1094  t187 = 2497.0/1536.0+t35-t37+t38-t39+t40-t41-t42+t43+t45+t46-t48+t50-t51;
1095  t189 = t117+t118+t119-t120-t121-t122-t123+t60-t61-t124-t125-t126+t127;
1096  t190 = 1125.0/256.0+t65-t129+t130-t131+t132-t133-t134+t135+t136+t137+t138-t78;
1097  t192 = 125.0/256.0+t141-t142-t143+t144+t145+t146+t147-t87+t88+t148-t149+t150;
1098  t193 = -t152-t92+t153-t154+t155-t156+t157-t158-t159-t160-t161-t162+t105;
1099  t195 = 125.0/256.0+t141-t142-t143-t144+t145-t146-t147+t87+t88-t148-t149+t150;
1100  t196 = t152+t92-t153-t154+t155+t156-t157-t158+t159+t160-t161-t162+t105;
1101  t198 = t117+t118+t119+t120-t121+t122+t123-t60-t61+t124-t125-t126-t127;
1102  t199 = 1125.0/256.0-t65+t129+t130-t131-t132+t133-t134-t135-t136+t137+t138-t78;
1103  t217 = 5625.0/64.0*t4;
1104  t218 = 18375.0/64.0*t22;
1105  t254 = 99225.0/32.0*t34;
1106  t256 = 416745.0/64.0*t20;
1107  t312 = 4455.0/64.0*t4;
1108 
1109  values[0] = t33+t52;
1110  values[1] = t66+t79;
1111  values[2] = t93+t106;
1112  values[3] = t108+t109;
1113  values[4] = t111+t112;
1114  values[5] = t114+t115;
1115  values[6] = t128+t139;
1116  values[7] = t151+t163;
1117  values[8] = t165+t166;
1118  values[9] = t168+t169;
1119  values[10] = t171+t172;
1120  values[11] = t174+t175;
1121  values[12] = t177+t178;
1122  values[13] = t180+t181;
1123  values[14] = t183+t184;
1124  values[15] = t186+t187;
1125  values[16] = t189+t190;
1126  values[17] = t192+t193;
1127  values[18] = t195+t196;
1128  values[19] = t198+t199;
1129  values[20] = 45.0/2.0*t4+3675.0/64.0*t22+6405.0/64.0*t1-453.0/64.0-105.0/4.0*t14+10395.0/64.0*t27-10395.0/64.0*t2+693.0/64.0*t47-315.0/2.0*t5;
1130  values[21] = 2205.0/64.0*t9-315.0/64.0*eta+19845.0/64.0*t18-11025.0/32.0*t44-2835.0/64.0*t15+1575.0/32.0*t24;
1131  values[22] = t217+t218+3675.0/64.0*t1-345.0/32.0-525.0/4.0*t14+3465.0/64.0*t47-11025.0/32.0*t5;
1132  values[23] = -33075.0/32.0*t44+19845.0/64.0*t9+4725.0/32.0*t24-2835.0/64.0*eta+46305.0/64.0*t18-6615.0/64.0*t15;
1133  values[24] = -4725.0/64.0*t11+1575.0/64.0*xi+33075.0/64.0*t49-14175.0/32.0*t12-11025.0/64.0*t29+4725.0/32.0*t7;
1134  values[25] = 19845.0/64.0*t31-6615.0/64.0*t36+178605.0/64.0*t20-99225.0/32.0*t25-59535.0/64.0*t16+33075.0/32.0*t34;
1135  values[26] = 33075.0/64.0*t11-99225.0/32.0*t12-11025.0/64.0*xi+33075.0/32.0*t7+165375.0/64.0*t49-55125.0/64.0*t29;
1136  values[27] = -297675.0/32.0*t25+178605.0/64.0*t31+t254-59535.0/64.0*t36+t256-138915.0/64.0*t16;
1137  values[28] = t217+t218+32025.0/64.0*t1-975.0/32.0-3675.0/64.0*t14+51975.0/64.0*t27-51975.0/64.0*t2-1575.0/2.0*t5;
1138  values[29] = 11025.0/64.0*t9-2205.0/64.0*eta-19845.0/64.0*t15+11025.0/32.0*t24+99225.0/64.0*t18-55125.0/32.0*t44;
1139  values[30] = 11025.0/32.0*t4+91875.0/64.0*t22+18375.0/64.0*t1-3675.0/64.0-18375.0/64.0*t14-55125.0/32.0*t5;
1140  values[31] = -46305.0/64.0*t15-165375.0/32.0*t44+99225.0/64.0*t9-19845.0/64.0*eta+33075.0/32.0*t24+231525.0/64.0*t18;
1141  values[32] = -11025.0/64.0*t11+4725.0/64.0*xi-33075.0/64.0*t29+14175.0/32.0*t7+77175.0/64.0*t49-33075.0/32.0*t12;
1142  values[33] = -178605.0/64.0*t16+t254-19845.0/64.0*t36+46305.0/64.0*t31+t256-231525.0/32.0*t25;
1143  values[34] = -165375.0/64.0*t29-231525.0/32.0*t12+77175.0/64.0*t11-33075.0/64.0*xi+99225.0/32.0*t7+385875.0/64.0*t49;
1144  values[35] = -416745.0/64.0*t16+972405.0/64.0*t20+416745.0/64.0*t31-694575.0/32.0*t25-178605.0/64.0*t36+297675.0/32.0*t34;
1145  values[36] = -297.0/64.0-10395.0/64.0*t14+t312+6237.0/64.0*t47;
1146  values[37] = 31185.0/32.0*t1-4455.0/64.0-31185.0/32.0*t5+t312-93555.0/64.0*t2+93555.0/64.0*t27;
1147 }
1148 
1149 // values of the derivatives in xi-eta direction
1150 static void C_Q_UL5_2D_DeriveXiEta(double xi, double eta, double *values)
1151 {
1152  double t1, t2, t3, t4, t5, t6, t7, t8, t9, t10, t12, t13, t14, t15;
1153  double t16, t17, t18, t19, t20, t21, t22, t23, t24, t25, t26, t28;
1154  double t29, t30, t31, t32, t33, t34, t35, t36, t37, t38, t39, t40;
1155  double t41, t42, t43, t44, t45, t46, t47, t48, t49, t50, t51, t52;
1156  double t53, t54, t56, t57, t58, t59, t60, t61, t62, t63, t64, t65;
1157  double t66, t67, t68, t69, t70, t71, t72, t73, t74, t75, t76, t77;
1158  double t78, t79, t80, t81, t82, t84, t85, t86, t87, t88, t89, t90;
1159  double t91, t92, t93, t94, t95, t96, t97, t98, t99, t100, t101, t102;
1160  double t103, t104, t105, t106, t107, t108, t109, t110, t112, t113;
1161  double t115, t116, t118, t119, t121, t122, t123, t124, t125, t126;
1162  double t127, t128, t129, t130, t131, t132, t133, t134, t135, t136;
1163  double t137, t138, t139, t140, t141, t142, t143, t145, t146, t147;
1164  double t148, t149, t150, t151, t152, t153, t154, t155, t156, t157;
1165  double t158, t159, t160, t161, t162, t163, t164, t165, t166, t167;
1166  double t169, t170, t172, t173, t175, t176, t178, t179, t181, t182;
1167  double t184, t185, t187, t188, t190, t191, t193, t194, t196, t197;
1168  double t199, t200, t202, t203, t219, t220, t259, t261, t318;
1169 
1170  t1 = xi*xi;
1171  t2 = 11505.0/1024.0*t1;
1172  t3 = eta*eta;
1173  t4 = xi*t3;
1174  t5 = 9475.0/256.0*t4;
1175  t6 = t1*xi;
1176  t7 = 2585.0/512.0*t6;
1177  t8 = t3*eta;
1178  t9 = 2585.0/512.0*t8;
1179  t10 = t1*t1;
1180  t12 = t10*xi*eta;
1181  t13 = 4389.0/512.0*t12;
1182  t14 = t6*eta;
1183  t15 = 6255.0/256.0*t14;
1184  t16 = xi*t8;
1185  t17 = 6255.0/256.0*t16;
1186  t18 = xi*eta;
1187  t19 = 3395.0/256.0*t18;
1188  t20 = t1*eta;
1189  t21 = 9475.0/256.0*t20;
1190  t22 = t6*t3;
1191  t23 = 16275.0/256.0*t22;
1192  t24 = t1*t8;
1193  t25 = 16275.0/256.0*t24;
1194  t26 = t3*t3;
1195  t28 = xi*t26*eta;
1196  t29 = 4389.0/512.0*t28;
1197  t30 = -1755.0/2048.0+t2+t5+t7+t9-t13+t15+t17-t19+t21-t23-t25-t29;
1198  t31 = t10*t26;
1199  t32 = 465675.0/2048.0*t31;
1200  t33 = t10*t8;
1201  t34 = 111475.0/1536.0*t33;
1202  t35 = t6*t26;
1203  t36 = 111475.0/1536.0*t35;
1204  t37 = t1*t26;
1205  t38 = 190225.0/1024.0*t37;
1206  t39 = t10*t3;
1207  t40 = 190225.0/1024.0*t39;
1208  t41 = xi*t26;
1209  t42 = 23425.0/512.0*t41;
1210  t43 = t10*eta;
1211  t44 = 23425.0/512.0*t43;
1212  t45 = t6*t8;
1213  t46 = 7525.0/384.0*t45;
1214  t47 = 1465.0/512.0*xi;
1215  t48 = 1465.0/512.0*eta;
1216  t49 = 11505.0/1024.0*t3;
1217  t50 = t1*t3;
1218  t51 = 75075.0/512.0*t50;
1219  t52 = 29515.0/2048.0*t10;
1220  t53 = 29515.0/2048.0*t26;
1221  t54 = -t32+t34+t36+t38+t40-t42-t44-t46-t47-t48+t49-t51-t52-t53;
1222  t56 = 24375.0/2048.0*t1;
1223  t57 = 34125.0/512.0*t4;
1224  t58 = 9375.0/1024.0*t6;
1225  t59 = 4375.0/3072.0*t8;
1226  t60 = 9375.0/256.0*t14;
1227  t61 = 2625.0/32.0*t16;
1228  t62 = 13875.0/512.0*t18;
1229  t63 = 24375.0/512.0*t20;
1230  t64 = 65625.0/512.0*t22;
1231  t65 = 56875.0/512.0*t24;
1232  t66 = 17325.0/512.0*t28;
1233  t67 = 1640625.0/4096.0*t31;
1234  t68 = 546875.0/3072.0*t33;
1235  t69 = t56+t57+t58+t59+t60+t61-t62+t63-t64-t65-t66-t67+t68;
1236  t70 = 196875.0/1024.0*t35;
1237  t71 = 511875.0/2048.0*t37;
1238  t72 = 546875.0/2048.0*t39;
1239  t73 = 102375.0/1024.0*t41;
1240  t74 = 78125.0/1024.0*t43;
1241  t75 = 21875.0/256.0*t45;
1242  t76 = 4875.0/1024.0*xi;
1243  t77 = 625.0/1024.0*eta;
1244  t78 = 4375.0/2048.0*t3;
1245  t79 = 170625.0/1024.0*t50;
1246  t80 = 78125.0/4096.0*t10;
1247  t81 = 13125.0/4096.0*t26;
1248  t82 = -625.0/4096.0+t70+t71+t72-t73-t74-t75-t76-t77+t78-t79-t80-t81;
1249  t84 = 31875.0/1024.0*t1;
1250  t85 = 14875.0/256.0*t4;
1251  t86 = 3125.0/512.0*t6;
1252  t87 = 13125.0/512.0*t8;
1253  t88 = 3125.0/128.0*t14;
1254  t89 = 875.0/64.0*t16;
1255  t90 = 2875.0/256.0*t18;
1256  t91 = 31875.0/256.0*t20;
1257  t92 = 21875.0/256.0*t22;
1258  t93 = 74375.0/256.0*t24;
1259  t94 = 5775.0/256.0*t28;
1260  t95 = 1640625.0/2048.0*t31;
1261  t96 = 546875.0/1536.0*t33;
1262  t97 = -t84-t85-t86-t87-t88-t89+t90-t91+t92+t93-t94+t95-t96;
1263  t98 = 65625.0/512.0*t35;
1264  t99 = 669375.0/1024.0*t37;
1265  t100 = 546875.0/1024.0*t39;
1266  t101 = 44625.0/512.0*t41;
1267  t102 = 78125.0/512.0*t43;
1268  t103 = 21875.0/384.0*t45;
1269  t104 = 2125.0/512.0*xi;
1270  t105 = 5625.0/512.0*eta;
1271  t106 = 39375.0/1024.0*t3;
1272  t107 = 223125.0/512.0*t50;
1273  t108 = 78125.0/2048.0*t10;
1274  t109 = 118125.0/2048.0*t26;
1275  t110 = 5625.0/2048.0-t98-t99-t100+t101+t102+t103+t104+t105-t106+t107+t108+t109;
1276  t112 = t84-t85-t86+t87-t88-t89+t90+t91+t92-t93-t94-t95+t96;
1277  t113 = -5625.0/2048.0-t98+t99+t100+t101-t102+t103+t104-t105+t106-t107-t108-t109;
1278  t115 = -t56+t57+t58-t59+t60+t61-t62-t63-t64+t65-t66+t67-t68;
1279  t116 = 625.0/4096.0+t70-t71-t72-t73+t74-t75-t76+t77-t78+t79+t80+t81;
1280  t118 = -t2+t5+t7-t9-t13+t15+t17-t19-t21-t23+t25-t29+t32;
1281  t119 = 1755.0/2048.0-t34+t36-t38-t40-t42+t44-t46-t47+t48-t49+t51+t52+t53;
1282  t121 = 4375.0/2048.0*t1;
1283  t122 = 24375.0/512.0*t4;
1284  t123 = 4375.0/3072.0*t6;
1285  t124 = 9375.0/1024.0*t8;
1286  t125 = 17325.0/512.0*t12;
1287  t126 = 2625.0/32.0*t14;
1288  t127 = 9375.0/256.0*t16;
1289  t128 = 34125.0/512.0*t20;
1290  t129 = 56875.0/512.0*t22;
1291  t130 = 65625.0/512.0*t24;
1292  t131 = 196875.0/1024.0*t33;
1293  t132 = -t121+t122+t123-t124-t125+t126+t127-t62-t128-t129+t130+t67-t131;
1294  t133 = 546875.0/3072.0*t35;
1295  t134 = 546875.0/2048.0*t37;
1296  t135 = 511875.0/2048.0*t39;
1297  t136 = 78125.0/1024.0*t41;
1298  t137 = 102375.0/1024.0*t43;
1299  t138 = 625.0/1024.0*xi;
1300  t139 = 4875.0/1024.0*eta;
1301  t140 = 24375.0/2048.0*t3;
1302  t141 = 13125.0/4096.0*t10;
1303  t142 = 78125.0/4096.0*t26;
1304  t143 = 625.0/4096.0+t133-t134-t135-t136+t137-t75-t138+t139-t140+t79+t141+t142;
1305  t145 = 39375.0/1024.0*t1;
1306  t146 = 31875.0/256.0*t4;
1307  t147 = 13125.0/512.0*t6;
1308  t148 = 3125.0/512.0*t8;
1309  t149 = 5775.0/256.0*t12;
1310  t150 = 875.0/64.0*t14;
1311  t151 = 3125.0/128.0*t16;
1312  t152 = 14875.0/256.0*t20;
1313  t153 = 74375.0/256.0*t22;
1314  t154 = 21875.0/256.0*t24;
1315  t155 = 65625.0/512.0*t33;
1316  t156 = t145-t146-t147+t148-t149-t150-t151+t90+t152+t153-t154-t95+t155;
1317  t157 = 546875.0/1536.0*t35;
1318  t158 = 546875.0/1024.0*t37;
1319  t159 = 669375.0/1024.0*t39;
1320  t160 = 78125.0/512.0*t41;
1321  t161 = 44625.0/512.0*t43;
1322  t162 = 5625.0/512.0*xi;
1323  t163 = 2125.0/512.0*eta;
1324  t164 = 31875.0/1024.0*t3;
1325  t165 = 118125.0/2048.0*t10;
1326  t166 = 78125.0/2048.0*t26;
1327  t167 = -5625.0/2048.0-t157+t158+t159+t160-t161+t103+t162-t163+t164-t107-t165-t166;
1328  t169 = -t145+t146+t147+t148-t149-t150-t151+t90+t152-t153-t154+t95+t155;
1329  t170 = 5625.0/2048.0+t157-t158-t159-t160-t161+t103-t162-t163-t164+t107+t165+t166;
1330  t172 = t121-t122-t123-t124-t125+t126+t127-t62-t128+t129+t130-t67-t131;
1331  t173 = -625.0/4096.0-t133+t134+t135+t136+t137-t75+t138+t139+t140-t79-t141-t142;
1332  t175 = -1755.0/2048.0+t2-t5-t7-t9-t13+t15+t17-t19-t21+t23+t25-t29;
1333  t176 = -t32-t34-t36+t38+t40+t42+t44-t46+t47+t48+t49-t51-t52-t53;
1334  t178 = t56-t57-t58-t59+t60+t61-t62-t63+t64+t65-t66-t67-t68;
1335  t179 = -625.0/4096.0-t70+t71+t72+t73+t74-t75+t76+t77+t78-t79-t80-t81;
1336  t181 = -t84+t85+t86+t87-t88-t89+t90+t91-t92-t93-t94+t95+t96;
1337  t182 = 5625.0/2048.0+t98-t99-t100-t101-t102+t103-t104-t105-t106+t107+t108+t109;
1338  t184 = t84+t85+t86-t87-t88-t89+t90-t91-t92+t93-t94-t95-t96;
1339  t185 = -5625.0/2048.0+t98+t99+t100-t101+t102+t103-t104+t105+t106-t107-t108-t109;
1340  t187 = -t56-t57-t58+t59+t60+t61-t62+t63+t64-t65-t66+t67+t68;
1341  t188 = 625.0/4096.0-t70-t71-t72+t73-t74-t75+t76-t77-t78+t79+t80+t81;
1342  t190 = -t2-t5-t7+t9-t13+t15+t17-t19+t21+t23-t25-t29+t32;
1343  t191 = 1755.0/2048.0+t34-t36-t38-t40+t42-t44-t46+t47-t48-t49+t51+t52+t53;
1344  t193 = -t121-t122-t123+t124-t125+t126+t127-t62+t128+t129-t130+t67+t131;
1345  t194 = 625.0/4096.0-t133-t134-t135+t136-t137-t75+t138-t139-t140+t79+t141+t142;
1346  t196 = t145+t146+t147-t148-t149-t150-t151+t90-t152-t153+t154-t95-t155;
1347  t197 = -5625.0/2048.0+t157+t158+t159-t160+t161+t103-t162+t163+t164-t107-t165-t166;
1348  t199 = -t145-t146-t147-t148-t149-t150-t151+t90-t152+t153+t154+t95-t155;
1349  t200 = 5625.0/2048.0-t157-t158-t159+t160+t161+t103+t162+t163-t164+t107+t165+t166;
1350  t202 = t121+t122+t123+t124-t125+t126+t127-t62+t128-t129-t130-t67+t131;
1351  t203 = -625.0/4096.0+t133+t134+t135-t136-t137-t75-t138-t139+t140-t79-t141-t142;
1352  t219 = 5625.0/32.0*t18;
1353  t220 = 6125.0/16.0*t45;
1354  t259 = 297675.0/64.0*t50;
1355  t261 = 2083725.0/256.0*t31;
1356  t318 = 4455.0/32.0*t18;
1357 
1358  values[0] = t30+t54;
1359  values[1] = t69+t82;
1360  values[2] = t97+t110;
1361  values[3] = t112+t113;
1362  values[4] = t115+t116;
1363  values[5] = t118+t119;
1364  values[6] = t132+t143;
1365  values[7] = t156+t167;
1366  values[8] = t169+t170;
1367  values[9] = t172+t173;
1368  values[10] = t175+t176;
1369  values[11] = t178+t179;
1370  values[12] = t181+t182;
1371  values[13] = t184+t185;
1372  values[14] = t187+t188;
1373  values[15] = t190+t191;
1374  values[16] = t193+t194;
1375  values[17] = t196+t197;
1376  values[18] = t199+t200;
1377  values[19] = t202+t203;
1378  values[20] = 45.0*t18+1225.0/16.0*t45-105.0*t16+2079.0/32.0*t12+2079.0/32.0*t28-105.0*t14;
1379  values[21] = 735.0/64.0*t6-315.0/64.0*xi+33075.0/64.0*t35-11025.0/32.0*t22-14175.0/64.0*t41+4725.0/32.0*t4;
1380  values[22] = t219+t220-525.0*t16+10395.0/32.0*t28-3675.0/16.0*t14;
1381  values[23] = -33075.0/32.0*t22+6615.0/64.0*t6+14175.0/32.0*t4-2835.0/64.0*xi+77175.0/64.0*t35-33075.0/64.0*t41;
1382  values[24] = 735.0/64.0*t8-315.0/64.0*eta+33075.0/64.0*t33-14175.0/64.0*t43-11025.0/32.0*t24+4725.0/32.0*t20;
1383  values[25] = 441.0/256.0+19845.0/256.0*t26-6615.0/128.0*t3+19845.0/256.0*t10-6615.0/128.0*t1+893025.0/256.0*t31-297675.0/128.0*t39-297675.0/128.0*t37+99225.0/64.0*t50;
1384  values[26] = -2205.0/64.0*eta+3675.0/64.0*t8-99225.0/64.0*t43+33075.0/32.0*t20+165375.0/64.0*t33-55125.0/32.0*t24;
1385  values[27] = 46305.0/256.0*t26-19845.0/128.0*t3+3969.0/256.0-893025.0/128.0*t39+178605.0/256.0*t10+t259-59535.0/128.0*t1+t261-694575.0/128.0*t37;
1386  values[28] = t219+t220-3675.0/16.0*t16+10395.0/32.0*t12-525.0*t14;
1387  values[29] = 3675.0/64.0*t6-2205.0/64.0*xi-99225.0/64.0*t41+33075.0/32.0*t4+165375.0/64.0*t35-55125.0/32.0*t22;
1388  values[30] = 11025.0/16.0*t18+30625.0/16.0*t45-18375.0/16.0*t16-18375.0/16.0*t14;
1389  values[31] = -231525.0/64.0*t41-165375.0/32.0*t22+33075.0/64.0*t6-19845.0/64.0*xi+99225.0/32.0*t4+385875.0/64.0*t35;
1390  values[32] = -33075.0/32.0*t24+14175.0/32.0*t20+6615.0/64.0*t8-2835.0/64.0*eta+77175.0/64.0*t33-33075.0/64.0*t43;
1391  values[33] = -893025.0/128.0*t37+t259-19845.0/128.0*t1+178605.0/256.0*t26-59535.0/128.0*t3+3969.0/256.0+46305.0/256.0*t10+t261-694575.0/128.0*t39;
1392  values[34] = 33075.0/64.0*t8-165375.0/32.0*t24-231525.0/64.0*t43+99225.0/32.0*t20-19845.0/64.0*eta+385875.0/64.0*t33;
1393  values[35] = 416745.0/256.0*t26+35721.0/256.0-2083725.0/128.0*t37+4862025.0/256.0*t31-178605.0/128.0*t3+416745.0/256.0*t10-2083725.0/128.0*t39-178605.0/128.0*t1+893025.0/64.0*t50;
1394  values[36] = -10395.0/16.0*t16+t318+18711.0/32.0*t28;
1395  values[37] = -10395.0/16.0*t14+t318+18711.0/32.0*t12;
1396 }
1397 
1398 // values of the derivatives in eta-eta direction
1399 static void C_Q_UL5_2D_DeriveEtaEta(double xi, double eta, double *values)
1400 {
1401  double t1, t2, t3, t4, t5, t6, t7, t8, t9, t10, t11, t12, t13, t14, t15;
1402  double t16, t17, t18, t19, t20, t21, t22, t23, t24, t25, t26, t27, t28;
1403  double t29, t30, t31, t32, t33, t34, t35, t36, t37, t38, t39, t40, t41;
1404  double t42, t43, t44, t45, t46, t47, t48, t49, t50, t51, t52, t54, t55;
1405  double t56, t57, t58, t59, t60, t61, t62, t63, t64, t65, t66, t67, t68;
1406  double t69, t70, t71, t72, t73, t74, t75, t76, t77, t78, t79, t80, t82;
1407  double t83, t84, t85, t86, t87, t88, t89, t90, t91, t92, t93, t94, t95;
1408  double t96, t97, t98, t99, t100, t101, t102, t103, t104, t105, t106, t107;
1409  double t108, t110, t111, t113, t114, t116, t117, t119, t120, t121, t122;
1410  double t123, t124, t125, t126, t127, t128, t129, t130, t131, t132, t133;
1411  double t134, t135, t136, t137, t138, t139, t140, t142, t143, t144, t145;
1412  double t146, t147, t148, t149, t150, t151, t152, t153, t154, t155, t156;
1413  double t157, t158, t159, t160, t161, t162, t163, t165, t166, t168, t169;
1414  double t171, t172, t174, t175, t177, t178, t180, t181, t183, t184, t186;
1415  double t187, t189, t190, t192, t193, t195, t196, t198, t199, t217, t218;
1416  double t256, t257, t313;
1417 
1418  t1 = 2175.0/512.0*eta;
1419  t2 = 1465.0/512.0*xi;
1420  t3 = xi*xi;
1421  t4 = t3*t3;
1422  t5 = 6255.0/1024.0*t4;
1423  t6 = t4*t3;
1424  t7 = 1463.0/1024.0*t6;
1425  t8 = t3*xi;
1426  t9 = 9475.0/768.0*t8;
1427  t10 = t4*xi;
1428  t11 = t10*eta;
1429  t12 = 38045.0/512.0*t11;
1430  t13 = t8*eta;
1431  t14 = 25025.0/256.0*t13;
1432  t15 = eta*eta;
1433  t16 = t15*eta;
1434  t17 = xi*t16;
1435  t18 = 29515.0/512.0*t17;
1436  t19 = xi*eta;
1437  t20 = 11505.0/512.0*t19;
1438  t21 = t3*t15;
1439  t22 = 18765.0/512.0*t21;
1440  t23 = xi*t15;
1441  t24 = 7755.0/512.0*t23;
1442  t25 = t3*eta;
1443  t26 = 9475.0/256.0*t25;
1444  t27 = 2497.0/1536.0-t1-t2+t5-t7+t9+t12-t14-t18+t20+t22+t24+t26;
1445  t28 = t8*t15;
1446  t29 = 16275.0/256.0*t28;
1447  t30 = t3*t16;
1448  t31 = 23425.0/256.0*t30;
1449  t32 = t4*t16;
1450  t33 = 111475.0/1536.0*t32;
1451  t34 = t10*t15;
1452  t35 = 22295.0/512.0*t34;
1453  t36 = t10*t16;
1454  t37 = 93135.0/512.0*t36;
1455  t38 = t15*t15;
1456  t39 = t3*t38;
1457  t40 = 21945.0/1024.0*t39;
1458  t41 = t8*t16;
1459  t42 = 190225.0/768.0*t41;
1460  t43 = t4*t15;
1461  t44 = 7525.0/512.0*t43;
1462  t45 = t4*eta;
1463  t46 = 16275.0/512.0*t45;
1464  t47 = 5525.0/512.0*t16;
1465  t48 = 4685.0/512.0*t10;
1466  t49 = 3395.0/512.0*t3;
1467  t50 = 2185.0/128.0*t15;
1468  t51 = 21945.0/1024.0*t38;
1469  t52 = -t29-t31+t33+t35-t37-t40+t42-t44-t46+t47-t48-t49-t50+t51;
1470  t54 = 2625.0/1024.0*eta;
1471  t55 = 625.0/1024.0*xi;
1472  t56 = 9375.0/1024.0*t4;
1473  t57 = 8125.0/512.0*t8;
1474  t58 = 109375.0/1024.0*t11;
1475  t59 = 56875.0/512.0*t13;
1476  t60 = 13125.0/1024.0*t17;
1477  t61 = 4375.0/1024.0*t19;
1478  t62 = 7875.0/64.0*t21;
1479  t63 = 4375.0/1024.0*t23;
1480  t64 = 34125.0/512.0*t25;
1481  t65 = 56875.0/512.0*t28;
1482  t66 = 102375.0/512.0*t30;
1483  t67 = -t54-t55+t56+t57+t58-t59-t60+t61+t62+t63+t64-t65-t66;
1484  t68 = 196875.0/1024.0*t32;
1485  t69 = 109375.0/1024.0*t34;
1486  t70 = 328125.0/1024.0*t36;
1487  t71 = 86625.0/1024.0*t39;
1488  t72 = 170625.0/512.0*t41;
1489  t73 = 65625.0/1024.0*t43;
1490  t74 = 65625.0/1024.0*t45;
1491  t75 = 7875.0/1024.0*t16;
1492  t76 = 15625.0/1024.0*t10;
1493  t77 = 13875.0/1024.0*t3;
1494  t78 = 60375.0/1024.0*t15;
1495  t79 = 86625.0/1024.0*t38;
1496  t80 = 1125.0/256.0+t68+t69-t70-t71+t72-t73-t74+t75-t76-t77-t78+t79;
1497  t82 = 7875.0/512.0*eta;
1498  t83 = 5625.0/512.0*xi;
1499  t84 = 3125.0/512.0*t4;
1500  t85 = 10625.0/256.0*t8;
1501  t86 = 109375.0/512.0*t11;
1502  t87 = 74375.0/256.0*t13;
1503  t88 = 118125.0/512.0*t17;
1504  t89 = 39375.0/512.0*t19;
1505  t90 = 2625.0/128.0*t21;
1506  t91 = 39375.0/512.0*t23;
1507  t92 = 14875.0/256.0*t25;
1508  t93 = 74375.0/256.0*t28;
1509  t94 = 125.0/256.0+t82+t83-t84-t85-t86+t87+t88-t89-t90-t91-t92+t93;
1510  t95 = 44625.0/256.0*t30;
1511  t96 = 65625.0/512.0*t32;
1512  t97 = 109375.0/512.0*t34;
1513  t98 = 328125.0/512.0*t36;
1514  t99 = 28875.0/512.0*t39;
1515  t100 = 223125.0/256.0*t41;
1516  t101 = 21875.0/512.0*t43;
1517  t102 = 21875.0/512.0*t45;
1518  t103 = 23625.0/512.0*t16;
1519  t104 = 15625.0/512.0*t10;
1520  t105 = 2875.0/512.0*t3;
1521  t106 = 11375.0/512.0*t15;
1522  t107 = 28875.0/512.0*t38;
1523  t108 = t95-t96-t97+t98-t99-t100+t101+t102-t103+t104+t105-t106+t107;
1524  t110 = 125.0/256.0+t82-t83-t84+t85+t86-t87-t88+t89-t90+t91-t92-t93;
1525  t111 = t95-t96+t97-t98-t99+t100+t101+t102-t103-t104+t105-t106+t107;
1526  t113 = -t54+t55+t56-t57-t58+t59+t60-t61+t62-t63+t64+t65-t66;
1527  t114 = 1125.0/256.0+t68-t69+t70-t71-t72-t73-t74+t75+t76-t77-t78+t79;
1528  t116 = 2497.0/1536.0-t1+t2+t5-t7-t9-t12+t14+t18-t20+t22-t24+t26;
1529  t117 = t29-t31+t33-t35+t37-t40-t42-t44-t46+t47+t48-t49-t50+t51;
1530  t119 = 4875.0/1024.0*eta;
1531  t120 = 4875.0/1024.0*xi;
1532  t121 = 2625.0/128.0*t4;
1533  t122 = 5775.0/1024.0*t6;
1534  t123 = 11375.0/512.0*t8;
1535  t124 = 102375.0/1024.0*t11;
1536  t125 = 78125.0/1024.0*t17;
1537  t126 = 24375.0/1024.0*t19;
1538  t127 = 28125.0/512.0*t21;
1539  t128 = 28125.0/1024.0*t23;
1540  t129 = 24375.0/512.0*t25;
1541  t130 = -t119+t120+t121-t122-t123-t124+t59+t125-t126+t127-t128+t129;
1542  t131 = 65625.0/512.0*t28;
1543  t132 = 78125.0/512.0*t30;
1544  t133 = 546875.0/3072.0*t32;
1545  t134 = 118125.0/1024.0*t34;
1546  t135 = 546875.0/1536.0*t41;
1547  t136 = 56875.0/1024.0*t45;
1548  t137 = 15625.0/1024.0*t16;
1549  t138 = 20475.0/1024.0*t10;
1550  t139 = 5625.0/1024.0*t15;
1551  t140 = 625.0/512.0+t131-t132+t133-t134+t70-t135-t73-t136+t137+t138-t77-t139;
1552  t142 = 6375.0/512.0*eta;
1553  t143 = 2125.0/512.0*xi;
1554  t144 = 875.0/256.0*t4;
1555  t145 = 1925.0/512.0*t6;
1556  t146 = 14875.0/768.0*t8;
1557  t147 = 133875.0/512.0*t11;
1558  t148 = 78125.0/512.0*t17;
1559  t149 = 31875.0/512.0*t19;
1560  t150 = 9375.0/256.0*t21;
1561  t151 = 9375.0/512.0*t23;
1562  t152 = -125.0/192.0+t142-t143-t144-t145+t146+t147-t87-t148+t149-t150+t151;
1563  t153 = 31875.0/256.0*t25;
1564  t154 = 21875.0/256.0*t28;
1565  t155 = 78125.0/256.0*t30;
1566  t156 = 546875.0/1536.0*t32;
1567  t157 = 39375.0/512.0*t34;
1568  t158 = 546875.0/768.0*t41;
1569  t159 = 74375.0/512.0*t45;
1570  t160 = 15625.0/512.0*t16;
1571  t161 = 8925.0/512.0*t10;
1572  t162 = 1875.0/512.0*t15;
1573  t163 = -t153-t154+t155-t156+t157-t98+t158+t101+t159-t160-t161+t105+t162;
1574  t165 = -125.0/192.0-t142-t143-t144-t145+t146-t147+t87+t148-t149-t150+t151;
1575  t166 = t153-t154-t155+t156+t157+t98-t158+t101-t159+t160-t161+t105+t162;
1576  t168 = t119+t120+t121-t122-t123+t124-t59-t125+t126+t127-t128-t129;
1577  t169 = 625.0/512.0+t131+t132-t133-t134-t70+t135-t73+t136-t137+t138-t77-t139;
1578  t171 = 2497.0/1536.0+t1+t2+t5-t7-t9+t12-t14-t18+t20+t22-t24-t26;
1579  t172 = t29+t31-t33-t35-t37-t40+t42-t44+t46-t47+t48-t49-t50+t51;
1580  t174 = t54+t55+t56-t57+t58-t59-t60+t61+t62-t63-t64+t65+t66;
1581  t175 = 1125.0/256.0-t68-t69-t70-t71+t72-t73+t74-t75+t76-t77-t78+t79;
1582  t177 = 125.0/256.0-t82-t83-t84+t85-t86+t87+t88-t89-t90+t91+t92-t93;
1583  t178 = -t95+t96+t97+t98-t99-t100+t101-t102+t103-t104+t105-t106+t107;
1584  t180 = 125.0/256.0-t82+t83-t84-t85+t86-t87-t88+t89-t90-t91+t92+t93;
1585  t181 = -t95+t96-t97-t98-t99+t100+t101-t102+t103+t104+t105-t106+t107;
1586  t183 = t54-t55+t56+t57-t58+t59+t60-t61+t62+t63-t64-t65+t66;
1587  t184 = 1125.0/256.0-t68+t69+t70-t71-t72-t73+t74-t75-t76-t77-t78+t79;
1588  t186 = 2497.0/1536.0+t1-t2+t5-t7+t9-t12+t14+t18-t20+t22+t24-t26;
1589  t187 = -t29+t31-t33+t35+t37-t40-t42-t44+t46-t47-t48-t49-t50+t51;
1590  t189 = t119-t120+t121-t122+t123-t124+t59+t125-t126+t127+t128-t129;
1591  t190 = 625.0/512.0-t131+t132-t133+t134+t70-t135-t73+t136-t137-t138-t77-t139;
1592  t192 = -125.0/192.0-t142+t143-t144-t145-t146+t147-t87-t148+t149-t150-t151;
1593  t193 = t153+t154-t155+t156-t157-t98+t158+t101-t159+t160+t161+t105+t162;
1594  t195 = -125.0/192.0+t142+t143-t144-t145-t146-t147+t87+t148-t149-t150-t151;
1595  t196 = -t153+t154+t155-t156-t157+t98-t158+t101+t159-t160+t161+t105+t162;
1596  t198 = -t119-t120+t121-t122+t123+t124-t59-t125+t126+t127+t128+t129;
1597  t199 = 625.0/512.0-t131-t132+t133+t134-t70+t135-t73-t136+t137-t138-t77-t139;
1598  t217 = 5625.0/64.0*t3;
1599  t218 = 18375.0/64.0*t43;
1600  t256 = 99225.0/32.0*t13;
1601  t257 = 416745.0/64.0*t36;
1602  t313 = 4455.0/64.0*t3;
1603 
1604  values[0] = t27+t52;
1605  values[1] = t67+t80;
1606  values[2] = t94+t108;
1607  values[3] = t110+t111;
1608  values[4] = t113+t114;
1609  values[5] = t116+t117;
1610  values[6] = t130+t140;
1611  values[7] = t152+t163;
1612  values[8] = t165+t166;
1613  values[9] = t168+t169;
1614  values[10] = t171+t172;
1615  values[11] = t174+t175;
1616  values[12] = t177+t178;
1617  values[13] = t180+t181;
1618  values[14] = t183+t184;
1619  values[15] = t186+t187;
1620  values[16] = t189+t190;
1621  values[17] = t192+t193;
1622  values[18] = t195+t196;
1623  values[19] = t198+t199;
1624  values[20] = 45.0/2.0*t3+3675.0/64.0*t43-453.0/64.0-315.0/2.0*t21+693.0/64.0*t6+6405.0/64.0*t15-10395.0/64.0*t38+10395.0/64.0*t39-105.0/4.0*t4;
1625  values[21] = -4725.0/64.0*t16+1575.0/64.0*eta+33075.0/64.0*t32-11025.0/64.0*t45-14175.0/32.0*t30+4725.0/32.0*t25;
1626  values[22] = t217+t218-975.0/32.0-1575.0/2.0*t21+32025.0/64.0*t15-51975.0/64.0*t38+51975.0/64.0*t39-3675.0/64.0*t4;
1627  values[23] = -11025.0/64.0*t16+4725.0/64.0*eta-33075.0/64.0*t45+14175.0/32.0*t25+77175.0/64.0*t32-33075.0/32.0*t30;
1628  values[24] = 2205.0/64.0*t23-315.0/64.0*xi+19845.0/64.0*t34-2835.0/64.0*t10-11025.0/32.0*t28+1575.0/32.0*t8;
1629  values[25] = 19845.0/64.0*t17-6615.0/64.0*t19+178605.0/64.0*t36-59535.0/64.0*t11-99225.0/32.0*t41+33075.0/32.0*t13;
1630  values[26] = -2205.0/64.0*xi+11025.0/64.0*t23-19845.0/64.0*t10+11025.0/32.0*t8+99225.0/64.0*t34-55125.0/32.0*t28;
1631  values[27] = 46305.0/64.0*t17-19845.0/64.0*t19-178605.0/64.0*t11+t256+t257-231525.0/32.0*t41;
1632  values[28] = t217+t218-345.0/32.0-11025.0/32.0*t21+3465.0/64.0*t6+3675.0/64.0*t15-525.0/4.0*t4;
1633  values[29] = 33075.0/64.0*t16-11025.0/64.0*eta-99225.0/32.0*t30+33075.0/32.0*t25+165375.0/64.0*t32-55125.0/64.0*t45;
1634  values[30] = 11025.0/32.0*t3+91875.0/64.0*t43-3675.0/64.0-55125.0/32.0*t21+18375.0/64.0*t15-18375.0/64.0*t4;
1635  values[31] = -231525.0/32.0*t30-165375.0/64.0*t45+77175.0/64.0*t16-33075.0/64.0*eta+99225.0/32.0*t25+385875.0/64.0*t32;
1636  values[32] = -33075.0/32.0*t28+4725.0/32.0*t8+19845.0/64.0*t23-2835.0/64.0*xi+46305.0/64.0*t34-6615.0/64.0*t10;
1637  values[33] = -297675.0/32.0*t41+t256+178605.0/64.0*t17-59535.0/64.0*t19+t257-138915.0/64.0*t11;
1638  values[34] = 99225.0/64.0*t23-165375.0/32.0*t28-46305.0/64.0*t10+33075.0/32.0*t8-19845.0/64.0*xi+231525.0/64.0*t34;
1639  values[35] = 416745.0/64.0*t17-694575.0/32.0*t41+972405.0/64.0*t36-178605.0/64.0*t19-416745.0/64.0*t11+297675.0/32.0*t13;
1640  values[36] = 31185.0/32.0*t15-4455.0/64.0-31185.0/32.0*t21+t313-93555.0/64.0*t38+93555.0/64.0*t39;
1641  values[37] = -297.0/64.0-10395.0/64.0*t4+t313+6237.0/64.0*t6;
1642 }
1643 
1644 
1645 // ***********************************************************************
1646 
1647 TBaseFunct2D *BF_C_Q_UL5_2D_Obj = new TBaseFunct2D
1648  (38, BF_C_Q_UL5_2D, BFUnitSquare,
1649  C_Q_UL5_2D_Funct, C_Q_UL5_2D_DeriveXi,
1650  C_Q_UL5_2D_DeriveEta, C_Q_UL5_2D_DeriveXiXi,
1651  C_Q_UL5_2D_DeriveXiEta, C_Q_UL5_2D_DeriveEtaEta, 6, 5,
1652  0, NULL);
TBaseFunct2D
Definition: BaseFunct2D.h:27
Generated on Tue May 10 2016 09:07:58 for ParMooN by   doxygen 1.8.5