ParMooN
 All Classes Functions Variables Friends Pages
BF_C_Q_UL6SE_2D.h
1 // ***********************************************************************
2 // UL6SE element, conforming, 2D
3 // ***********************************************************************
4 
5 // base function values
6 static void C_Q_UL6SE_2D_Funct(double xi, double eta, double *values)
7 {
8  double t2 = xi*xi;
9  double t3 = 1.0-t2;
10  double t5 = t2*t2;
11  double t8 = 3.0/8.0+35.0/8.0*t5-15.0/4.0*t2;
12  double t9 = (1.0-eta)*t3*t8;
13  double t10 = 81.0/700.0*t9;
14  double t11 = 421161.0/11468800.0*eta;
15  double t12 = t5*xi;
16  double t13 = 4651.0/327680.0*t12;
17  double t14 = eta*eta;
18  double t15 = 1564309.0/5734400.0*t14;
19  double t16 = t14*eta;
20  double t17 = 6163.0/819200.0*t16;
21  double t18 = t14*t14;
22  double t19 = 111833.0/458752.0*t18;
23  double t20 = t18*eta;
24  double t21 = 4651.0/327680.0*t20;
25  double t23 = 1.0-t14;
26  double t27 = 3.0/8.0+35.0/8.0*t18-15.0/4.0*t14;
27  double t28 = (1.0-xi)*t23*t27;
28  double t29 = 81.0/700.0*t28;
29  double t30 = 1564309.0/5734400.0*t2;
30  double t31 = 111833.0/458752.0*t5;
31  double t32 = t2*xi;
32  double t33 = 6163.0/819200.0*t32;
33  double t34 = 421161.0/11468800.0*xi;
34  double t35 = t2*t16;
35  double t36 = 727871.0/3686400.0*t35;
36  double t37 = xi*eta;
37  double t38 = 339.0/327680.0*t37;
38  double t39 = t2*eta;
39  double t40 = 4891279.0/17203200.0*t39;
40  double t41 = t32*eta;
41  double t42 = 1873.0/98304.0*t41;
42  double t43 = t5*eta;
43  double t44 = 1875847.0/6881280.0*t43;
44  double t45 = t12*eta;
45  double t46 = 42289.0/983040.0*t45;
46  double t47 = xi*t14;
47  double t48 = 4891279.0/17203200.0*t47;
48  double t49 = -t10-t11-t13-t15+t17+t19-t21-t29-t30+t31+t33-t34-t36+t38+t40-t42-t44+t46+t48;
49  double t50 = t2*t14;
50  double t51 = 496367.0/3686400.0*t50;
51  double t52 = t32*t14;
52  double t53 = 727871.0/3686400.0*t52;
53  double t54 = t5*t14;
54  double t55 = 463271.0/1474560.0*t54;
55  double t56 = t12*t14;
56  double t57 = 538871.0/1474560.0*t56;
57  double t58 = xi*t16;
58  double t59 = 1873.0/98304.0*t58;
59  double t60 = t32*t16;
60  double t61 = 38375.0/147456.0*t60;
61  double t62 = t5*t16;
62  double t63 = 138955.0/294912.0*t62;
63  double t64 = t12*t16;
64  double t65 = 154075.0/294912.0*t64;
65  double t66 = xi*t18;
66  double t67 = 1875847.0/6881280.0*t66;
67  double t68 = t2*t18;
68  double t69 = 463271.0/1474560.0*t68;
69  double t70 = t32*t18;
70  double t71 = 138955.0/294912.0*t70;
71  double t72 = t5*t18;
72  double t73 = 430175.0/589824.0*t72;
73  double t74 = t12*t18;
74  double t75 = 505775.0/589824.0*t74;
75  double t76 = xi*t20;
76  double t77 = 42289.0/983040.0*t76;
77  double t78 = t2*t20;
78  double t79 = 538871.0/1474560.0*t78;
79  double t80 = t32*t20;
80  double t81 = 154075.0/294912.0*t80;
81  double t82 = t5*t20;
82  double t83 = 505775.0/589824.0*t82;
83  double t84 = t12*t20;
84  double t85 = 581375.0/589824.0*t84;
85  double t86 = 830457.0/11468800.0+t51-t53-t55+t57-t59+t61+t63-t65-t67-t69+t71+t73-t75+t77+t79-t81-t83+t85;
86  double t88 = 243.0/350.0*t9;
87  double t89 = 78003.0/358400.0*eta;
88  double t90 = 243.0/5120.0*t12;
89  double t91 = 729.0/25600.0*t14;
90  double t92 = 729.0/25600.0*t16;
91  double t93 = 729.0/10240.0*t18;
92  double t94 = 729.0/10240.0*t20;
93  double t95 = 291519.0/179200.0*t2;
94  double t96 = 101007.0/71680.0*t5;
95  double t97 = 27.0/512.0*t32;
96  double t98 = 27.0/5120.0*xi;
97  double t99 = 7677.0/12800.0*t35;
98  double t100 = 27.0/5120.0*t37;
99  double t101 = 291519.0/179200.0*t39;
100  double t102 = 27.0/512.0*t41;
101  double t103 = 101007.0/71680.0*t43;
102  double t104 = 243.0/5120.0*t45;
103  double t105 = 75.0/512.0*t47;
104  double t106 = t88+t89+t90+t91-t92-t93+t94+t95-t96-t97+t98+t99-t100-t101+t102+t103-t104-t105;
105  double t107 = 7677.0/12800.0*t50;
106  double t108 = 375.0/256.0*t52;
107  double t109 = 585.0/1024.0*t54;
108  double t110 = 675.0/512.0*t56;
109  double t111 = 75.0/512.0*t58;
110  double t112 = 375.0/256.0*t60;
111  double t113 = 585.0/1024.0*t62;
112  double t114 = 675.0/512.0*t64;
113  double t115 = 375.0/1024.0*t66;
114  double t116 = 7677.0/5120.0*t68;
115  double t117 = 1875.0/512.0*t70;
116  double t118 = 2925.0/2048.0*t72;
117  double t119 = 3375.0/1024.0*t74;
118  double t120 = 375.0/1024.0*t76;
119  double t121 = 7677.0/5120.0*t78;
120  double t122 = 1875.0/512.0*t80;
121  double t123 = 2925.0/2048.0*t82;
122  double t124 = 3375.0/1024.0*t84;
123  double t125 = -78003.0/358400.0-t107+t108+t109-t110+t111-t112-t113+t114+t115+t116-t117-t118+t119-t120-t121+t122+t123-t124;
124  double t127 = 243.0/140.0*t9;
125  double t128 = 78003.0/143360.0*eta;
126  double t129 = 243.0/4096.0*t12;
127  double t130 = 729.0/10240.0*t14;
128  double t131 = 729.0/10240.0*t16;
129  double t132 = 729.0/4096.0*t18;
130  double t133 = 729.0/4096.0*t20;
131  double t134 = 574533.0/143360.0*t2;
132  double t135 = 49653.0/14336.0*t5;
133  double t136 = 351.0/4096.0*t32;
134  double t137 = 27.0/1024.0*xi;
135  double t138 = 1521.0/10240.0*t35;
136  double t139 = 27.0/1024.0*t37;
137  double t140 = 574533.0/143360.0*t39;
138  double t141 = 351.0/4096.0*t41;
139  double t142 = 49653.0/14336.0*t43;
140  double t143 = 243.0/4096.0*t45;
141  double t144 = 375.0/512.0*t47;
142  double t145 = -t127-t128-t129-t130+t131+t132-t133-t134+t135+t136-t137+t138+t139+t140-t141-t142+t143+t144;
143  double t146 = 1521.0/10240.0*t50;
144  double t147 = 4875.0/2048.0*t52;
145  double t148 = 225.0/1024.0*t54;
146  double t149 = 3375.0/2048.0*t56;
147  double t150 = 375.0/512.0*t58;
148  double t151 = 4875.0/2048.0*t60;
149  double t152 = 225.0/1024.0*t62;
150  double t153 = 3375.0/2048.0*t64;
151  double t154 = 1875.0/1024.0*t66;
152  double t155 = 1521.0/4096.0*t68;
153  double t156 = 24375.0/4096.0*t70;
154  double t157 = 1125.0/2048.0*t72;
155  double t158 = 16875.0/4096.0*t74;
156  double t159 = 1875.0/1024.0*t76;
157  double t160 = 1521.0/4096.0*t78;
158  double t161 = 24375.0/4096.0*t80;
159  double t162 = 1125.0/2048.0*t82;
160  double t163 = 16875.0/4096.0*t84;
161  double t164 = 78003.0/143360.0-t146-t147+t148+t149-t150+t151-t152-t153-t154+t155+t156-t157-t158+t159-t160-t161+t162+t163;
162  double t166 = 25581.0/35840.0*eta;
163  double t167 = 95073.0/17920.0*t39;
164  double t168 = 32913.0/7168.0*t43;
165  double t169 = 4823.0/3840.0*t50;
166  double t170 = 525.0/512.0*t54;
167  double t171 = 4823.0/3840.0*t35;
168  double t172 = 525.0/512.0*t62;
169  double t173 = 95073.0/17920.0*t2;
170  double t174 = 32913.0/7168.0*t5;
171  double t175 = 1771.0/7680.0*t14;
172  double t176 = 1771.0/7680.0*t16;
173  double t177 = 1771.0/3072.0*t18;
174  double t178 = 1771.0/3072.0*t20;
175  double t179 = 4823.0/1536.0*t68;
176  double t180 = 2625.0/1024.0*t72;
177  double t181 = 4823.0/1536.0*t78;
178  double t182 = 2625.0/1024.0*t82;
179  double t184 = t166-t167+t168+t169-t170-t171+t172+t173-t174-t175+t176+t177-t178-t179+t180+t181-t182+81.0/35.0*t9-25581.0/35840.0;
180  double t185 = -t127-t128+t129-t130+t131+t132-t133-t134+t135-t136+t137+t138-t139+t140+t141-t142-t143-t144;
181  double t186 = 78003.0/143360.0-t146+t147+t148-t149+t150-t151-t152+t153+t154+t155-t156-t157+t158-t159-t160+t161+t162-t163;
182  double t188 = t88+t89-t90+t91-t92-t93+t94+t95-t96+t97-t98+t99+t100-t101-t102+t103+t104+t105;
183  double t189 = -78003.0/358400.0-t107-t108+t109+t110-t111+t112-t113-t114-t115+t116+t117-t118-t119+t120-t121-t122+t123+t124;
184  double t193 = (1.0+xi)*t23*t27;
185  double t194 = 81.0/700.0*t193;
186  double t195 = -t194-t10-t11+t13-t15+t17+t19-t21-t30+t31-t33+t34-t36-t38+t40+t42-t44-t46-t48;
187  double t196 = 830457.0/11468800.0+t51+t53-t55-t57+t59-t61+t63+t65+t67-t69-t71+t73+t75-t77+t79+t81-t83-t85;
188  double t198 = 243.0/350.0*t193;
189  double t199 = 27.0/5120.0*eta;
190  double t200 = 729.0/10240.0*t12;
191  double t201 = 291519.0/179200.0*t14;
192  double t202 = 27.0/512.0*t16;
193  double t203 = 101007.0/71680.0*t18;
194  double t204 = 243.0/5120.0*t20;
195  double t205 = 729.0/25600.0*t2;
196  double t206 = 729.0/10240.0*t5;
197  double t207 = 729.0/25600.0*t32;
198  double t208 = 78003.0/358400.0*xi;
199  double t209 = 375.0/256.0*t35;
200  double t210 = 75.0/512.0*t39;
201  double t211 = 75.0/512.0*t41;
202  double t212 = 375.0/1024.0*t43;
203  double t213 = 375.0/1024.0*t45;
204  double t214 = 291519.0/179200.0*t47;
205  double t215 = t198+t199-t200+t201-t202-t203+t204+t205-t206+t207-t208+t209+t100-t210-t211+t212+t213+t214;
206  double t216 = 7677.0/12800.0*t52;
207  double t217 = 7677.0/5120.0*t54;
208  double t218 = 7677.0/5120.0*t56;
209  double t219 = 27.0/512.0*t58;
210  double t220 = 1875.0/512.0*t62;
211  double t221 = 1875.0/512.0*t64;
212  double t222 = 101007.0/71680.0*t66;
213  double t223 = 585.0/1024.0*t68;
214  double t224 = 585.0/1024.0*t70;
215  double t225 = 2925.0/2048.0*t74;
216  double t226 = 243.0/5120.0*t76;
217  double t227 = 675.0/512.0*t78;
218  double t228 = 675.0/512.0*t80;
219  double t229 = 3375.0/1024.0*t82;
220  double t230 = -78003.0/358400.0-t107-t216+t217+t218-t219+t112-t220-t221-t222+t223+t224-t118-t225+t226-t227-t228+t229+t124;
221  double t232 = 243.0/140.0*t193;
222  double t233 = 27.0/1024.0*eta;
223  double t234 = 729.0/4096.0*t12;
224  double t235 = 574533.0/143360.0*t14;
225  double t236 = 351.0/4096.0*t16;
226  double t237 = 49653.0/14336.0*t18;
227  double t238 = 243.0/4096.0*t20;
228  double t239 = 729.0/10240.0*t2;
229  double t240 = 729.0/4096.0*t5;
230  double t241 = 729.0/10240.0*t32;
231  double t242 = 78003.0/143360.0*xi;
232  double t243 = 4875.0/2048.0*t35;
233  double t244 = 375.0/512.0*t39;
234  double t245 = 375.0/512.0*t41;
235  double t246 = 1875.0/1024.0*t43;
236  double t247 = 1875.0/1024.0*t45;
237  double t248 = 574533.0/143360.0*t47;
238  double t249 = -t232-t233+t234-t235+t236+t237-t238-t239+t240-t241+t242-t243-t139+t244+t245-t246-t247-t248;
239  double t250 = 1521.0/10240.0*t52;
240  double t251 = 1521.0/4096.0*t54;
241  double t252 = 1521.0/4096.0*t56;
242  double t253 = 351.0/4096.0*t58;
243  double t254 = 24375.0/4096.0*t62;
244  double t255 = 24375.0/4096.0*t64;
245  double t256 = 49653.0/14336.0*t66;
246  double t257 = 225.0/1024.0*t68;
247  double t258 = 225.0/1024.0*t70;
248  double t259 = 1125.0/2048.0*t74;
249  double t260 = 243.0/4096.0*t76;
250  double t261 = 3375.0/2048.0*t78;
251  double t262 = 3375.0/2048.0*t80;
252  double t263 = 16875.0/4096.0*t82;
253  double t264 = 78003.0/143360.0-t146-t250+t251+t252+t253-t151+t254+t255+t256+t257+t258-t157-t259-t260+t261+t262-t263-t163;
254  double t266 = 25581.0/35840.0*xi;
255  double t267 = 95073.0/17920.0*t47;
256  double t268 = 4823.0/3840.0*t52;
257  double t269 = 4823.0/1536.0*t54;
258  double t270 = 4823.0/1536.0*t56;
259  double t271 = 1771.0/7680.0*t2;
260  double t272 = 1771.0/7680.0*t32;
261  double t273 = 1771.0/3072.0*t5;
262  double t274 = 1771.0/3072.0*t12;
263  double t275 = 95073.0/17920.0*t14;
264  double t276 = 32913.0/7168.0*t18;
265  double t277 = 32913.0/7168.0*t66;
266  double t278 = 525.0/512.0*t68;
267  double t279 = 525.0/512.0*t70;
268  double t280 = 2625.0/1024.0*t74;
269  double t282 = -t266+t267+t169+t268-t269-t270-t271-t272+t273+t274+t275-t276-t277-t278-t279+t180+t280+81.0/35.0*t193-25581.0/35840.0;
270  double t283 = -t232+t233+t234-t235-t236+t237+t238-t239+t240-t241+t242+t243+t139-t244-t245+t246+t247-t248;
271  double t284 = 78003.0/143360.0-t146-t250+t251+t252-t253+t151-t254-t255+t256+t257+t258-t157-t259+t260-t261-t262+t263+t163;
272  double t286 = t198-t199-t200+t201+t202-t203-t204+t205-t206+t207-t208-t209-t100+t210+t211-t212-t213+t214;
273  double t287 = -78003.0/358400.0-t107-t216+t217+t218+t219-t112+t220+t221-t222+t223+t224-t118-t225-t226+t227+t228-t229-t124;
274  double t291 = (1.0+eta)*t3*t8;
275  double t292 = 81.0/700.0*t291;
276  double t293 = -t194+t11-t292+t13-t15-t17+t19+t21-t30+t31-t33+t34+t36+t38-t40-t42+t44+t46-t48;
277  double t294 = 830457.0/11468800.0+t51+t53-t55-t57-t59+t61-t63-t65+t67-t69-t71+t73+t75+t77-t79-t81+t83+t85;
278  double t296 = 243.0/350.0*t291;
279  double t297 = -t89+t296-t90+t91+t92-t93-t94+t95-t96+t97-t98-t99-t100+t101+t102-t103-t104+t105;
280  double t298 = -78003.0/358400.0-t107-t108+t109+t110+t111-t112+t113+t114-t115+t116+t117-t118-t119-t120+t121+t122-t123-t124;
281  double t300 = 243.0/140.0*t291;
282  double t301 = t128-t300+t129-t130-t131+t132+t133-t134+t135-t136+t137-t138+t139-t140-t141+t142+t143-t144;
283  double t302 = 78003.0/143360.0-t146+t147+t148-t149-t150+t151+t152-t153+t154+t155-t156-t157+t158+t159+t160-t161-t162+t163;
284  double t305 = -t166+t167-t168+t169-t170+t171-t172+t173-t174-t175-t176+t177+t178-t179+t180-t181+t182+81.0/35.0*t291-25581.0/35840.0;
285  double t306 = t128-t300-t129-t130-t131+t132+t133-t134+t135+t136-t137-t138-t139-t140+t141+t142-t143+t144;
286  double t307 = 78003.0/143360.0-t146-t147+t148+t149+t150-t151+t152+t153-t154+t155+t156-t157-t158-t159+t160+t161-t162-t163;
287  double t309 = -t89+t296+t90+t91+t92-t93-t94+t95-t96-t97+t98-t99+t100+t101-t102-t103+t104-t105;
288  double t310 = -78003.0/358400.0-t107+t108+t109-t110-t111+t112+t113-t114+t115+t116-t117-t118+t119+t120+t121-t122-t123+t124;
289  double t312 = t11-t292-t13-t15-t17+t19+t21-t29-t30+t31+t33-t34+t36-t38-t40+t42+t44-t46+t48;
290  double t313 = 830457.0/11468800.0+t51-t53-t55+t57+t59-t61-t63+t65-t67-t69+t71+t73-t75-t77-t79+t81+t83-t85;
291  double t315 = 243.0/350.0*t28;
292  double t316 = -t199+t200+t201+t202-t203-t204+t315+t205-t206-t207+t208-t209+t100+t210-t211-t212+t213-t214;
293  double t317 = -78003.0/358400.0-t107+t216+t217-t218-t219+t112+t220-t221+t222+t223-t224-t118+t225+t226+t227-t228-t229+t124;
294  double t319 = 243.0/140.0*t28;
295  double t320 = t233-t234-t235-t236+t237+t238-t319-t239+t240+t241-t242+t243-t139-t244+t245+t246-t247+t248;
296  double t321 = 78003.0/143360.0-t146+t250+t251-t252+t253-t151-t254+t255-t256+t257-t258-t157+t259-t260-t261+t262+t263-t163;
297  double t324 = t266-t267+t169-t268-t269+t270-t271+t272+t273-t274+t275-t276+t277-t278+t279+t180-t280+81.0/35.0*t28-25581.0/35840.0;
298  double t325 = -t233-t234-t235+t236+t237-t238-t319-t239+t240+t241-t242-t243+t139+t244-t245-t246+t247+t248;
299  double t326 = 78003.0/143360.0-t146+t250+t251-t252-t253+t151+t254-t255-t256+t257-t258-t157+t259+t260+t261-t262-t263+t163;
300  double t328 = t199+t200+t201-t202-t203+t204+t315+t205-t206-t207+t208+t209-t100-t210+t211+t212-t213-t214;
301  double t329 = -78003.0/358400.0-t107+t216+t217-t218+t219-t112-t220+t221+t222+t223-t224-t118+t225-t226-t227+t228+t229-t124;
302  double t331 = 3125.0/196608.0*eta;
303  double t332 = 78125.0/196608.0*t12;
304  double t333 = 8125.0/32768.0*t14;
305  double t334 = 40625.0/98304.0*t16;
306  double t335 = 15625.0/65536.0*t18;
307  double t336 = 78125.0/196608.0*t20;
308  double t337 = 8125.0/32768.0*t2;
309  double t338 = 15625.0/65536.0*t5;
310  double t339 = 40625.0/98304.0*t32;
311  double t340 = 3125.0/196608.0*xi;
312  double t341 = 528125.0/49152.0*t35;
313  double t342 = 15625.0/589824.0*t37;
314  double t343 = 40625.0/98304.0*t39;
315  double t344 = 203125.0/294912.0*t41;
316  double t345 = 78125.0/196608.0*t43;
317  double t346 = 390625.0/589824.0*t45;
318  double t347 = 40625.0/98304.0*t47;
319  double t348 = 625.0/65536.0-t331-t332-t333+t334+t335-t336-t337+t338+t339-t340-t341+t342+t343-t344-t345+t346+t347;
320  double t349 = 105625.0/16384.0*t50;
321  double t350 = 528125.0/49152.0*t52;
322  double t351 = 203125.0/32768.0*t54;
323  double t352 = 1015625.0/98304.0*t56;
324  double t353 = 203125.0/294912.0*t58;
325  double t354 = 2640625.0/147456.0*t60;
326  double t355 = 1015625.0/98304.0*t62;
327  double t356 = 5078125.0/294912.0*t64;
328  double t357 = 78125.0/196608.0*t66;
329  double t358 = 203125.0/32768.0*t68;
330  double t359 = 1015625.0/98304.0*t70;
331  double t360 = 390625.0/65536.0*t72;
332  double t361 = 1953125.0/196608.0*t74;
333  double t362 = 390625.0/589824.0*t76;
334  double t363 = 1015625.0/98304.0*t78;
335  double t364 = 5078125.0/294912.0*t80;
336  double t365 = 1953125.0/196608.0*t82;
337  double t366 = 9765625.0/589824.0*t84;
338  double t367 = t349-t350-t351+t352-t353+t354+t355-t356-t357-t358+t359+t360-t361+t362+t363-t364-t365+t366;
339  double t369 = 9375.0/32768.0*eta;
340  double t370 = 78125.0/32768.0*t12;
341  double t371 = 10625.0/49152.0*t14;
342  double t372 = 53125.0/49152.0*t16;
343  double t373 = 15625.0/98304.0*t18;
344  double t374 = 78125.0/98304.0*t20;
345  double t375 = 24375.0/16384.0*t2;
346  double t376 = 46875.0/32768.0*t5;
347  double t377 = 40625.0/16384.0*t32;
348  double t378 = 3125.0/32768.0*xi;
349  double t379 = 690625.0/24576.0*t35;
350  double t380 = 15625.0/32768.0*t37;
351  double t381 = 121875.0/16384.0*t39;
352  double t382 = 203125.0/16384.0*t41;
353  double t383 = 234375.0/32768.0*t43;
354  double t384 = 390625.0/32768.0*t45;
355  double t385 = 53125.0/147456.0*t47;
356  double t386 = -1875.0/32768.0+t369+t370+t371-t372-t373+t374+t375-t376-t377+t378+t379-t380-t381+t382+t383-t384-t385;
357  double t387 = 138125.0/24576.0*t50;
358  double t388 = 690625.0/73728.0*t52;
359  double t389 = 265625.0/49152.0*t54;
360  double t390 = 1328125.0/147456.0*t56;
361  double t391 = 265625.0/147456.0*t58;
362  double t392 = 3453125.0/73728.0*t60;
363  double t393 = 1328125.0/49152.0*t62;
364  double t394 = 6640625.0/147456.0*t64;
365  double t395 = 78125.0/294912.0*t66;
366  double t396 = 203125.0/49152.0*t68;
367  double t397 = 1015625.0/147456.0*t70;
368  double t398 = 390625.0/98304.0*t72;
369  double t399 = 1953125.0/294912.0*t74;
370  double t400 = 390625.0/294912.0*t76;
371  double t401 = 1015625.0/49152.0*t78;
372  double t402 = 5078125.0/147456.0*t80;
373  double t403 = 1953125.0/98304.0*t82;
374  double t404 = 9765625.0/294912.0*t84;
375  double t405 = -t387+t388+t389-t390+t391-t392-t393+t394+t395+t396-t397-t398+t399-t400-t401+t402+t403-t404;
376  double t407 = -1875.0/32768.0-t369+t370+t371+t372-t373-t374+t375-t376-t377+t378-t379+t380+t381-t382-t383+t384-t385;
377  double t408 = -t387+t388+t389-t390-t391+t392+t393-t394+t395+t396-t397-t398+t399+t400+t401-t402-t403+t404;
378  double t410 = 625.0/65536.0+t331-t332-t333-t334+t335+t336-t337+t338+t339-t340+t341-t342-t343+t344+t345-t346+t347;
379  double t411 = t349-t350-t351+t352+t353-t354-t355+t356-t357-t358+t359+t360-t361-t362-t363+t364+t365-t366;
380  double t413 = 3125.0/32768.0*eta;
381  double t414 = 78125.0/98304.0*t12;
382  double t415 = 24375.0/16384.0*t14;
383  double t416 = 40625.0/16384.0*t16;
384  double t417 = 46875.0/32768.0*t18;
385  double t418 = 78125.0/32768.0*t20;
386  double t419 = 10625.0/49152.0*t2;
387  double t420 = 15625.0/98304.0*t5;
388  double t421 = 53125.0/49152.0*t32;
389  double t422 = 9375.0/32768.0*xi;
390  double t423 = 690625.0/73728.0*t35;
391  double t424 = 53125.0/147456.0*t39;
392  double t425 = 265625.0/147456.0*t41;
393  double t426 = 78125.0/294912.0*t43;
394  double t427 = 390625.0/294912.0*t45;
395  double t428 = 121875.0/16384.0*t47;
396  double t429 = -1875.0/32768.0+t413+t414+t415-t416-t417+t418+t419-t420-t421+t422+t423-t380-t424+t425+t426-t427-t428;
397  double t430 = 690625.0/24576.0*t52;
398  double t431 = 203125.0/49152.0*t54;
399  double t432 = 1015625.0/49152.0*t56;
400  double t433 = 203125.0/16384.0*t58;
401  double t434 = 1015625.0/147456.0*t62;
402  double t435 = 5078125.0/147456.0*t64;
403  double t436 = 234375.0/32768.0*t66;
404  double t437 = 265625.0/49152.0*t68;
405  double t438 = 1328125.0/49152.0*t70;
406  double t439 = 1953125.0/98304.0*t74;
407  double t440 = 390625.0/32768.0*t76;
408  double t441 = 1328125.0/147456.0*t78;
409  double t442 = 6640625.0/147456.0*t80;
410  double t443 = 1953125.0/294912.0*t82;
411  double t444 = -t387+t430+t431-t432+t433-t392-t434+t435+t436+t437-t438-t398+t439-t440-t441+t442+t443-t404;
412  double t446 = 28125.0/16384.0*eta;
413  double t447 = 78125.0/16384.0*t12;
414  double t448 = 10625.0/8192.0*t14;
415  double t449 = 53125.0/8192.0*t16;
416  double t450 = 15625.0/16384.0*t18;
417  double t451 = 78125.0/16384.0*t20;
418  double t452 = 10625.0/8192.0*t2;
419  double t453 = 15625.0/16384.0*t5;
420  double t454 = 53125.0/8192.0*t32;
421  double t455 = 28125.0/16384.0*xi;
422  double t456 = 903125.0/36864.0*t35;
423  double t457 = 140625.0/16384.0*t37;
424  double t458 = 53125.0/8192.0*t39;
425  double t459 = 265625.0/8192.0*t41;
426  double t460 = 78125.0/16384.0*t43;
427  double t461 = 390625.0/16384.0*t45;
428  double t462 = 53125.0/8192.0*t47;
429  double t463 = 5625.0/16384.0-t446-t447-t448+t449+t450-t451-t452+t453+t454-t455-t456+t457+t458-t459-t460+t461+t462;
430  double t464 = 180625.0/36864.0*t50;
431  double t465 = 903125.0/36864.0*t52;
432  double t466 = 265625.0/73728.0*t54;
433  double t467 = 1328125.0/73728.0*t56;
434  double t468 = 265625.0/8192.0*t58;
435  double t469 = 4515625.0/36864.0*t60;
436  double t470 = 1328125.0/73728.0*t62;
437  double t471 = 6640625.0/73728.0*t64;
438  double t472 = 78125.0/16384.0*t66;
439  double t473 = 265625.0/73728.0*t68;
440  double t474 = 1328125.0/73728.0*t70;
441  double t475 = 390625.0/147456.0*t72;
442  double t476 = 1953125.0/147456.0*t74;
443  double t477 = 390625.0/16384.0*t76;
444  double t478 = 1328125.0/73728.0*t78;
445  double t479 = 6640625.0/73728.0*t80;
446  double t480 = 1953125.0/147456.0*t82;
447  double t481 = 9765625.0/147456.0*t84;
448  double t482 = t464-t465-t466+t467-t468+t469+t470-t471-t472-t473+t474+t475-t476+t477+t478-t479-t480+t481;
449  double t484 = 5625.0/16384.0+t446-t447-t448-t449+t450+t451-t452+t453+t454-t455+t456-t457-t458+t459+t460-t461+t462;
450  double t485 = t464-t465-t466+t467+t468-t469-t470+t471-t472-t473+t474+t475-t476-t477-t478+t479+t480-t481;
451  double t487 = -1875.0/32768.0-t413+t414+t415+t416-t417-t418+t419-t420-t421+t422-t423+t380+t424-t425-t426+t427-t428;
452  double t488 = -t387+t430+t431-t432-t433+t392+t434-t435+t436+t437-t438-t398+t439+t440+t441-t442-t443+t404;
453  double t490 = -1875.0/32768.0+t413-t414+t415-t416-t417+t418+t419-t420+t421-t422+t423+t380-t424-t425+t426+t427+t428;
454  double t491 = -t387-t430+t431+t432-t433+t392-t434-t435-t436+t437+t438-t398-t439+t440-t441-t442+t443+t404;
455  double t493 = 5625.0/16384.0-t446+t447-t448+t449+t450-t451-t452+t453-t454+t455-t456-t457+t458+t459-t460-t461-t462;
456  double t494 = t464+t465-t466-t467+t468-t469+t470+t471+t472-t473-t474+t475+t476-t477+t478+t479-t480-t481;
457  double t496 = 5625.0/16384.0+t446+t447-t448-t449+t450+t451-t452+t453-t454+t455+t456+t457-t458-t459+t460+t461-t462;
458  double t497 = t464+t465-t466-t467-t468+t469-t470-t471+t472-t473-t474+t475+t476+t477-t478-t479+t480+t481;
459  double t499 = -1875.0/32768.0-t413-t414+t415+t416-t417-t418+t419-t420+t421-t422-t423-t380+t424+t425-t426-t427+t428;
460  double t500 = -t387-t430+t431+t432+t433-t392+t434+t435-t436+t437+t438-t398-t439-t440+t441+t442-t443-t404;
461  double t502 = 625.0/65536.0-t331+t332-t333+t334+t335-t336-t337+t338-t339+t340-t341-t342+t343+t344-t345-t346-t347;
462  double t503 = t349+t350-t351-t352+t353-t354+t355+t356+t357-t358-t359+t360+t361-t362+t363+t364-t365-t366;
463  double t505 = -1875.0/32768.0+t369-t370+t371-t372-t373+t374+t375-t376+t377-t378+t379+t380-t381-t382+t383+t384+t385;
464  double t506 = -t387-t388+t389+t390-t391+t392-t393-t394-t395+t396+t397-t398-t399+t400-t401-t402+t403+t404;
465  double t508 = -1875.0/32768.0-t369-t370+t371+t372-t373-t374+t375-t376+t377-t378-t379-t380+t381+t382-t383-t384+t385;
466  double t509 = -t387-t388+t389+t390+t391-t392+t393+t394-t395+t396+t397-t398-t399-t400+t401+t402-t403-t404;
467  double t511 = 625.0/65536.0+t331+t332-t333-t334+t335+t336-t337+t338-t339+t340+t341+t342-t343-t344+t345+t346-t347;
468  double t512 = t349+t350-t351-t352-t353+t354-t355-t356+t357-t358-t359+t360+t361+t362-t363-t364+t365+t366;
469 
470  values[0] = t49+t86;
471  values[1] = t106+t125;
472  values[2] = t145+t164;
473  values[3] = t184;
474  values[4] = t185+t186;
475  values[5] = t188+t189;
476  values[6] = t195+t196;
477  values[7] = t215+t230;
478  values[8] = t249+t264;
479  values[9] = t282;
480  values[10] = t283+t284;
481  values[11] = t286+t287;
482  values[12] = t293+t294;
483  values[13] = t297+t298;
484  values[14] = t301+t302;
485  values[15] = t305;
486  values[16] = t306+t307;
487  values[17] = t309+t310;
488  values[18] = t312+t313;
489  values[19] = t316+t317;
490  values[20] = t320+t321;
491  values[21] = t324;
492  values[22] = t325+t326;
493  values[23] = t328+t329;
494  values[24] = t348+t367;
495  values[25] = t386+t405;
496  values[26] = t407+t408;
497  values[27] = t410+t411;
498  values[28] = t429+t444;
499  values[29] = t463+t482;
500  values[30] = t484+t485;
501  values[31] = t487+t488;
502  values[32] = t490+t491;
503  values[33] = t493+t494;
504  values[34] = t496+t497;
505  values[35] = t499+t500;
506  values[36] = t502+t503;
507  values[37] = t505+t506;
508  values[38] = t508+t509;
509  values[39] = t511+t512;
510 }
511 
512 // values of the derivatives in xi direction
513 static void C_Q_UL6SE_2D_DeriveXi(double xi, double eta, double *values)
514 {
515  double t1 = xi*xi;
516  double t2 = eta*eta;
517  double t3 = t2*eta;
518  double t4 = t1*t3;
519  double t5 = 38375.0/49152.0*t4;
520  double t6 = xi*eta;
521  double t7 = 4891279.0/8601600.0*t6;
522  double t8 = t1*eta;
523  double t9 = 1873.0/32768.0*t8;
524  double t10 = t1*xi;
525  double t11 = t10*eta;
526  double t12 = 1875847.0/1720320.0*t11;
527  double t13 = t1*t1;
528  double t14 = t13*eta;
529  double t15 = 42289.0/196608.0*t14;
530  double t16 = t2*xi;
531  double t17 = 496367.0/1843200.0*t16;
532  double t18 = t1*t2;
533  double t19 = 727871.0/1228800.0*t18;
534  double t20 = t10*t2;
535  double t21 = 463271.0/368640.0*t20;
536  double t22 = t13*t2;
537  double t23 = 538871.0/294912.0*t22;
538  double t24 = xi*t3;
539  double t25 = 727871.0/1843200.0*t24;
540  double t26 = t10*t3;
541  double t27 = 138955.0/73728.0*t26;
542  double t28 = t13*t3;
543  double t29 = 770375.0/294912.0*t28;
544  double t30 = t2*t2;
545  double t31 = xi*t30;
546  double t32 = 463271.0/737280.0*t31;
547  double t33 = t1*t30;
548  double t34 = 138955.0/98304.0*t33;
549  double t35 = t10*t30;
550  double t36 = 430175.0/147456.0*t35;
551  double t37 = t13*t30;
552  double t38 = 2528875.0/589824.0*t37;
553  double t39 = t5+t7-t9-t12+t15+t17-t19-t21+t23-t25+t27-t29-t32+t34+t36-t38;
554  double t40 = t30*eta;
555  double t41 = xi*t40;
556  double t42 = 538871.0/737280.0*t41;
557  double t43 = t1*t40;
558  double t44 = 154075.0/98304.0*t43;
559  double t45 = t10*t40;
560  double t46 = 505775.0/147456.0*t45;
561  double t47 = t13*t40;
562  double t48 = 2906875.0/589824.0*t47;
563  double t53 = (1.0-t2)*(3.0/8.0+35.0/8.0*t30-15.0/4.0*t2);
564  double t54 = 81.0/700.0*t53;
565  double t55 = 18489.0/819200.0*t1;
566  double t56 = 111833.0/114688.0*t10;
567  double t57 = 4651.0/65536.0*t13;
568  double t58 = 1564309.0/2867200.0*xi;
569  double t59 = 4891279.0/17203200.0*t2;
570  double t60 = 1873.0/98304.0*t3;
571  double t61 = 339.0/327680.0*eta;
572  double t62 = 1875847.0/6881280.0*t30;
573  double t63 = 42289.0/983040.0*t40;
574  double t64 = 1.0-eta;
575  double t65 = 1.0-t1;
576  double t69 = 35.0/2.0*t10-15.0/2.0*xi;
577  double t70 = t64*t65*t69;
578  double t71 = 81.0/700.0*t70;
579  double t75 = 3.0/8.0+35.0/8.0*t13-15.0/4.0*t1;
580  double t76 = t64*xi*t75;
581  double t77 = 81.0/350.0*t76;
582  double t78 = -421161.0/11468800.0+t42-t44-t46+t48+t54+t55+t56-t57-t58+t59-t60+t61-t62+t63-t71+t77;
583  double t80 = 1125.0/256.0*t4;
584  double t81 = 291519.0/89600.0*t6;
585  double t82 = 81.0/512.0*t8;
586  double t83 = 101007.0/17920.0*t11;
587  double t84 = 243.0/1024.0*t14;
588  double t85 = 7677.0/6400.0*t16;
589  double t86 = 1125.0/256.0*t18;
590  double t87 = 585.0/256.0*t20;
591  double t88 = 3375.0/512.0*t22;
592  double t89 = 7677.0/6400.0*t24;
593  double t90 = 585.0/256.0*t26;
594  double t91 = 3375.0/512.0*t28;
595  double t92 = 7677.0/2560.0*t31;
596  double t93 = 5625.0/512.0*t33;
597  double t94 = 2925.0/512.0*t35;
598  double t95 = 16875.0/1024.0*t37;
599  double t96 = -t80-t81+t82+t83-t84-t85+t86+t87-t88+t89-t90+t91+t92-t93-t94+t95;
600  double t97 = 7677.0/2560.0*t41;
601  double t98 = 5625.0/512.0*t43;
602  double t99 = 2925.0/512.0*t45;
603  double t100 = 16875.0/1024.0*t47;
604  double t101 = 81.0/512.0*t1;
605  double t102 = 101007.0/17920.0*t10;
606  double t103 = 243.0/1024.0*t13;
607  double t104 = 291519.0/89600.0*xi;
608  double t105 = 75.0/512.0*t2;
609  double t106 = 75.0/512.0*t3;
610  double t107 = 27.0/5120.0*eta;
611  double t108 = 375.0/1024.0*t30;
612  double t109 = 375.0/1024.0*t40;
613  double t110 = 243.0/350.0*t70;
614  double t111 = 243.0/175.0*t76;
615  double t112 = 27.0/5120.0-t97+t98+t99-t100-t101-t102+t103+t104-t105+t106-t107+t108-t109+t110-t111;
616  double t114 = 14625.0/2048.0*t4;
617  double t115 = 574533.0/71680.0*t6;
618  double t116 = 1053.0/4096.0*t8;
619  double t117 = 49653.0/3584.0*t11;
620  double t118 = 1215.0/4096.0*t14;
621  double t119 = 1521.0/5120.0*t16;
622  double t120 = 14625.0/2048.0*t18;
623  double t121 = 225.0/256.0*t20;
624  double t122 = 16875.0/2048.0*t22;
625  double t123 = 1521.0/5120.0*t24;
626  double t124 = 225.0/256.0*t26;
627  double t125 = 16875.0/2048.0*t28;
628  double t126 = 1521.0/2048.0*t31;
629  double t127 = 73125.0/4096.0*t33;
630  double t128 = 1125.0/512.0*t35;
631  double t129 = 84375.0/4096.0*t37;
632  double t130 = t114+t115-t116-t117+t118-t119-t120+t121+t122+t123-t124-t125+t126+t127-t128-t129;
633  double t131 = 1521.0/2048.0*t41;
634  double t132 = 73125.0/4096.0*t43;
635  double t133 = 1125.0/512.0*t45;
636  double t134 = 84375.0/4096.0*t47;
637  double t135 = 1053.0/4096.0*t1;
638  double t136 = 49653.0/3584.0*t10;
639  double t137 = 1215.0/4096.0*t13;
640  double t138 = 574533.0/71680.0*xi;
641  double t139 = 375.0/512.0*t2;
642  double t140 = 375.0/512.0*t3;
643  double t141 = 27.0/1024.0*eta;
644  double t142 = 1875.0/1024.0*t30;
645  double t143 = 1875.0/1024.0*t40;
646  double t144 = 243.0/140.0*t70;
647  double t145 = 243.0/70.0*t76;
648  double t146 = -27.0/1024.0-t131-t132+t133+t134+t135+t136-t137-t138+t139-t140+t141-t142+t143-t144+t145;
649  double t148 = 95073.0/8960.0*t6;
650  double t149 = 32913.0/1792.0*t11;
651  double t150 = 4823.0/1920.0*t16;
652  double t151 = 525.0/128.0*t20;
653  double t152 = 4823.0/1920.0*t24;
654  double t153 = 525.0/128.0*t26;
655  double t154 = 95073.0/8960.0*xi;
656  double t155 = 32913.0/1792.0*t10;
657  double t156 = 4823.0/768.0*t31;
658  double t157 = 2625.0/256.0*t35;
659  double t158 = 4823.0/768.0*t41;
660  double t159 = 2625.0/256.0*t45;
661  double t162 = -t148+t149+t150-t151-t152+t153+t154-t155-t156+t157+t158-t159-162.0/35.0*t76+81.0/35.0*t70;
662  double t163 = -t114+t115+t116-t117-t118-t119+t120+t121-t122+t123-t124+t125+t126-t127-t128+t129;
663  double t164 = 27.0/1024.0-t131+t132+t133-t134-t135+t136+t137-t138-t139+t140-t141+t142-t143-t144+t145;
664  double t166 = t80-t81-t82+t83+t84-t85-t86+t87+t88+t89-t90-t91+t92+t93-t94-t95;
665  double t167 = -27.0/5120.0-t97-t98+t99+t100+t101-t102-t103+t104+t105-t106+t107-t108+t109+t110-t111;
666  double t169 = -t5+t7+t9-t12-t15+t17+t19-t21-t23-t25+t27+t29-t32-t34+t36+t38;
667  double t170 = 421161.0/11468800.0+t42+t44-t46-t48-t54-t55+t56+t57-t58-t59+t60-t61+t62-t63-t71+t77;
668  double t172 = 75.0/256.0*t6;
669  double t173 = 225.0/512.0*t8;
670  double t174 = 375.0/256.0*t11;
671  double t175 = 1875.0/1024.0*t14;
672  double t176 = 23031.0/12800.0*t18;
673  double t177 = 7677.0/1280.0*t20;
674  double t178 = 7677.0/1024.0*t22;
675  double t179 = 375.0/128.0*t24;
676  double t180 = 1875.0/128.0*t26;
677  double t181 = 9375.0/512.0*t28;
678  double t182 = 585.0/512.0*t31;
679  double t183 = 1755.0/1024.0*t33;
680  double t184 = t80-t172-t173+t174+t175-t85-t176+t177+t178+t179-t180-t181+t182+t183-t94;
681  double t185 = 14625.0/2048.0*t37;
682  double t186 = 675.0/256.0*t41;
683  double t187 = 2025.0/512.0*t43;
684  double t188 = 3375.0/256.0*t45;
685  double t189 = 243.0/350.0*t53;
686  double t190 = 2187.0/25600.0*t1;
687  double t191 = 729.0/2560.0*t10;
688  double t192 = 729.0/2048.0*t13;
689  double t193 = 729.0/12800.0*xi;
690  double t194 = 291519.0/179200.0*t2;
691  double t195 = 27.0/512.0*t3;
692  double t196 = 101007.0/71680.0*t30;
693  double t197 = 243.0/5120.0*t40;
694  double t198 = -78003.0/358400.0-t185-t186-t187+t188+t100+t189+t190-t191-t192+t193+t194-t195+t107-t196+t197;
695  double t200 = 375.0/256.0*t6;
696  double t201 = 1125.0/512.0*t8;
697  double t202 = 1875.0/256.0*t11;
698  double t203 = 9375.0/1024.0*t14;
699  double t204 = 4563.0/10240.0*t18;
700  double t205 = 1521.0/1024.0*t20;
701  double t206 = 7605.0/4096.0*t22;
702  double t207 = 4875.0/1024.0*t24;
703  double t208 = 24375.0/1024.0*t26;
704  double t209 = 121875.0/4096.0*t28;
705  double t210 = 225.0/512.0*t31;
706  double t211 = 675.0/1024.0*t33;
707  double t212 = -t114+t200+t201-t202-t203-t119-t204+t205+t206-t207+t208+t209+t210+t211-t128;
708  double t213 = 5625.0/2048.0*t37;
709  double t214 = 3375.0/1024.0*t41;
710  double t215 = 10125.0/2048.0*t43;
711  double t216 = 16875.0/1024.0*t45;
712  double t217 = 243.0/140.0*t53;
713  double t218 = 2187.0/10240.0*t1;
714  double t219 = 729.0/1024.0*t10;
715  double t220 = 3645.0/4096.0*t13;
716  double t221 = 729.0/5120.0*xi;
717  double t222 = 574533.0/143360.0*t2;
718  double t223 = 351.0/4096.0*t3;
719  double t224 = 49653.0/14336.0*t30;
720  double t225 = 243.0/4096.0*t40;
721  double t226 = 78003.0/143360.0-t213+t214+t215-t216-t134-t217-t218+t219+t220-t221-t222+t223-t141+t224-t225;
722  double t228 = 95073.0/17920.0*t2;
723  double t229 = 4823.0/1280.0*t18;
724  double t230 = 4823.0/384.0*t20;
725  double t231 = 24115.0/1536.0*t22;
726  double t232 = 1771.0/3840.0*xi;
727  double t233 = 1771.0/2560.0*t1;
728  double t234 = 1771.0/768.0*t10;
729  double t235 = 8855.0/3072.0*t13;
730  double t236 = 32913.0/7168.0*t30;
731  double t237 = 525.0/256.0*t31;
732  double t238 = 1575.0/512.0*t33;
733  double t239 = 13125.0/1024.0*t37;
734  double t240 = 81.0/35.0*t53;
735  double t241 = -25581.0/35840.0+t228+t150+t229-t230-t231-t232-t233+t234+t235-t236-t237-t238+t157+t239+t240;
736  double t242 = t114-t200-t201+t202+t203-t119-t204+t205+t206+t207-t208-t209+t210+t211-t128;
737  double t243 = 78003.0/143360.0-t213-t214-t215+t216+t134-t217-t218+t219+t220-t221-t222-t223+t141+t224+t225;
738  double t245 = -t80+t172+t173-t174-t175-t85-t176+t177+t178-t179+t180+t181+t182+t183-t94;
739  double t246 = -78003.0/358400.0-t185+t186+t187-t188-t100+t189+t190-t191-t192+t193+t194+t195-t107-t196-t197;
740  double t248 = 1.0+eta;
741  double t250 = t248*xi*t75;
742  double t251 = 81.0/350.0*t250;
743  double t253 = t248*t65*t69;
744  double t254 = 81.0/700.0*t253;
745  double t255 = t251-t254+t5-t7-t9+t12+t15+t17+t19-t21-t23+t25-t27-t29-t32-t34;
746  double t256 = 421161.0/11468800.0+t36+t38-t42-t44+t46+t48-t54-t55+t56+t57-t58-t59-t60+t61+t62+t63;
747  double t258 = 243.0/175.0*t250;
748  double t259 = 243.0/350.0*t253;
749  double t260 = -t258+t259-t80+t81+t82-t83-t84-t85-t86+t87+t88-t89+t90+t91+t92+t93;
750  double t261 = -27.0/5120.0-t94-t95+t97+t98-t99-t100+t101-t102-t103+t104+t105+t106-t107-t108-t109;
751  double t263 = 243.0/70.0*t250;
752  double t264 = 243.0/140.0*t253;
753  double t265 = t263-t264+t114-t115-t116+t117+t118-t119+t120+t121-t122-t123+t124-t125+t126-t127;
754  double t266 = 27.0/1024.0-t128+t129+t131-t132-t133+t134-t135+t136+t137-t138-t139-t140+t141+t142+t143;
755  double t270 = t148-t149+t150-t151+t152-t153+t154-t155-t156+t157-t158+t159-162.0/35.0*t250+81.0/35.0*t253;
756  double t271 = t263-t264-t114-t115+t116+t117-t118-t119-t120+t121+t122-t123+t124+t125+t126+t127;
757  double t272 = -27.0/1024.0-t128-t129+t131+t132-t133-t134+t135+t136-t137-t138+t139+t140-t141-t142-t143;
758  double t274 = -t258+t259+t80+t81-t82-t83+t84-t85+t86+t87-t88-t89+t90-t91+t92-t93;
759  double t275 = 27.0/5120.0-t94+t95+t97-t98-t99+t100-t101-t102+t103+t104-t105-t106+t107+t108+t109;
760  double t277 = t251-t254-t5-t7+t9+t12-t15+t17-t19-t21+t23+t25-t27+t29-t32+t34;
761  double t278 = -421161.0/11468800.0+t36-t38-t42+t44+t46-t48+t54+t55+t56-t57-t58+t59+t60-t61-t62-t63;
762  double t280 = t80+t172-t173-t174+t175-t85+t176+t177-t178-t179+t180-t181+t182-t183-t94;
763  double t281 = 78003.0/358400.0+t185+t186-t187-t188+t100-t189-t190-t191+t192+t193-t194-t195+t107+t196+t197;
764  double t283 = -t114-t200+t201+t202-t203-t119+t204+t205-t206+t207-t208+t209+t210-t211-t128;
765  double t284 = -78003.0/143360.0+t213-t214+t215+t216-t134+t217+t218+t219-t220-t221+t222+t223-t141-t224-t225;
766  double t286 = 25581.0/35840.0-t228+t150-t229-t230+t231-t232+t233+t234-t235+t236-t237+t238+t157-t239-t240;
767  double t287 = t114+t200-t201-t202+t203-t119+t204+t205-t206-t207+t208-t209+t210-t211-t128;
768  double t288 = -78003.0/143360.0+t213+t214-t215-t216+t134+t217+t218+t219-t220-t221+t222-t223+t141-t224+t225;
769  double t290 = -t80-t172+t173+t174-t175-t85+t176+t177-t178+t179-t180+t181+t182-t183-t94;
770  double t291 = 78003.0/358400.0+t185-t186+t187+t188-t100-t189-t190-t191+t192+t193-t194+t195-t107+t196-t197;
771  double t293 = 2640625.0/49152.0*t4;
772  double t294 = 40625.0/49152.0*t6;
773  double t295 = 203125.0/98304.0*t8;
774  double t296 = 78125.0/49152.0*t11;
775  double t297 = 1953125.0/589824.0*t14;
776  double t298 = 105625.0/8192.0*t16;
777  double t299 = 528125.0/16384.0*t18;
778  double t300 = 203125.0/8192.0*t20;
779  double t301 = 5078125.0/98304.0*t22;
780  double t302 = 528125.0/24576.0*t24;
781  double t303 = 1015625.0/24576.0*t26;
782  double t304 = 25390625.0/294912.0*t28;
783  double t305 = 203125.0/16384.0*t31;
784  double t306 = 1015625.0/32768.0*t33;
785  double t307 = 390625.0/16384.0*t35;
786  double t308 = t293+t294-t295-t296+t297+t298-t299-t300+t301-t302+t303-t304-t305+t306+t307;
787  double t309 = 9765625.0/196608.0*t37;
788  double t310 = 1015625.0/49152.0*t41;
789  double t311 = 5078125.0/98304.0*t43;
790  double t312 = 1953125.0/49152.0*t45;
791  double t313 = 48828125.0/589824.0*t47;
792  double t314 = 40625.0/32768.0*t1;
793  double t315 = 15625.0/16384.0*t10;
794  double t316 = 390625.0/196608.0*t13;
795  double t317 = 8125.0/16384.0*xi;
796  double t318 = 40625.0/98304.0*t2;
797  double t319 = 203125.0/294912.0*t3;
798  double t320 = 15625.0/589824.0*eta;
799  double t321 = 78125.0/196608.0*t30;
800  double t322 = 390625.0/589824.0*t40;
801  double t323 = -3125.0/196608.0-t309+t310-t311-t312+t313+t314+t315-t316-t317+t318-t319+t320-t321+t322;
802  double t325 = 3453125.0/24576.0*t4;
803  double t326 = 121875.0/8192.0*t6;
804  double t327 = 609375.0/16384.0*t8;
805  double t328 = 234375.0/8192.0*t11;
806  double t329 = 1953125.0/32768.0*t14;
807  double t330 = 138125.0/12288.0*t16;
808  double t331 = 690625.0/24576.0*t18;
809  double t332 = 265625.0/12288.0*t20;
810  double t333 = 6640625.0/147456.0*t22;
811  double t334 = 690625.0/12288.0*t24;
812  double t335 = 1328125.0/12288.0*t26;
813  double t336 = 33203125.0/147456.0*t28;
814  double t337 = 203125.0/24576.0*t31;
815  double t338 = 1015625.0/49152.0*t33;
816  double t339 = 390625.0/24576.0*t35;
817  double t340 = -t325-t326+t327+t328-t329-t330+t331+t332-t333+t334-t335+t336+t337-t338-t339;
818  double t341 = 9765625.0/294912.0*t37;
819  double t342 = 1015625.0/24576.0*t41;
820  double t343 = 5078125.0/49152.0*t43;
821  double t344 = 1953125.0/24576.0*t45;
822  double t345 = 48828125.0/294912.0*t47;
823  double t346 = 121875.0/16384.0*t1;
824  double t347 = 46875.0/8192.0*t10;
825  double t348 = 390625.0/32768.0*t13;
826  double t349 = 24375.0/8192.0*xi;
827  double t350 = 53125.0/147456.0*t2;
828  double t351 = 265625.0/147456.0*t3;
829  double t352 = 15625.0/32768.0*eta;
830  double t353 = 78125.0/294912.0*t30;
831  double t354 = 390625.0/294912.0*t40;
832  double t355 = 3125.0/32768.0+t341-t342+t343+t344-t345-t346-t347+t348+t349-t350+t351-t352+t353-t354;
833  double t357 = t325+t326-t327-t328+t329-t330+t331+t332-t333-t334+t335-t336+t337-t338-t339;
834  double t358 = 3125.0/32768.0+t341+t342-t343-t344+t345-t346-t347+t348+t349-t350-t351+t352+t353+t354;
835  double t360 = -t293-t294+t295+t296-t297+t298-t299-t300+t301+t302-t303+t304-t305+t306+t307;
836  double t361 = -3125.0/196608.0-t309-t310+t311+t312-t313+t314+t315-t316-t317+t318+t319-t320-t321-t322;
837  double t363 = 53125.0/73728.0*t6;
838  double t364 = 265625.0/49152.0*t8;
839  double t365 = 78125.0/73728.0*t11;
840  double t366 = 1953125.0/294912.0*t14;
841  double t367 = 690625.0/8192.0*t18;
842  double t368 = 203125.0/12288.0*t20;
843  double t369 = 5078125.0/49152.0*t22;
844  double t370 = 690625.0/36864.0*t24;
845  double t371 = 1015625.0/36864.0*t26;
846  double t372 = 25390625.0/147456.0*t28;
847  double t373 = 265625.0/24576.0*t31;
848  double t374 = 1328125.0/16384.0*t33;
849  double t375 = -t325-t363+t364+t365-t366-t330+t367+t368-t369+t370-t371+t372+t373-t374-t339;
850  double t376 = 9765625.0/98304.0*t37;
851  double t377 = 1328125.0/73728.0*t41;
852  double t378 = 6640625.0/49152.0*t43;
853  double t379 = 1953125.0/73728.0*t45;
854  double t380 = 53125.0/16384.0*t1;
855  double t381 = 15625.0/24576.0*t10;
856  double t382 = 390625.0/98304.0*t13;
857  double t383 = 10625.0/24576.0*xi;
858  double t384 = 121875.0/16384.0*t2;
859  double t385 = 203125.0/16384.0*t3;
860  double t386 = 234375.0/32768.0*t30;
861  double t387 = 390625.0/32768.0*t40;
862  double t388 = 9375.0/32768.0+t376-t377+t378+t379-t345-t380-t381+t382+t383-t384+t385-t352+t386-t387;
863  double t390 = 4515625.0/12288.0*t4;
864  double t391 = 53125.0/4096.0*t6;
865  double t392 = 796875.0/8192.0*t8;
866  double t393 = 78125.0/4096.0*t11;
867  double t394 = 1953125.0/16384.0*t14;
868  double t395 = 180625.0/18432.0*t16;
869  double t396 = 903125.0/12288.0*t18;
870  double t397 = 265625.0/18432.0*t20;
871  double t398 = 6640625.0/73728.0*t22;
872  double t399 = 903125.0/18432.0*t24;
873  double t400 = 1328125.0/18432.0*t26;
874  double t401 = 33203125.0/73728.0*t28;
875  double t402 = 265625.0/36864.0*t31;
876  double t403 = 1328125.0/24576.0*t33;
877  double t404 = 390625.0/36864.0*t35;
878  double t405 = t390+t391-t392-t393+t394+t395-t396-t397+t398-t399+t400-t401-t402+t403+t404;
879  double t406 = 9765625.0/147456.0*t37;
880  double t407 = 1328125.0/36864.0*t41;
881  double t408 = 6640625.0/24576.0*t43;
882  double t409 = 1953125.0/36864.0*t45;
883  double t410 = 48828125.0/147456.0*t47;
884  double t411 = 159375.0/8192.0*t1;
885  double t412 = 15625.0/4096.0*t10;
886  double t413 = 390625.0/16384.0*t13;
887  double t414 = 10625.0/4096.0*xi;
888  double t415 = 53125.0/8192.0*t2;
889  double t416 = 265625.0/8192.0*t3;
890  double t417 = 140625.0/16384.0*eta;
891  double t418 = 78125.0/16384.0*t30;
892  double t419 = 390625.0/16384.0*t40;
893  double t420 = -28125.0/16384.0-t406+t407-t408-t409+t410+t411+t412-t413-t414+t415-t416+t417-t418+t419;
894  double t422 = -t390-t391+t392+t393-t394+t395-t396-t397+t398+t399-t400+t401-t402+t403+t404;
895  double t423 = -28125.0/16384.0-t406-t407+t408+t409-t410+t411+t412-t413-t414+t415+t416-t417-t418-t419;
896  double t425 = t325+t363-t364-t365+t366-t330+t367+t368-t369-t370+t371-t372+t373-t374-t339;
897  double t426 = 9375.0/32768.0+t376+t377-t378-t379+t345-t380-t381+t382+t383-t384-t385+t352+t386+t387;
898  double t428 = -9375.0/32768.0+t325-t363-t364+t365+t366-t330-t367+t368+t369+t370-t371-t372+t373+t374;
899  double t429 = -t339-t376-t377-t378+t379+t345+t380-t381-t382+t383+t384-t385+t352-t386+t387;
900  double t431 = -t390+t391+t392-t393-t394+t395+t396-t397-t398-t399+t400+t401-t402-t403+t404;
901  double t432 = 28125.0/16384.0+t406+t407+t408-t409-t410-t411+t412+t413-t414-t415+t416-t417+t418-t419;
902  double t434 = t390-t391-t392+t393+t394+t395+t396-t397-t398+t399-t400-t401-t402-t403+t404;
903  double t435 = 28125.0/16384.0+t406-t407-t408+t409+t410-t411+t412+t413-t414-t415-t416+t417+t418+t419;
904  double t437 = -9375.0/32768.0-t325+t363+t364-t365-t366-t330-t367+t368+t369-t370+t371+t372+t373+t374;
905  double t438 = -t339-t376+t377+t378-t379-t345+t380-t381-t382+t383+t384+t385-t352-t386-t387;
906  double t440 = 3125.0/196608.0-t293+t294+t295-t296-t297+t298+t299-t300-t301-t302+t303+t304-t305-t306;
907  double t441 = t307+t309+t310+t311-t312-t313-t314+t315+t316-t317-t318+t319-t320+t321-t322;
908  double t443 = -3125.0/32768.0+t325-t326-t327+t328+t329-t330-t331+t332+t333+t334-t335-t336+t337+t338;
909  double t444 = -t339-t341-t342-t343+t344+t345+t346-t347-t348+t349+t350-t351+t352-t353+t354;
910  double t446 = -3125.0/32768.0-t325+t326+t327-t328-t329-t330-t331+t332+t333-t334+t335+t336+t337+t338;
911  double t447 = -t339-t341+t342+t343-t344-t345+t346-t347-t348+t349+t350+t351-t352-t353-t354;
912  double t449 = 3125.0/196608.0+t293-t294-t295+t296+t297+t298+t299-t300-t301+t302-t303-t304-t305-t306;
913  double t450 = t307+t309-t310-t311+t312+t313-t314+t315+t316-t317-t318-t319+t320+t321+t322;
914 
915  values[0] = t39+t78;
916  values[1] = t96+t112;
917  values[2] = t130+t146;
918  values[3] = t162;
919  values[4] = t163+t164;
920  values[5] = t166+t167;
921  values[6] = t169+t170;
922  values[7] = t184+t198;
923  values[8] = t212+t226;
924  values[9] = t241;
925  values[10] = t242+t243;
926  values[11] = t245+t246;
927  values[12] = t255+t256;
928  values[13] = t260+t261;
929  values[14] = t265+t266;
930  values[15] = t270;
931  values[16] = t271+t272;
932  values[17] = t274+t275;
933  values[18] = t277+t278;
934  values[19] = t280+t281;
935  values[20] = t283+t284;
936  values[21] = t286;
937  values[22] = t287+t288;
938  values[23] = t290+t291;
939  values[24] = t308+t323;
940  values[25] = t340+t355;
941  values[26] = t357+t358;
942  values[27] = t360+t361;
943  values[28] = t375+t388;
944  values[29] = t405+t420;
945  values[30] = t422+t423;
946  values[31] = t425+t426;
947  values[32] = t428+t429;
948  values[33] = t431+t432;
949  values[34] = t434+t435;
950  values[35] = t437+t438;
951  values[36] = t440+t441;
952  values[37] = t443+t444;
953  values[38] = t446+t447;
954  values[39] = t449+t450;
955 }
956 
957 // values of the derivatives in eta direction
958 static void C_Q_UL6SE_2D_DeriveEta(double xi, double eta, double *values)
959 {
960  double t1 = 1.0-xi;
961  double t3 = eta*eta;
962  double t4 = t3*t3;
963  double t7 = 3.0/8.0+35.0/8.0*t4-15.0/4.0*t3;
964  double t8 = t1*eta*t7;
965  double t9 = 81.0/350.0*t8;
966  double t10 = xi*xi;
967  double t12 = t10*t10;
968  double t16 = (1.0-t10)*(3.0/8.0+35.0/8.0*t12-15.0/4.0*t10);
969  double t17 = 81.0/700.0*t16;
970  double t18 = t3*eta;
971  double t19 = t10*t18;
972  double t20 = 463271.0/368640.0*t19;
973  double t21 = xi*eta;
974  double t22 = 4891279.0/8601600.0*t21;
975  double t23 = t10*eta;
976  double t24 = 496367.0/1843200.0*t23;
977  double t25 = t10*xi;
978  double t26 = t25*eta;
979  double t27 = 727871.0/1843200.0*t26;
980  double t28 = t12*eta;
981  double t29 = 463271.0/737280.0*t28;
982  double t30 = t12*xi;
983  double t31 = t30*eta;
984  double t32 = 538871.0/737280.0*t31;
985  double t33 = xi*t3;
986  double t34 = 1873.0/32768.0*t33;
987  double t35 = t10*t3;
988  double t36 = 727871.0/1228800.0*t35;
989  double t37 = t25*t3;
990  double t38 = 38375.0/49152.0*t37;
991  double t39 = t12*t3;
992  double t40 = 138955.0/98304.0*t39;
993  double t41 = t30*t3;
994  double t42 = 154075.0/98304.0*t41;
995  double t43 = xi*t18;
996  double t44 = 1875847.0/1720320.0*t43;
997  double t45 = t25*t18;
998  double t46 = 138955.0/73728.0*t45;
999  double t47 = t12*t18;
1000  double t48 = 430175.0/147456.0*t47;
1001  double t49 = t9+t17-t20+t22+t24-t27-t29+t32-t34-t36+t38+t40-t42-t44+t46+t48;
1002  double t50 = t30*t18;
1003  double t51 = 505775.0/147456.0*t50;
1004  double t52 = xi*t4;
1005  double t53 = 42289.0/196608.0*t52;
1006  double t54 = t10*t4;
1007  double t55 = 538871.0/294912.0*t54;
1008  double t56 = t25*t4;
1009  double t57 = 770375.0/294912.0*t56;
1010  double t58 = 4891279.0/17203200.0*t10;
1011  double t59 = 1873.0/98304.0*t25;
1012  double t60 = 1875847.0/6881280.0*t12;
1013  double t61 = 42289.0/983040.0*t30;
1014  double t62 = 18489.0/819200.0*t3;
1015  double t63 = 111833.0/114688.0*t18;
1016  double t64 = 4651.0/65536.0*t4;
1017  double t65 = t30*t4;
1018  double t66 = 2906875.0/589824.0*t65;
1019  double t67 = 339.0/327680.0*xi;
1020  double t68 = 1564309.0/2867200.0*eta;
1021  double t69 = t12*t4;
1022  double t70 = 2528875.0/589824.0*t69;
1023  double t71 = 1.0-t3;
1024  double t75 = 35.0/2.0*t18-15.0/2.0*eta;
1025  double t76 = t1*t71*t75;
1026  double t77 = 81.0/700.0*t76;
1027  double t78 = -421161.0/11468800.0-t51+t53+t55-t57+t58-t59-t60+t61+t62+t63-t64+t66+t67-t68-t70-t77;
1028  double t80 = 243.0/350.0*t16;
1029  double t81 = 7677.0/1280.0*t19;
1030  double t82 = 75.0/256.0*t21;
1031  double t83 = 7677.0/6400.0*t23;
1032  double t84 = 375.0/128.0*t26;
1033  double t85 = 585.0/512.0*t28;
1034  double t86 = 675.0/256.0*t31;
1035  double t87 = 225.0/512.0*t33;
1036  double t88 = 23031.0/12800.0*t35;
1037  double t89 = 1125.0/256.0*t37;
1038  double t90 = 1755.0/1024.0*t39;
1039  double t91 = 2025.0/512.0*t41;
1040  double t92 = 375.0/256.0*t43;
1041  double t93 = 1875.0/128.0*t45;
1042  double t94 = 2925.0/512.0*t47;
1043  double t95 = -t80+t81-t82-t83+t84+t85-t86+t87+t88-t89-t90+t91+t92-t93-t94;
1044  double t96 = 3375.0/256.0*t50;
1045  double t97 = 1875.0/1024.0*t52;
1046  double t98 = 7677.0/1024.0*t54;
1047  double t99 = 9375.0/512.0*t56;
1048  double t100 = 291519.0/179200.0*t10;
1049  double t101 = 27.0/512.0*t25;
1050  double t102 = 101007.0/71680.0*t12;
1051  double t103 = 243.0/5120.0*t30;
1052  double t104 = 2187.0/25600.0*t3;
1053  double t105 = 729.0/2560.0*t18;
1054  double t106 = 729.0/2048.0*t4;
1055  double t107 = 16875.0/1024.0*t65;
1056  double t108 = 27.0/5120.0*xi;
1057  double t109 = 729.0/12800.0*eta;
1058  double t110 = 14625.0/2048.0*t69;
1059  double t111 = 78003.0/358400.0+t96-t97-t98+t99-t100+t101+t102-t103-t104-t105+t106-t107-t108+t109+t110;
1060  double t113 = 243.0/140.0*t16;
1061  double t114 = 1521.0/1024.0*t19;
1062  double t115 = 375.0/256.0*t21;
1063  double t116 = 1521.0/5120.0*t23;
1064  double t117 = 4875.0/1024.0*t26;
1065  double t118 = 225.0/512.0*t28;
1066  double t119 = 3375.0/1024.0*t31;
1067  double t120 = 1125.0/512.0*t33;
1068  double t121 = 4563.0/10240.0*t35;
1069  double t122 = 14625.0/2048.0*t37;
1070  double t123 = 675.0/1024.0*t39;
1071  double t124 = 10125.0/2048.0*t41;
1072  double t125 = 1875.0/256.0*t43;
1073  double t126 = 24375.0/1024.0*t45;
1074  double t127 = 1125.0/512.0*t47;
1075  double t128 = t113+t114+t115-t116-t117+t118+t119-t120+t121+t122-t123-t124-t125+t126-t127;
1076  double t129 = 16875.0/1024.0*t50;
1077  double t130 = 9375.0/1024.0*t52;
1078  double t131 = 7605.0/4096.0*t54;
1079  double t132 = 121875.0/4096.0*t56;
1080  double t133 = 574533.0/143360.0*t10;
1081  double t134 = 351.0/4096.0*t25;
1082  double t135 = 49653.0/14336.0*t12;
1083  double t136 = 243.0/4096.0*t30;
1084  double t137 = 2187.0/10240.0*t3;
1085  double t138 = 729.0/1024.0*t18;
1086  double t139 = 3645.0/4096.0*t4;
1087  double t140 = 84375.0/4096.0*t65;
1088  double t141 = 27.0/1024.0*xi;
1089  double t142 = 729.0/5120.0*eta;
1090  double t143 = 5625.0/2048.0*t69;
1091  double t144 = -78003.0/143360.0-t129+t130-t131-t132+t133-t134-t135+t136+t137+t138-t139+t140+t141-t142+t143;
1092  double t146 = 95073.0/17920.0*t10;
1093  double t147 = 32913.0/7168.0*t12;
1094  double t148 = 4823.0/1920.0*t23;
1095  double t149 = 525.0/256.0*t28;
1096  double t150 = 4823.0/1280.0*t35;
1097  double t151 = 1575.0/512.0*t39;
1098  double t152 = 1771.0/3840.0*eta;
1099  double t153 = 1771.0/2560.0*t3;
1100  double t154 = 1771.0/768.0*t18;
1101  double t155 = 8855.0/3072.0*t4;
1102  double t156 = 4823.0/384.0*t19;
1103  double t157 = 2625.0/256.0*t47;
1104  double t158 = 24115.0/1536.0*t54;
1105  double t159 = 13125.0/1024.0*t69;
1106  double t160 = 81.0/35.0*t16;
1107  double t161 = 25581.0/35840.0-t146+t147+t148-t149-t150+t151-t152+t153+t154-t155-t156+t157+t158-t159-t160;
1108  double t162 = t113+t114-t115-t116+t117+t118-t119+t120+t121-t122-t123+t124+t125-t126-t127;
1109  double t163 = -78003.0/143360.0+t129-t130-t131+t132+t133+t134-t135-t136+t137+t138-t139-t140-t141-t142+t143;
1110  double t165 = -t80+t81+t82-t83-t84+t85+t86-t87+t88+t89-t90-t91-t92+t93-t94;
1111  double t166 = 78003.0/358400.0-t96+t97-t98-t99-t100-t101+t102+t103-t104-t105+t106+t107+t108+t109+t110;
1112  double t168 = 1.0+xi;
1113  double t170 = t168*t71*t75;
1114  double t171 = 81.0/700.0*t170;
1115  double t173 = t168*eta*t7;
1116  double t174 = 81.0/350.0*t173;
1117  double t175 = -t171+t174+t17-t20-t22+t24+t27-t29-t32+t34-t36-t38+t40+t42+t44-t46;
1118  double t176 = -421161.0/11468800.0+t48+t51-t53+t55+t57+t58+t59-t60-t61+t62+t63-t64-t66-t67-t68-t70;
1119  double t178 = 243.0/350.0*t170;
1120  double t179 = 243.0/175.0*t173;
1121  double t180 = 585.0/256.0*t19;
1122  double t181 = 291519.0/89600.0*t21;
1123  double t182 = 7677.0/6400.0*t26;
1124  double t183 = 7677.0/2560.0*t28;
1125  double t184 = 7677.0/2560.0*t31;
1126  double t185 = 81.0/512.0*t33;
1127  double t186 = 1125.0/256.0*t35;
1128  double t187 = 5625.0/512.0*t39;
1129  double t188 = 5625.0/512.0*t41;
1130  double t189 = 101007.0/17920.0*t43;
1131  double t190 = 585.0/256.0*t45;
1132  double t191 = t178-t179+t180+t181-t83-t182+t183+t184-t185+t186+t89-t187-t188-t189+t190-t94;
1133  double t192 = 2925.0/512.0*t50;
1134  double t193 = 243.0/1024.0*t52;
1135  double t194 = 3375.0/512.0*t54;
1136  double t195 = 3375.0/512.0*t56;
1137  double t196 = 75.0/512.0*t10;
1138  double t197 = 75.0/512.0*t25;
1139  double t198 = 375.0/1024.0*t12;
1140  double t199 = 375.0/1024.0*t30;
1141  double t200 = 81.0/512.0*t3;
1142  double t201 = 101007.0/17920.0*t18;
1143  double t202 = 243.0/1024.0*t4;
1144  double t203 = 291519.0/89600.0*eta;
1145  double t204 = 16875.0/1024.0*t69;
1146  double t205 = 27.0/5120.0-t192+t193-t194-t195-t196-t197+t198+t199-t200-t201+t202+t107+t108+t203+t204;
1147  double t207 = 243.0/140.0*t170;
1148  double t208 = 243.0/70.0*t173;
1149  double t209 = 225.0/256.0*t19;
1150  double t210 = 574533.0/71680.0*t21;
1151  double t211 = 1521.0/5120.0*t26;
1152  double t212 = 1521.0/2048.0*t28;
1153  double t213 = 1521.0/2048.0*t31;
1154  double t214 = 1053.0/4096.0*t33;
1155  double t215 = 14625.0/2048.0*t35;
1156  double t216 = 73125.0/4096.0*t39;
1157  double t217 = 73125.0/4096.0*t41;
1158  double t218 = 49653.0/3584.0*t43;
1159  double t219 = 225.0/256.0*t45;
1160  double t220 = -t207+t208+t209-t210-t116-t211+t212+t213+t214-t215-t122+t216+t217+t218+t219-t127;
1161  double t221 = 1125.0/512.0*t50;
1162  double t222 = 1215.0/4096.0*t52;
1163  double t223 = 16875.0/2048.0*t54;
1164  double t224 = 16875.0/2048.0*t56;
1165  double t225 = 375.0/512.0*t10;
1166  double t226 = 375.0/512.0*t25;
1167  double t227 = 1875.0/1024.0*t12;
1168  double t228 = 1875.0/1024.0*t30;
1169  double t229 = 1053.0/4096.0*t3;
1170  double t230 = 49653.0/3584.0*t18;
1171  double t231 = 1215.0/4096.0*t4;
1172  double t232 = 574533.0/71680.0*eta;
1173  double t233 = 84375.0/4096.0*t69;
1174  double t234 = -27.0/1024.0-t221-t222+t223+t224+t225+t226-t227-t228+t229+t230-t231-t140-t141-t232-t233;
1175  double t236 = 95073.0/8960.0*t21;
1176  double t237 = 4823.0/1920.0*t26;
1177  double t238 = 4823.0/768.0*t28;
1178  double t239 = 4823.0/768.0*t31;
1179  double t240 = 95073.0/8960.0*eta;
1180  double t241 = 32913.0/1792.0*t18;
1181  double t242 = 32913.0/1792.0*t43;
1182  double t243 = 525.0/128.0*t19;
1183  double t244 = 525.0/128.0*t45;
1184  double t245 = 2625.0/256.0*t50;
1185  double t248 = t236+t148+t237-t238-t239+t240-t241-t242-t243-t244+t157+t245-162.0/35.0*t173+81.0/35.0*t170;
1186  double t249 = -t207+t208+t209-t210-t116-t211+t212+t213-t214+t215+t122-t216-t217+t218+t219-t127;
1187  double t250 = 27.0/1024.0-t221+t222-t223-t224-t225-t226+t227+t228-t229+t230+t231+t140+t141-t232+t233;
1188  double t252 = t178-t179+t180+t181-t83-t182+t183+t184+t185-t186-t89+t187+t188-t189+t190-t94;
1189  double t253 = -27.0/5120.0-t192-t193+t194+t195+t196+t197-t198-t199+t200-t201-t202-t107-t108+t203-t204;
1190  double t255 = -t171+t174-t17-t20-t22+t24+t27-t29-t32-t34+t36+t38-t40-t42+t44-t46;
1191  double t256 = 421161.0/11468800.0+t48+t51+t53-t55-t57-t58-t59+t60+t61-t62+t63+t64+t66+t67-t68+t70;
1192  double t258 = t80+t81+t82-t83-t84+t85+t86+t87-t88-t89+t90+t91-t92+t93-t94;
1193  double t259 = -78003.0/358400.0-t96-t97+t98+t99+t100+t101-t102-t103+t104-t105-t106-t107-t108+t109-t110;
1194  double t261 = -t113+t114-t115-t116+t117+t118-t119-t120-t121+t122+t123-t124+t125-t126-t127;
1195  double t262 = 78003.0/143360.0+t129+t130+t131-t132-t133-t134+t135+t136-t137+t138+t139+t140+t141-t142-t143;
1196  double t264 = -25581.0/35840.0+t146-t147+t148-t149+t150-t151-t152-t153+t154+t155-t156+t157-t158+t159+t160;
1197  double t265 = -t113+t114+t115-t116-t117+t118+t119+t120-t121-t122+t123+t124-t125+t126-t127;
1198  double t266 = 78003.0/143360.0-t129-t130+t131+t132-t133+t134+t135-t136-t137+t138+t139-t140-t141-t142-t143;
1199  double t268 = t80+t81-t82-t83+t84+t85-t86-t87-t88+t89+t90-t91+t92-t93-t94;
1200  double t269 = -78003.0/358400.0+t96+t97+t98-t99+t100-t101-t102+t103+t104-t105-t106+t107+t108+t109-t110;
1201  double t271 = t9-t17-t20+t22+t24-t27-t29+t32+t34+t36-t38-t40+t42-t44+t46+t48;
1202  double t272 = 421161.0/11468800.0-t51-t53-t55+t57-t58+t59+t60-t61-t62+t63+t64-t66-t67-t68+t70-t77;
1203  double t274 = 243.0/175.0*t8;
1204  double t275 = -t274+t180-t181-t83+t182+t183-t184-t185-t186+t89+t187-t188+t189-t190-t94+t192;
1205  double t276 = 243.0/350.0*t76;
1206  double t277 = -27.0/5120.0+t193+t194-t195+t196-t197-t198+t199+t200-t201-t202+t107+t108+t203-t204+t276;
1207  double t279 = 243.0/70.0*t8;
1208  double t280 = t279+t209+t210-t116+t211+t212-t213+t214+t215-t122-t216+t217-t218-t219-t127+t221;
1209  double t281 = 243.0/140.0*t76;
1210  double t282 = 27.0/1024.0-t222-t223+t224-t225+t226+t227-t228-t229+t230+t231-t140-t141-t232+t233-t281;
1211  double t286 = -t236+t148-t237-t238+t239+t240-t241+t242-t243+t244+t157-t245-162.0/35.0*t8+81.0/35.0*t76;
1212  double t287 = t279+t209+t210-t116+t211+t212-t213-t214-t215+t122+t216-t217-t218-t219-t127+t221;
1213  double t288 = -27.0/1024.0+t222+t223-t224+t225-t226-t227+t228+t229+t230-t231+t140+t141-t232-t233-t281;
1214  double t290 = -t274+t180-t181-t83+t182+t183-t184+t185+t186-t89-t187+t188+t189-t190-t94+t192;
1215  double t291 = 27.0/5120.0-t193-t194+t195-t196+t197+t198-t199-t200-t201+t202-t107-t108+t203+t204+t276;
1216  double t293 = 203125.0/8192.0*t19;
1217  double t294 = 40625.0/49152.0*t21;
1218  double t295 = 105625.0/8192.0*t23;
1219  double t296 = 528125.0/24576.0*t26;
1220  double t297 = 203125.0/16384.0*t28;
1221  double t298 = 1015625.0/49152.0*t31;
1222  double t299 = 203125.0/98304.0*t33;
1223  double t300 = 528125.0/16384.0*t35;
1224  double t301 = 2640625.0/49152.0*t37;
1225  double t302 = 1015625.0/32768.0*t39;
1226  double t303 = 5078125.0/98304.0*t41;
1227  double t304 = 78125.0/49152.0*t43;
1228  double t305 = 1015625.0/24576.0*t45;
1229  double t306 = 390625.0/16384.0*t47;
1230  double t307 = 1953125.0/49152.0*t50;
1231  double t308 = -t293+t294+t295-t296-t297+t298-t299-t300+t301+t302-t303-t304+t305+t306-t307;
1232  double t309 = 1953125.0/589824.0*t52;
1233  double t310 = 5078125.0/98304.0*t54;
1234  double t311 = 25390625.0/294912.0*t56;
1235  double t312 = 40625.0/98304.0*t10;
1236  double t313 = 203125.0/294912.0*t25;
1237  double t314 = 78125.0/196608.0*t12;
1238  double t315 = 390625.0/589824.0*t30;
1239  double t316 = 40625.0/32768.0*t3;
1240  double t317 = 15625.0/16384.0*t18;
1241  double t318 = 390625.0/196608.0*t4;
1242  double t319 = 48828125.0/589824.0*t65;
1243  double t320 = 15625.0/589824.0*xi;
1244  double t321 = 8125.0/16384.0*eta;
1245  double t322 = 9765625.0/196608.0*t69;
1246  double t323 = -3125.0/196608.0+t309+t310-t311+t312-t313-t314+t315+t316+t317-t318+t319+t320-t321-t322;
1247  double t325 = 203125.0/12288.0*t19;
1248  double t326 = 53125.0/73728.0*t21;
1249  double t327 = 138125.0/12288.0*t23;
1250  double t328 = 690625.0/36864.0*t26;
1251  double t329 = 265625.0/24576.0*t28;
1252  double t330 = 1328125.0/73728.0*t31;
1253  double t331 = 265625.0/49152.0*t33;
1254  double t332 = 690625.0/8192.0*t35;
1255  double t333 = 3453125.0/24576.0*t37;
1256  double t334 = 1328125.0/16384.0*t39;
1257  double t335 = 6640625.0/49152.0*t41;
1258  double t336 = 78125.0/73728.0*t43;
1259  double t337 = 1015625.0/36864.0*t45;
1260  double t338 = 390625.0/24576.0*t47;
1261  double t339 = 1953125.0/73728.0*t50;
1262  double t340 = t325-t326-t327+t328+t329-t330+t331+t332-t333-t334+t335+t336-t337-t338+t339;
1263  double t341 = 1953125.0/294912.0*t52;
1264  double t342 = 5078125.0/49152.0*t54;
1265  double t343 = 25390625.0/147456.0*t56;
1266  double t344 = 121875.0/16384.0*t10;
1267  double t345 = 203125.0/16384.0*t25;
1268  double t346 = 234375.0/32768.0*t12;
1269  double t347 = 390625.0/32768.0*t30;
1270  double t348 = 53125.0/16384.0*t3;
1271  double t349 = 15625.0/24576.0*t18;
1272  double t350 = 390625.0/98304.0*t4;
1273  double t351 = 48828125.0/294912.0*t65;
1274  double t352 = 15625.0/32768.0*xi;
1275  double t353 = 10625.0/24576.0*eta;
1276  double t354 = 9765625.0/98304.0*t69;
1277  double t355 = 9375.0/32768.0-t341-t342+t343-t344+t345+t346-t347-t348-t349+t350-t351-t352+t353+t354;
1278  double t357 = -9375.0/32768.0+t325-t326-t327+t328+t329-t330-t331-t332+t333+t334-t335+t336-t337-t338;
1279  double t358 = t339+t341+t342-t343+t344-t345-t346+t347+t348-t349-t350+t351+t352+t353-t354;
1280  double t360 = 3125.0/196608.0-t293+t294+t295-t296-t297+t298+t299+t300-t301-t302+t303-t304+t305+t306;
1281  double t361 = -t307-t309-t310+t311-t312+t313+t314-t315-t316+t317+t318-t319-t320-t321+t322;
1282  double t363 = 265625.0/12288.0*t19;
1283  double t364 = 121875.0/8192.0*t21;
1284  double t365 = 690625.0/12288.0*t26;
1285  double t366 = 203125.0/24576.0*t28;
1286  double t367 = 1015625.0/24576.0*t31;
1287  double t368 = 609375.0/16384.0*t33;
1288  double t369 = 690625.0/24576.0*t35;
1289  double t370 = 1015625.0/49152.0*t39;
1290  double t371 = 5078125.0/49152.0*t41;
1291  double t372 = 234375.0/8192.0*t43;
1292  double t373 = 1328125.0/12288.0*t45;
1293  double t374 = 1953125.0/24576.0*t50;
1294  double t375 = t363-t364-t327+t365+t366-t367+t368+t369-t333-t370+t371+t372-t373-t338+t374;
1295  double t376 = 1953125.0/32768.0*t52;
1296  double t377 = 6640625.0/147456.0*t54;
1297  double t378 = 33203125.0/147456.0*t56;
1298  double t379 = 53125.0/147456.0*t10;
1299  double t380 = 265625.0/147456.0*t25;
1300  double t381 = 78125.0/294912.0*t12;
1301  double t382 = 390625.0/294912.0*t30;
1302  double t383 = 121875.0/16384.0*t3;
1303  double t384 = 46875.0/8192.0*t18;
1304  double t385 = 390625.0/32768.0*t4;
1305  double t386 = 24375.0/8192.0*eta;
1306  double t387 = 9765625.0/294912.0*t69;
1307  double t388 = 3125.0/32768.0-t376-t377+t378-t379+t380+t381-t382-t383-t384+t385-t351-t352+t386+t387;
1308  double t390 = 265625.0/18432.0*t19;
1309  double t391 = 53125.0/4096.0*t21;
1310  double t392 = 180625.0/18432.0*t23;
1311  double t393 = 903125.0/18432.0*t26;
1312  double t394 = 265625.0/36864.0*t28;
1313  double t395 = 1328125.0/36864.0*t31;
1314  double t396 = 796875.0/8192.0*t33;
1315  double t397 = 903125.0/12288.0*t35;
1316  double t398 = 4515625.0/12288.0*t37;
1317  double t399 = 1328125.0/24576.0*t39;
1318  double t400 = 6640625.0/24576.0*t41;
1319  double t401 = 78125.0/4096.0*t43;
1320  double t402 = 1328125.0/18432.0*t45;
1321  double t403 = 390625.0/36864.0*t47;
1322  double t404 = 1953125.0/36864.0*t50;
1323  double t405 = -t390+t391+t392-t393-t394+t395-t396-t397+t398+t399-t400-t401+t402+t403-t404;
1324  double t406 = 1953125.0/16384.0*t52;
1325  double t407 = 6640625.0/73728.0*t54;
1326  double t408 = 33203125.0/73728.0*t56;
1327  double t409 = 53125.0/8192.0*t10;
1328  double t410 = 265625.0/8192.0*t25;
1329  double t411 = 78125.0/16384.0*t12;
1330  double t412 = 390625.0/16384.0*t30;
1331  double t413 = 159375.0/8192.0*t3;
1332  double t414 = 15625.0/4096.0*t18;
1333  double t415 = 390625.0/16384.0*t4;
1334  double t416 = 48828125.0/147456.0*t65;
1335  double t417 = 140625.0/16384.0*xi;
1336  double t418 = 10625.0/4096.0*eta;
1337  double t419 = 9765625.0/147456.0*t69;
1338  double t420 = -28125.0/16384.0+t406+t407-t408+t409-t410-t411+t412+t413+t414-t415+t416+t417-t418-t419;
1339  double t422 = -t390+t391+t392-t393-t394+t395+t396+t397-t398-t399+t400-t401+t402+t403-t404;
1340  double t423 = 28125.0/16384.0-t406-t407+t408-t409+t410+t411-t412-t413+t414+t415-t416-t417-t418+t419;
1341  double t425 = -3125.0/32768.0+t363-t364-t327+t365+t366-t367-t368-t369+t333+t370-t371+t372-t373-t338;
1342  double t426 = t374+t376+t377-t378+t379-t380-t381+t382+t383-t384-t385+t351+t352+t386-t387;
1343  double t428 = t363+t364-t327-t365+t366+t367-t368+t369+t333-t370-t371-t372+t373-t338-t374;
1344  double t429 = 3125.0/32768.0+t376-t377-t378-t379-t380+t381+t382-t383-t384+t385+t351+t352+t386+t387;
1345  double t431 = -t390-t391+t392+t393-t394-t395+t396-t397-t398+t399+t400+t401-t402+t403+t404;
1346  double t432 = -28125.0/16384.0-t406+t407+t408+t409+t410-t411-t412+t413+t414-t415-t416-t417-t418-t419;
1347  double t434 = -t390-t391+t392+t393-t394-t395-t396+t397+t398-t399-t400+t401-t402+t403+t404;
1348  double t435 = 28125.0/16384.0+t406-t407-t408-t409-t410+t411+t412-t413+t414+t415+t416+t417-t418+t419;
1349  double t437 = -3125.0/32768.0+t363+t364-t327-t365+t366+t367+t368-t369-t333+t370+t371-t372+t373-t338;
1350  double t438 = -t374-t376+t377+t378+t379+t380-t381-t382+t383-t384-t385-t351-t352+t386-t387;
1351  double t440 = -t293-t294+t295+t296-t297-t298+t299-t300-t301+t302+t303+t304-t305+t306+t307;
1352  double t441 = -3125.0/196608.0-t309+t310+t311+t312+t313-t314-t315+t316+t317-t318-t319-t320-t321-t322;
1353  double t443 = t325+t326-t327-t328+t329+t330-t331+t332+t333-t334-t335-t336+t337-t338-t339;
1354  double t444 = 9375.0/32768.0+t341-t342-t343-t344-t345+t346+t347-t348-t349+t350+t351+t352+t353+t354;
1355  double t446 = -9375.0/32768.0+t325+t326-t327-t328+t329+t330+t331-t332-t333+t334+t335-t336+t337-t338;
1356  double t447 = -t339-t341+t342+t343+t344+t345-t346-t347+t348-t349-t350-t351-t352+t353-t354;
1357  double t449 = 3125.0/196608.0-t293-t294+t295+t296-t297-t298-t299+t300+t301-t302-t303+t304-t305+t306;
1358  double t450 = t307+t309-t310-t311-t312-t313+t314+t315-t316+t317+t318+t319+t320-t321+t322;
1359 
1360  values[0] = t49+t78;
1361  values[1] = t95+t111;
1362  values[2] = t128+t144;
1363  values[3] = t161;
1364  values[4] = t162+t163;
1365  values[5] = t165+t166;
1366  values[6] = t175+t176;
1367  values[7] = t191+t205;
1368  values[8] = t220+t234;
1369  values[9] = t248;
1370  values[10] = t249+t250;
1371  values[11] = t252+t253;
1372  values[12] = t255+t256;
1373  values[13] = t258+t259;
1374  values[14] = t261+t262;
1375  values[15] = t264;
1376  values[16] = t265+t266;
1377  values[17] = t268+t269;
1378  values[18] = t271+t272;
1379  values[19] = t275+t277;
1380  values[20] = t280+t282;
1381  values[21] = t286;
1382  values[22] = t287+t288;
1383  values[23] = t290+t291;
1384  values[24] = t308+t323;
1385  values[25] = t340+t355;
1386  values[26] = t357+t358;
1387  values[27] = t360+t361;
1388  values[28] = t375+t388;
1389  values[29] = t405+t420;
1390  values[30] = t422+t423;
1391  values[31] = t425+t426;
1392  values[32] = t428+t429;
1393  values[33] = t431+t432;
1394  values[34] = t434+t435;
1395  values[35] = t437+t438;
1396  values[36] = t440+t441;
1397  values[37] = t443+t444;
1398  values[38] = t446+t447;
1399  values[39] = t449+t450;
1400 }
1401 
1402 // values of the derivatives in xi-xi direction
1403 static void C_Q_UL6SE_2D_DeriveXiXi(double xi, double eta, double *values)
1404 {
1405  double t1 = xi*xi;
1406  double t2 = t1*xi;
1407  double t3 = 4651.0/16384.0*t2;
1408  double t4 = eta*eta;
1409  double t5 = 496367.0/1843200.0*t4;
1410  double t6 = t4*eta;
1411  double t7 = 727871.0/1843200.0*t6;
1412  double t8 = t4*t4;
1413  double t9 = 463271.0/737280.0*t8;
1414  double t10 = t8*eta;
1415  double t11 = 538871.0/737280.0*t10;
1416  double t12 = 18489.0/409600.0*xi;
1417  double t13 = 4891279.0/8601600.0*eta;
1418  double t14 = 335499.0/114688.0*t1;
1419  double t15 = t1*t6;
1420  double t16 = 138955.0/24576.0*t15;
1421  double t17 = xi*eta;
1422  double t18 = 1873.0/16384.0*t17;
1423  double t19 = t1*eta;
1424  double t20 = 1875847.0/573440.0*t19;
1425  double t21 = t2*eta;
1426  double t22 = 42289.0/49152.0*t21;
1427  double t23 = xi*t4;
1428  double t24 = 727871.0/614400.0*t23;
1429  double t25 = -t3+t5-t7-t9+t11+t12+t13+t14+t16-t18-t20+t22-t24;
1430  double t26 = t1*t4;
1431  double t27 = 463271.0/122880.0*t26;
1432  double t28 = t2*t4;
1433  double t29 = 538871.0/73728.0*t28;
1434  double t30 = xi*t6;
1435  double t31 = 38375.0/24576.0*t30;
1436  double t32 = t2*t6;
1437  double t33 = 770375.0/73728.0*t32;
1438  double t34 = xi*t8;
1439  double t35 = 138955.0/49152.0*t34;
1440  double t36 = t1*t8;
1441  double t37 = 430175.0/49152.0*t36;
1442  double t38 = t2*t8;
1443  double t39 = 2528875.0/147456.0*t38;
1444  double t40 = t10*xi;
1445  double t41 = 154075.0/49152.0*t40;
1446  double t42 = t1*t10;
1447  double t43 = 505775.0/49152.0*t42;
1448  double t44 = t2*t10;
1449  double t45 = 2906875.0/147456.0*t44;
1450  double t46 = 1.0-eta;
1451  double t47 = t1*t1;
1452  double t50 = 3.0/8.0+35.0/8.0*t47-15.0/4.0*t1;
1453  double t51 = t46*t50;
1454  double t52 = 81.0/350.0*t51;
1455  double t56 = 35.0/2.0*t2-15.0/2.0*xi;
1456  double t57 = t46*xi*t56;
1457  double t58 = 81.0/175.0*t57;
1458  double t59 = 1.0-t1;
1459  double t62 = 105.0/2.0*t1-15.0/2.0;
1460  double t63 = t46*t59*t62;
1461  double t64 = 81.0/700.0*t63;
1462  double t65 = -1564309.0/2867200.0-t27+t29+t31-t33+t35+t37-t39-t41-t43+t45+t52+t58-t64;
1463  double t67 = 243.0/256.0*t2;
1464  double t68 = 7677.0/6400.0*t4;
1465  double t69 = 7677.0/6400.0*t6;
1466  double t70 = 7677.0/2560.0*t8;
1467  double t71 = 7677.0/2560.0*t10;
1468  double t72 = 81.0/256.0*xi;
1469  double t73 = 291519.0/89600.0*eta;
1470  double t74 = 303021.0/17920.0*t1;
1471  double t75 = 1755.0/256.0*t15;
1472  double t76 = 81.0/256.0*t17;
1473  double t77 = 303021.0/17920.0*t19;
1474  double t78 = 243.0/256.0*t21;
1475  double t79 = 1125.0/128.0*t23;
1476  double t80 = t67-t68+t69+t70-t71-t72-t73-t74-t75+t76+t77-t78+t79;
1477  double t81 = 1755.0/256.0*t26;
1478  double t82 = 3375.0/128.0*t28;
1479  double t83 = 1125.0/128.0*t30;
1480  double t84 = 3375.0/128.0*t32;
1481  double t85 = 5625.0/256.0*t34;
1482  double t86 = 8775.0/512.0*t36;
1483  double t87 = 16875.0/256.0*t38;
1484  double t88 = 5625.0/256.0*t40;
1485  double t89 = 8775.0/512.0*t42;
1486  double t90 = 16875.0/256.0*t44;
1487  double t91 = 243.0/175.0*t51;
1488  double t92 = 486.0/175.0*t57;
1489  double t93 = 243.0/350.0*t63;
1490  double t94 = 291519.0/89600.0+t81-t82-t83+t84-t85-t86+t87+t88+t89-t90-t91-t92+t93;
1491  double t96 = 1215.0/1024.0*t2;
1492  double t97 = 1521.0/5120.0*t4;
1493  double t98 = 1521.0/5120.0*t6;
1494  double t99 = 1521.0/2048.0*t8;
1495  double t100 = 1521.0/2048.0*t10;
1496  double t101 = 1053.0/2048.0*xi;
1497  double t102 = 574533.0/71680.0*eta;
1498  double t103 = 148959.0/3584.0*t1;
1499  double t104 = 675.0/256.0*t15;
1500  double t105 = 1053.0/2048.0*t17;
1501  double t106 = 148959.0/3584.0*t19;
1502  double t107 = 1215.0/1024.0*t21;
1503  double t108 = -574533.0/71680.0-t96-t97+t98+t99-t100+t101+t102+t103-t104-t105-t106+t107;
1504  double t109 = 14625.0/1024.0*t23;
1505  double t110 = 675.0/256.0*t26;
1506  double t111 = 16875.0/512.0*t28;
1507  double t112 = 14625.0/1024.0*t30;
1508  double t113 = 16875.0/512.0*t32;
1509  double t114 = 73125.0/2048.0*t34;
1510  double t115 = 3375.0/512.0*t36;
1511  double t116 = 84375.0/1024.0*t38;
1512  double t117 = 73125.0/2048.0*t40;
1513  double t118 = 3375.0/512.0*t42;
1514  double t119 = 84375.0/1024.0*t44;
1515  double t120 = 243.0/70.0*t51;
1516  double t121 = 243.0/35.0*t57;
1517  double t122 = 243.0/140.0*t63;
1518  double t123 = -t109+t110+t111+t112-t113+t114-t115-t116-t117+t118+t119+t120+t121-t122;
1519  double t125 = 95073.0/8960.0*eta;
1520  double t126 = 98739.0/1792.0*t19;
1521  double t127 = 4823.0/1920.0*t4;
1522  double t128 = 1575.0/128.0*t26;
1523  double t129 = 4823.0/1920.0*t6;
1524  double t130 = 1575.0/128.0*t15;
1525  double t131 = 98739.0/1792.0*t1;
1526  double t132 = 4823.0/768.0*t8;
1527  double t133 = 7875.0/256.0*t36;
1528  double t134 = 4823.0/768.0*t10;
1529  double t135 = 7875.0/256.0*t42;
1530  double t139 = -t125+t126+t127-t128-t129+t130+95073.0/8960.0-t131-t132+t133+t134-t135-162.0/35.0*t51-324.0/35.0*t57+81.0/35.0*t63;
1531  double t140 = -574533.0/71680.0+t96-t97+t98+t99-t100-t101+t102+t103-t104+t105-t106-t107;
1532  double t141 = t109+t110-t111-t112+t113-t114-t115+t116+t117+t118-t119+t120+t121-t122;
1533  double t143 = -t67-t68+t69+t70-t71+t72-t73-t74-t75-t76+t77+t78-t79;
1534  double t144 = 291519.0/89600.0+t81+t82+t83-t84+t85-t86-t87-t88+t89+t90-t91-t92+t93;
1535  double t146 = t3+t5-t7-t9+t11-t12+t13+t14+t16+t18-t20-t22+t24;
1536  double t147 = -1564309.0/2867200.0-t27-t29-t31+t33-t35+t37+t39+t41-t43-t45+t52+t58-t64;
1537  double t149 = 2187.0/12800.0*xi;
1538  double t150 = 2187.0/2560.0*t1;
1539  double t151 = 729.0/512.0*t2;
1540  double t152 = 75.0/256.0*eta;
1541  double t153 = 225.0/256.0*t17;
1542  double t154 = 1125.0/256.0*t19;
1543  double t155 = 1875.0/256.0*t21;
1544  double t156 = 23031.0/6400.0*t23;
1545  double t157 = 23031.0/1280.0*t26;
1546  double t158 = 7677.0/256.0*t28;
1547  double t159 = t149-t150-t151+t83-t152-t153+t154+t155-t68-t156+t157+t158;
1548  double t160 = 375.0/128.0*t6;
1549  double t161 = 5625.0/128.0*t15;
1550  double t162 = 9375.0/128.0*t32;
1551  double t163 = 585.0/512.0*t8;
1552  double t164 = 1755.0/512.0*t34;
1553  double t165 = 14625.0/512.0*t38;
1554  double t166 = 675.0/256.0*t10;
1555  double t167 = 2025.0/256.0*t40;
1556  double t168 = 10125.0/256.0*t42;
1557  double t169 = 729.0/12800.0+t160-t161-t162+t163+t164-t86-t165-t166-t167+t168+t90;
1558  double t171 = 2187.0/5120.0*xi;
1559  double t172 = 2187.0/1024.0*t1;
1560  double t173 = 3645.0/1024.0*t2;
1561  double t174 = 375.0/256.0*eta;
1562  double t175 = 1125.0/256.0*t17;
1563  double t176 = 5625.0/256.0*t19;
1564  double t177 = 9375.0/256.0*t21;
1565  double t178 = 4563.0/5120.0*t23;
1566  double t179 = 4563.0/1024.0*t26;
1567  double t180 = 7605.0/1024.0*t28;
1568  double t181 = -t171+t172+t173-t112+t174+t175-t176-t177-t97-t178+t179+t180;
1569  double t182 = 4875.0/1024.0*t6;
1570  double t183 = 73125.0/1024.0*t15;
1571  double t184 = 121875.0/1024.0*t32;
1572  double t185 = 225.0/512.0*t8;
1573  double t186 = 675.0/512.0*t34;
1574  double t187 = 5625.0/512.0*t38;
1575  double t188 = 3375.0/1024.0*t10;
1576  double t189 = 10125.0/1024.0*t40;
1577  double t190 = 50625.0/1024.0*t42;
1578  double t191 = -729.0/5120.0-t182+t183+t184+t185+t186-t115-t187+t188+t189-t190-t119;
1579  double t193 = 4823.0/640.0*t23;
1580  double t194 = 4823.0/128.0*t26;
1581  double t195 = 24115.0/384.0*t28;
1582  double t196 = 1771.0/1280.0*xi;
1583  double t197 = 1771.0/256.0*t1;
1584  double t198 = 8855.0/768.0*t2;
1585  double t199 = 525.0/256.0*t8;
1586  double t200 = 1575.0/256.0*t34;
1587  double t201 = 13125.0/256.0*t38;
1588  double t202 = t127+t193-t194-t195-1771.0/3840.0-t196+t197+t198-t199-t200+t133+t201;
1589  double t203 = -t171+t172+t173+t112-t174-t175+t176+t177-t97-t178+t179+t180;
1590  double t204 = -729.0/5120.0+t182-t183-t184+t185+t186-t115-t187-t188-t189+t190+t119;
1591  double t206 = t149-t150-t151-t83+t152+t153-t154-t155-t68-t156+t157+t158;
1592  double t207 = 729.0/12800.0-t160+t161+t162+t163+t164-t86-t165+t166+t167-t168-t90;
1593  double t209 = 1.0+eta;
1594  double t210 = t209*t50;
1595  double t211 = 81.0/350.0*t210;
1596  double t212 = t3+t5+t7-t9-t11-t12-t13+t14+t211-t16-t18+t20+t22;
1597  double t214 = t209*t59*t62;
1598  double t215 = 81.0/700.0*t214;
1599  double t217 = t209*xi*t56;
1600  double t218 = 81.0/175.0*t217;
1601  double t219 = -1564309.0/2867200.0+t24-t27-t29+t31-t33-t35+t37+t39-t41+t43+t45-t215+t218;
1602  double t221 = 243.0/175.0*t210;
1603  double t222 = -t67-t68-t69+t70+t71+t72+t73-t74-t221+t75+t76-t77-t78;
1604  double t223 = 243.0/350.0*t214;
1605  double t224 = 486.0/175.0*t217;
1606  double t225 = 291519.0/89600.0-t79+t81+t82-t83+t84+t85-t86-t87+t88-t89-t90+t223-t224;
1607  double t227 = 243.0/70.0*t210;
1608  double t228 = -574533.0/71680.0+t96-t97-t98+t99+t100-t101-t102+t103+t227+t104-t105+t106;
1609  double t229 = 243.0/140.0*t214;
1610  double t230 = 243.0/35.0*t217;
1611  double t231 = t107+t109+t110-t111+t112-t113-t114-t115+t116-t117-t118+t119-t229+t230;
1612  double t236 = t125-t126+t127-t128+t129-t130+95073.0/8960.0-t131-t132+t133-t134+t135-162.0/35.0*t210-324.0/35.0*t217+81.0/35.0*t214;
1613  double t237 = -574533.0/71680.0-t96-t97-t98+t99+t100+t101-t102+t103+t227+t104+t105+t106;
1614  double t238 = -t107-t109+t110+t111-t112+t113+t114-t115-t116+t117-t118-t119-t229+t230;
1615  double t240 = t67-t68-t69+t70+t71-t72+t73-t74-t221+t75-t76-t77+t78;
1616  double t241 = 291519.0/89600.0+t79+t81-t82+t83-t84-t85-t86+t87-t88-t89+t90+t223-t224;
1617  double t243 = -t3+t5+t7-t9-t11+t12-t13+t14+t211-t16+t18+t20-t22;
1618  double t244 = -1564309.0/2867200.0-t24-t27+t29-t31+t33+t35+t37-t39+t41+t43-t45-t215+t218;
1619  double t246 = -t149-t150+t151+t83+t152-t153-t154+t155-t68+t156+t157-t158;
1620  double t247 = 729.0/12800.0-t160+t161-t162+t163-t164-t86+t165+t166-t167-t168+t90;
1621  double t249 = t171+t172-t173-t112-t174+t175+t176-t177-t97+t178+t179-t180;
1622  double t250 = -729.0/5120.0+t182-t183+t184+t185-t186-t115+t187-t188+t189+t190-t119;
1623  double t252 = t127-t193-t194+t195-1771.0/3840.0+t196+t197-t198-t199+t200+t133-t201;
1624  double t253 = t171+t172-t173+t112+t174-t175-t176+t177-t97+t178+t179-t180;
1625  double t254 = -729.0/5120.0-t182+t183-t184+t185-t186-t115+t187+t188-t189-t190+t119;
1626  double t256 = -t149-t150+t151-t83-t152+t153+t154-t155-t68+t156+t157-t158;
1627  double t257 = 729.0/12800.0+t160-t161+t162+t163-t164-t86+t165-t166+t167+t168-t90;
1628  double t259 = 40625.0/16384.0*xi;
1629  double t260 = 46875.0/16384.0*t1;
1630  double t261 = 390625.0/49152.0*t2;
1631  double t262 = 2640625.0/24576.0*t30;
1632  double t263 = 40625.0/49152.0*eta;
1633  double t264 = 203125.0/49152.0*t17;
1634  double t265 = 78125.0/16384.0*t19;
1635  double t266 = 1953125.0/147456.0*t21;
1636  double t267 = 105625.0/8192.0*t4;
1637  double t268 = 528125.0/8192.0*t23;
1638  double t269 = 609375.0/8192.0*t26;
1639  double t270 = 5078125.0/24576.0*t28;
1640  double t271 = t259+t260-t261+t262+t263-t264-t265+t266+t267-t268-t269+t270;
1641  double t272 = 528125.0/24576.0*t6;
1642  double t273 = 1015625.0/8192.0*t15;
1643  double t274 = 25390625.0/73728.0*t32;
1644  double t275 = 203125.0/16384.0*t8;
1645  double t276 = 1015625.0/16384.0*t34;
1646  double t277 = 1171875.0/16384.0*t36;
1647  double t278 = 9765625.0/49152.0*t38;
1648  double t279 = 1015625.0/49152.0*t10;
1649  double t280 = 5078125.0/49152.0*t40;
1650  double t281 = 1953125.0/16384.0*t42;
1651  double t282 = 48828125.0/147456.0*t44;
1652  double t283 = -8125.0/16384.0-t272+t273-t274-t275+t276+t277-t278+t279-t280-t281+t282;
1653  double t285 = 121875.0/8192.0*xi;
1654  double t286 = 140625.0/8192.0*t1;
1655  double t287 = 390625.0/8192.0*t2;
1656  double t288 = 3453125.0/12288.0*t30;
1657  double t289 = 121875.0/8192.0*eta;
1658  double t290 = 609375.0/8192.0*t17;
1659  double t291 = 703125.0/8192.0*t19;
1660  double t292 = 1953125.0/8192.0*t21;
1661  double t293 = 138125.0/12288.0*t4;
1662  double t294 = 690625.0/12288.0*t23;
1663  double t295 = 265625.0/4096.0*t26;
1664  double t296 = 6640625.0/36864.0*t28;
1665  double t297 = -t285-t286+t287-t288-t289+t290+t291-t292-t293+t294+t295-t296;
1666  double t298 = 690625.0/12288.0*t6;
1667  double t299 = 1328125.0/4096.0*t15;
1668  double t300 = 33203125.0/36864.0*t32;
1669  double t301 = 203125.0/24576.0*t8;
1670  double t302 = 1015625.0/24576.0*t34;
1671  double t303 = 390625.0/8192.0*t36;
1672  double t304 = 9765625.0/73728.0*t38;
1673  double t305 = 1015625.0/24576.0*t10;
1674  double t306 = 5078125.0/24576.0*t40;
1675  double t307 = 1953125.0/8192.0*t42;
1676  double t308 = 48828125.0/73728.0*t44;
1677  double t309 = 24375.0/8192.0+t298-t299+t300+t301-t302-t303+t304-t305+t306+t307-t308;
1678  double t311 = -t285-t286+t287+t288+t289-t290-t291+t292-t293+t294+t295-t296;
1679  double t312 = 24375.0/8192.0-t298+t299-t300+t301-t302-t303+t304+t305-t306-t307+t308;
1680  double t314 = t259+t260-t261-t262-t263+t264+t265-t266+t267-t268-t269+t270;
1681  double t315 = -8125.0/16384.0+t272-t273+t274-t275+t276+t277-t278-t279+t280+t281-t282;
1682  double t317 = 53125.0/8192.0*xi;
1683  double t318 = 15625.0/8192.0*t1;
1684  double t319 = 390625.0/24576.0*t2;
1685  double t320 = 53125.0/73728.0*eta;
1686  double t321 = 265625.0/24576.0*t17;
1687  double t322 = 78125.0/24576.0*t19;
1688  double t323 = 1953125.0/73728.0*t21;
1689  double t324 = 690625.0/4096.0*t23;
1690  double t325 = 203125.0/4096.0*t26;
1691  double t326 = 5078125.0/12288.0*t28;
1692  double t327 = -t317-t318+t319-t288-t320+t321+t322-t323-t293+t324+t325-t326;
1693  double t328 = 690625.0/36864.0*t6;
1694  double t329 = 1015625.0/12288.0*t15;
1695  double t330 = 25390625.0/36864.0*t32;
1696  double t331 = 265625.0/24576.0*t8;
1697  double t332 = 1328125.0/8192.0*t34;
1698  double t333 = 9765625.0/24576.0*t38;
1699  double t334 = 1328125.0/73728.0*t10;
1700  double t335 = 6640625.0/24576.0*t40;
1701  double t336 = 1953125.0/24576.0*t42;
1702  double t337 = 10625.0/24576.0+t328-t329+t330+t331-t332-t303+t333-t334+t335+t336-t308;
1703  double t339 = 159375.0/4096.0*xi;
1704  double t340 = 46875.0/4096.0*t1;
1705  double t341 = 390625.0/4096.0*t2;
1706  double t342 = 4515625.0/6144.0*t30;
1707  double t343 = 53125.0/4096.0*eta;
1708  double t344 = 796875.0/4096.0*t17;
1709  double t345 = 234375.0/4096.0*t19;
1710  double t346 = 1953125.0/4096.0*t21;
1711  double t347 = 180625.0/18432.0*t4;
1712  double t348 = 903125.0/6144.0*t23;
1713  double t349 = 265625.0/6144.0*t26;
1714  double t350 = 6640625.0/18432.0*t28;
1715  double t351 = t339+t340-t341+t342+t343-t344-t345+t346+t347-t348-t349+t350;
1716  double t352 = 903125.0/18432.0*t6;
1717  double t353 = 1328125.0/6144.0*t15;
1718  double t354 = 33203125.0/18432.0*t32;
1719  double t355 = 265625.0/36864.0*t8;
1720  double t356 = 1328125.0/12288.0*t34;
1721  double t357 = 390625.0/12288.0*t36;
1722  double t358 = 9765625.0/36864.0*t38;
1723  double t359 = 1328125.0/36864.0*t10;
1724  double t360 = 6640625.0/12288.0*t40;
1725  double t361 = 1953125.0/12288.0*t42;
1726  double t362 = 48828125.0/36864.0*t44;
1727  double t363 = -10625.0/4096.0-t352+t353-t354-t355+t356+t357-t358+t359-t360-t361+t362;
1728  double t365 = t339+t340-t341-t342-t343+t344+t345-t346+t347-t348-t349+t350;
1729  double t366 = -10625.0/4096.0+t352-t353+t354-t355+t356+t357-t358-t359+t360+t361-t362;
1730  double t368 = -t317-t318+t319+t288+t320-t321-t322+t323-t293+t324+t325-t326;
1731  double t369 = 10625.0/24576.0-t328+t329-t330+t331-t332-t303+t333+t334-t335-t336+t308;
1732  double t371 = t317-t318-t319+t288-t320-t321+t322+t323-t293-t324+t325+t326;
1733  double t372 = 10625.0/24576.0+t328-t329-t330+t331+t332-t303-t333-t334-t335+t336+t308;
1734  double t374 = -t339+t340+t341-t342+t343+t344-t345-t346+t347+t348-t349-t350;
1735  double t375 = -10625.0/4096.0-t352+t353+t354-t355-t356+t357+t358+t359+t360-t361-t362;
1736  double t377 = -t339+t340+t341+t342-t343-t344+t345+t346+t347+t348-t349-t350;
1737  double t378 = -10625.0/4096.0+t352-t353-t354-t355-t356+t357+t358-t359-t360+t361+t362;
1738  double t380 = t317-t318-t319-t288+t320+t321-t322-t323-t293-t324+t325+t326;
1739  double t381 = 10625.0/24576.0-t328+t329+t330+t331+t332-t303-t333+t334+t335-t336-t308;
1740  double t383 = -t259+t260+t261-t262+t263+t264-t265-t266+t267+t268-t269-t270;
1741  double t384 = -8125.0/16384.0-t272+t273+t274-t275-t276+t277+t278+t279+t280-t281-t282;
1742  double t386 = t285-t286-t287+t288-t289-t290+t291+t292-t293-t294+t295+t296;
1743  double t387 = 24375.0/8192.0+t298-t299-t300+t301+t302-t303-t304-t305-t306+t307+t308;
1744  double t389 = t285-t286-t287-t288+t289+t290-t291-t292-t293-t294+t295+t296;
1745  double t390 = 24375.0/8192.0-t298+t299+t300+t301+t302-t303-t304+t305+t306-t307-t308;
1746  double t392 = -t259+t260+t261+t262-t263-t264+t265+t266+t267+t268-t269-t270;
1747  double t393 = -8125.0/16384.0+t272-t273-t274-t275-t276+t277+t278-t279-t280+t281+t282;
1748 
1749  values[0] = t25+t65;
1750  values[1] = t80+t94;
1751  values[2] = t108+t123;
1752  values[3] = t139;
1753  values[4] = t140+t141;
1754  values[5] = t143+t144;
1755  values[6] = t146+t147;
1756  values[7] = t159+t169;
1757  values[8] = t181+t191;
1758  values[9] = t202;
1759  values[10] = t203+t204;
1760  values[11] = t206+t207;
1761  values[12] = t212+t219;
1762  values[13] = t222+t225;
1763  values[14] = t228+t231;
1764  values[15] = t236;
1765  values[16] = t237+t238;
1766  values[17] = t240+t241;
1767  values[18] = t243+t244;
1768  values[19] = t246+t247;
1769  values[20] = t249+t250;
1770  values[21] = t252;
1771  values[22] = t253+t254;
1772  values[23] = t256+t257;
1773  values[24] = t271+t283;
1774  values[25] = t297+t309;
1775  values[26] = t311+t312;
1776  values[27] = t314+t315;
1777  values[28] = t327+t337;
1778  values[29] = t351+t363;
1779  values[30] = t365+t366;
1780  values[31] = t368+t369;
1781  values[32] = t371+t372;
1782  values[33] = t374+t375;
1783  values[34] = t377+t378;
1784  values[35] = t380+t381;
1785  values[36] = t383+t384;
1786  values[37] = t386+t387;
1787  values[38] = t389+t390;
1788  values[39] = t392+t393;
1789 }
1790 
1791 // values of the derivatives in xi-eta direction
1792 static void C_Q_UL6SE_2D_DeriveXiEta(double xi, double eta, double *values)
1793 {
1794  double t1 = xi*xi;
1795  double t2 = 1873.0/32768.0*t1;
1796  double t3 = t1*xi;
1797  double t4 = 1875847.0/1720320.0*t3;
1798  double t5 = t1*t1;
1799  double t6 = 42289.0/196608.0*t5;
1800  double t7 = eta*eta;
1801  double t8 = 1873.0/32768.0*t7;
1802  double t9 = t7*eta;
1803  double t10 = 1875847.0/1720320.0*t9;
1804  double t11 = 4891279.0/8601600.0*xi;
1805  double t12 = 4891279.0/8601600.0*eta;
1806  double t13 = t7*t7;
1807  double t14 = 42289.0/196608.0*t13;
1808  double t15 = t1*t9;
1809  double t16 = 138955.0/24576.0*t15;
1810  double t17 = xi*eta;
1811  double t18 = 496367.0/921600.0*t17;
1812  double t19 = t1*eta;
1813  double t20 = 727871.0/614400.0*t19;
1814  double t21 = t3*eta;
1815  double t22 = 463271.0/184320.0*t21;
1816  double t23 = t5*eta;
1817  double t24 = 538871.0/147456.0*t23;
1818  double t25 = xi*t7;
1819  double t26 = 727871.0/614400.0*t25;
1820  double t27 = -t2-t4+t6-t8-t10+t11+t12+t14+t16+t18-t20-t22+t24-t26;
1821  double t28 = t1*t7;
1822  double t29 = 38375.0/16384.0*t28;
1823  double t30 = t3*t7;
1824  double t31 = 138955.0/24576.0*t30;
1825  double t32 = t5*t7;
1826  double t33 = 770375.0/98304.0*t32;
1827  double t34 = xi*t9;
1828  double t35 = 463271.0/184320.0*t34;
1829  double t36 = t3*t9;
1830  double t37 = 430175.0/36864.0*t36;
1831  double t38 = t5*t9;
1832  double t39 = 2528875.0/147456.0*t38;
1833  double t40 = xi*t13;
1834  double t41 = 538871.0/147456.0*t40;
1835  double t42 = t1*t13;
1836  double t43 = 770375.0/98304.0*t42;
1837  double t44 = t3*t13;
1838  double t45 = 2528875.0/147456.0*t44;
1839  double t46 = t5*t13;
1840  double t47 = 14534375.0/589824.0*t46;
1841  double t52 = (1.0-t7)*(35.0/2.0*t9-15.0/2.0*eta);
1842  double t53 = 81.0/700.0*t52;
1843  double t57 = eta*(3.0/8.0+35.0/8.0*t13-15.0/4.0*t7);
1844  double t58 = 81.0/350.0*t57;
1845  double t63 = (1.0-t1)*(35.0/2.0*t3-15.0/2.0*xi);
1846  double t64 = 81.0/700.0*t63;
1847  double t68 = xi*(3.0/8.0+35.0/8.0*t5-15.0/4.0*t1);
1848  double t69 = 81.0/350.0*t68;
1849  double t70 = 339.0/327680.0+t29+t31-t33-t35+t37-t39+t41-t43-t45+t47+t53-t58+t64-t69;
1850  double t72 = 81.0/512.0*t1;
1851  double t73 = 101007.0/17920.0*t3;
1852  double t74 = 243.0/1024.0*t5;
1853  double t75 = 225.0/512.0*t7;
1854  double t76 = 375.0/256.0*t9;
1855  double t77 = 291519.0/89600.0*xi;
1856  double t78 = 75.0/256.0*eta;
1857  double t79 = 1875.0/1024.0*t13;
1858  double t80 = 5625.0/128.0*t15;
1859  double t81 = 7677.0/3200.0*t17;
1860  double t82 = 1125.0/128.0*t19;
1861  double t83 = 585.0/128.0*t21;
1862  double t84 = -27.0/5120.0+t72+t73-t74+t75+t76-t77-t78-t79-t80-t81+t82+t83;
1863  double t85 = 3375.0/256.0*t23;
1864  double t86 = 23031.0/6400.0*t25;
1865  double t87 = 3375.0/256.0*t28;
1866  double t88 = 1755.0/256.0*t30;
1867  double t89 = 10125.0/512.0*t32;
1868  double t90 = 7677.0/640.0*t34;
1869  double t91 = 2925.0/128.0*t36;
1870  double t92 = 16875.0/256.0*t38;
1871  double t93 = 7677.0/512.0*t40;
1872  double t94 = 28125.0/512.0*t42;
1873  double t95 = 14625.0/512.0*t44;
1874  double t96 = 84375.0/1024.0*t46;
1875  double t97 = 243.0/350.0*t63;
1876  double t98 = 243.0/175.0*t68;
1877  double t99 = -t85+t86-t87-t88+t89+t90-t91+t92-t93+t94+t95-t96-t97+t98;
1878  double t101 = 1053.0/4096.0*t1;
1879  double t102 = 49653.0/3584.0*t3;
1880  double t103 = 1215.0/4096.0*t5;
1881  double t104 = 1125.0/512.0*t7;
1882  double t105 = 1875.0/256.0*t9;
1883  double t106 = 574533.0/71680.0*xi;
1884  double t107 = 375.0/256.0*eta;
1885  double t108 = 9375.0/1024.0*t13;
1886  double t109 = 73125.0/1024.0*t15;
1887  double t110 = 1521.0/2560.0*t17;
1888  double t111 = 14625.0/1024.0*t19;
1889  double t112 = 225.0/128.0*t21;
1890  double t113 = 16875.0/1024.0*t23;
1891  double t114 = -t101-t102+t103-t104-t105+t106+t107+t108+t109-t110-t111+t112+t113;
1892  double t115 = 4563.0/5120.0*t25;
1893  double t116 = 43875.0/2048.0*t28;
1894  double t117 = 675.0/256.0*t30;
1895  double t118 = 50625.0/2048.0*t32;
1896  double t119 = 1521.0/512.0*t34;
1897  double t120 = 1125.0/128.0*t36;
1898  double t121 = 84375.0/1024.0*t38;
1899  double t122 = 7605.0/2048.0*t40;
1900  double t123 = 365625.0/4096.0*t42;
1901  double t124 = 5625.0/512.0*t44;
1902  double t125 = 421875.0/4096.0*t46;
1903  double t126 = 243.0/140.0*t63;
1904  double t127 = 243.0/70.0*t68;
1905  double t128 = 27.0/1024.0+t115+t116-t117-t118+t119-t120-t121-t122-t123+t124+t125+t126-t127;
1906  double t130 = 95073.0/8960.0*xi;
1907  double t131 = 32913.0/1792.0*t3;
1908  double t132 = 4823.0/960.0*t17;
1909  double t133 = 525.0/64.0*t21;
1910  double t134 = 4823.0/640.0*t25;
1911  double t135 = 1575.0/128.0*t30;
1912  double t136 = 4823.0/192.0*t34;
1913  double t137 = 2625.0/64.0*t36;
1914  double t138 = 24115.0/768.0*t40;
1915  double t139 = 13125.0/256.0*t44;
1916  double t140 = 162.0/35.0*t68;
1917  double t141 = 81.0/35.0*t63;
1918  double t142 = -t130+t131+t132-t133-t134+t135-t136+t137+t138-t139+t140-t141;
1919  double t143 = t101-t102-t103+t104+t105+t106-t107-t108-t109-t110+t111+t112-t113;
1920  double t144 = -27.0/1024.0+t115-t116-t117+t118+t119-t120+t121-t122+t123+t124-t125+t126-t127;
1921  double t146 = 27.0/5120.0-t72+t73+t74-t75-t76-t77+t78+t79+t80-t81-t82+t83;
1922  double t147 = t85+t86+t87-t88-t89+t90-t91-t92-t93-t94+t95+t96-t97+t98;
1923  double t149 = t2-t4-t6+t8+t10+t11-t12-t14-t16+t18+t20-t22-t24-t26;
1924  double t150 = -339.0/327680.0-t29+t31+t33-t35+t37+t39+t41+t43-t45-t47-t53+t58+t64-t69;
1925  double t152 = 225.0/512.0*t1;
1926  double t153 = 375.0/256.0*t3;
1927  double t154 = 1875.0/1024.0*t5;
1928  double t155 = 81.0/512.0*t7;
1929  double t156 = 101007.0/17920.0*t9;
1930  double t157 = 75.0/256.0*xi;
1931  double t158 = 291519.0/89600.0*eta;
1932  double t159 = 243.0/1024.0*t13;
1933  double t160 = 1755.0/256.0*t15;
1934  double t161 = 23031.0/6400.0*t19;
1935  double t162 = 7677.0/640.0*t21;
1936  double t163 = 27.0/5120.0-t152+t153+t154-t155-t156-t157+t158+t159+t160-t81-t161+t162;
1937  double t164 = 7677.0/512.0*t23;
1938  double t165 = 1125.0/128.0*t25;
1939  double t166 = 5625.0/128.0*t30;
1940  double t167 = 28125.0/512.0*t32;
1941  double t168 = 585.0/128.0*t34;
1942  double t169 = 14625.0/512.0*t38;
1943  double t170 = 3375.0/256.0*t40;
1944  double t171 = 10125.0/512.0*t42;
1945  double t172 = 16875.0/256.0*t44;
1946  double t173 = 243.0/350.0*t52;
1947  double t174 = 243.0/175.0*t57;
1948  double t175 = t164+t165+t87-t166-t167+t168-t91-t169-t170-t171+t172+t96+t173-t174;
1949  double t177 = 1125.0/512.0*t1;
1950  double t178 = 1875.0/256.0*t3;
1951  double t179 = 9375.0/1024.0*t5;
1952  double t180 = 1053.0/4096.0*t7;
1953  double t181 = 49653.0/3584.0*t9;
1954  double t182 = 375.0/256.0*xi;
1955  double t183 = 574533.0/71680.0*eta;
1956  double t184 = 1215.0/4096.0*t13;
1957  double t185 = 675.0/256.0*t15;
1958  double t186 = 4563.0/5120.0*t19;
1959  double t187 = 1521.0/512.0*t21;
1960  double t188 = 7605.0/2048.0*t23;
1961  double t189 = t177-t178-t179+t180+t181+t182-t183-t184+t185-t110-t186+t187+t188;
1962  double t190 = 14625.0/1024.0*t25;
1963  double t191 = 73125.0/1024.0*t30;
1964  double t192 = 365625.0/4096.0*t32;
1965  double t193 = 225.0/128.0*t34;
1966  double t194 = 5625.0/512.0*t38;
1967  double t195 = 16875.0/1024.0*t40;
1968  double t196 = 50625.0/2048.0*t42;
1969  double t197 = 84375.0/1024.0*t44;
1970  double t198 = 243.0/140.0*t52;
1971  double t199 = 243.0/70.0*t57;
1972  double t200 = -27.0/1024.0-t190-t116+t191+t192+t193-t120-t194+t195+t196-t197-t125-t198+t199;
1973  double t202 = 95073.0/8960.0*eta;
1974  double t203 = 4823.0/640.0*t19;
1975  double t204 = 4823.0/192.0*t21;
1976  double t205 = 24115.0/768.0*t23;
1977  double t206 = 32913.0/1792.0*t9;
1978  double t207 = 525.0/64.0*t34;
1979  double t208 = 1575.0/128.0*t15;
1980  double t209 = 13125.0/256.0*t38;
1981  double t210 = 162.0/35.0*t57;
1982  double t211 = 81.0/35.0*t52;
1983  double t212 = t202+t132+t203-t204-t205-t206-t207-t208+t137+t209-t210+t211;
1984  double t213 = -t177+t178+t179-t180+t181-t182-t183+t184+t185-t110-t186+t187+t188;
1985  double t214 = 27.0/1024.0+t190+t116-t191-t192+t193-t120-t194-t195-t196+t197+t125-t198+t199;
1986  double t216 = -27.0/5120.0+t152-t153-t154+t155-t156+t157+t158-t159+t160-t81-t161+t162;
1987  double t217 = t164-t165-t87+t166+t167+t168-t91-t169+t170+t171-t172-t96+t173-t174;
1988  double t219 = -t2+t4+t6-t8+t10-t11-t12+t14-t16+t18+t20-t22-t24+t26;
1989  double t220 = 339.0/327680.0+t29-t31-t33-t35+t37+t39-t41-t43+t45+t47-t53+t58-t64+t69;
1990  double t222 = -27.0/5120.0+t72-t73-t74+t75-t76+t77+t78-t79+t80-t81-t82+t83;
1991  double t223 = t85-t86-t87+t88+t89+t90-t91-t92+t93+t94-t95-t96+t97-t98;
1992  double t225 = -t101+t102+t103-t104+t105-t106-t107+t108-t109-t110+t111+t112-t113;
1993  double t226 = 27.0/1024.0-t115+t116+t117-t118+t119-t120+t121+t122-t123-t124+t125-t126+t127;
1994  double t228 = t130-t131+t132-t133+t134-t135-t136+t137-t138+t139-t140+t141;
1995  double t229 = t101+t102-t103+t104-t105-t106+t107-t108+t109-t110-t111+t112+t113;
1996  double t230 = -27.0/1024.0-t115-t116+t117+t118+t119-t120-t121+t122+t123-t124-t125-t126+t127;
1997  double t232 = 27.0/5120.0-t72-t73+t74-t75+t76+t77-t78+t79-t80-t81+t82+t83;
1998  double t233 = -t85-t86+t87+t88-t89+t90-t91+t92+t93-t94-t95+t96+t97-t98;
1999  double t235 = t2+t4-t6+t8-t10-t11+t12-t14+t16+t18-t20-t22+t24+t26;
2000  double t236 = -339.0/327680.0-t29-t31+t33-t35+t37-t39-t41+t43+t45-t47+t53-t58-t64+t69;
2001  double t238 = 27.0/5120.0-t152-t153+t154-t155+t156+t157-t158+t159-t160-t81+t161+t162;
2002  double t239 = -t164-t165+t87+t166-t167+t168-t91+t169+t170-t171-t172+t96-t173+t174;
2003  double t241 = t177+t178-t179+t180-t181-t182+t183-t184-t185-t110+t186+t187-t188;
2004  double t242 = -27.0/1024.0+t190-t116-t191+t192+t193-t120+t194-t195+t196+t197-t125+t198-t199;
2005  double t244 = -t202+t132-t203-t204+t205+t206-t207+t208+t137-t209+t210-t211;
2006  double t245 = -t177-t178+t179-t180-t181+t182+t183+t184-t185-t110+t186+t187-t188;
2007  double t246 = 27.0/1024.0-t190+t116+t191-t192+t193-t120+t194+t195-t196-t197+t125+t198-t199;
2008  double t248 = -27.0/5120.0+t152+t153-t154+t155+t156-t157-t158-t159-t160-t81+t161+t162;
2009  double t249 = -t164+t165-t87-t166+t167+t168-t91+t169-t170+t171+t172-t96-t173+t174;
2010  double t251 = 203125.0/98304.0*t1;
2011  double t252 = 78125.0/49152.0*t3;
2012  double t253 = 1953125.0/589824.0*t5;
2013  double t254 = 203125.0/98304.0*t7;
2014  double t255 = 78125.0/49152.0*t9;
2015  double t256 = 40625.0/49152.0*xi;
2016  double t257 = 40625.0/49152.0*eta;
2017  double t258 = 1953125.0/589824.0*t13;
2018  double t259 = 1015625.0/8192.0*t15;
2019  double t260 = 105625.0/4096.0*t17;
2020  double t261 = 528125.0/8192.0*t19;
2021  double t262 = 15625.0/589824.0-t251-t252+t253-t254-t255+t256+t257+t258+t259+t260-t261;
2022  double t263 = 203125.0/4096.0*t21;
2023  double t264 = 5078125.0/49152.0*t23;
2024  double t265 = 528125.0/8192.0*t25;
2025  double t266 = 2640625.0/16384.0*t28;
2026  double t267 = 1015625.0/8192.0*t30;
2027  double t268 = 25390625.0/98304.0*t32;
2028  double t269 = 203125.0/4096.0*t34;
2029  double t270 = 390625.0/4096.0*t36;
2030  double t271 = 9765625.0/49152.0*t38;
2031  double t272 = 5078125.0/49152.0*t40;
2032  double t273 = 25390625.0/98304.0*t42;
2033  double t274 = 9765625.0/49152.0*t44;
2034  double t275 = 244140625.0/589824.0*t46;
2035  double t276 = -t263+t264-t265+t266+t267-t268-t269+t270-t271+t272-t273-t274+t275;
2036  double t278 = 609375.0/16384.0*t1;
2037  double t279 = 234375.0/8192.0*t3;
2038  double t280 = 1953125.0/32768.0*t5;
2039  double t281 = 265625.0/49152.0*t7;
2040  double t282 = 78125.0/73728.0*t9;
2041  double t283 = 121875.0/8192.0*xi;
2042  double t284 = 53125.0/73728.0*eta;
2043  double t285 = 1953125.0/294912.0*t13;
2044  double t286 = 1015625.0/12288.0*t15;
2045  double t287 = 138125.0/6144.0*t17;
2046  double t288 = 690625.0/12288.0*t19;
2047  double t289 = -15625.0/32768.0+t278+t279-t280+t281+t282-t283-t284-t285-t286-t287+t288;
2048  double t290 = 265625.0/6144.0*t21;
2049  double t291 = 6640625.0/73728.0*t23;
2050  double t292 = 690625.0/4096.0*t25;
2051  double t293 = 3453125.0/8192.0*t28;
2052  double t294 = 1328125.0/4096.0*t30;
2053  double t295 = 33203125.0/49152.0*t32;
2054  double t296 = 203125.0/6144.0*t34;
2055  double t297 = 390625.0/6144.0*t36;
2056  double t298 = 9765625.0/73728.0*t38;
2057  double t299 = 5078125.0/24576.0*t40;
2058  double t300 = 25390625.0/49152.0*t42;
2059  double t301 = 9765625.0/24576.0*t44;
2060  double t302 = 244140625.0/294912.0*t46;
2061  double t303 = t290-t291+t292-t293-t294+t295+t296-t297+t298-t299+t300+t301-t302;
2062  double t305 = 15625.0/32768.0-t278-t279+t280-t281+t282+t283-t284+t285-t286-t287+t288;
2063  double t306 = t290-t291-t292+t293+t294-t295+t296-t297+t298+t299-t300-t301+t302;
2064  double t308 = t251+t252-t253+t254-t255-t256+t257-t258+t259+t260-t261-t263;
2065  double t309 = -15625.0/589824.0+t264+t265-t266-t267+t268-t269+t270-t271-t272+t273+t274-t275;
2066  double t311 = 265625.0/49152.0*t1;
2067  double t312 = 78125.0/73728.0*t3;
2068  double t313 = 1953125.0/294912.0*t5;
2069  double t314 = 609375.0/16384.0*t7;
2070  double t315 = 234375.0/8192.0*t9;
2071  double t316 = 53125.0/73728.0*xi;
2072  double t317 = 121875.0/8192.0*eta;
2073  double t318 = 1953125.0/32768.0*t13;
2074  double t319 = 1328125.0/4096.0*t15;
2075  double t320 = 690625.0/4096.0*t19;
2076  double t321 = -15625.0/32768.0+t311+t312-t313+t314+t315-t316-t317-t318-t319-t287+t320;
2077  double t322 = 203125.0/6144.0*t21;
2078  double t323 = 5078125.0/24576.0*t23;
2079  double t324 = 690625.0/12288.0*t25;
2080  double t325 = 1015625.0/12288.0*t30;
2081  double t326 = 25390625.0/49152.0*t32;
2082  double t327 = 265625.0/6144.0*t34;
2083  double t328 = 9765625.0/24576.0*t38;
2084  double t329 = 6640625.0/73728.0*t40;
2085  double t330 = 33203125.0/49152.0*t42;
2086  double t331 = 9765625.0/73728.0*t44;
2087  double t332 = t322-t323+t324-t293-t325+t326+t327-t297+t328-t329+t330+t331-t302;
2088  double t334 = 796875.0/8192.0*t1;
2089  double t335 = 78125.0/4096.0*t3;
2090  double t336 = 1953125.0/16384.0*t5;
2091  double t337 = 796875.0/8192.0*t7;
2092  double t338 = 78125.0/4096.0*t9;
2093  double t339 = 53125.0/4096.0*xi;
2094  double t340 = 53125.0/4096.0*eta;
2095  double t341 = 1953125.0/16384.0*t13;
2096  double t342 = 1328125.0/6144.0*t15;
2097  double t343 = 180625.0/9216.0*t17;
2098  double t344 = 903125.0/6144.0*t19;
2099  double t345 = 140625.0/16384.0-t334-t335+t336-t337-t338+t339+t340+t341+t342+t343-t344;
2100  double t346 = 265625.0/9216.0*t21;
2101  double t347 = 6640625.0/36864.0*t23;
2102  double t348 = 903125.0/6144.0*t25;
2103  double t349 = 4515625.0/4096.0*t28;
2104  double t350 = 1328125.0/6144.0*t30;
2105  double t351 = 33203125.0/24576.0*t32;
2106  double t352 = 265625.0/9216.0*t34;
2107  double t353 = 390625.0/9216.0*t36;
2108  double t354 = 9765625.0/36864.0*t38;
2109  double t355 = 6640625.0/36864.0*t40;
2110  double t356 = 33203125.0/24576.0*t42;
2111  double t357 = 9765625.0/36864.0*t44;
2112  double t358 = 244140625.0/147456.0*t46;
2113  double t359 = -t346+t347-t348+t349+t350-t351-t352+t353-t354+t355-t356-t357+t358;
2114  double t361 = t334+t335-t336+t337-t338-t339+t340-t341+t342+t343-t344-t346;
2115  double t362 = -140625.0/16384.0+t347+t348-t349-t350+t351-t352+t353-t354-t355+t356+t357-t358;
2116  double t364 = 15625.0/32768.0-t311-t312+t313-t314+t315+t316-t317+t318-t319-t287+t320;
2117  double t365 = t322-t323-t324+t293+t325-t326+t327-t297+t328+t329-t330-t331+t302;
2118  double t367 = 15625.0/32768.0-t311+t312+t313-t314-t315-t316+t317+t318+t319-t287-t320;
2119  double t368 = t322+t323+t324+t293-t325-t326+t327-t297-t328-t329-t330+t331+t302;
2120  double t370 = t334-t335-t336+t337+t338+t339-t340-t341-t342+t343+t344-t346;
2121  double t371 = -140625.0/16384.0-t347-t348-t349+t350+t351-t352+t353+t354+t355+t356-t357-t358;
2122  double t373 = 140625.0/16384.0-t334+t335+t336-t337+t338-t339-t340+t341-t342+t343+t344;
2123  double t374 = -t346-t347+t348+t349-t350-t351-t352+t353+t354-t355-t356+t357+t358;
2124  double t376 = -15625.0/32768.0+t311-t312-t313+t314-t315+t316+t317-t318+t319-t287-t320;
2125  double t377 = t322+t323-t324-t293+t325+t326+t327-t297-t328+t329+t330-t331-t302;
2126  double t379 = t251-t252-t253+t254+t255+t256-t257-t258-t259+t260+t261-t263;
2127  double t380 = -15625.0/589824.0-t264-t265-t266+t267+t268-t269+t270+t271+t272+t273-t274-t275;
2128  double t382 = 15625.0/32768.0-t278+t279+t280-t281-t282-t283+t284+t285+t286-t287-t288;
2129  double t383 = t290+t291+t292+t293-t294-t295+t296-t297-t298-t299-t300+t301+t302;
2130  double t385 = -15625.0/32768.0+t278-t279-t280+t281-t282+t283+t284-t285+t286-t287-t288;
2131  double t386 = t290+t291-t292-t293+t294+t295+t296-t297-t298+t299+t300-t301-t302;
2132  double t388 = 15625.0/589824.0-t251+t252+t253-t254+t255-t256-t257+t258-t259+t260+t261;
2133  double t389 = -t263-t264+t265+t266-t267-t268-t269+t270+t271-t272-t273+t274+t275;
2134 
2135  values[0] = t27+t70;
2136  values[1] = t84+t99;
2137  values[2] = t114+t128;
2138  values[3] = t142;
2139  values[4] = t143+t144;
2140  values[5] = t146+t147;
2141  values[6] = t149+t150;
2142  values[7] = t163+t175;
2143  values[8] = t189+t200;
2144  values[9] = t212;
2145  values[10] = t213+t214;
2146  values[11] = t216+t217;
2147  values[12] = t219+t220;
2148  values[13] = t222+t223;
2149  values[14] = t225+t226;
2150  values[15] = t228;
2151  values[16] = t229+t230;
2152  values[17] = t232+t233;
2153  values[18] = t235+t236;
2154  values[19] = t238+t239;
2155  values[20] = t241+t242;
2156  values[21] = t244;
2157  values[22] = t245+t246;
2158  values[23] = t248+t249;
2159  values[24] = t262+t276;
2160  values[25] = t289+t303;
2161  values[26] = t305+t306;
2162  values[27] = t308+t309;
2163  values[28] = t321+t332;
2164  values[29] = t345+t359;
2165  values[30] = t361+t362;
2166  values[31] = t364+t365;
2167  values[32] = t367+t368;
2168  values[33] = t370+t371;
2169  values[34] = t373+t374;
2170  values[35] = t376+t377;
2171  values[36] = t379+t380;
2172  values[37] = t382+t383;
2173  values[38] = t385+t386;
2174  values[39] = t388+t389;
2175 }
2176 
2177 // values of the derivatives in eta-eta direction
2178 static void C_Q_UL6SE_2D_DeriveEtaEta(double xi, double eta, double *values)
2179 {
2180  double t1 = 1.0-xi;
2181  double t3 = eta*eta;
2182  double t4 = t3*eta;
2183  double t7 = 35.0/2.0*t4-15.0/2.0*eta;
2184  double t8 = t1*eta*t7;
2185  double t9 = 81.0/175.0*t8;
2186  double t10 = xi*xi;
2187  double t11 = t10*t10;
2188  double t12 = t11*xi;
2189  double t13 = 538871.0/737280.0*t12;
2190  double t14 = 335499.0/114688.0*t3;
2191  double t15 = 1.0-t3;
2192  double t18 = 105.0/2.0*t3-15.0/2.0;
2193  double t19 = t1*t15*t18;
2194  double t20 = 81.0/700.0*t19;
2195  double t21 = t10*t4;
2196  double t22 = 538871.0/73728.0*t21;
2197  double t23 = xi*eta;
2198  double t24 = 1873.0/16384.0*t23;
2199  double t25 = t10*eta;
2200  double t26 = 727871.0/614400.0*t25;
2201  double t27 = t10*xi;
2202  double t28 = t27*eta;
2203  double t29 = 38375.0/24576.0*t28;
2204  double t30 = t11*eta;
2205  double t31 = 138955.0/49152.0*t30;
2206  double t32 = t12*eta;
2207  double t33 = 154075.0/49152.0*t32;
2208  double t34 = xi*t3;
2209  double t35 = 1875847.0/573440.0*t34;
2210  double t36 = t10*t3;
2211  double t37 = 463271.0/122880.0*t36;
2212  double t38 = t27*t3;
2213  double t39 = 138955.0/24576.0*t38;
2214  double t40 = t9+t13+t14-t20+t22-t24-t26+t29+t31-t33-t35-t37+t39;
2215  double t41 = t11*t3;
2216  double t42 = 430175.0/49152.0*t41;
2217  double t43 = t12*t3;
2218  double t44 = 505775.0/49152.0*t43;
2219  double t45 = xi*t4;
2220  double t46 = 42289.0/49152.0*t45;
2221  double t47 = t27*t4;
2222  double t48 = 770375.0/73728.0*t47;
2223  double t49 = t11*t4;
2224  double t50 = 2528875.0/147456.0*t49;
2225  double t51 = t12*t4;
2226  double t52 = 2906875.0/147456.0*t51;
2227  double t53 = t3*t3;
2228  double t56 = 3.0/8.0+35.0/8.0*t53-15.0/4.0*t3;
2229  double t57 = t1*t56;
2230  double t58 = 81.0/350.0*t57;
2231  double t59 = 4651.0/16384.0*t4;
2232  double t60 = 496367.0/1843200.0*t10;
2233  double t61 = 727871.0/1843200.0*t27;
2234  double t62 = 463271.0/737280.0*t11;
2235  double t63 = 4891279.0/8601600.0*xi;
2236  double t64 = 18489.0/409600.0*eta;
2237  double t65 = -1564309.0/2867200.0+t42-t44+t46-t48-t50+t52+t58-t59+t60-t61-t62+t63+t64;
2238  double t67 = 2187.0/12800.0*eta;
2239  double t68 = 2187.0/2560.0*t3;
2240  double t69 = 729.0/512.0*t4;
2241  double t70 = 23031.0/1280.0*t36;
2242  double t71 = 75.0/256.0*xi;
2243  double t72 = 7677.0/6400.0*t10;
2244  double t73 = 375.0/128.0*t27;
2245  double t74 = 585.0/512.0*t11;
2246  double t75 = 675.0/256.0*t12;
2247  double t76 = 225.0/256.0*t23;
2248  double t77 = 23031.0/6400.0*t25;
2249  double t78 = 1125.0/128.0*t28;
2250  double t79 = -t67-t68+t69+t70-t71-t72+t73+t74-t75+t76+t77-t78;
2251  double t80 = 1755.0/512.0*t30;
2252  double t81 = 2025.0/256.0*t32;
2253  double t82 = 1125.0/256.0*t34;
2254  double t83 = 5625.0/128.0*t38;
2255  double t84 = 8775.0/512.0*t41;
2256  double t85 = 10125.0/256.0*t43;
2257  double t86 = 1875.0/256.0*t45;
2258  double t87 = 7677.0/256.0*t21;
2259  double t88 = 9375.0/128.0*t47;
2260  double t89 = 14625.0/512.0*t49;
2261  double t90 = 16875.0/256.0*t51;
2262  double t91 = 729.0/12800.0-t80+t81+t82-t83-t84+t85-t86-t87+t88+t89-t90;
2263  double t93 = 2187.0/5120.0*eta;
2264  double t94 = 2187.0/1024.0*t3;
2265  double t95 = 3645.0/1024.0*t4;
2266  double t96 = 4563.0/1024.0*t36;
2267  double t97 = 375.0/256.0*xi;
2268  double t98 = 1521.0/5120.0*t10;
2269  double t99 = 4875.0/1024.0*t27;
2270  double t100 = 225.0/512.0*t11;
2271  double t101 = 3375.0/1024.0*t12;
2272  double t102 = 1125.0/256.0*t23;
2273  double t103 = 4563.0/5120.0*t25;
2274  double t104 = 14625.0/1024.0*t28;
2275  double t105 = t93+t94-t95+t96+t97-t98-t99+t100+t101-t102+t103+t104;
2276  double t106 = 675.0/512.0*t30;
2277  double t107 = 10125.0/1024.0*t32;
2278  double t108 = 5625.0/256.0*t34;
2279  double t109 = 73125.0/1024.0*t38;
2280  double t110 = 3375.0/512.0*t41;
2281  double t111 = 50625.0/1024.0*t43;
2282  double t112 = 9375.0/256.0*t45;
2283  double t113 = 7605.0/1024.0*t21;
2284  double t114 = 121875.0/1024.0*t47;
2285  double t115 = 5625.0/512.0*t49;
2286  double t116 = 84375.0/1024.0*t51;
2287  double t117 = -729.0/5120.0-t106-t107-t108+t109-t110-t111+t112-t113-t114+t115+t116;
2288  double t119 = 4823.0/1920.0*t10;
2289  double t120 = 525.0/256.0*t11;
2290  double t121 = 4823.0/640.0*t25;
2291  double t122 = 1575.0/256.0*t30;
2292  double t123 = 1771.0/1280.0*eta;
2293  double t124 = 1771.0/256.0*t3;
2294  double t125 = 8855.0/768.0*t4;
2295  double t126 = 4823.0/128.0*t36;
2296  double t127 = 7875.0/256.0*t41;
2297  double t128 = 24115.0/384.0*t21;
2298  double t129 = 13125.0/256.0*t49;
2299  double t130 = t119-t120-t121+t122-1771.0/3840.0+t123+t124-t125-t126+t127+t128-t129;
2300  double t131 = t93+t94-t95+t96-t97-t98+t99+t100-t101+t102+t103-t104;
2301  double t132 = -729.0/5120.0-t106+t107+t108-t109-t110+t111-t112-t113+t114+t115-t116;
2302  double t134 = -t67-t68+t69+t70+t71-t72-t73+t74+t75-t76+t77+t78;
2303  double t135 = 729.0/12800.0-t80-t81-t82+t83-t84-t85+t86-t87-t88+t89+t90;
2304  double t137 = -t13+t14+t22+t24-t26-t29+t31+t33+t35-t37-t39+t42+t44;
2305  double t138 = 1.0+xi;
2306  double t140 = t138*t15*t18;
2307  double t141 = 81.0/700.0*t140;
2308  double t143 = t138*eta*t7;
2309  double t144 = 81.0/175.0*t143;
2310  double t145 = t138*t56;
2311  double t146 = 81.0/350.0*t145;
2312  double t147 = -1564309.0/2867200.0-t46+t48-t50-t52-t141+t144-t59+t60+t61-t62-t63+t64+t146;
2313  double t149 = 7677.0/2560.0*t12;
2314  double t150 = 303021.0/17920.0*t3;
2315  double t151 = 3375.0/128.0*t21;
2316  double t152 = 81.0/256.0*t23;
2317  double t153 = 1125.0/128.0*t25;
2318  double t154 = 5625.0/256.0*t30;
2319  double t155 = 5625.0/256.0*t32;
2320  double t156 = 303021.0/17920.0*t34;
2321  double t157 = 1755.0/256.0*t36;
2322  double t158 = 1755.0/256.0*t38;
2323  double t159 = 8775.0/512.0*t43;
2324  double t160 = t149-t150-t151-t152+t153+t78-t154-t155-t156+t157+t158-t84-t159;
2325  double t161 = 243.0/256.0*t45;
2326  double t162 = 3375.0/128.0*t47;
2327  double t163 = 16875.0/256.0*t49;
2328  double t164 = 243.0/350.0*t140;
2329  double t165 = 486.0/175.0*t143;
2330  double t166 = 243.0/256.0*t4;
2331  double t167 = 7677.0/6400.0*t27;
2332  double t168 = 7677.0/2560.0*t11;
2333  double t169 = 291519.0/89600.0*xi;
2334  double t170 = 81.0/256.0*eta;
2335  double t171 = 243.0/175.0*t145;
2336  double t172 = 291519.0/89600.0+t161-t162+t163+t90+t164-t165+t166-t72-t167+t168+t169-t170-t171;
2337  double t174 = 1521.0/2048.0*t12;
2338  double t175 = 148959.0/3584.0*t3;
2339  double t176 = 16875.0/512.0*t21;
2340  double t177 = 1053.0/2048.0*t23;
2341  double t178 = 14625.0/1024.0*t25;
2342  double t179 = 73125.0/2048.0*t30;
2343  double t180 = 73125.0/2048.0*t32;
2344  double t181 = 148959.0/3584.0*t34;
2345  double t182 = 675.0/256.0*t36;
2346  double t183 = 675.0/256.0*t38;
2347  double t184 = 3375.0/512.0*t43;
2348  double t185 = t174+t175+t176+t177-t178-t104+t179+t180+t181+t182+t183-t110-t184;
2349  double t186 = 1215.0/1024.0*t45;
2350  double t187 = 16875.0/512.0*t47;
2351  double t188 = 84375.0/1024.0*t49;
2352  double t189 = 243.0/140.0*t140;
2353  double t190 = 243.0/35.0*t143;
2354  double t191 = 1215.0/1024.0*t4;
2355  double t192 = 1521.0/5120.0*t27;
2356  double t193 = 1521.0/2048.0*t11;
2357  double t194 = 574533.0/71680.0*xi;
2358  double t195 = 1053.0/2048.0*eta;
2359  double t196 = 243.0/70.0*t145;
2360  double t197 = -574533.0/71680.0-t186+t187-t188-t116-t189+t190-t191-t98-t192+t193-t194+t195+t196;
2361  double t199 = 95073.0/8960.0*xi;
2362  double t200 = 4823.0/1920.0*t27;
2363  double t201 = 4823.0/768.0*t11;
2364  double t202 = 4823.0/768.0*t12;
2365  double t203 = 98739.0/1792.0*t3;
2366  double t204 = 98739.0/1792.0*t34;
2367  double t205 = 1575.0/128.0*t36;
2368  double t206 = 1575.0/128.0*t38;
2369  double t207 = 7875.0/256.0*t43;
2370  double t211 = t199+t119+t200-t201-t202+95073.0/8960.0-t203-t204-t205-t206+t127+t207-162.0/35.0*t145-324.0/35.0*t143+81.0/35.0*t140;
2371  double t212 = t174+t175-t176-t177+t178+t104-t179-t180+t181+t182+t183-t110-t184;
2372  double t213 = -574533.0/71680.0+t186-t187+t188+t116-t189+t190+t191-t98-t192+t193-t194-t195+t196;
2373  double t215 = t149-t150+t151+t152-t153-t78+t154+t155-t156+t157+t158-t84-t159;
2374  double t216 = 291519.0/89600.0-t161+t162-t163-t90+t164-t165-t166-t72-t167+t168+t169+t170-t171;
2375  double t218 = -t13+t14-t22-t24+t26+t29-t31-t33+t35-t37-t39+t42+t44;
2376  double t219 = -1564309.0/2867200.0+t46-t48+t50+t52-t141+t144+t59+t60+t61-t62-t63-t64+t146;
2377  double t221 = t67-t68-t69+t70+t71-t72-t73+t74+t75+t76-t77-t78;
2378  double t222 = 729.0/12800.0+t80+t81-t82+t83-t84-t85-t86+t87+t88-t89-t90;
2379  double t224 = -t93+t94+t95+t96-t97-t98+t99+t100-t101-t102-t103+t104;
2380  double t225 = -729.0/5120.0+t106-t107+t108-t109-t110+t111+t112+t113-t114-t115+t116;
2381  double t227 = t119-t120+t121-t122-1771.0/3840.0-t123+t124+t125-t126+t127-t128+t129;
2382  double t228 = -t93+t94+t95+t96+t97-t98-t99+t100+t101+t102-t103-t104;
2383  double t229 = -729.0/5120.0+t106+t107-t108+t109-t110-t111-t112+t113+t114-t115-t116;
2384  double t231 = t67-t68-t69+t70-t71-t72+t73+t74-t75-t76-t77+t78;
2385  double t232 = 729.0/12800.0+t80-t81+t82-t83-t84+t85+t86+t87-t88-t89+t90;
2386  double t234 = t9+t13+t14-t20-t22+t24+t26-t29-t31+t33-t35-t37+t39;
2387  double t235 = -1564309.0/2867200.0+t42-t44-t46+t48+t50-t52+t58+t59+t60-t61-t62+t63-t64;
2388  double t237 = 486.0/175.0*t8;
2389  double t238 = 243.0/350.0*t19;
2390  double t239 = -t237-t149-t150+t238+t151-t152-t153+t78+t154-t155+t156+t157-t158;
2391  double t240 = 243.0/175.0*t57;
2392  double t241 = 291519.0/89600.0-t84+t159+t161-t162-t163+t90-t240-t166-t72+t167+t168-t169+t170;
2393  double t243 = 243.0/35.0*t8;
2394  double t244 = 243.0/140.0*t19;
2395  double t245 = t243-t174+t175-t244-t176+t177+t178-t104-t179+t180-t181+t182-t183;
2396  double t246 = 243.0/70.0*t57;
2397  double t247 = -574533.0/71680.0-t110+t184-t186+t187+t188-t116+t246+t191-t98+t192+t193+t194-t195;
2398  double t252 = -t199+t119-t200-t201+t202+95073.0/8960.0-t203+t204-t205+t206+t127-t207-162.0/35.0*t57-324.0/35.0*t8+81.0/35.0*t19;
2399  double t253 = t243-t174+t175-t244+t176-t177-t178+t104+t179-t180-t181+t182-t183;
2400  double t254 = -574533.0/71680.0-t110+t184+t186-t187-t188+t116+t246-t191-t98+t192+t193+t194+t195;
2401  double t256 = -t237-t149-t150+t238-t151+t152+t153-t78-t154+t155+t156+t157-t158;
2402  double t257 = 291519.0/89600.0-t84+t159-t161+t162+t163-t90-t240+t166-t72+t167+t168-t169-t170;
2403  double t259 = 40625.0/16384.0*eta;
2404  double t260 = 46875.0/16384.0*t3;
2405  double t261 = 390625.0/49152.0*t4;
2406  double t262 = 609375.0/8192.0*t36;
2407  double t263 = 40625.0/49152.0*xi;
2408  double t264 = 105625.0/8192.0*t10;
2409  double t265 = 528125.0/24576.0*t27;
2410  double t266 = 203125.0/16384.0*t11;
2411  double t267 = 1015625.0/49152.0*t12;
2412  double t268 = 203125.0/49152.0*t23;
2413  double t269 = 528125.0/8192.0*t25;
2414  double t270 = 2640625.0/24576.0*t28;
2415  double t271 = t259+t260-t261-t262+t263+t264-t265-t266+t267-t268-t269+t270;
2416  double t272 = 1015625.0/16384.0*t30;
2417  double t273 = 5078125.0/49152.0*t32;
2418  double t274 = 78125.0/16384.0*t34;
2419  double t275 = 1015625.0/8192.0*t38;
2420  double t276 = 1171875.0/16384.0*t41;
2421  double t277 = 1953125.0/16384.0*t43;
2422  double t278 = 1953125.0/147456.0*t45;
2423  double t279 = 5078125.0/24576.0*t21;
2424  double t280 = 25390625.0/73728.0*t47;
2425  double t281 = 9765625.0/49152.0*t49;
2426  double t282 = 48828125.0/147456.0*t51;
2427  double t283 = -8125.0/16384.0+t272-t273-t274+t275+t276-t277+t278+t279-t280-t281+t282;
2428  double t285 = 53125.0/8192.0*eta;
2429  double t286 = 15625.0/8192.0*t3;
2430  double t287 = 390625.0/24576.0*t4;
2431  double t288 = 203125.0/4096.0*t36;
2432  double t289 = 53125.0/73728.0*xi;
2433  double t290 = 138125.0/12288.0*t10;
2434  double t291 = 690625.0/36864.0*t27;
2435  double t292 = 265625.0/24576.0*t11;
2436  double t293 = 1328125.0/73728.0*t12;
2437  double t294 = 265625.0/24576.0*t23;
2438  double t295 = 690625.0/4096.0*t25;
2439  double t296 = 3453125.0/12288.0*t28;
2440  double t297 = -t285-t286+t287+t288-t289-t290+t291+t292-t293+t294+t295-t296;
2441  double t298 = 1328125.0/8192.0*t30;
2442  double t299 = 6640625.0/24576.0*t32;
2443  double t300 = 78125.0/24576.0*t34;
2444  double t301 = 1015625.0/12288.0*t38;
2445  double t302 = 390625.0/8192.0*t41;
2446  double t303 = 1953125.0/24576.0*t43;
2447  double t304 = 1953125.0/73728.0*t45;
2448  double t305 = 5078125.0/12288.0*t21;
2449  double t306 = 25390625.0/36864.0*t47;
2450  double t307 = 9765625.0/24576.0*t49;
2451  double t308 = 48828125.0/73728.0*t51;
2452  double t309 = 10625.0/24576.0-t298+t299+t300-t301-t302+t303-t304-t305+t306+t307-t308;
2453  double t311 = t285-t286-t287+t288-t289-t290+t291+t292-t293-t294-t295+t296;
2454  double t312 = 10625.0/24576.0+t298-t299+t300-t301-t302+t303+t304+t305-t306-t307+t308;
2455  double t314 = -t259+t260+t261-t262+t263+t264-t265-t266+t267+t268+t269-t270;
2456  double t315 = -8125.0/16384.0-t272+t273-t274+t275+t276-t277-t278-t279+t280+t281-t282;
2457  double t317 = 121875.0/8192.0*eta;
2458  double t318 = 140625.0/8192.0*t3;
2459  double t319 = 390625.0/8192.0*t4;
2460  double t320 = 265625.0/4096.0*t36;
2461  double t321 = 121875.0/8192.0*xi;
2462  double t322 = 690625.0/12288.0*t27;
2463  double t323 = 203125.0/24576.0*t11;
2464  double t324 = 1015625.0/24576.0*t12;
2465  double t325 = 609375.0/8192.0*t23;
2466  double t326 = 690625.0/12288.0*t25;
2467  double t327 = -t317-t318+t319+t320-t321-t290+t322+t323-t324+t325+t326-t296;
2468  double t328 = 1015625.0/24576.0*t30;
2469  double t329 = 5078125.0/24576.0*t32;
2470  double t330 = 703125.0/8192.0*t34;
2471  double t331 = 1328125.0/4096.0*t38;
2472  double t332 = 1953125.0/8192.0*t43;
2473  double t333 = 1953125.0/8192.0*t45;
2474  double t334 = 6640625.0/36864.0*t21;
2475  double t335 = 33203125.0/36864.0*t47;
2476  double t336 = 9765625.0/73728.0*t49;
2477  double t337 = 24375.0/8192.0-t328+t329+t330-t331-t302+t332-t333-t334+t335+t336-t308;
2478  double t339 = 159375.0/4096.0*eta;
2479  double t340 = 46875.0/4096.0*t3;
2480  double t341 = 390625.0/4096.0*t4;
2481  double t342 = 265625.0/6144.0*t36;
2482  double t343 = 53125.0/4096.0*xi;
2483  double t344 = 180625.0/18432.0*t10;
2484  double t345 = 903125.0/18432.0*t27;
2485  double t346 = 265625.0/36864.0*t11;
2486  double t347 = 1328125.0/36864.0*t12;
2487  double t348 = 796875.0/4096.0*t23;
2488  double t349 = 903125.0/6144.0*t25;
2489  double t350 = 4515625.0/6144.0*t28;
2490  double t351 = t339+t340-t341-t342+t343+t344-t345-t346+t347-t348-t349+t350;
2491  double t352 = 1328125.0/12288.0*t30;
2492  double t353 = 6640625.0/12288.0*t32;
2493  double t354 = 234375.0/4096.0*t34;
2494  double t355 = 1328125.0/6144.0*t38;
2495  double t356 = 390625.0/12288.0*t41;
2496  double t357 = 1953125.0/12288.0*t43;
2497  double t358 = 1953125.0/4096.0*t45;
2498  double t359 = 6640625.0/18432.0*t21;
2499  double t360 = 33203125.0/18432.0*t47;
2500  double t361 = 9765625.0/36864.0*t49;
2501  double t362 = 48828125.0/36864.0*t51;
2502  double t363 = -10625.0/4096.0+t352-t353-t354+t355+t356-t357+t358+t359-t360-t361+t362;
2503  double t365 = -t339+t340+t341-t342+t343+t344-t345-t346+t347+t348+t349-t350;
2504  double t366 = -10625.0/4096.0-t352+t353-t354+t355+t356-t357-t358-t359+t360+t361-t362;
2505  double t368 = t317-t318-t319+t320-t321-t290+t322+t323-t324-t325-t326+t296;
2506  double t369 = 24375.0/8192.0+t328-t329+t330-t331-t302+t332+t333+t334-t335-t336+t308;
2507  double t371 = -t317-t318+t319+t320+t321-t290-t322+t323+t324-t325+t326+t296;
2508  double t372 = 24375.0/8192.0-t328-t329-t330+t331-t302-t332+t333-t334-t335+t336+t308;
2509  double t374 = t339+t340-t341-t342-t343+t344+t345-t346-t347+t348-t349-t350;
2510  double t375 = -10625.0/4096.0+t352+t353+t354-t355+t356+t357-t358+t359+t360-t361-t362;
2511  double t377 = -t339+t340+t341-t342-t343+t344+t345-t346-t347-t348+t349+t350;
2512  double t378 = -10625.0/4096.0-t352-t353+t354-t355+t356+t357+t358-t359-t360+t361+t362;
2513  double t380 = t317-t318-t319+t320+t321-t290-t322+t323+t324+t325-t326-t296;
2514  double t381 = 24375.0/8192.0+t328+t329-t330+t331-t302-t332-t333+t334+t335-t336-t308;
2515  double t383 = t259+t260-t261-t262-t263+t264+t265-t266-t267+t268-t269-t270;
2516  double t384 = -8125.0/16384.0+t272+t273+t274-t275+t276+t277-t278+t279+t280-t281-t282;
2517  double t386 = -t285-t286+t287+t288+t289-t290-t291+t292+t293-t294+t295+t296;
2518  double t387 = 10625.0/24576.0-t298-t299-t300+t301-t302-t303+t304-t305-t306+t307+t308;
2519  double t389 = t285-t286-t287+t288+t289-t290-t291+t292+t293+t294-t295-t296;
2520  double t390 = 10625.0/24576.0+t298+t299-t300+t301-t302-t303-t304+t305+t306-t307-t308;
2521  double t392 = -t259+t260+t261-t262-t263+t264+t265-t266-t267-t268+t269+t270;
2522  double t393 = -8125.0/16384.0-t272-t273+t274-t275+t276+t277+t278-t279-t280+t281+t282;
2523 
2524  values[0] = t40+t65;
2525  values[1] = t79+t91;
2526  values[2] = t105+t117;
2527  values[3] = t130;
2528  values[4] = t131+t132;
2529  values[5] = t134+t135;
2530  values[6] = t137+t147;
2531  values[7] = t160+t172;
2532  values[8] = t185+t197;
2533  values[9] = t211;
2534  values[10] = t212+t213;
2535  values[11] = t215+t216;
2536  values[12] = t218+t219;
2537  values[13] = t221+t222;
2538  values[14] = t224+t225;
2539  values[15] = t227;
2540  values[16] = t228+t229;
2541  values[17] = t231+t232;
2542  values[18] = t234+t235;
2543  values[19] = t239+t241;
2544  values[20] = t245+t247;
2545  values[21] = t252;
2546  values[22] = t253+t254;
2547  values[23] = t256+t257;
2548  values[24] = t271+t283;
2549  values[25] = t297+t309;
2550  values[26] = t311+t312;
2551  values[27] = t314+t315;
2552  values[28] = t327+t337;
2553  values[29] = t351+t363;
2554  values[30] = t365+t366;
2555  values[31] = t368+t369;
2556  values[32] = t371+t372;
2557  values[33] = t374+t375;
2558  values[34] = t377+t378;
2559  values[35] = t380+t381;
2560  values[36] = t383+t384;
2561  values[37] = t386+t387;
2562  values[38] = t389+t390;
2563  values[39] = t392+t393;
2564 }
2565 
2566 // ***********************************************************************
2567 
2568 TBaseFunct2D *BF_C_Q_UL6SE_2D_Obj = new TBaseFunct2D
2569  (40, BF_C_Q_UL6SE_2D, BFUnitSquare,
2570  C_Q_UL6SE_2D_Funct, C_Q_UL6SE_2D_DeriveXi,
2571  C_Q_UL6SE_2D_DeriveEta, C_Q_UL6SE_2D_DeriveXiXi,
2572  C_Q_UL6SE_2D_DeriveXiEta, C_Q_UL6SE_2D_DeriveEtaEta, 6, 6,
2573  0, NULL);
TBaseFunct2D
Definition: BaseFunct2D.h:27
Generated on Tue May 10 2016 09:07:58 for ParMooN by   doxygen 1.8.5