ParMooN
 All Classes Functions Variables Friends Pages
BF_C_T_UL5_2D.h
1 // ***********************************************************************
2 // P5 element with bubbles, conforming, 2D
3 // ***********************************************************************
4 
5 // base function values
6 static void C_T_UL5_2D_Funct(double xi, double eta, double *values)
7 {
8  double t1, t2, t5, t7, t9, t11, t12, t13, t14, t16, t18, t19, t20, t22, t23, t25;
9  double t26, t29, t30, t32, t33, t34, t35, t36, t39, t42, t45, t47, t49;
10  double t51, t53, t59, t61, t64, t66, t68, t75, t77, t78, t80, t85, t87;
11  double t88, t90, t93, t100, t102, t103, t105, t110, t112, t113, t115, t119;
12  double t126, t128, t130, t136, t138, t142, t149, t151, t153, t159, t161;
13  double t164, t171, t174, t183, t185, t188, t192, t196, t203, t205, t207;
14  double t211, t215, t222, t224, t232, t240, t248, t256, t266, t276, t287;
15  double t289, t294, t297, t298, t309, t311, t316, t319, t320, t322, t331;
16  double t339, t349, t357, t378, t380, t384, t396, t398, t415, t417, t421;
17  double t432, t433, t435, t452, t453, t469, t471, t475, t486, t488, t493;
18  double t497, t498, t504, t505, t506, t507, t520, t536, t538, t542, t553;
19  double t554, t556, t561, t572, t574, t593, t609, t610, t626;
20 
21  t1 = xi*xi;
22  t2 = t1*xi;
23  t5 = eta*eta;
24  t7 = t1*t1;
25  t9 = t5*t5;
26  t11 = t7*xi;
27  t12 = 625.0/24.0*t11;
28  t13 = t9*eta;
29  t14 = 625.0/24.0*t13;
30  t16 = t5*eta;
31  t18 = t2*t5;
32  t19 = 1.0-xi-eta;
33  t20 = t18*t19;
34  t22 = t7*eta;
35  t23 = t22*t19;
36  t25 = xi*t9;
37  t26 = t25*t19;
38  t29 = t11*eta*t19;
39  t30 = 43927.0/12.0*t29;
40  t32 = xi*t13*t19;
41  t33 = 43927.0/12.0*t32;
42  t34 = 1.0-2125.0/24.0*t2-137.0/12.0*eta+375.0/8.0*t5+625.0/8.0*t7+625.0/8.0*t9-t12-t14+375.0/8.0*t1-2125.0/24.0*t16+154495.0/6.0*t20+53185.0/6.0*t23+53185.0/6.0*t26-t30-t33;
43  t35 = t1*t16;
44  t36 = t35*t19;
45  t39 = t7*t5*t19;
46  t42 = t1*t9*t19;
47  t45 = t2*t16*t19;
48  t47 = xi*t5;
49  t49 = t2*eta;
50  t51 = t1*t5;
51  t53 = xi*t16;
52  t59 = xi*eta;
53  t61 = t1*eta;
54  t64 = 154495.0/6.0*t36-333385.0/24.0*t39-333385.0/24.0*t42-62335.0/3.0*t45+23365.0/8.0*t47-61375.0/6.0*t49-20350.0*t51-61375.0/6.0*t53+182075.0/24.0*t22+22625.0*t18+22625.0*t35+182075.0/24.0*t25-535.0/2.0*t59+23365.0/8.0*t61-137.0/12.0*xi;
55  t66 = 25.0*xi;
56  t68 = 5275.0/6.0*t59;
57  t75 = 34375.0*t51;
58  t77 = t66+110625.0/4.0*t20-t68+5000.0*t61+59125.0/8.0*t47-1925.0/12.0*t1+8875.0/24.0*t2-4375.0/12.0*t7-31250.0/3.0*t49-t75-75625.0/3.0*t53;
59  t78 = 3125.0*t23;
60  t80 = 227825.0/24.0*t32;
61  t85 = 3125.0/24.0*t11;
62  t87 = 40625.0/4.0*t39;
63  t88 = 105625.0/4.0*t45;
64  t90 = t78+536875.0/24.0*t26-t80+75625.0/12.0*t22+30000.0*t18+42500.0*t35+149375.0/8.0*t25+t85+270625.0/6.0*t36-t87-t88-641875.0/24.0*t42;
65  t93 = 4525.0/6.0*t59;
66  t100 = 38750.0*t51;
67  t102 = -t66-45000.0*t20+t93-25375.0/4.0*t61-10625.0/2.0*t47+2675.0/12.0*t1-7375.0/12.0*t2+8125.0/12.0*t7+49375.0/3.0*t49+t100+15000.0*t53;
68  t103 = 6250.0*t23;
69  t105 = 19175.0/6.0*t32;
70  t110 = 3125.0/12.0*t11;
71  t112 = 40625.0/2.0*t39;
72  t113 = 113750.0/3.0*t45;
73  t115 = -t103-29375.0/3.0*t26+t105-130625.0/12.0*t22-39375.0*t18-117500.0/3.0*t35-31250.0/3.0*t25-t110-140000.0/3.0*t36+t112+t113+284375.0/12.0*t42;
74  t119 = 475.0*t59;
75  t126 = 23125.0*t51;
76  t128 = 50.0/3.0*xi+69375.0/2.0*t20-t119+4250.0*t61+12875.0/4.0*t47-325.0/2.0*t1+6125.0/12.0*t2-625.0*t7-38750.0/3.0*t49-t126-27500.0/3.0*t53;
77  t130 = 34775.0/12.0*t32;
78  t136 = 138125.0/6.0*t45;
79  t138 = t103+83125.0/12.0*t26-t130+27500.0/3.0*t22+25625.0*t18+23125.0*t35+25625.0/4.0*t25+t110+64375.0/3.0*t36-t112-t136-105625.0/12.0*t42;
80  t142 = 2075.0/12.0*t59;
81  t149 = 8125.0*t51;
82  t151 = -25.0/4.0*xi-24375.0/2.0*t20+t142-12625.0/8.0*t61-4625.0/4.0*t47+1525.0/24.0*t1-5125.0/24.0*t2+6875.0/24.0*t7+30625.0/6.0*t49+t149+10000.0/3.0*t53;
83  t153 = 12025.0/12.0*t32;
84  t159 = 8125.0/2.0*t45;
85  t161 = -t78-15625.0/6.0*t26+t153-89375.0/24.0*t22-9375.0*t18-8125.0*t35-9375.0/4.0*t25-t85-21875.0/3.0*t36+t87+t159+105625.0/24.0*t42;
86  t164 = 148.0/3.0*t59;
87  t171 = 3325.0/2.0*t51;
88  t174 = xi+33853.0/24.0*t20-t164+1415.0/2.0*t61+1775.0/8.0*t47-125.0/12.0*t1+875.0/24.0*t2-625.0/12.0*t7-22475.0/6.0*t49-t171-1325.0/2.0*t53+34669.0/6.0*t23;
89  t183 = 25649.0/24.0*t45;
90  t185 = 3809.0/6.0*t26-t30-3809.0/12.0*t32+9263.0/3.0*t22+94067.0/24.0*t18+5287.0/4.0*t35+11759.0/24.0*t25+t12+15529.0/12.0*t36-18031.0/24.0*t39-t183-3809.0/6.0*t42;
91  t188 = 175.0/2.0*t59;
92  t192 = 10125.0/2.0*t51;
93  t196 = 227825.0/24.0*t29;
94  t203 = 25025.0/8.0*t45;
95  t205 = 140075.0/12.0*t20-t188+7125.0/4.0*t61+375.0/2.0*t47-20375.0/2.0*t49-t192+875.0/2.0*t53+374425.0/24.0*t23-8425.0/6.0*t26-t196+t153+34075.0/4.0*t22+140425.0/12.0*t18+66475.0/24.0*t35-2175.0/4.0*t25+40775.0/24.0*t36-269425.0/24.0*t39-t203-9425.0/24.0*t42;
96  t207 = 175.0/3.0*t59;
97  t211 = 7125.0*t51;
98  t215 = 19175.0/6.0*t29;
99  t222 = 168025.0/12.0*t45;
100  t224 = 185425.0/12.0*t20-t207+375.0*t61+1875.0/2.0*t47-125.0*t49-t211-3875.0*t53-8975.0/3.0*t23+55975.0/12.0*t26+t215-t130-650.0/3.0*t22+70325.0/12.0*t18+52975.0/6.0*t35+6025.0/2.0*t25+54275.0/6.0*t36-130975.0/12.0*t39-t222-33475.0/12.0*t42;
101  t232 = 34775.0/12.0*t29;
102  t240 = 54275.0/6.0*t20-t207+1875.0/2.0*t61+375.0*t47-3875.0*t49-t211-125.0*t53+55975.0/12.0*t23-8975.0/3.0*t26-t232+t105+6025.0/2.0*t22+52975.0/6.0*t18+70325.0/12.0*t35-650.0/3.0*t25+185425.0/12.0*t36-33475.0/12.0*t39-t222-130975.0/12.0*t42;
103  t248 = 12025.0/12.0*t29;
104  t256 = 40775.0/24.0*t20-t188+375.0/2.0*t61+7125.0/4.0*t47+875.0/2.0*t49-t192-20375.0/2.0*t53-8425.0/6.0*t23+374425.0/24.0*t26+t248-t80-2175.0/4.0*t22+66475.0/24.0*t18+140425.0/12.0*t35+34075.0/4.0*t25+140075.0/12.0*t36-9425.0/24.0*t39-t203-269425.0/24.0*t42;
105  t266 = eta+15529.0/12.0*t20-t164+1775.0/8.0*t61+1415.0/2.0*t47-125.0/12.0*t5+875.0/24.0*t16-1325.0/2.0*t49-t171-22475.0/6.0*t53-625.0/12.0*t9+3809.0/6.0*t23;
106  t276 = 34669.0/6.0*t26-3809.0/12.0*t29-t33+11759.0/24.0*t22+5287.0/4.0*t18+94067.0/24.0*t35+9263.0/3.0*t25+t14+33853.0/24.0*t36-3809.0/6.0*t39-t183-18031.0/24.0*t42;
107  t287 = -25.0/4.0*eta-21875.0/3.0*t20+t142-4625.0/4.0*t61-12625.0/8.0*t47+1525.0/24.0*t5-5125.0/24.0*t16+10000.0/3.0*t49+t149+30625.0/6.0*t53+6875.0/24.0*t9;
108  t289 = 3125.0*t26;
109  t294 = 3125.0/24.0*t13;
110  t297 = 40625.0/4.0*t42;
111  t298 = -15625.0/6.0*t23-t289+t248-9375.0/4.0*t22-8125.0*t18-9375.0*t35-89375.0/24.0*t25-t294-24375.0/2.0*t36+105625.0/24.0*t39+t159+t297;
112  t309 = 50.0/3.0*eta+64375.0/3.0*t20-t119+12875.0/4.0*t61+4250.0*t47-325.0/2.0*t5+6125.0/12.0*t16-27500.0/3.0*t49-t126-38750.0/3.0*t53-625.0*t9;
113  t311 = 6250.0*t26;
114  t316 = 3125.0/12.0*t13;
115  t319 = 40625.0/2.0*t42;
116  t320 = 83125.0/12.0*t23+t311-t232+25625.0/4.0*t22+23125.0*t18+25625.0*t35+27500.0/3.0*t25+t316+69375.0/2.0*t36-105625.0/12.0*t39-t136-t319;
117  t322 = 25.0*eta;
118  t331 = -t322-140000.0/3.0*t20+t93-10625.0/2.0*t61-25375.0/4.0*t47+2675.0/12.0*t5-7375.0/12.0*t16+15000.0*t49+t100+49375.0/3.0*t53+8125.0/12.0*t9;
119  t339 = -29375.0/3.0*t23-t311+t215-31250.0/3.0*t22-117500.0/3.0*t18-39375.0*t35-130625.0/12.0*t25-t316-45000.0*t36+284375.0/12.0*t39+t113+t319;
120  t349 = t322+270625.0/6.0*t20-t68+59125.0/8.0*t61+5000.0*t47-1925.0/12.0*t5+8875.0/24.0*t16-75625.0/3.0*t49-t75-31250.0/3.0*t53-4375.0/12.0*t9;
121  t357 = 536875.0/24.0*t23+t289-t196+149375.0/8.0*t22+42500.0*t18+30000.0*t35+75625.0/12.0*t25+t294+110625.0/4.0*t36-641875.0/24.0*t39-t88-t297;
122  t378 = -84084.0*t20+1512.0*t59-12600.0*t61-12600.0*t47+46200.0*t49+69300.0*t51+46200.0*t53-48048.0*t23-48048.0*t26+24024.0*t29+24024.0*t32-35112.0*t22-82236.0*t18-82236.0*t35-35112.0*t25-84084.0*t36+48048.0*t39+72072.0*t45+48048.0*t42;
123  t380 = 6912.0*t59;
124  t384 = 252000.0*t51;
125  t396 = 123552.0*t45;
126  t398 = 130944.0*t20-t380+40320.0*t61+57600.0*t47-67200.0*t49-t384-211200.0*t53-33792.0*t23+219648.0*t26+54912.0*t29-109824.0*t32+33792.0*t22+201696.0*t18+328416.0*t35+160512.0*t25+336864.0*t36-13728.0*t39-t396-219648.0*t42;
127  t415 = 336864.0*t20-t380+57600.0*t61+40320.0*t47-211200.0*t49-t384-67200.0*t53+219648.0*t23-33792.0*t26-109824.0*t29+54912.0*t32+160512.0*t22+328416.0*t18+201696.0*t35+33792.0*t25+130944.0*t36-219648.0*t39-t396-13728.0*t42;
128  t417 = 924.0*t59;
129  t421 = 48510.0*t51;
130  t432 = 36036.0*t39;
131  t433 = 45474.0*t45;
132  t435 = -62073.0*t20+t417-9405.0*t61-6435.0*t47+34980.0*t49+t421+18480.0*t53-36036.0*t23-13101.0*t26+18018.0*t29+4719.0*t32-26499.0*t22-61017.0*t18-47487.0*t35-12969.0*t25-53493.0*t36+t432+t433+27456.0*t42;
133  t452 = 36036.0*t42;
134  t453 = -63558.0*t20+t417-7920.0*t61-7920.0*t47+25740.0*t49+51480.0*t51+25740.0*t53-21681.0*t23-21681.0*t26+9009.0*t29+9009.0*t32-18744.0*t22-56232.0*t18-56232.0*t35-18744.0*t25-63558.0*t36+t432+49764.0*t45+t452;
135  t469 = -53493.0*t20+t417-6435.0*t61-9405.0*t47+18480.0*t49+t421+34980.0*t53-13101.0*t23-36036.0*t26+4719.0*t29+18018.0*t32-12969.0*t22-47487.0*t18-61017.0*t35-26499.0*t25-62073.0*t36+27456.0*t39+t433+t452;
136  t471 = 14784.0*t59;
137  t475 = 855360.0*t51;
138  t486 = 864864.0*t45;
139  t488 = 1175328.0*t20-t471+147840.0*t61+100320.0*t47-443520.0*t49-t475-253440.0*t53+206976.0*t23+138336.0*t26-41184.0*t32+310464.0*t22+947232.0*t18+804672.0*t35+167904.0*t25+874368.0*t36-576576.0*t39-t486-370656.0*t42;
140  t493 = 834240.0*t51;
141  t497 = 64064.0*t29;
142  t498 = 64064.0*t32;
143  t504 = 416416.0*t39;
144  t505 = 727584.0*t45;
145  t506 = 416416.0*t42;
146  t507 = 962368.0*t20-t471+132000.0*t61+116160.0*t47-387200.0*t49-t493-323840.0*t53+231616.0*t23+199936.0*t26-t497-t498+269984.0*t22+885632.0*t18+838112.0*t35+222464.0*t25+914848.0*t36-t504-t505-t506;
147  t520 = 914848.0*t20-t471+116160.0*t61+132000.0*t47-323840.0*t49-t493-387200.0*t53+199936.0*t23+231616.0*t26-t497-t498+222464.0*t22+838112.0*t18+885632.0*t35+269984.0*t25+962368.0*t36-t504-t505-t506;
148  t536 = 874368.0*t20-t471+100320.0*t61+147840.0*t47-253440.0*t49-t475-443520.0*t53+138336.0*t23+206976.0*t26-41184.0*t29+167904.0*t22+804672.0*t18+947232.0*t35+310464.0*t25+1175328.0*t36-370656.0*t39-t486-576576.0*t42;
149  t538 = 17160.0*t59;
150  t542 = 1029600.0*t51;
151  t553 = 1441440.0*t39;
152  t554 = 960960.0*t45;
153  t556 = -1921920.0*t20+t538-257400.0*t61-85800.0*t47+1201200.0*t49+t542+171600.0*t53-1441440.0*t23-68640.0*t26+720720.0*t29+17160.0*t32-960960.0*t22-1681680.0*t18-823680.0*t35-102960.0*t25-720720.0*t36+t553+t554+240240.0*t42;
154  t561 = 1287000.0*t51;
155  t572 = 1561560.0*t45;
156  t574 = -2222220.0*t20+t538-214500.0*t61-128700.0*t47+858000.0*t49+t561+343200.0*t53-840840.0*t23-197340.0*t26+360360.0*t29+60060.0*t32-660660.0*t22-1638780.0*t18-1209780.0*t35-231660.0*t25-1364220.0*t36+t553+t572+583440.0*t42;
157  t593 = -1950520.0*t20+t538-171600.0*t61-171600.0*t47+572000.0*t49+1372800.0*t51+572000.0*t53-440440.0*t23-440440.0*t26+160160.0*t29+160160.0*t32-417560.0*t22-1481480.0*t18-1481480.0*t35-417560.0*t25-1950520.0*t36+1041040.0*t39+1818960.0*t45+1041040.0*t42;
158  t609 = 1441440.0*t42;
159  t610 = -1364220.0*t20+t538-128700.0*t61-214500.0*t47+343200.0*t49+t561+858000.0*t53-197340.0*t23-840840.0*t26+60060.0*t29+360360.0*t32-231660.0*t22-1209780.0*t18-1638780.0*t35-660660.0*t25-2222220.0*t36+583440.0*t39+t572+t609;
160  t626 = -720720.0*t20+t538-85800.0*t61-257400.0*t47+171600.0*t49+t542+1201200.0*t53-68640.0*t23-1441440.0*t26+17160.0*t29+720720.0*t32-102960.0*t22-823680.0*t18-1681680.0*t35-960960.0*t25-1921920.0*t36+240240.0*t39+t554+t609;
161 
162  values[0] = t34+t64;
163  values[1] = t77+t90;
164  values[2] = t102+t115;
165  values[3] = t128+t138;
166  values[4] = t151+t161;
167  values[5] = t174+t185;
168  values[6] = t205;
169  values[7] = t224;
170  values[8] = t240;
171  values[9] = t256;
172  values[10] = t266+t276;
173  values[11] = t287+t298;
174  values[12] = t309+t320;
175  values[13] = t331+t339;
176  values[14] = t349+t357;
177  values[15] = t378;
178  values[16] = t398;
179  values[17] = t415;
180  values[18] = t435;
181  values[19] = t453;
182  values[20] = t469;
183  values[21] = t488;
184  values[22] = t507;
185  values[23] = t520;
186  values[24] = t536;
187  values[25] = t556;
188  values[26] = t574;
189  values[27] = t593;
190  values[28] = t610;
191  values[29] = t626;
192 }
193 
194 // values of the derivatives in xi direction
195 static void C_T_UL5_2D_DeriveXi(double xi, double eta, double *values)
196 {
197  double t1, t2, t3, t4, t5, t7, t8, t9, t11, t12, t13, t15, t16, t18, t19;
198  double t20, t22, t23, t24, t25, t26, t27, t30, t31, t32, t33, t34, t37;
199  double t42, t44, t47, t49, t51, t53, t55, t62, t64, t65, t67, t70, t72;
200  double t74, t75, t78, t81, t82, t83, t89, t91, t96, t98, t100, t103, t105;
201  double t107, t108, t111, t114, t115, t116, t122, t124, t129, t131, t133;
202  double t137, t139, t144, t145, t151, t153, t159, t161, t165, t167, t172;
203  double t174, t179, t181, t187, t189, t194, t204, t205, t211, t213, t221;
204  double t226, t228, t230, t236, t238, t244, t246, t250, t255, t257, t259;
205  double t265, t267, t273, t275, t279, t285, t287, t293, t304, t310, t312;
206  double t318, t329, t343, t354, t359, t361, t366, t372, t373, t377, t382;
207  double t384, t389, t395, t396, t400, t410, t419, t429, t438, t454, t469;
208  double t474, t485, t487, t493, t495, t500, t515, t527, t532, t535, t540;
209  double t543, t545, t551, t553, t558, t564, t572, t581, t586, t600, t611;
210  double t616, t626, t627, t633, t635, t640, t645, t646, t647, t648, t649;
211  double t650, t653, t656, t658, t664, t665, t668, t670, t679, t689, t703;
212  double t714, t719, t722, t727, t730, t732, t738, t740, t745, t750, t759;
213  double t761, t767, t773, t789, t803, t808, t818, t825, t830, t844, t855;
214 
215  t1 = xi*xi;
216  t2 = eta*eta;
217  t3 = t1*t2;
218  t4 = 1.0-xi-eta;
219  t5 = t3*t4;
220  t7 = t1*xi;
221  t8 = t7*eta;
222  t9 = t8*t4;
223  t11 = t2*eta;
224  t12 = xi*t11;
225  t13 = t12*t4;
226  t15 = t1*t11;
227  t16 = t15*t4;
228  t18 = t2*t2;
229  t19 = xi*t18;
230  t20 = t19*t4;
231  t22 = t1*t1;
232  t23 = t22*eta;
233  t24 = t23*t4;
234  t25 = 219635.0/12.0*t24;
235  t26 = t7*t2;
236  t27 = t26*t4;
237  t30 = t22*xi*eta;
238  t31 = 43927.0/12.0*t30;
239  t32 = t18*eta;
240  t33 = xi*t32;
241  t34 = 43927.0/12.0*t33;
242  t37 = t22*t2;
243  t42 = xi*t2;
244  t44 = 154495.0/2.0*t5+106370.0/3.0*t9+154495.0/3.0*t13-62335.0*t16-333385.0/12.0*t20-t25-333385.0/6.0*t27+t31+t34-53185.0/6.0*t19-154495.0/6.0*t26+333385.0/24.0*t37-154495.0/6.0*t15+45250.0*t12-53185.0/6.0*t23-40700.0*t42;
245  t47 = xi*eta;
246  t49 = t1*eta;
247  t51 = t18*t4;
248  t53 = t7*t11;
249  t55 = t1*t18;
250  t62 = 3125.0/24.0*t22;
251  t64 = t32*t4;
252  t65 = 43927.0/12.0*t64;
253  t67 = -137.0/12.0+182075.0/6.0*t8+67875.0*t3+23365.0/4.0*t47-61375.0/2.0*t49+53185.0/6.0*t51+62335.0/3.0*t53+333385.0/24.0*t55-535.0/2.0*eta+182075.0/24.0*t18+375.0/4.0*xi+23365.0/8.0*t2+625.0/2.0*t7-t62-2125.0/8.0*t1-t65-61375.0/6.0*t11;
254  t70 = 12500.0*t9;
255  t72 = 316875.0/4.0*t16;
256  t74 = 40625.0*t27;
257  t75 = 227825.0/24.0*t33;
258  t78 = 40625.0/4.0*t37;
259  t81 = 3125.0*t23;
260  t82 = 68750.0*t42;
261  t83 = 25.0+331875.0/4.0*t5+t70+270625.0/3.0*t13-t72-641875.0/12.0*t20-t74+t75-536875.0/24.0*t19-110625.0/4.0*t26+t78-270625.0/6.0*t15+85000.0*t12-t81-t82;
262  t89 = 105625.0/4.0*t53;
263  t91 = 5275.0/6.0*eta;
264  t96 = 15625.0/24.0*t22;
265  t98 = 227825.0/24.0*t64;
266  t100 = 75625.0/3.0*t8+90000.0*t3+10000.0*t47-31250.0*t49+536875.0/24.0*t51+t89+641875.0/24.0*t55-t91+149375.0/8.0*t18-1925.0/6.0*xi+59125.0/8.0*t2-4375.0/3.0*t7+t96+8875.0/8.0*t1-t98-75625.0/3.0*t11;
267  t103 = 25000.0*t9;
268  t105 = 113750.0*t16;
269  t107 = 81250.0*t27;
270  t108 = 19175.0/6.0*t33;
271  t111 = 40625.0/2.0*t37;
272  t114 = 6250.0*t23;
273  t115 = 77500.0*t42;
274  t116 = -25.0-135000.0*t5-t103-280000.0/3.0*t13+t105+284375.0/6.0*t20+t107-t108+29375.0/3.0*t19+45000.0*t26-t111+140000.0/3.0*t15-235000.0/3.0*t12+t114+t115;
275  t122 = 113750.0/3.0*t53;
276  t124 = 4525.0/6.0*eta;
277  t129 = 15625.0/12.0*t22;
278  t131 = 19175.0/6.0*t64;
279  t133 = -130625.0/3.0*t8-118125.0*t3-25375.0/2.0*t47+49375.0*t49-29375.0/3.0*t51-t122-284375.0/12.0*t55+t124-31250.0/3.0*t18+2675.0/6.0*xi-10625.0/2.0*t2+8125.0/3.0*t7-t129-7375.0/4.0*t1+t131+15000.0*t11;
280  t137 = 138125.0/2.0*t16;
281  t139 = 34775.0/12.0*t33;
282  t144 = 46250.0*t42;
283  t145 = 50.0/3.0+208125.0/2.0*t5+t103+128750.0/3.0*t13-t137-105625.0/6.0*t20-t107+t139-83125.0/12.0*t19-69375.0/2.0*t26+t111-64375.0/3.0*t15+46250.0*t12-t114-t144;
284  t151 = 138125.0/6.0*t53;
285  t153 = 475.0*eta;
286  t159 = 34775.0/12.0*t64;
287  t161 = 110000.0/3.0*t8+76875.0*t3+8500.0*t47-38750.0*t49+83125.0/12.0*t51+t151+105625.0/12.0*t55-t153+25625.0/4.0*t18-325.0*xi+12875.0/4.0*t2-2500.0*t7+t129+6125.0/4.0*t1-t159-27500.0/3.0*t11;
288  t165 = 24375.0/2.0*t16;
289  t167 = 12025.0/12.0*t33;
290  t172 = 16250.0*t42;
291  t174 = -73125.0/2.0*t5-t70-43750.0/3.0*t13+t165+105625.0/12.0*t20+t74-t167+15625.0/6.0*t19+24375.0/2.0*t26-t78+21875.0/3.0*t15-16250.0*t12+t81+t172-89375.0/6.0*t8;
292  t179 = 8125.0/2.0*t53;
293  t181 = 2075.0/12.0*eta;
294  t187 = 12025.0/12.0*t64;
295  t189 = -25.0/4.0-28125.0*t3-12625.0/4.0*t47+30625.0/2.0*t49-15625.0/6.0*t51-t179-105625.0/24.0*t55+t181-9375.0/4.0*t18+1525.0/12.0*xi-4625.0/4.0*t2+6875.0/6.0*t7-t96-5125.0/8.0*t1+t187+10000.0/3.0*t11;
296  t194 = 25649.0/8.0*t16;
297  t204 = 1.0+33853.0/8.0*t5+69338.0/3.0*t9+15529.0/6.0*t13-t194-3809.0/3.0*t20-t25-18031.0/6.0*t27+t31+3809.0/12.0*t33-3809.0/6.0*t19-33853.0/24.0*t26+18031.0/24.0*t37-15529.0/12.0*t15+5287.0/2.0*t12-34669.0/6.0*t23;
298  t205 = 3325.0*t42;
299  t211 = 25649.0/24.0*t53;
300  t213 = 148.0/3.0*eta;
301  t221 = -t205+37052.0/3.0*t8+94067.0/8.0*t3+1415.0*t47-22475.0/2.0*t49+3809.0/6.0*t51+t211+3809.0/6.0*t55-t213+11759.0/24.0*t18-125.0/6.0*xi+1775.0/8.0*t2-625.0/3.0*t7+t62+875.0/8.0*t1-3809.0/12.0*t64-1325.0/2.0*t11;
302  t226 = 75075.0/8.0*t16;
303  t228 = 1139125.0/24.0*t24;
304  t230 = 227825.0/24.0*t30;
305  t236 = 140075.0/4.0*t5+374425.0/6.0*t9+40775.0/12.0*t13-t226-9425.0/12.0*t20-t228-269425.0/6.0*t27+t230-t167+8425.0/6.0*t19-140075.0/12.0*t26+269425.0/24.0*t37-40775.0/24.0*t15+66475.0/12.0*t12;
306  t238 = 10125.0*t42;
307  t244 = 25025.0/8.0*t53;
308  t246 = 175.0/2.0*eta;
309  t250 = -374425.0/24.0*t23-t238+34075.0*t8+140425.0/4.0*t3+7125.0/2.0*t47-61125.0/2.0*t49-8425.0/6.0*t51+t244+9425.0/24.0*t55-t246-2175.0/4.0*t18+375.0/2.0*t2+t187+875.0/2.0*t11;
310  t255 = 168025.0/4.0*t16;
311  t257 = 95875.0/6.0*t24;
312  t259 = 19175.0/6.0*t30;
313  t265 = 185425.0/4.0*t5-35900.0/3.0*t9+54275.0/3.0*t13-t255-33475.0/6.0*t20+t257-130975.0/3.0*t27-t259+t139-55975.0/12.0*t19-185425.0/12.0*t26+130975.0/12.0*t37-54275.0/6.0*t15+52975.0/3.0*t12;
314  t267 = 14250.0*t42;
315  t273 = 168025.0/12.0*t53;
316  t275 = 175.0/3.0*eta;
317  t279 = 8975.0/3.0*t23-t267-2600.0/3.0*t8+70325.0/4.0*t3+750.0*t47-375.0*t49+55975.0/12.0*t51+t273+33475.0/12.0*t55-t275+6025.0/2.0*t18+1875.0/2.0*t2-t159-3875.0*t11;
318  t285 = 173875.0/12.0*t24;
319  t287 = 34775.0/12.0*t30;
320  t293 = 54275.0/2.0*t5+55975.0/3.0*t9+185425.0/6.0*t13-t255-130975.0/6.0*t20-t285-33475.0/3.0*t27+t287-t108+8975.0/3.0*t19-54275.0/6.0*t26+33475.0/12.0*t37-185425.0/12.0*t15+70325.0/6.0*t12;
321  t304 = -55975.0/12.0*t23-t267+12050.0*t8+52975.0/2.0*t3+1875.0*t47-11625.0*t49-8975.0/3.0*t51+t273+130975.0/12.0*t55-t275-650.0/3.0*t18+375.0*t2+t131-125.0*t11;
322  t310 = 60125.0/12.0*t24;
323  t312 = 12025.0/12.0*t30;
324  t318 = 40775.0/8.0*t5-16850.0/3.0*t9+140075.0/6.0*t13-t226-269425.0/12.0*t20+t310-9425.0/6.0*t27-t312+t75-374425.0/24.0*t19-40775.0/24.0*t26+9425.0/24.0*t37-140075.0/12.0*t15+140425.0/6.0*t12;
325  t329 = 8425.0/6.0*t23-t238-2175.0*t8+66475.0/8.0*t3+375.0*t47+2625.0/2.0*t49+374425.0/24.0*t51+t244+269425.0/24.0*t55-t246+34075.0/4.0*t18+7125.0/4.0*t2-t98-20375.0/2.0*t11;
326  t343 = 15529.0/4.0*t5+7618.0/3.0*t9+33853.0/12.0*t13-t194-18031.0/12.0*t20-19045.0/12.0*t24-7618.0/3.0*t27+3809.0/12.0*t30+t34-34669.0/6.0*t19-15529.0/12.0*t26+3809.0/6.0*t37-33853.0/24.0*t15+94067.0/12.0*t12;
327  t354 = -3809.0/6.0*t23-t205+11759.0/6.0*t8+15861.0/4.0*t3+1775.0/4.0*t47-3975.0/2.0*t49+34669.0/6.0*t51+t211+18031.0/24.0*t55-t213+9263.0/3.0*t18+1415.0/2.0*t2-t65-22475.0/6.0*t11;
328  t359 = 40625.0/2.0*t20;
329  t361 = 3125.0*t19;
330  t366 = -21875.0*t5-31250.0/3.0*t9-24375.0*t13+t165+t359+t310+105625.0/6.0*t27-t312+t361+21875.0/3.0*t26-105625.0/24.0*t37+24375.0/2.0*t15-18750.0*t12;
331  t372 = 3125.0*t51;
332  t373 = 40625.0/4.0*t55;
333  t377 = 15625.0/6.0*t23+t172-9375.0*t8-24375.0*t3-4625.0/2.0*t47+10000.0*t49-t372-t179-t373+t181-89375.0/24.0*t18-12625.0/8.0*t2+30625.0/6.0*t11;
334  t382 = 40625.0*t20;
335  t384 = 6250.0*t19;
336  t389 = 64375.0*t5+83125.0/3.0*t9+69375.0*t13-t137-t382-t285-105625.0/3.0*t27+t287-t384-64375.0/3.0*t26+105625.0/12.0*t37-69375.0/2.0*t15+51250.0*t12;
337  t395 = 6250.0*t51;
338  t396 = 40625.0/2.0*t55;
339  t400 = -83125.0/12.0*t23-t144+25625.0*t8+69375.0*t3+12875.0/2.0*t47-27500.0*t49+t395+t151+t396-t153+27500.0/3.0*t18+4250.0*t2-38750.0/3.0*t11;
340  t410 = -140000.0*t5-117500.0/3.0*t9-90000.0*t13+t105+t382+t257+284375.0/3.0*t27-t259+t384+140000.0/3.0*t26-284375.0/12.0*t37+45000.0*t15-78750.0*t12;
341  t419 = 29375.0/3.0*t23+t115-125000.0/3.0*t8-117500.0*t3-10625.0*t47+45000.0*t49-t395-t122-t396+t124-130625.0/12.0*t18-25375.0/4.0*t2+49375.0/3.0*t11;
342  t429 = 270625.0/2.0*t5+536875.0/6.0*t9+110625.0/2.0*t13-t72-t359-t228-641875.0/6.0*t27+t230-t361-270625.0/6.0*t26+641875.0/24.0*t37-110625.0/4.0*t15+60000.0*t12;
343  t438 = -536875.0/24.0*t23-t82+149375.0/2.0*t8+127500.0*t3+59125.0/4.0*t47-75625.0*t49+t372+t89+t373-t91+75625.0/12.0*t18+5000.0*t2-31250.0/3.0*t11;
344  t454 = -252252.0*t5-192192.0*t9-168168.0*t13+216216.0*t16+96096.0*t20+120120.0*t24+192192.0*t27-24024.0*t30-24024.0*t33+48048.0*t19+84084.0*t26-48048.0*t37+84084.0*t15-164472.0*t12;
345  t469 = 48048.0*t23+138600.0*t42-140448.0*t8-246708.0*t3-25200.0*t47+138600.0*t49-48048.0*t51-72072.0*t53-48048.0*t55+1512.0*eta-35112.0*t18-12600.0*t2+24024.0*t64+46200.0*t11;
346  t474 = 370656.0*t16;
347  t485 = 392832.0*t5-135168.0*t9+673728.0*t13-t474-439296.0*t20+274560.0*t24-54912.0*t27-54912.0*t30+109824.0*t33-219648.0*t19-130944.0*t26+13728.0*t37-336864.0*t15+656832.0*t12;
348  t487 = 504000.0*t42;
349  t493 = 123552.0*t53;
350  t495 = 6912.0*eta;
351  t500 = 33792.0*t23-t487+135168.0*t8+605088.0*t3+80640.0*t47-201600.0*t49+219648.0*t51+t493+219648.0*t55-t495+160512.0*t18+57600.0*t2-109824.0*t64-211200.0*t11;
352  t515 = 1010592.0*t5+878592.0*t9+261888.0*t13-t474-27456.0*t20-549120.0*t24-878592.0*t27+109824.0*t30-54912.0*t33+33792.0*t19-336864.0*t26+219648.0*t37-130944.0*t15+403392.0*t12;
353  t527 = -219648.0*t23-t487+642048.0*t8+985248.0*t3+115200.0*t47-633600.0*t49-33792.0*t51+t493+13728.0*t55-t495+33792.0*t18+40320.0*t2+54912.0*t64-67200.0*t11;
354  t532 = 136422.0*t16;
355  t535 = 144144.0*t27;
356  t540 = 36036.0*t37;
357  t543 = -186219.0*t5-144144.0*t9-106986.0*t13+t532+54912.0*t20+90090.0*t24+t535-18018.0*t30-4719.0*t33+13101.0*t19+62073.0*t26-t540+53493.0*t15-94974.0*t12;
358  t545 = 97020.0*t42;
359  t551 = 45474.0*t53;
360  t553 = 924.0*eta;
361  t558 = 36036.0*t23+t545-105996.0*t8-183051.0*t3-18810.0*t47+104940.0*t49-13101.0*t51-t551-27456.0*t55+t553-12969.0*t18-6435.0*t2+4719.0*t64+18480.0*t11;
362  t564 = 72072.0*t20;
363  t572 = -190674.0*t5-86724.0*t9-127116.0*t13+149292.0*t16+t564+45045.0*t24+t535-9009.0*t30-9009.0*t33+21681.0*t19+63558.0*t26-t540+63558.0*t15-112464.0*t12;
364  t581 = 36036.0*t55;
365  t586 = 21681.0*t23+102960.0*t42-74976.0*t8-168696.0*t3-15840.0*t47+77220.0*t49-21681.0*t51-49764.0*t53-t581+t553-18744.0*t18-7920.0*t2+9009.0*t64+25740.0*t11;
366  t600 = -160479.0*t5-52404.0*t9-124146.0*t13+t532+t564+23595.0*t24+109824.0*t27-4719.0*t30-18018.0*t33+36036.0*t19+53493.0*t26-27456.0*t37+62073.0*t15-122034.0*t12;
367  t611 = 13101.0*t23+t545-51876.0*t8-142461.0*t3-12870.0*t47+55440.0*t49-36036.0*t51-t551-t581+t553-26499.0*t18-9405.0*t2+18018.0*t64+34980.0*t11;
368  t616 = 2594592.0*t16;
369  t626 = 3525984.0*t5+827904.0*t9+1748736.0*t13-t616-741312.0*t20-2306304.0*t27+41184.0*t33-138336.0*t19-1175328.0*t26+576576.0*t37-874368.0*t15+1609344.0*t12-206976.0*t23;
370  t627 = 1710720.0*t42;
371  t633 = 864864.0*t53;
372  t635 = 14784.0*eta;
373  t640 = -t627+1241856.0*t8+2841696.0*t3+295680.0*t47-1330560.0*t49+138336.0*t51+t633+370656.0*t55-t635+167904.0*t18+100320.0*t2-41184.0*t64-253440.0*t11;
374  t645 = 2182752.0*t16;
375  t646 = 832832.0*t20;
376  t647 = 320320.0*t24;
377  t648 = 1665664.0*t27;
378  t649 = 64064.0*t30;
379  t650 = 64064.0*t33;
380  t653 = 416416.0*t37;
381  t656 = 2887104.0*t5+926464.0*t9+1829696.0*t13-t645-t646-t647-t648+t649+t650-199936.0*t19-962368.0*t26+t653-914848.0*t15+1676224.0*t12;
382  t658 = 1668480.0*t42;
383  t664 = 727584.0*t53;
384  t665 = 416416.0*t55;
385  t668 = 64064.0*t64;
386  t670 = -231616.0*t23-t658+1079936.0*t8+2656896.0*t3+264000.0*t47-1161600.0*t49+199936.0*t51+t664+t665-t635+222464.0*t18+116160.0*t2-t668-323840.0*t11;
387  t679 = 2744544.0*t5+799744.0*t9+1924736.0*t13-t645-t646-t647-t648+t649+t650-231616.0*t19-914848.0*t26+t653-962368.0*t15+1771264.0*t12;
388  t689 = -199936.0*t23-t658+889856.0*t8+2514336.0*t3+232320.0*t47-971520.0*t49+231616.0*t51+t664+t665-t635+269984.0*t18+132000.0*t2-t668-387200.0*t11;
389  t703 = 2623104.0*t5+553344.0*t9+2350656.0*t13-t616-1153152.0*t20-205920.0*t24-1482624.0*t27+41184.0*t30-206976.0*t19-874368.0*t26+370656.0*t37-1175328.0*t15+1894464.0*t12;
390  t714 = -138336.0*t23-t627+671616.0*t8+2414016.0*t3+200640.0*t47-760320.0*t49+206976.0*t51+t633+576576.0*t55-t635+310464.0*t18+147840.0*t2-443520.0*t11;
391  t719 = 2882880.0*t16;
392  t722 = 5765760.0*t27;
393  t727 = 1441440.0*t37;
394  t730 = -5765760.0*t5-5765760.0*t9-1441440.0*t13+t719+480480.0*t20+3603600.0*t24+t722-720720.0*t30-17160.0*t33+68640.0*t19+1921920.0*t26-t727+720720.0*t15-1647360.0*t12;
395  t732 = 2059200.0*t42;
396  t738 = 960960.0*t53;
397  t740 = 17160.0*eta;
398  t745 = 1441440.0*t23+t732-3843840.0*t8-5045040.0*t3-514800.0*t47+3603600.0*t49-68640.0*t51-t738-240240.0*t55+t740-102960.0*t18-85800.0*t2+17160.0*t64+171600.0*t11;
399  t750 = 4684680.0*t16;
400  t759 = -6666660.0*t5-3363360.0*t9-2728440.0*t13+t750+1166880.0*t20+1801800.0*t24+t722-360360.0*t30-60060.0*t33+197340.0*t19+2222220.0*t26-t727+1364220.0*t15-2419560.0*t12;
401  t761 = 2574000.0*t42;
402  t767 = 1561560.0*t53;
403  t773 = 840840.0*t23+t761-2642640.0*t8-4916340.0*t3-429000.0*t47+2574000.0*t49-197340.0*t51-t767-583440.0*t55+t740-231660.0*t18-128700.0*t2+60060.0*t64+343200.0*t11;
404  t789 = -5851560.0*t5-1761760.0*t9-3901040.0*t13+5456880.0*t16+2082080.0*t20+800800.0*t24+4164160.0*t27-160160.0*t30-160160.0*t33+440440.0*t19+1950520.0*t26-1041040.0*t37+1950520.0*t15-2962960.0*t12;
405  t803 = 440440.0*t23+2745600.0*t42-1670240.0*t8-4444440.0*t3-343200.0*t47+1716000.0*t49-440440.0*t51-1818960.0*t53-1041040.0*t55+t740-417560.0*t18-171600.0*t2+160160.0*t64+572000.0*t11;
406  t808 = 2882880.0*t20;
407  t818 = -4092660.0*t5-789360.0*t9-4444440.0*t13+t750+t808+300300.0*t24+2333760.0*t27-60060.0*t30-360360.0*t33+840840.0*t19+1364220.0*t26-583440.0*t37+2222220.0*t15-3277560.0*t12;
408  t825 = 1441440.0*t55;
409  t830 = 197340.0*t23+t761-926640.0*t8-3629340.0*t3-257400.0*t47+1029600.0*t49-840840.0*t51-t767-t825+t740-660660.0*t18-214500.0*t2+360360.0*t64+858000.0*t11;
410  t844 = -2162160.0*t5-274560.0*t9-3843840.0*t13+t719+t808+85800.0*t24+960960.0*t27-17160.0*t30-720720.0*t33+1441440.0*t19+720720.0*t26-240240.0*t37+1921920.0*t15-3363360.0*t12;
411  t855 = 68640.0*t23+t732-411840.0*t8-2471040.0*t3-171600.0*t47+514800.0*t49-1441440.0*t51-t738-t825+t740-960960.0*t18-257400.0*t2+720720.0*t64+1201200.0*t11;
412 
413  values[0] = t44+t67;
414  values[1] = t83+t100;
415  values[2] = t116+t133;
416  values[3] = t145+t161;
417  values[4] = t174+t189;
418  values[5] = t204+t221;
419  values[6] = t236+t250;
420  values[7] = t265+t279;
421  values[8] = t293+t304;
422  values[9] = t318+t329;
423  values[10] = t343+t354;
424  values[11] = t366+t377;
425  values[12] = t389+t400;
426  values[13] = t410+t419;
427  values[14] = t429+t438;
428  values[15] = t454+t469;
429  values[16] = t485+t500;
430  values[17] = t515+t527;
431  values[18] = t543+t558;
432  values[19] = t572+t586;
433  values[20] = t600+t611;
434  values[21] = t626+t640;
435  values[22] = t656+t670;
436  values[23] = t679+t689;
437  values[24] = t703+t714;
438  values[25] = t730+t745;
439  values[26] = t759+t773;
440  values[27] = t789+t803;
441  values[28] = t818+t830;
442  values[29] = t844+t855;
443 }
444 
445 // values of the derivatives in eta direction
446 static void C_T_UL5_2D_DeriveEta(double xi, double eta, double *values)
447 {
448  double t1, t2, t3, t4, t5, t7, t8, t9, t11, t12, t13, t15, t16, t18, t19;
449  double t20, t21, t22, t23, t24, t26, t27, t30, t31, t32, t34, t35, t37;
450  double t38, t39, t41, t43, t45, t46, t48, t55, t61, t64, t67, t73, t74;
451  double t75, t76, t78, t79, t80, t82, t83, t90, t93, t96, t102, t103, t104;
452  double t105, t107, t108, t109, t111, t112, t119, t122, t125, t131, t132;
453  double t134, t135, t137, t138, t145, t150, t156, t157, t159, t160, t162;
454  double t163, t170, t175, t183, t188, t189, t191, t198, t204, t211, t213;
455  double t215, t217, t218, t220, t227, t233, t240, t242, t244, t246, t247;
456  double t249, t256, t262, t270, t272, t274, t287, t295, t297, t299, t312;
457  double t325, t340, t344, t345, t348, t350, t351, t358, t367, t371, t372;
458  double t375, t377, t378, t385, t394, t401, t416, t423, t438, t454, t469;
459  double t477, t484, t485, t487, t494, t500, t514, t527, t534, t535, t541;
460  double t542, t543, t545, t552, t558, t563, t572, t573, t586, t599, t611;
461  double t619, t624, t626, t627, t634, t640, t645, t646, t647, t648, t650;
462  double t652, t653, t654, t655, t656, t657, t658, t670, t677, t689, t702;
463  double t714, t721, t722, t728, t729, t730, t732, t739, t745, t752, t758;
464  double t759, t761, t773, t789, t803, t808, t817, t818, t830, t843, t855;
465 
466  t1 = xi*xi;
467  t2 = eta*eta;
468  t3 = t1*t2;
469  t4 = 1.0-xi-eta;
470  t5 = t3*t4;
471  t7 = t1*xi;
472  t8 = t7*eta;
473  t9 = t8*t4;
474  t11 = t2*eta;
475  t12 = xi*t11;
476  t13 = t12*t4;
477  t15 = t1*t11;
478  t16 = t15*t4;
479  t18 = t2*t2;
480  t19 = xi*t18;
481  t20 = t19*t4;
482  t21 = 219635.0/12.0*t20;
483  t22 = t1*t1;
484  t23 = t22*eta;
485  t24 = t23*t4;
486  t26 = t7*t2;
487  t27 = t26*t4;
488  t30 = t22*xi;
489  t31 = t30*t4;
490  t32 = 43927.0/12.0*t31;
491  t34 = t30*eta;
492  t35 = 43927.0/12.0*t34;
493  t37 = xi*t18*eta;
494  t38 = 43927.0/12.0*t37;
495  t39 = t22*t2;
496  t41 = t7*t11;
497  t43 = t1*t18;
498  t45 = -137.0/12.0+154495.0/2.0*t5+154495.0/3.0*t9+106370.0/3.0*t13-333385.0/6.0*t16-t21-333385.0/12.0*t24-62335.0*t27-53185.0/6.0*t23-t32-53185.0/6.0*t19+t35+t38+333385.0/24.0*t39+62335.0/3.0*t41+333385.0/24.0*t43;
499  t46 = t1*eta;
500  t48 = xi*t2;
501  t55 = 3125.0/24.0*t18;
502  t61 = t22*t4;
503  t64 = xi*eta;
504  t67 = -40700.0*t46-61375.0/2.0*t48+45250.0*t8+67875.0*t3+182075.0/6.0*t12-154495.0/6.0*t26-154495.0/6.0*t15-t55-535.0/2.0*xi+23365.0/8.0*t1+625.0/2.0*t11-61375.0/6.0*t7-2125.0/8.0*t2+53185.0/6.0*t61+375.0/4.0*eta+23365.0/4.0*t64+182075.0/24.0*t22;
505  t73 = 1139125.0/24.0*t20;
506  t74 = 40625.0/2.0*t24;
507  t75 = 316875.0/4.0*t27;
508  t76 = 3125.0*t23;
509  t78 = 227825.0/24.0*t37;
510  t79 = 40625.0/4.0*t39;
511  t80 = 105625.0/4.0*t41;
512  t82 = 270625.0/2.0*t5+110625.0/2.0*t9+536875.0/6.0*t13-641875.0/6.0*t16-t73-t74-t75-t76-536875.0/24.0*t19+t78+t79+t80+641875.0/24.0*t43;
513  t83 = 68750.0*t46;
514  t90 = 5275.0/6.0*xi;
515  t93 = 3125.0*t61;
516  t96 = -t83-75625.0*t48+60000.0*t8+127500.0*t3+149375.0/2.0*t12-110625.0/4.0*t26-270625.0/6.0*t15-t90+5000.0*t1-31250.0/3.0*t7+t93+59125.0/4.0*t64+75625.0/12.0*t22;
517  t102 = 95875.0/6.0*t20;
518  t103 = 40625.0*t24;
519  t104 = 113750.0*t27;
520  t105 = 6250.0*t23;
521  t107 = 19175.0/6.0*t37;
522  t108 = 40625.0/2.0*t39;
523  t109 = 113750.0/3.0*t41;
524  t111 = -140000.0*t5-90000.0*t9-117500.0/3.0*t13+284375.0/3.0*t16+t102+t103+t104+t105+29375.0/3.0*t19-t107-t108-t109-284375.0/12.0*t43;
525  t112 = 77500.0*t46;
526  t119 = 4525.0/6.0*xi;
527  t122 = 6250.0*t61;
528  t125 = t112+45000.0*t48-78750.0*t8-117500.0*t3-125000.0/3.0*t12+45000.0*t26+140000.0/3.0*t15+t119-25375.0/4.0*t1+49375.0/3.0*t7-t122-10625.0*t64-130625.0/12.0*t22;
529  t131 = 173875.0/12.0*t20;
530  t132 = 138125.0/2.0*t27;
531  t134 = 34775.0/12.0*t37;
532  t135 = 138125.0/6.0*t41;
533  t137 = 64375.0*t5+69375.0*t9+83125.0/3.0*t13-105625.0/3.0*t16-t131-t103-t132-t105-83125.0/12.0*t19+t134+t108+t135+105625.0/12.0*t43;
534  t138 = 46250.0*t46;
535  t145 = 475.0*xi;
536  t150 = -t138-27500.0*t48+51250.0*t8+69375.0*t3+25625.0*t12-69375.0/2.0*t26-64375.0/3.0*t15-t145+4250.0*t1-38750.0/3.0*t7+t122+12875.0/2.0*t64+27500.0/3.0*t22;
537  t156 = 60125.0/12.0*t20;
538  t157 = 24375.0/2.0*t27;
539  t159 = 12025.0/12.0*t37;
540  t160 = 8125.0/2.0*t41;
541  t162 = -21875.0*t5-24375.0*t9-31250.0/3.0*t13+105625.0/6.0*t16+t156+t74+t157+t76+15625.0/6.0*t19-t159-t79-t160-105625.0/24.0*t43;
542  t163 = 16250.0*t46;
543  t170 = 2075.0/12.0*xi;
544  t175 = t163+10000.0*t48-18750.0*t8-24375.0*t3-9375.0*t12+24375.0/2.0*t26+21875.0/3.0*t15+t170-12625.0/8.0*t1+30625.0/6.0*t7-t93-4625.0/2.0*t64-89375.0/24.0*t22;
545  t183 = 25649.0/8.0*t27;
546  t188 = 25649.0/24.0*t41;
547  t189 = 15529.0/4.0*t5+33853.0/12.0*t9+7618.0/3.0*t13-7618.0/3.0*t16-19045.0/12.0*t20-18031.0/12.0*t24-t183-34669.0/6.0*t23-t32-3809.0/6.0*t19+t35+3809.0/12.0*t37+18031.0/24.0*t39+t188;
548  t191 = 3325.0*t46;
549  t198 = 148.0/3.0*xi;
550  t204 = 3809.0/6.0*t43-t191-3975.0/2.0*t48+94067.0/12.0*t8+15861.0/4.0*t3+11759.0/6.0*t12-33853.0/24.0*t26-15529.0/12.0*t15-t198+1415.0/2.0*t1-22475.0/6.0*t7+34669.0/6.0*t61+1775.0/4.0*t64+9263.0/3.0*t22;
551  t211 = 75075.0/8.0*t27;
552  t213 = 227825.0/24.0*t31;
553  t215 = 227825.0/24.0*t34;
554  t217 = 25025.0/8.0*t41;
555  t218 = 40775.0/8.0*t5+140075.0/6.0*t9-16850.0/3.0*t13-9425.0/6.0*t16+t156-269425.0/12.0*t24-t211-374425.0/24.0*t23-t213+8425.0/6.0*t19+t215-t159+269425.0/24.0*t39+t217;
556  t220 = 10125.0*t46;
557  t227 = 175.0/2.0*xi;
558  t233 = 9425.0/24.0*t43-t220+2625.0/2.0*t48+140425.0/6.0*t8+66475.0/8.0*t3-2175.0*t12-140075.0/12.0*t26-40775.0/24.0*t15-t227+7125.0/4.0*t1-20375.0/2.0*t7+374425.0/24.0*t61+375.0*t64+34075.0/4.0*t22;
559  t240 = 168025.0/4.0*t27;
560  t242 = 19175.0/6.0*t31;
561  t244 = 19175.0/6.0*t34;
562  t246 = 168025.0/12.0*t41;
563  t247 = 54275.0/2.0*t5+185425.0/6.0*t9+55975.0/3.0*t13-33475.0/3.0*t16-t131-130975.0/6.0*t24-t240+8975.0/3.0*t23+t242-55975.0/12.0*t19-t244+t134+130975.0/12.0*t39+t246;
564  t249 = 14250.0*t46;
565  t256 = 175.0/3.0*xi;
566  t262 = 33475.0/12.0*t43-t249-11625.0*t48+70325.0/6.0*t8+52975.0/2.0*t3+12050.0*t12-185425.0/12.0*t26-54275.0/6.0*t15-t256+375.0*t1-125.0*t7-8975.0/3.0*t61+1875.0*t64-650.0/3.0*t22;
567  t270 = 34775.0/12.0*t31;
568  t272 = 34775.0/12.0*t34;
569  t274 = 185425.0/4.0*t5+54275.0/3.0*t9-35900.0/3.0*t13-130975.0/3.0*t16+t102-33475.0/6.0*t24-t240-55975.0/12.0*t23-t270+8975.0/3.0*t19+t272-t107+33475.0/12.0*t39+t246;
570  t287 = 130975.0/12.0*t43-t249-375.0*t48+52975.0/3.0*t8+70325.0/4.0*t3-2600.0/3.0*t12-54275.0/6.0*t26-185425.0/12.0*t15-t256+1875.0/2.0*t1-3875.0*t7+55975.0/12.0*t61+750.0*t64+6025.0/2.0*t22;
571  t295 = 12025.0/12.0*t31;
572  t297 = 12025.0/12.0*t34;
573  t299 = 140075.0/4.0*t5+40775.0/12.0*t9+374425.0/6.0*t13-269425.0/6.0*t16-t73-9425.0/12.0*t24-t211+8425.0/6.0*t23+t295-374425.0/24.0*t19-t297+t78+9425.0/24.0*t39+t217;
574  t312 = 269425.0/24.0*t43-t220-61125.0/2.0*t48+66475.0/12.0*t8+140425.0/4.0*t3+34075.0*t12-40775.0/24.0*t26-140075.0/12.0*t15-t227+375.0/2.0*t1+875.0/2.0*t7-8425.0/6.0*t61+7125.0/2.0*t64-2175.0/4.0*t22;
575  t325 = 1.0+33853.0/8.0*t5+15529.0/6.0*t9+69338.0/3.0*t13-18031.0/6.0*t16-t21-3809.0/3.0*t24-t183-3809.0/6.0*t23-3809.0/12.0*t31-34669.0/6.0*t19+3809.0/12.0*t34+t38+3809.0/6.0*t39+t188+18031.0/24.0*t43;
576  t340 = -t191-22475.0/2.0*t48+5287.0/2.0*t8+94067.0/8.0*t3+37052.0/3.0*t12-15529.0/12.0*t26-33853.0/24.0*t15+t55-t198+1775.0/8.0*t1-625.0/3.0*t11-1325.0/2.0*t7+875.0/8.0*t2+3809.0/6.0*t61-125.0/6.0*eta+1415.0*t64+11759.0/24.0*t22;
577  t344 = 12500.0*t13;
578  t345 = 40625.0*t16;
579  t348 = 3125.0*t19;
580  t350 = 40625.0/4.0*t43;
581  t351 = -25.0/4.0-73125.0/2.0*t5-43750.0/3.0*t9-t344+t345+105625.0/12.0*t24+t157+15625.0/6.0*t23+t295+t348-t297-105625.0/24.0*t39-t160-t350+t163;
582  t358 = 15625.0/24.0*t18;
583  t367 = 30625.0/2.0*t48-16250.0*t8-28125.0*t3-89375.0/6.0*t12+21875.0/3.0*t26+24375.0/2.0*t15-t358+t170-4625.0/4.0*t1+6875.0/6.0*t11+10000.0/3.0*t7-5125.0/8.0*t2-15625.0/6.0*t61+1525.0/12.0*eta-12625.0/4.0*t64-9375.0/4.0*t22;
584  t371 = 25000.0*t13;
585  t372 = 81250.0*t16;
586  t375 = 6250.0*t19;
587  t377 = 40625.0/2.0*t43;
588  t378 = 50.0/3.0+208125.0/2.0*t5+128750.0/3.0*t9+t371-t372-105625.0/6.0*t24-t132-83125.0/12.0*t23-t270-t375+t272+105625.0/12.0*t39+t135+t377-t138;
589  t385 = 15625.0/12.0*t18;
590  t394 = -38750.0*t48+46250.0*t8+76875.0*t3+110000.0/3.0*t12-64375.0/3.0*t26-69375.0/2.0*t15+t385-t145+12875.0/4.0*t1-2500.0*t11-27500.0/3.0*t7+6125.0/4.0*t2+83125.0/12.0*t61-325.0*eta+8500.0*t64+25625.0/4.0*t22;
591  t401 = -25.0-135000.0*t5-280000.0/3.0*t9-t371+t372+284375.0/6.0*t24+t104+29375.0/3.0*t23+t242+t375-t244-284375.0/12.0*t39-t109-t377+t112;
592  t416 = 49375.0*t48-235000.0/3.0*t8-118125.0*t3-130625.0/3.0*t12+140000.0/3.0*t26+45000.0*t15-t385+t119-10625.0/2.0*t1+8125.0/3.0*t11+15000.0*t7-7375.0/4.0*t2-29375.0/3.0*t61+2675.0/6.0*eta-25375.0/2.0*t64-31250.0/3.0*t22;
593  t423 = 25.0+331875.0/4.0*t5+270625.0/3.0*t9+t344-t345-641875.0/12.0*t24-t75-536875.0/24.0*t23-t213-t348+t215+641875.0/24.0*t39+t80+t350-t83;
594  t438 = -31250.0*t48+85000.0*t8+90000.0*t3+75625.0/3.0*t12-270625.0/6.0*t26-110625.0/4.0*t15+t358-t90+59125.0/8.0*t1-4375.0/3.0*t11-75625.0/3.0*t7+8875.0/8.0*t2+536875.0/24.0*t61-1925.0/6.0*eta+10000.0*t64+149375.0/8.0*t22;
595  t454 = -252252.0*t5-168168.0*t9-192192.0*t13+192192.0*t16+120120.0*t20+96096.0*t24+216216.0*t27+48048.0*t23+24024.0*t31+48048.0*t19-24024.0*t34-24024.0*t37-48048.0*t39-72072.0*t41;
596  t469 = -48048.0*t43+138600.0*t46+138600.0*t48-164472.0*t8-246708.0*t3-140448.0*t12+84084.0*t26+84084.0*t15+1512.0*xi-12600.0*t1+46200.0*t7-48048.0*t61-25200.0*t64-35112.0*t22;
597  t477 = 370656.0*t27;
598  t484 = 123552.0*t41;
599  t485 = 1010592.0*t5+261888.0*t9+878592.0*t13-878592.0*t16-549120.0*t20-27456.0*t24-t477+33792.0*t23+54912.0*t31-219648.0*t19-54912.0*t34+109824.0*t37+13728.0*t39+t484;
600  t487 = 504000.0*t46;
601  t494 = 6912.0*xi;
602  t500 = 219648.0*t43-t487-633600.0*t48+403392.0*t8+985248.0*t3+642048.0*t12-130944.0*t26-336864.0*t15-t494+40320.0*t1-67200.0*t7-33792.0*t61+115200.0*t64+33792.0*t22;
603  t514 = 392832.0*t5+673728.0*t9-135168.0*t13-54912.0*t16+274560.0*t20-439296.0*t24-t477-219648.0*t23-109824.0*t31+33792.0*t19+109824.0*t34-54912.0*t37+219648.0*t39+t484;
604  t527 = 13728.0*t43-t487-201600.0*t48+656832.0*t8+605088.0*t3+135168.0*t12-336864.0*t26-130944.0*t15-t494+57600.0*t1-211200.0*t7+219648.0*t61+80640.0*t64+160512.0*t22;
605  t534 = 72072.0*t24;
606  t535 = 136422.0*t27;
607  t541 = 36036.0*t39;
608  t542 = 45474.0*t41;
609  t543 = -160479.0*t5-124146.0*t9-52404.0*t13+109824.0*t16+23595.0*t20+t534+t535+36036.0*t23+18018.0*t31+13101.0*t19-18018.0*t34-4719.0*t37-t541-t542;
610  t545 = 97020.0*t46;
611  t552 = 924.0*xi;
612  t558 = -27456.0*t43+t545+55440.0*t48-122034.0*t8-142461.0*t3-51876.0*t12+62073.0*t26+53493.0*t15+t552-9405.0*t1+34980.0*t7-36036.0*t61-12870.0*t64-26499.0*t22;
613  t563 = 144144.0*t16;
614  t572 = -190674.0*t5-127116.0*t9-86724.0*t13+t563+45045.0*t20+t534+149292.0*t27+21681.0*t23+9009.0*t31+21681.0*t19-9009.0*t34-9009.0*t37-t541-49764.0*t41;
615  t573 = 36036.0*t43;
616  t586 = -t573+102960.0*t46+77220.0*t48-112464.0*t8-168696.0*t3-74976.0*t12+63558.0*t26+63558.0*t15+t552-7920.0*t1+25740.0*t7-21681.0*t61-15840.0*t64-18744.0*t22;
617  t599 = -186219.0*t5-106986.0*t9-144144.0*t13+t563+90090.0*t20+54912.0*t24+t535+13101.0*t23+4719.0*t31+36036.0*t19-4719.0*t34-18018.0*t37-27456.0*t39-t542;
618  t611 = -t573+t545+104940.0*t48-94974.0*t8-183051.0*t3-105996.0*t12+53493.0*t26+62073.0*t15+t552-6435.0*t1+18480.0*t7-13101.0*t61-18810.0*t64-12969.0*t22;
619  t619 = 2594592.0*t27;
620  t624 = 864864.0*t41;
621  t626 = 2623104.0*t5+2350656.0*t9+553344.0*t13-1482624.0*t16-205920.0*t20-1153152.0*t24-t619-206976.0*t23-138336.0*t19+41184.0*t37+576576.0*t39+t624+370656.0*t43;
622  t627 = 1710720.0*t46;
623  t634 = 14784.0*xi;
624  t640 = -t627-760320.0*t48+1894464.0*t8+2414016.0*t3+671616.0*t12-1175328.0*t26-874368.0*t15-t634+147840.0*t1-443520.0*t7+206976.0*t61+200640.0*t64+310464.0*t22;
625  t645 = 1665664.0*t16;
626  t646 = 320320.0*t20;
627  t647 = 832832.0*t24;
628  t648 = 2182752.0*t27;
629  t650 = 64064.0*t31;
630  t652 = 64064.0*t34;
631  t653 = 64064.0*t37;
632  t654 = 416416.0*t39;
633  t655 = 727584.0*t41;
634  t656 = 2744544.0*t5+1924736.0*t9+799744.0*t13-t645-t646-t647-t648-231616.0*t23-t650-199936.0*t19+t652+t653+t654+t655;
635  t657 = 416416.0*t43;
636  t658 = 1668480.0*t46;
637  t670 = t657-t658-971520.0*t48+1771264.0*t8+2514336.0*t3+889856.0*t12-962368.0*t26-914848.0*t15-t634+132000.0*t1-387200.0*t7+231616.0*t61+232320.0*t64+269984.0*t22;
638  t677 = 2887104.0*t5+1829696.0*t9+926464.0*t13-t645-t646-t647-t648-199936.0*t23-t650-231616.0*t19+t652+t653+t654+t655;
639  t689 = t657-t658-1161600.0*t48+1676224.0*t8+2656896.0*t3+1079936.0*t12-914848.0*t26-962368.0*t15-t634+116160.0*t1-323840.0*t7+199936.0*t61+264000.0*t64+222464.0*t22;
640  t702 = 3525984.0*t5+1748736.0*t9+827904.0*t13-2306304.0*t16-741312.0*t24-t619-138336.0*t23-41184.0*t31-206976.0*t19+41184.0*t34+370656.0*t39+t624+576576.0*t43;
641  t714 = -t627-1330560.0*t48+1609344.0*t8+2841696.0*t3+1241856.0*t12-874368.0*t26-1175328.0*t15-t634+100320.0*t1-253440.0*t7+138336.0*t61+295680.0*t64+167904.0*t22;
642  t721 = 2882880.0*t24;
643  t722 = 2882880.0*t27;
644  t728 = 1441440.0*t39;
645  t729 = 960960.0*t41;
646  t730 = -2162160.0*t5-3843840.0*t9-274560.0*t13+960960.0*t16+85800.0*t20+t721+t722+1441440.0*t23+720720.0*t31+68640.0*t19-720720.0*t34-17160.0*t37-t728-t729;
647  t732 = 2059200.0*t46;
648  t739 = 17160.0*xi;
649  t745 = -240240.0*t43+t732+514800.0*t48-3363360.0*t8-2471040.0*t3-411840.0*t12+1921920.0*t26+720720.0*t15+t739-257400.0*t1+1201200.0*t7-1441440.0*t61-171600.0*t64-960960.0*t22;
650  t752 = 4684680.0*t27;
651  t758 = 1561560.0*t41;
652  t759 = -4092660.0*t5-4444440.0*t9-789360.0*t13+2333760.0*t16+300300.0*t20+t721+t752+840840.0*t23+360360.0*t31+197340.0*t19-360360.0*t34-60060.0*t37-t728-t758;
653  t761 = 2574000.0*t46;
654  t773 = -583440.0*t43+t761+1029600.0*t48-3277560.0*t8-3629340.0*t3-926640.0*t12+2222220.0*t26+1364220.0*t15+t739-214500.0*t1+858000.0*t7-840840.0*t61-257400.0*t64-660660.0*t22;
655  t789 = -5851560.0*t5-3901040.0*t9-1761760.0*t13+4164160.0*t16+800800.0*t20+2082080.0*t24+5456880.0*t27+440440.0*t23+160160.0*t31+440440.0*t19-160160.0*t34-160160.0*t37-1041040.0*t39-1818960.0*t41;
656  t803 = -1041040.0*t43+2745600.0*t46+1716000.0*t48-2962960.0*t8-4444440.0*t3-1670240.0*t12+1950520.0*t26+1950520.0*t15+t739-171600.0*t1+572000.0*t7-440440.0*t61-343200.0*t64-417560.0*t22;
657  t808 = 5765760.0*t16;
658  t817 = -6666660.0*t5-2728440.0*t9-3363360.0*t13+t808+1801800.0*t20+1166880.0*t24+t752+197340.0*t23+60060.0*t31+840840.0*t19-60060.0*t34-360360.0*t37-583440.0*t39-t758;
659  t818 = 1441440.0*t43;
660  t830 = -t818+t761+2574000.0*t48-2419560.0*t8-4916340.0*t3-2642640.0*t12+1364220.0*t26+2222220.0*t15+t739-128700.0*t1+343200.0*t7-197340.0*t61-429000.0*t64-231660.0*t22;
661  t843 = -5765760.0*t5-1441440.0*t9-5765760.0*t13+t808+3603600.0*t20+480480.0*t24+t722+68640.0*t23+17160.0*t31+1441440.0*t19-17160.0*t34-720720.0*t37-240240.0*t39-t729;
662  t855 = -t818+t732+3603600.0*t48-1647360.0*t8-5045040.0*t3-3843840.0*t12+720720.0*t26+1921920.0*t15+t739-85800.0*t1+171600.0*t7-68640.0*t61-514800.0*t64-102960.0*t22;
663 
664  values[0] = t45+t67;
665  values[1] = t82+t96;
666  values[2] = t111+t125;
667  values[3] = t137+t150;
668  values[4] = t162+t175;
669  values[5] = t189+t204;
670  values[6] = t218+t233;
671  values[7] = t247+t262;
672  values[8] = t274+t287;
673  values[9] = t299+t312;
674  values[10] = t325+t340;
675  values[11] = t351+t367;
676  values[12] = t378+t394;
677  values[13] = t401+t416;
678  values[14] = t423+t438;
679  values[15] = t454+t469;
680  values[16] = t485+t500;
681  values[17] = t514+t527;
682  values[18] = t543+t558;
683  values[19] = t572+t586;
684  values[20] = t599+t611;
685  values[21] = t626+t640;
686  values[22] = t656+t670;
687  values[23] = t677+t689;
688  values[24] = t702+t714;
689  values[25] = t730+t745;
690  values[26] = t759+t773;
691  values[27] = t789+t803;
692  values[28] = t817+t830;
693  values[29] = t843+t855;
694 }
695 
696 // values of the derivatives in xi-xi direction
697 static void C_T_UL5_2D_DeriveXiXi(double xi, double eta, double *values)
698 {
699  double t1, t2, t3, t5, t6, t7, t8, t9, t11, t12, t14, t15, t17, t18, t20;
700  double t21, t22, t25, t30, t32, t33, t34, t35, t37, t39, t40, t42, t47;
701  double t49, t50, t52, t56, t57, t59, t64, t66, t68, t69, t72, t73, t76;
702  double t78, t82, t83, t85, t90, t91, t94, t95, t98, t99, t102, t104, t108;
703  double t114, t117, t120, t123, t125, t129, t136, t138, t141, t144, t146;
704  double t149, t153, t158, t162, t171, t175, t176, t184, t187, t189, t191;
705  double t195, t199, t200, t208, t211, t213, t215, t219, t223, t231, t234;
706  double t240, t244, t252, t255, t261, t273, t282, t284, t297, t299, t301;
707  double t303, t316, t318, t320, t336, t352, t365, t377, t381, t390, t395;
708  double t397, t402, t414, t424, t428, t430, t437, t441, t442, t444, t449;
709  double t451, t461, t466, t472, t483, t492, t496, t507, t509, t514, t515;
710  double t517, t518, t519, t526, t529, t530, t531, t532, t533, t536, t538;
711  double t547, t553, t574, t578, t580, t587, t591, t592, t594, t599, t603;
712  double t611, t615, t617, t622, t635, t647, t649, t659, t664, t669, t680;
713  double t689;
714 
715  t1 = xi*xi;
716  t2 = t1*xi;
717  t3 = 3125.0/6.0*t2;
718  t5 = eta*eta;
719  t6 = t5*eta;
720  t7 = xi*t6;
721  t8 = 1.0-xi-eta;
722  t9 = t7*t8;
723  t11 = t1*eta;
724  t12 = t11*t8;
725  t14 = xi*t5;
726  t15 = t14*t8;
727  t17 = t1*t5;
728  t18 = t17*t8;
729  t20 = t2*eta;
730  t21 = t20*t8;
731  t22 = 219635.0/3.0*t21;
732  t25 = t6*t8;
733  t30 = -t3+1875.0/2.0*t1-124670.0*t9+106370.0*t12+154495.0*t15-333385.0/2.0*t18-t22-40700.0*t5-2125.0/4.0*xi+154495.0/3.0*t25+182075.0/2.0*t11+135750.0*t14-154495.0*t17;
734  t32 = t1*t1;
735  t33 = t32*eta;
736  t34 = 219635.0/6.0*t33;
737  t35 = t2*t5;
738  t37 = t1*t6;
739  t39 = t5*t5;
740  t40 = xi*t39;
741  t42 = xi*eta;
742  t47 = t39*t8;
743  t49 = t39*eta;
744  t50 = 43927.0/6.0*t49;
745  t52 = 375.0/4.0-308990.0/3.0*t7+t34+333385.0/3.0*t35+124670.0*t37+333385.0/6.0*t40-61375.0*t42-212740.0/3.0*t20+45250.0*t6-53185.0/3.0*t39-333385.0/12.0*t47+t50+23365.0/4.0*eta;
746  t56 = 316875.0/2.0*t9;
747  t57 = 121875.0*t18;
748  t59 = 37500.0*t12;
749  t64 = 25000.0*t20;
750  t66 = -641875.0/12.0*t47+270625.0/3.0*t25-t56-t57+331875.0/2.0*t15+t59+8875.0/4.0*xi-62500.0*t42+75625.0*t11+180000.0*t14-t64-331875.0/2.0*t17;
751  t68 = 81250.0*t35;
752  t69 = 316875.0/2.0*t37;
753  t72 = 15625.0/6.0*t2;
754  t73 = 68750.0*t5;
755  t76 = 227825.0/12.0*t49;
756  t78 = -1925.0/6.0-541250.0/3.0*t7+t68+t69+641875.0/6.0*t40-4375.0*t1+t72-t73+85000.0*t6-536875.0/12.0*t39+t76+10000.0*eta;
757  t82 = 227500.0*t9;
758  t83 = 243750.0*t18;
759  t85 = 75000.0*t12;
760  t90 = 50000.0*t20;
761  t91 = 2675.0/6.0+284375.0/6.0*t47-280000.0/3.0*t25+t82+t83-270000.0*t15-t85-7375.0/2.0*xi+98750.0*t42-130625.0*t11-236250.0*t14+t90;
762  t94 = 162500.0*t35;
763  t95 = 227500.0*t37;
764  t98 = 15625.0/3.0*t2;
765  t99 = 77500.0*t5;
766  t102 = 19175.0/3.0*t49;
767  t104 = 270000.0*t17+560000.0/3.0*t7-t94-t95-284375.0/3.0*t40+8125.0*t1-t98+t99-235000.0/3.0*t6+58750.0/3.0*t39-t102-25375.0/2.0*eta;
768  t108 = 138125.0*t9;
769  t114 = -325.0-105625.0/6.0*t47+128750.0/3.0*t25-t108-t83+208125.0*t15+t85+6125.0/2.0*xi-77500.0*t42+110000.0*t11+153750.0*t14-t90;
770  t117 = 138125.0*t37;
771  t120 = 46250.0*t5;
772  t123 = 34775.0/6.0*t49;
773  t125 = -208125.0*t17-257500.0/3.0*t7+t94+t117+105625.0/3.0*t40-7500.0*t1+t98-t120+46250.0*t6-83125.0/6.0*t39+t123+8500.0*eta;
774  t129 = 24375.0*t9;
775  t136 = 105625.0/12.0*t47-43750.0/3.0*t25+t129+t57-73125.0*t15-t59-5125.0/4.0*xi+30625.0*t42-89375.0/2.0*t11-56250.0*t14+t64+73125.0*t17;
776  t138 = 24375.0*t37;
777  t141 = 16250.0*t5;
778  t144 = 12025.0/6.0*t49;
779  t146 = 1525.0/12.0+87500.0/3.0*t7-t68-t138-105625.0/6.0*t40+6875.0/2.0*t1-t72+t141-16250.0*t6+15625.0/3.0*t39-t144-12625.0/4.0*eta;
780  t149 = 25649.0/4.0*t9;
781  t153 = 3325.0*t5;
782  t158 = -125.0/6.0+t3-625.0*t1-t149+69338.0*t12+33853.0/4.0*t15-18031.0/2.0*t18-t22-t153+875.0/4.0*xi+15529.0/6.0*t25+37052.0*t11+94067.0/4.0*t14;
783  t162 = 25649.0/4.0*t37;
784  t171 = -33853.0/4.0*t17-15529.0/3.0*t7+t34+18031.0/3.0*t35+t162+7618.0/3.0*t40-22475.0*t42-138676.0/3.0*t20+5287.0/2.0*t6-3809.0/3.0*t39-3809.0/3.0*t47+3809.0/6.0*t49+1415.0*eta;
785  t175 = 75075.0/4.0*t9;
786  t176 = 1139125.0/6.0*t21;
787  t184 = -9425.0/12.0*t47+40775.0/12.0*t25-t175-t176-269425.0/2.0*t18+140075.0/2.0*t15+374425.0/2.0*t12-61125.0*t42+102225.0*t11+140425.0/2.0*t14-374425.0/3.0*t20;
788  t187 = 1139125.0/12.0*t33;
789  t189 = 75075.0/4.0*t37;
790  t191 = 10125.0*t5;
791  t195 = -140075.0/2.0*t17-40775.0/6.0*t7+t187+269425.0/3.0*t35+t189+9425.0/6.0*t40-t191+66475.0/12.0*t6+8425.0/3.0*t39-t144+7125.0/2.0*eta;
792  t199 = 168025.0/2.0*t9;
793  t200 = 191750.0/3.0*t21;
794  t208 = -33475.0/6.0*t47+54275.0/3.0*t25-t199+t200-130975.0*t18+185425.0/2.0*t15-35900.0*t12-750.0*t42-2600.0*t11+70325.0/2.0*t14+71800.0/3.0*t20;
795  t211 = 95875.0/3.0*t33;
796  t213 = 168025.0/2.0*t37;
797  t215 = 14250.0*t5;
798  t219 = -185425.0/2.0*t17-108550.0/3.0*t7-t211+261950.0/3.0*t35+t213+33475.0/3.0*t40-t215+52975.0/3.0*t6-55975.0/6.0*t39+t123+750.0*eta;
799  t223 = 173875.0/3.0*t21;
800  t231 = -130975.0/6.0*t47+185425.0/6.0*t25-t199-t223-33475.0*t18+54275.0*t15+55975.0*t12-23250.0*t42+36150.0*t11+52975.0*t14-111950.0/3.0*t20;
801  t234 = 173875.0/6.0*t33;
802  t240 = -54275.0*t17-185425.0/3.0*t7+t234+66950.0/3.0*t35+t213+130975.0/3.0*t40-t215+70325.0/6.0*t6+17950.0/3.0*t39-t102+1875.0*eta;
803  t244 = 60125.0/3.0*t21;
804  t252 = -269425.0/12.0*t47+140075.0/6.0*t25-t175+t244-9425.0/2.0*t18+40775.0/4.0*t15-16850.0*t12+2625.0*t42-6525.0*t11+66475.0/4.0*t14+33700.0/3.0*t20;
805  t255 = 60125.0/6.0*t33;
806  t261 = -40775.0/4.0*t17-140075.0/3.0*t7-t255+9425.0/3.0*t35+t189+269425.0/6.0*t40-t191+140425.0/6.0*t6-374425.0/12.0*t39+t76+375.0*eta;
807  t273 = -18031.0/12.0*t47+33853.0/12.0*t25-t149-19045.0/3.0*t21-7618.0*t18+15529.0/2.0*t15+7618.0*t12-3975.0*t42+11759.0/2.0*t11+15861.0/2.0*t14-15236.0/3.0*t20;
808  t282 = -15529.0/2.0*t17-33853.0/6.0*t7+19045.0/6.0*t33+15236.0/3.0*t35+t162+18031.0/6.0*t40-t153+94067.0/12.0*t6-34669.0/3.0*t39+t50+1775.0/4.0*eta;
809  t284 = 40625.0/2.0*t47;
810  t297 = 40625.0*t40;
811  t299 = 6250.0*t39;
812  t301 = 62500.0/3.0*t20+43750.0*t17+48750.0*t7-t255-105625.0/3.0*t35-t138-t297+t141-18750.0*t6+t299-4625.0/2.0*eta;
813  t303 = 40625.0*t47;
814  t316 = 81250.0*t40;
815  t318 = 12500.0*t39;
816  t320 = -166250.0/3.0*t20-128750.0*t17-138750.0*t7+t234+211250.0/3.0*t35+t117+t316-t120+51250.0*t6-t318+12875.0/2.0*eta;
817  t336 = 235000.0/3.0*t20+280000.0*t17+180000.0*t7-t211-568750.0/3.0*t35-t95-t316+t99-78750.0*t6+t318-10625.0*eta;
818  t352 = -536875.0/3.0*t20-270625.0*t17-110625.0*t7+t187+641875.0/3.0*t35+t69+t297-t73+60000.0*t6-t299+59125.0/4.0*eta;
819  t365 = 96096.0*t47-168168.0*t25+432432.0*t9+480480.0*t21+576576.0*t18-504504.0*t15-576576.0*t12+277200.0*t42-421344.0*t11-493416.0*t14+384384.0*t20;
820  t377 = 504504.0*t17+336336.0*t7-240240.0*t33-384384.0*t35-432432.0*t37-192192.0*t40+138600.0*t5-164472.0*t6+96096.0*t39-48048.0*t49-25200.0*eta;
821  t381 = 741312.0*t9;
822  t390 = -439296.0*t47+673728.0*t25-t381+1098240.0*t21-164736.0*t18+785664.0*t15-405504.0*t12-403200.0*t42+405504.0*t11+1210176.0*t14+270336.0*t20;
823  t395 = 741312.0*t37;
824  t397 = 504000.0*t5;
825  t402 = -785664.0*t17-1347456.0*t7-549120.0*t33+109824.0*t35+t395+878592.0*t40-t397+656832.0*t6-439296.0*t39+219648.0*t49+80640.0*eta;
826  t414 = -27456.0*t47+261888.0*t25-t381-2196480.0*t21-2635776.0*t18+2021184.0*t15+2635776.0*t12-1267200.0*t42+1926144.0*t11+1970496.0*t14-1757184.0*t20;
827  t424 = -2021184.0*t17-523776.0*t7+1098240.0*t33+1757184.0*t35+t395+54912.0*t40-t397+403392.0*t6+67584.0*t39-109824.0*t49+115200.0*eta;
828  t428 = 272844.0*t9;
829  t430 = 432432.0*t18;
830  t437 = 54912.0*t47-106986.0*t25+t428+360360.0*t21+t430-372438.0*t15-432432.0*t12+209880.0*t42-317988.0*t11-366102.0*t14+288288.0*t20;
831  t441 = 288288.0*t35;
832  t442 = 272844.0*t37;
833  t444 = 97020.0*t5;
834  t449 = 372438.0*t17+213972.0*t7-180180.0*t33-t441-t442-109824.0*t40+t444-94974.0*t6+26202.0*t39-9438.0*t49-18810.0*eta;
835  t451 = 72072.0*t47;
836  t461 = t451-127116.0*t25+298584.0*t9+180180.0*t21+t430-381348.0*t15-260172.0*t12+154440.0*t42-224928.0*t11-337392.0*t14+173448.0*t20;
837  t466 = 144144.0*t40;
838  t472 = 381348.0*t17+254232.0*t7-90090.0*t33-t441-298584.0*t37-t466+102960.0*t5-112464.0*t6+43362.0*t39-18018.0*t49-15840.0*eta;
839  t483 = t451-124146.0*t25+t428+94380.0*t21+329472.0*t18-320958.0*t15-157212.0*t12+110880.0*t42-155628.0*t11-284922.0*t14+104808.0*t20;
840  t492 = 320958.0*t17+248292.0*t7-47190.0*t33-219648.0*t35-t442-t466+t444-122034.0*t6+72072.0*t39-36036.0*t49-12870.0*eta;
841  t496 = 5189184.0*t9;
842  t507 = 5189184.0*t37;
843  t509 = 1710720.0*t5;
844  t514 = -741312.0*t47+1748736.0*t25-t496-6918912.0*t18+7051968.0*t15+2483712.0*t12-2661120.0*t42+3725568.0*t11+5683392.0*t14-1655808.0*t20-7051968.0*t17-3497472.0*t7+4612608.0*t35+t507+1482624.0*t40-t509+1609344.0*t6-276672.0*t39+82368.0*t49+295680.0*eta;
845  t515 = 832832.0*t47;
846  t517 = 4365504.0*t9;
847  t518 = 1281280.0*t21;
848  t519 = 4996992.0*t18;
849  t526 = -t515+1829696.0*t25-t517-t518-t519+5774208.0*t15+2779392.0*t12-2323200.0*t42+3239808.0*t11+5313792.0*t14-1852928.0*t20;
850  t529 = 640640.0*t33;
851  t530 = 3331328.0*t35;
852  t531 = 4365504.0*t37;
853  t532 = 1665664.0*t40;
854  t533 = 1668480.0*t5;
855  t536 = 128128.0*t49;
856  t538 = -5774208.0*t17-3659392.0*t7+t529+t530+t531+t532-t533+1676224.0*t6-399872.0*t39+t536+264000.0*eta;
857  t547 = -t515+1924736.0*t25-t517-t518-t519+5489088.0*t15+2399232.0*t12-1943040.0*t42+2669568.0*t11+5028672.0*t14-1599488.0*t20;
858  t553 = -5489088.0*t17-3849472.0*t7+t529+t530+t531+t532-t533+1771264.0*t6-463232.0*t39+t536+232320.0*eta;
859  t574 = -1106688.0*t20-5246208.0*t17-4701312.0*t7+411840.0*t33+2965248.0*t35+t507+2306304.0*t40-t509+1894464.0*t6-413952.0*t39+200640.0*eta;
860  t578 = 5765760.0*t9;
861  t580 = 17297280.0*t18;
862  t587 = 480480.0*t47-1441440.0*t25+t578+14414400.0*t21+t580-11531520.0*t15-17297280.0*t12+7207200.0*t42-11531520.0*t11-10090080.0*t14+11531520.0*t20;
863  t591 = 11531520.0*t35;
864  t592 = 5765760.0*t37;
865  t594 = 2059200.0*t5;
866  t599 = 11531520.0*t17+2882880.0*t7-7207200.0*t33-t591-t592-960960.0*t40+t594-1647360.0*t6+137280.0*t39-34320.0*t49-514800.0*eta;
867  t603 = 9369360.0*t9;
868  t611 = 1166880.0*t47-2728440.0*t25+t603+7207200.0*t21+t580-13333320.0*t15-10090080.0*t12+5148000.0*t42-7927920.0*t11-9832680.0*t14+6726720.0*t20;
869  t615 = 9369360.0*t37;
870  t617 = 2574000.0*t5;
871  t622 = 13333320.0*t17+5456880.0*t7-3603600.0*t33-t591-t615-2333760.0*t40+t617-2419560.0*t6+394680.0*t39-120120.0*t49-429000.0*eta;
872  t635 = 2082080.0*t47-3901040.0*t25+10913760.0*t9+3203200.0*t21+12492480.0*t18-11703120.0*t15-5285280.0*t12+3432000.0*t42-5010720.0*t11-8888880.0*t14+3523520.0*t20;
873  t647 = 11703120.0*t17+7802080.0*t7-1601600.0*t33-8328320.0*t35-10913760.0*t37-4164160.0*t40+2745600.0*t5-2962960.0*t6+880880.0*t39-320320.0*t49-343200.0*eta;
874  t649 = 2882880.0*t47;
875  t659 = t649-4444440.0*t25+t603+1201200.0*t21+7001280.0*t18-8185320.0*t15-2368080.0*t12+2059200.0*t42-2779920.0*t11-7258680.0*t14+1578720.0*t20;
876  t664 = 5765760.0*t40;
877  t669 = 8185320.0*t17+8888880.0*t7-600600.0*t33-4667520.0*t35-t615-t664+t617-3277560.0*t6+1681680.0*t39-720720.0*t49-257400.0*eta;
878  t680 = t649-3843840.0*t25+t578+343200.0*t21+2882880.0*t18-4324320.0*t15-823680.0*t12+1029600.0*t42-1235520.0*t11-4942080.0*t14+549120.0*t20;
879  t689 = 4324320.0*t17+7687680.0*t7-171600.0*t33-1921920.0*t35-t592-t664+t594-3363360.0*t6+2882880.0*t39-1441440.0*t49-171600.0*eta;
880 
881  values[0] = t30+t52;
882  values[1] = t66+t78;
883  values[2] = t91+t104;
884  values[3] = t114+t125;
885  values[4] = t136+t146;
886  values[5] = t158+t171;
887  values[6] = t184+t195;
888  values[7] = t208+t219;
889  values[8] = t231+t240;
890  values[9] = t252+t261;
891  values[10] = t273+t282;
892  values[11] = t284-24375.0*t25+t129+t244+105625.0/2.0*t18-43750.0*t15-31250.0*t12+20000.0*t42-28125.0*t11-48750.0*t14+t301;
893  values[12] = -t303+69375.0*t25-t108-t223-105625.0*t18+128750.0*t15+83125.0*t12-55000.0*t42+76875.0*t11+138750.0*t14+t320;
894  values[13] = t303-90000.0*t25+t82+t200+284375.0*t18-280000.0*t15-117500.0*t12+90000.0*t42-125000.0*t11-235000.0*t14+t336;
895  values[14] = -t284+110625.0/2.0*t25-t56-t176-641875.0/2.0*t18+270625.0*t15+536875.0/2.0*t12-151250.0*t42+448125.0/2.0*t11+255000.0*t14+t352;
896  values[15] = t365+t377;
897  values[16] = t390+t402;
898  values[17] = t414+t424;
899  values[18] = t437+t449;
900  values[19] = t461+t472;
901  values[20] = t483+t492;
902  values[21] = t514;
903  values[22] = t526+t538;
904  values[23] = t547+t553;
905  values[24] = -1153152.0*t47+2350656.0*t25-t496-823680.0*t21-4447872.0*t18+5246208.0*t15+1660032.0*t12-1520640.0*t42+2014848.0*t11+4828032.0*t14+t574;
906  values[25] = t587+t599;
907  values[26] = t611+t622;
908  values[27] = t635+t647;
909  values[28] = t659+t669;
910  values[29] = t680+t689;
911 }
912 
913 // values of the derivatives in xi-eta direction
914 static void C_T_UL5_2D_DeriveXiEta(double xi, double eta, double *values)
915 {
916  double t1, t2, t3, t5, t6, t7, t8, t9, t10, t11, t12, t13, t14, t16, t17;
917  double t18, t20, t21, t23, t24, t26, t27, t29, t30, t31, t32, t33, t35;
918  double t37, t40, t45, t47, t49, t59, t62, t63, t67, t68, t69, t71, t74;
919  double t76, t78, t84, t89, t92, t93, t97, t98, t99, t101, t104, t106, t108;
920  double t114, t119, t122, t126, t127, t129, t132, t134, t145, t148, t152;
921  double t153, t155, t158, t160, t171, t179, t183, t185, t188, t201, t204;
922  double t209, t210, t213, t215, t218, t231, t234, t239, t240, t243, t245;
923  double t248, t261, t264, t269, t273, t288, t291, t296, t300, t315, t327;
924  double t342, t344, t346, t350, t353, t361, t367, t369, t371, t375, t378;
925  double t386, t392, t400, t413, t421, t434, t450, t466, t475, t480, t482;
926  double t485, t498, t512, t527, t532, t536, t541, t543, t546, t559, t564;
927  double t574, t590, t603, t618, t626, t630, t633, t635, t637, t648, t650;
928  double t651, t652, t653, t654, t657, t658, t659, t662, t664, t667, t669;
929  double t670, t679, t686, t698, t710, t724, t729, t733, t738, t740, t743;
930  double t756, t764, t769, t771, t774, t787, t803, t819, t825, t833, t848;
931  double t861, t876;
932 
933  t1 = eta*eta;
934  t2 = t1*t1;
935  t3 = xi*t2;
936  t5 = xi*xi;
937  t6 = t5*t5;
938  t7 = t6*xi;
939  t8 = 43927.0/12.0*t7;
940  t9 = t2*eta;
941  t10 = 43927.0/12.0*t9;
942  t11 = t5*xi;
943  t12 = t11*eta;
944  t13 = 1.0-xi-eta;
945  t14 = t12*t13;
946  t16 = t1*eta;
947  t17 = xi*t16;
948  t18 = t17*t13;
949  t20 = t5*eta;
950  t21 = t20*t13;
951  t23 = xi*t1;
952  t24 = t23*t13;
953  t26 = t5*t1;
954  t27 = t26*t13;
955  t29 = t2*t13;
956  t30 = 219635.0/12.0*t29;
957  t31 = t6*t13;
958  t32 = 219635.0/12.0*t31;
959  t33 = t11*t13;
960  t35 = t16*t13;
961  t37 = xi*eta;
962  t40 = -535.0/2.0+46085.0*t3+t8+t10-333385.0/3.0*t14-333385.0/3.0*t18+154495.0*t21+154495.0*t24-187005.0*t27-t30-t32+106370.0/3.0*t33+106370.0/3.0*t35-81400.0*t37+135750.0*t20;
963  t45 = t6*eta;
964  t47 = t11*t1;
965  t49 = t5*t16;
966  t59 = 135750.0*t23-86955.0*t12-154495.0*t26-86955.0*t17+46085.0*t45+707395.0/6.0*t47+707395.0/6.0*t49+23365.0/4.0*eta+23365.0/4.0*xi-53185.0/6.0*t2-53185.0/6.0*t6+182075.0/6.0*t11+182075.0/6.0*t16-61375.0/2.0*t5-61375.0/2.0*t1;
967  t62 = 227825.0/24.0*t9;
968  t63 = 81250.0*t14;
969  t67 = 950625.0/4.0*t27;
970  t68 = 1139125.0/24.0*t29;
971  t69 = 12500.0*t33;
972  t71 = 137500.0*t37;
973  t74 = -5275.0/6.0+807625.0/8.0*t3+t62-t63-641875.0/3.0*t18+331875.0/2.0*t21+270625.0*t24-t67-t68+t69+536875.0/6.0*t35-t71+180000.0*t20+255000.0*t23;
974  t76 = 873125.0/4.0*t26;
975  t78 = 40625.0/2.0*t45;
976  t84 = 3125.0*t6;
977  t89 = -135625.0/2.0*t12-t76-359375.0/2.0*t17+t78+479375.0/4.0*t47+2234375.0/12.0*t49+59125.0/4.0*eta+10000.0*xi-536875.0/24.0*t2-t84+75625.0/3.0*t11+149375.0/2.0*t16-31250.0*t5-75625.0*t1;
978  t92 = 19175.0/6.0*t9;
979  t93 = 162500.0*t14;
980  t97 = 341250.0*t27;
981  t98 = 95875.0/6.0*t29;
982  t99 = 25000.0*t33;
983  t101 = 155000.0*t37;
984  t104 = 4525.0/6.0-63375.0*t3-t92+t93+568750.0/3.0*t18-270000.0*t21-280000.0*t24+t97+t98-t99-117500.0/3.0*t35+t101-236250.0*t20-235000.0*t23;
985  t106 = 275000.0*t26;
986  t108 = 40625.0*t45;
987  t114 = 6250.0*t6;
988  t119 = 115000.0*t12+t106+132500.0*t17-t108-195000.0*t47-625625.0/3.0*t49-10625.0*eta-25375.0/2.0*xi+29375.0/3.0*t2+t114-130625.0/3.0*t11-125000.0/3.0*t16+49375.0*t5+45000.0*t1;
989  t122 = 34775.0/12.0*t9;
990  t126 = 414375.0/2.0*t27;
991  t127 = 173875.0/12.0*t29;
992  t129 = 92500.0*t37;
993  t132 = -475.0+128375.0/4.0*t3+t122-t93-211250.0/3.0*t18+208125.0*t21+128750.0*t24-t126-t127+t99+83125.0/3.0*t35-t129+153750.0*t20+138750.0*t23;
994  t134 = 336875.0/2.0*t26;
995  t145 = -94375.0*t12-t134-70625.0*t17+t108+300625.0/2.0*t47+625625.0/6.0*t49+12875.0/2.0*eta+8500.0*xi-83125.0/12.0*t2-t114+110000.0/3.0*t11+25625.0*t16-38750.0*t5-27500.0*t1;
996  t148 = 12025.0/12.0*t9;
997  t152 = 73125.0/2.0*t27;
998  t153 = 60125.0/12.0*t29;
999  t155 = 32500.0*t37;
1000  t158 = 2075.0/12.0-27625.0/2.0*t3-t148+t63+105625.0/3.0*t18-73125.0*t21-43750.0*t24+t152+t153-t69-31250.0/3.0*t35+t155-56250.0*t20-48750.0*t23;
1001  t160 = 116875.0/2.0*t26;
1002  t171 = 36875.0*t12+t160+25000.0*t17-t78-105625.0/2.0*t47-89375.0/3.0*t49-4625.0/2.0*eta-12625.0/4.0*xi+15625.0/6.0*t2+t84-89375.0/6.0*t11-9375.0*t16+30625.0/2.0*t5+10000.0*t1;
1003  t179 = 76947.0/8.0*t27;
1004  t183 = 6650.0*t37;
1005  t185 = -148.0/3.0+11427.0/4.0*t3+t8+3809.0/12.0*t9-18031.0/3.0*t14-15236.0/3.0*t18+33853.0/4.0*t21+15529.0/2.0*t24-t179-19045.0/12.0*t29-t32+69338.0/3.0*t33+7618.0/3.0*t35-t183+94067.0/4.0*t20;
1006  t188 = 64911.0/8.0*t26;
1007  t201 = 15861.0/2.0*t23-103735.0/4.0*t12-t188-10255.0/2.0*t17+39611.0/2.0*t45+149071.0/24.0*t47+137891.0/24.0*t49+1775.0/4.0*eta+1415.0*xi-3809.0/6.0*t2-34669.0/6.0*t6+37052.0/3.0*t11+11759.0/6.0*t16-22475.0/2.0*t5-3975.0/2.0*t1;
1008  t204 = 227825.0/24.0*t7;
1009  t209 = 225225.0/8.0*t27;
1010  t210 = 1139125.0/24.0*t31;
1011  t213 = 20250.0*t37;
1012  t215 = -175.0/2.0-4225.0*t3+t204-t148-269425.0/3.0*t14-9425.0/3.0*t18+140075.0/2.0*t21+40775.0/4.0*t24-t209+t153-t210+374425.0/6.0*t33-16850.0/3.0*t35-t213+140425.0/2.0*t20;
1013  t218 = 320925.0/8.0*t26;
1014  t231 = 66475.0/4.0*t23-85750.0*t12-t218+8875.0/4.0*t17+559325.0/8.0*t45+1302925.0/24.0*t47+262925.0/24.0*t49+375.0*eta+7125.0/2.0*xi+8425.0/6.0*t2-374425.0/24.0*t6+34075.0*t11-2175.0*t16-61125.0/2.0*t5+2625.0/2.0*t1;
1015  t234 = 19175.0/6.0*t7;
1016  t239 = 504075.0/4.0*t27;
1017  t240 = 95875.0/6.0*t31;
1018  t243 = 28500.0*t37;
1019  t245 = -175.0/3.0+80275.0/4.0*t3-t234+t122-261950.0/3.0*t14-66950.0/3.0*t18+185425.0/2.0*t21+54275.0*t24-t239-t127+t240-35900.0/3.0*t33+55975.0/3.0*t35-t243+70325.0/2.0*t20;
1020  t248 = 293975.0/4.0*t26;
1021  t261 = 52975.0*t23-37875.0/2.0*t12-t248-36750.0*t17+5850.0*t45+1027975.0/12.0*t47+637975.0/12.0*t49+1875.0*eta+750.0*xi-55975.0/12.0*t2+8975.0/3.0*t6-2600.0/3.0*t11+12050.0*t16-375.0*t5-11625.0*t1;
1022  t264 = 34775.0/12.0*t7;
1023  t269 = 173875.0/12.0*t31;
1024  t273 = -175.0/3.0+5850.0*t3+t264-t92-66950.0/3.0*t14-261950.0/3.0*t18+54275.0*t21+185425.0/2.0*t24-t239+t98-t269+55975.0/3.0*t33-35900.0/3.0*t35-t243+52975.0*t20;
1025  t288 = 70325.0/2.0*t23-36750.0*t12-t248-37875.0/2.0*t17+80275.0/4.0*t45+637975.0/12.0*t47+1027975.0/12.0*t49+750.0*eta+1875.0*xi+8975.0/3.0*t2-55975.0/12.0*t6+12050.0*t11-2600.0/3.0*t16-11625.0*t5-375.0*t1;
1026  t291 = 12025.0/12.0*t7;
1027  t296 = 60125.0/12.0*t31;
1028  t300 = -175.0/2.0+559325.0/8.0*t3-t291+t62-9425.0/3.0*t14-269425.0/3.0*t18+40775.0/4.0*t21+140075.0/2.0*t24-t209-t68+t296-16850.0/3.0*t33+374425.0/6.0*t35-t213+66475.0/4.0*t20;
1029  t315 = 140425.0/2.0*t23+8875.0/4.0*t12-t218-85750.0*t17-4225.0*t45+262925.0/24.0*t47+1302925.0/24.0*t49+7125.0/2.0*eta+375.0*xi-374425.0/24.0*t2+8425.0/6.0*t6-2175.0*t11+34075.0*t16+2625.0/2.0*t5-61125.0/2.0*t1;
1030  t327 = -148.0/3.0+39611.0/2.0*t3+3809.0/12.0*t7+t10-15236.0/3.0*t14-18031.0/3.0*t18+15529.0/2.0*t21+33853.0/4.0*t24-t179-t30-19045.0/12.0*t31+7618.0/3.0*t33+69338.0/3.0*t35-t183+15861.0/2.0*t20;
1031  t342 = 94067.0/4.0*t23-10255.0/2.0*t12-t188-103735.0/4.0*t17+11427.0/4.0*t45+137891.0/24.0*t47+149071.0/24.0*t49+1415.0*eta+1775.0/4.0*xi-34669.0/6.0*t2-3809.0/6.0*t6+11759.0/6.0*t11+37052.0/3.0*t16-3975.0/2.0*t5-22475.0/2.0*t1;
1032  t344 = 40625.0/2.0*t3;
1033  t346 = 81250.0*t18;
1034  t350 = 12500.0*t35;
1035  t353 = 2075.0/12.0-t344-t291+105625.0/3.0*t14+t346-43750.0*t21-73125.0*t24+t152+t296-31250.0/3.0*t33-t350+t155-48750.0*t20-56250.0*t23;
1036  t361 = 3125.0*t2;
1037  t367 = 25000.0*t12+t160+36875.0*t17-27625.0/2.0*t45-89375.0/3.0*t47-105625.0/2.0*t49-12625.0/4.0*eta-4625.0/2.0*xi+t361+15625.0/6.0*t6-9375.0*t11-89375.0/6.0*t16+10000.0*t5+30625.0/2.0*t1;
1038  t369 = 40625.0*t3;
1039  t371 = 162500.0*t18;
1040  t375 = 25000.0*t35;
1041  t378 = -475.0+t369+t264-211250.0/3.0*t14-t371+128750.0*t21+208125.0*t24-t126-t269+83125.0/3.0*t33+t375-t129+138750.0*t20+153750.0*t23;
1042  t386 = 6250.0*t2;
1043  t392 = -70625.0*t12-t134-94375.0*t17+128375.0/4.0*t45+625625.0/6.0*t47+300625.0/2.0*t49+8500.0*eta+12875.0/2.0*xi-t386-83125.0/12.0*t6+25625.0*t11+110000.0/3.0*t16-27500.0*t5-38750.0*t1;
1044  t400 = 4525.0/6.0-t369-t234+568750.0/3.0*t14+t371-280000.0*t21-270000.0*t24+t97+t240-117500.0/3.0*t33-t375+t101-235000.0*t20-236250.0*t23;
1045  t413 = 132500.0*t12+t106+115000.0*t17-63375.0*t45-625625.0/3.0*t47-195000.0*t49-25375.0/2.0*eta-10625.0*xi+t386+29375.0/3.0*t6-125000.0/3.0*t11-130625.0/3.0*t16+45000.0*t5+49375.0*t1;
1046  t421 = -5275.0/6.0+t344+t204-641875.0/3.0*t14-t346+270625.0*t21+331875.0/2.0*t24-t67-t210+536875.0/6.0*t33+t350-t71+255000.0*t20+180000.0*t23;
1047  t434 = -359375.0/2.0*t12-t76-135625.0/2.0*t17+807625.0/8.0*t45+2234375.0/12.0*t47+479375.0/4.0*t49+10000.0*eta+59125.0/4.0*xi-t361-536875.0/24.0*t6+149375.0/2.0*t11+75625.0/3.0*t16-75625.0*t5-31250.0*t1;
1048  t450 = 1512.0-216216.0*t3-24024.0*t7-24024.0*t9+384384.0*t14+384384.0*t18-504504.0*t21-504504.0*t24+648648.0*t27+120120.0*t29+120120.0*t31-192192.0*t33-192192.0*t35+277200.0*t37-493416.0*t20;
1049  t466 = -493416.0*t23+360360.0*t12+504504.0*t26+360360.0*t17-216216.0*t45-408408.0*t47-408408.0*t49-25200.0*eta-25200.0*xi+48048.0*t2+48048.0*t6-140448.0*t11-140448.0*t16+138600.0*t5+138600.0*t1;
1050  t475 = 1111968.0*t27;
1051  t480 = 1008000.0*t37;
1052  t482 = -6912.0+988416.0*t3-54912.0*t7+109824.0*t9-109824.0*t14-1757184.0*t18+785664.0*t21+2021184.0*t24-t475-549120.0*t29+274560.0*t31-135168.0*t33+878592.0*t35-t480+1210176.0*t20;
1053  t485 = 1403424.0*t26;
1054  t498 = 1970496.0*t23-126720.0*t12-t485-1552320.0*t17-247104.0*t45+425568.0*t47+1249248.0*t49+115200.0*eta+80640.0*xi-219648.0*t2+33792.0*t6+135168.0*t11+642048.0*t16-201600.0*t5-633600.0*t1;
1055  t512 = -6912.0-247104.0*t3+109824.0*t7-54912.0*t9-1757184.0*t14-109824.0*t18+2021184.0*t21+785664.0*t24-t475+274560.0*t29-549120.0*t31+878592.0*t33-135168.0*t35-t480+1970496.0*t20;
1056  t527 = 1210176.0*t23-1552320.0*t12-t485-126720.0*t17+988416.0*t45+1249248.0*t47+425568.0*t49+80640.0*eta+115200.0*xi+33792.0*t2-219648.0*t6+642048.0*t11+135168.0*t16-633600.0*t5-201600.0*t1;
1057  t532 = 288288.0*t14;
1058  t536 = 409266.0*t27;
1059  t541 = 194040.0*t37;
1060  t543 = 924.0-78507.0*t3-18018.0*t7-4719.0*t9+t532+219648.0*t18-372438.0*t21-320958.0*t24+t536+23595.0*t29+90090.0*t31-144144.0*t33-52404.0*t35+t541-366102.0*t20;
1061  t546 = 346698.0*t26;
1062  t559 = -284922.0*t23+268290.0*t12+t546+159390.0*t17-162162.0*t45-280566.0*t47-246246.0*t49-12870.0*eta-18810.0*xi+13101.0*t2+36036.0*t6-105996.0*t11-51876.0*t16+104940.0*t5+55440.0*t1;
1063  t564 = 288288.0*t18;
1064  t574 = 924.0-117117.0*t3-9009.0*t7-9009.0*t9+t532+t564-381348.0*t21-381348.0*t24+447876.0*t27+45045.0*t29+45045.0*t31-86724.0*t33-86724.0*t35+205920.0*t37-337392.0*t20;
1065  t590 = -337392.0*t23+213840.0*t12+381348.0*t26+213840.0*t17-117117.0*t45-293436.0*t47-293436.0*t49-15840.0*eta-15840.0*xi+21681.0*t2+21681.0*t6-74976.0*t11-74976.0*t16+77220.0*t5+77220.0*t1;
1066  t603 = 924.0-162162.0*t3-4719.0*t7-18018.0*t9+219648.0*t14+t564-320958.0*t21-372438.0*t24+t536+90090.0*t29+23595.0*t31-52404.0*t33-144144.0*t35+t541-284922.0*t20;
1067  t618 = -366102.0*t23+159390.0*t12+t546+268290.0*t17-78507.0*t45-246246.0*t47-280566.0*t49-18810.0*eta-12870.0*xi+36036.0*t2+13101.0*t6-51876.0*t11-105996.0*t16+55440.0*t5+104940.0*t1;
1068  t626 = 7783776.0*t27;
1069  t630 = 3421440.0*t37;
1070  t633 = -14784.0+947232.0*t3+41184.0*t9-4612608.0*t14-2965248.0*t18+7051968.0*t21+5246208.0*t24-t626-205920.0*t29+827904.0*t33+553344.0*t35-t630+5683392.0*t20+4828032.0*t23;
1071  t635 = 6149088.0*t26;
1072  t637 = 1153152.0*t45;
1073  t648 = -3178560.0*t12-t635-2302080.0*t17+t637+4900896.0*t47+4077216.0*t49+200640.0*eta+295680.0*xi-138336.0*t2-206976.0*t6+1241856.0*t11+671616.0*t16-1330560.0*t5-760320.0*t1;
1074  t650 = 1153152.0*t3;
1075  t651 = 64064.0*t7;
1076  t652 = 64064.0*t9;
1077  t653 = 3331328.0*t14;
1078  t654 = 3331328.0*t18;
1079  t657 = 6548256.0*t27;
1080  t658 = 320320.0*t29;
1081  t659 = 320320.0*t31;
1082  t662 = 3336960.0*t37;
1083  t664 = -14784.0+t650+t651+t652-t653-t654+5774208.0*t21+5489088.0*t24-t657-t658-t659+926464.0*t33+799744.0*t35-t662+5313792.0*t20;
1084  t667 = 5631648.0*t26;
1085  t669 = 3848416.0*t47;
1086  t670 = 3848416.0*t49;
1087  t679 = 5028672.0*t23-2851200.0*t12-t667-2629440.0*t17+t637+t669+t670+232320.0*eta+264000.0*xi-199936.0*t2-231616.0*t6+1079936.0*t11+889856.0*t16-1161600.0*t5-971520.0*t1;
1088  t686 = -14784.0+t650+t651+t652-t653-t654+5489088.0*t21+5774208.0*t24-t657-t658-t659+799744.0*t33+926464.0*t35-t662+5028672.0*t20;
1089  t698 = 5313792.0*t23-2629440.0*t12-t667-2851200.0*t17+t637+t669+t670+264000.0*eta+232320.0*xi-231616.0*t2-199936.0*t6+889856.0*t11+1079936.0*t16-971520.0*t5-1161600.0*t1;
1090  t710 = -14784.0+t650+41184.0*t7-2965248.0*t14-4612608.0*t18+5246208.0*t21+7051968.0*t24-t626-205920.0*t31+553344.0*t33+827904.0*t35-t630+4828032.0*t20+5683392.0*t23;
1091  t724 = -2302080.0*t12-t635-3178560.0*t17+947232.0*t45+4077216.0*t47+4900896.0*t49+295680.0*eta+200640.0*xi-206976.0*t2-138336.0*t6+671616.0*t11+1241856.0*t16-760320.0*t5-1330560.0*t1;
1092  t729 = 11531520.0*t14;
1093  t733 = 8648640.0*t27;
1094  t738 = 4118400.0*t37;
1095  t740 = 17160.0-566280.0*t3-720720.0*t7-17160.0*t9+t729+1921920.0*t18-11531520.0*t21-4324320.0*t24+t733+85800.0*t29+3603600.0*t31-5765760.0*t33-274560.0*t35+t738-10090080.0*t20;
1096  t743 = 7927920.0*t26;
1097  t756 = -4942080.0*t23+9609600.0*t12+t743+1716000.0*t17-6486480.0*t45-8648640.0*t47-3843840.0*t49-171600.0*eta-514800.0*xi+68640.0*t2+1441440.0*t6-3843840.0*t11-411840.0*t16+3603600.0*t5+514800.0*t1;
1098  t764 = 14054040.0*t27;
1099  t769 = 5148000.0*t37;
1100  t771 = 17160.0-1467180.0*t3-360360.0*t7-60060.0*t9+t729+4667520.0*t18-13333320.0*t21-8185320.0*t24+t764+300300.0*t29+1801800.0*t31-3363360.0*t33-789360.0*t35+t769-9832680.0*t20;
1101  t774 = 10759320.0*t26;
1102  t787 = -7258680.0*t23+7807800.0*t12+t774+3517800.0*t17-4684680.0*t45-10450440.0*t47-7018440.0*t49-257400.0*eta-429000.0*xi+197340.0*t2+840840.0*t6-2642640.0*t11-926640.0*t16+2574000.0*t5+1029600.0*t1;
1103  t803 = 17160.0-2882880.0*t3-160160.0*t7-160160.0*t9+8328320.0*t14+8328320.0*t18-11703120.0*t21-11703120.0*t24+16370640.0*t27+800800.0*t29+800800.0*t31-1761760.0*t33-1761760.0*t35+5491200.0*t37-8888880.0*t20;
1104  t819 = -8888880.0*t23+5662800.0*t12+11703120.0*t26+5662800.0*t17-2882880.0*t45-9621040.0*t47-9621040.0*t49-343200.0*eta-343200.0*xi+440440.0*t2+440440.0*t6-1670240.0*t11-1670240.0*t16+1716000.0*t5+1716000.0*t1;
1105  t825 = 11531520.0*t18;
1106  t833 = 17160.0-4684680.0*t3-60060.0*t7-360360.0*t9+4667520.0*t14+t825-8185320.0*t21-13333320.0*t24+t764+1801800.0*t29+300300.0*t31-789360.0*t33-3363360.0*t35+t769-7258680.0*t20;
1107  t848 = -9832680.0*t23+3517800.0*t12+t774+7807800.0*t17-1467180.0*t45-7018440.0*t47-10450440.0*t49-429000.0*eta-257400.0*xi+840840.0*t2+197340.0*t6-926640.0*t11-2642640.0*t16+1029600.0*t5+2574000.0*t1;
1108  t861 = 17160.0-6486480.0*t3-17160.0*t7-720720.0*t9+1921920.0*t14+t825-4324320.0*t21-11531520.0*t24+t733+3603600.0*t29+85800.0*t31-274560.0*t33-5765760.0*t35+t738-4942080.0*t20;
1109  t876 = -10090080.0*t23+1716000.0*t12+t743+9609600.0*t17-566280.0*t45-3843840.0*t47-8648640.0*t49-514800.0*eta-171600.0*xi+1441440.0*t2+68640.0*t6-411840.0*t11-3843840.0*t16+514800.0*t5+3603600.0*t1;
1110 
1111  values[0] = t40+t59;
1112  values[1] = t74+t89;
1113  values[2] = t104+t119;
1114  values[3] = t132+t145;
1115  values[4] = t158+t171;
1116  values[5] = t185+t201;
1117  values[6] = t215+t231;
1118  values[7] = t245+t261;
1119  values[8] = t273+t288;
1120  values[9] = t300+t315;
1121  values[10] = t327+t342;
1122  values[11] = t353+t367;
1123  values[12] = t378+t392;
1124  values[13] = t400+t413;
1125  values[14] = t421+t434;
1126  values[15] = t450+t466;
1127  values[16] = t482+t498;
1128  values[17] = t512+t527;
1129  values[18] = t543+t559;
1130  values[19] = t574+t590;
1131  values[20] = t603+t618;
1132  values[21] = t633+t648;
1133  values[22] = t664+t679;
1134  values[23] = t686+t698;
1135  values[24] = t710+t724;
1136  values[25] = t740+t756;
1137  values[26] = t771+t787;
1138  values[27] = t803+t819;
1139  values[28] = t833+t848;
1140  values[29] = t861+t876;
1141 }
1142 
1143 // values of the derivatives in eta-eta direction
1144 static void C_T_UL5_2D_DeriveEtaEta(double xi, double eta, double *values)
1145 {
1146  double t1, t3, t5, t6, t8, t9, t11, t13, t14, t16, t17, t18, t20, t21;
1147  double t22, t23, t24, t29, t31, t32, t33, t35, t36, t37, t38, t40, t41;
1148  double t44, t47, t51, t52, t54, t56, t57, t69, t70, t72, t73, t75, t76;
1149  double t78, t80, t81, t93, t94, t96, t97, t99, t100, t103, t104, t116;
1150  double t118, t119, t121, t124, t125, t137, t139, t140, t142, t147, t155;
1151  double t160, t163, t166, t170, t178, t183, t185, t188, t189, t193, t201;
1152  double t206, t208, t211, t212, t223, t231, t232, t243, t251, t252, t264;
1153  double t274, t278, t279, t287, t288, t290, t294, t296, t300, t301, t308;
1154  double t310, t312, t316, t318, t328, t335, t346, t352, t365, t377, t382;
1155  double t390, t395, t398, t402, t414, t424, t426, t429, t437, t441, t442;
1156  double t445, t449, t454, t461, t466, t472, t483, t492, t497, t510, t513;
1157  double t516, t518, t520, t521, t522, t529, t533, t534, t535, t536, t537;
1158  double t540, t541, t550, t556, t575, t576, t579, t587, t591, t592, t595;
1159  double t599, t603, t611, t615, t618, t622, t635, t647, t652, t659, t664;
1160  double t669, t680, t689;
1161 
1162  t1 = eta*eta;
1163  t3 = xi*t1;
1164  t5 = xi*xi;
1165  t6 = t5*eta;
1166  t8 = 1.0-xi-eta;
1167  t9 = t6*t8;
1168  t11 = t3*t8;
1169  t13 = t5*t1;
1170  t14 = t13*t8;
1171  t16 = t5*xi;
1172  t17 = t16*eta;
1173  t18 = t17*t8;
1174  t20 = t1*eta;
1175  t21 = xi*t20;
1176  t22 = t21*t8;
1177  t23 = 219635.0/3.0*t22;
1178  t24 = xi*eta;
1179  t29 = t16*t1;
1180  t31 = 1875.0/2.0*t1+182075.0/2.0*t3+135750.0*t6+154495.0*t9+106370.0*t11-333385.0/2.0*t14-124670.0*t18-t23-61375.0*t24-308990.0/3.0*t17-212740.0/3.0*t21-154495.0*t13+124670.0*t29;
1181  t32 = t5*t5;
1182  t33 = t32*eta;
1183  t35 = t1*t1;
1184  t36 = xi*t35;
1185  t37 = 219635.0/6.0*t36;
1186  t38 = t5*t20;
1187  t40 = t32*xi;
1188  t41 = 43927.0/6.0*t40;
1189  t44 = t16*t8;
1190  t47 = t32*t8;
1191  t51 = 3125.0/6.0*t20;
1192  t52 = 375.0/4.0+333385.0/6.0*t33+t37+333385.0/3.0*t38+t41-53185.0/3.0*t32+23365.0/4.0*xi+154495.0/3.0*t44-2125.0/4.0*eta-333385.0/12.0*t47+45250.0*t16-40700.0*t5-t51;
1193  t54 = 40625.0/2.0*t47;
1194  t56 = 1139125.0/6.0*t22;
1195  t57 = 316875.0/2.0*t18;
1196  t69 = 40625.0*t33;
1197  t70 = 316875.0/2.0*t29;
1198  t72 = 1139125.0/12.0*t36;
1199  t73 = 68750.0*t5;
1200  t75 = 6250.0*t32;
1201  t76 = 448125.0/2.0*t3-110625.0*t17-270625.0*t13-536875.0/3.0*t21+t69+t70+641875.0/3.0*t38+t72-t73+60000.0*t16-t75;
1202  t78 = 40625.0*t47;
1203  t80 = 191750.0/3.0*t22;
1204  t81 = 227500.0*t18;
1205  t93 = 81250.0*t33;
1206  t94 = 227500.0*t29;
1207  t96 = 95875.0/3.0*t36;
1208  t97 = 77500.0*t5;
1209  t99 = 12500.0*t32;
1210  t100 = -125000.0*t3+180000.0*t17+280000.0*t13+235000.0/3.0*t21-t93-t94-568750.0/3.0*t38-t96+t97-78750.0*t16+t99;
1211  t103 = 173875.0/3.0*t22;
1212  t104 = 138125.0*t18;
1213  t116 = 138125.0*t29;
1214  t118 = 173875.0/6.0*t36;
1215  t119 = 46250.0*t5;
1216  t121 = 76875.0*t3-138750.0*t17-128750.0*t13-166250.0/3.0*t21+t93+t116+211250.0/3.0*t38+t118-t119+51250.0*t16-t99;
1217  t124 = 60125.0/3.0*t22;
1218  t125 = 24375.0*t18;
1219  t137 = 24375.0*t29;
1220  t139 = 60125.0/6.0*t36;
1221  t140 = 16250.0*t5;
1222  t142 = -28125.0*t3+48750.0*t17+43750.0*t13+62500.0/3.0*t21-t69-t137-105625.0/3.0*t38-t139+t140-18750.0*t16+t75;
1223  t147 = 25649.0/4.0*t18;
1224  t155 = -18031.0/12.0*t47+33853.0/12.0*t44-19045.0/3.0*t22-t147-7618.0*t14+7618.0*t11+15529.0/2.0*t9+1775.0/4.0*xi-3975.0*t24+15861.0/2.0*t6+11759.0/2.0*t3;
1225  t160 = 25649.0/4.0*t29;
1226  t163 = 3325.0*t5;
1227  t166 = -33853.0/6.0*t17-15529.0/2.0*t13-15236.0/3.0*t21+18031.0/6.0*t33+t160+15236.0/3.0*t38+19045.0/6.0*t36-t163+94067.0/12.0*t16-34669.0/3.0*t32+t41;
1228  t170 = 75075.0/4.0*t18;
1229  t178 = -269425.0/12.0*t47+140075.0/6.0*t44+t124-t170-9425.0/2.0*t14-16850.0*t11+40775.0/4.0*t9+375.0*xi+2625.0*t24+66475.0/4.0*t6-6525.0*t3;
1230  t183 = 75075.0/4.0*t29;
1231  t185 = 10125.0*t5;
1232  t188 = 227825.0/12.0*t40;
1233  t189 = -140075.0/3.0*t17-40775.0/4.0*t13+33700.0/3.0*t21+269425.0/6.0*t33+t183+9425.0/3.0*t38-t139-t185+140425.0/6.0*t16-374425.0/12.0*t32+t188;
1234  t193 = 168025.0/2.0*t18;
1235  t201 = -130975.0/6.0*t47+185425.0/6.0*t44-t103-t193-33475.0*t14+55975.0*t11+54275.0*t9+1875.0*xi-23250.0*t24+52975.0*t6+36150.0*t3;
1236  t206 = 168025.0/2.0*t29;
1237  t208 = 14250.0*t5;
1238  t211 = 19175.0/3.0*t40;
1239  t212 = -185425.0/3.0*t17-54275.0*t13-111950.0/3.0*t21+130975.0/3.0*t33+t206+66950.0/3.0*t38+t118-t208+70325.0/6.0*t16+17950.0/3.0*t32-t211;
1240  t223 = -33475.0/6.0*t47+54275.0/3.0*t44+t80-t193-130975.0*t14-35900.0*t11+185425.0/2.0*t9+750.0*xi-750.0*t24+70325.0/2.0*t6-2600.0*t3;
1241  t231 = 34775.0/6.0*t40;
1242  t232 = -108550.0/3.0*t17-185425.0/2.0*t13+71800.0/3.0*t21+33475.0/3.0*t33+t206+261950.0/3.0*t38-t96-t208+52975.0/3.0*t16-55975.0/6.0*t32+t231;
1243  t243 = -9425.0/12.0*t47+40775.0/12.0*t44-t56-t170-269425.0/2.0*t14+374425.0/2.0*t11+140075.0/2.0*t9+7125.0/2.0*xi-61125.0*t24+140425.0/2.0*t6+102225.0*t3;
1244  t251 = 12025.0/6.0*t40;
1245  t252 = -40775.0/6.0*t17-140075.0/2.0*t13-374425.0/3.0*t21+9425.0/6.0*t33+t183+269425.0/3.0*t38+t72-t185+66475.0/12.0*t16+8425.0/3.0*t32-t251;
1246  t264 = -125.0/6.0-625.0*t1+37052.0*t3+94067.0/4.0*t6+33853.0/4.0*t9+69338.0*t11-18031.0/2.0*t14-t147-t23-22475.0*t24-15529.0/3.0*t17-138676.0/3.0*t21-33853.0/4.0*t13;
1247  t274 = t160+7618.0/3.0*t33+t37+18031.0/3.0*t38+3809.0/6.0*t40-3809.0/3.0*t32+1415.0*xi+15529.0/6.0*t44+875.0/4.0*eta-3809.0/3.0*t47+5287.0/2.0*t16-t163+t51;
1248  t278 = 121875.0*t14;
1249  t279 = 37500.0*t11;
1250  t287 = 105625.0/12.0*t47-43750.0/3.0*t44+t125+t278-t279-73125.0*t9-12625.0/4.0*xi+30625.0*t24-56250.0*t6-89375.0/2.0*t3+87500.0/3.0*t17+73125.0*t13;
1251  t288 = 25000.0*t21;
1252  t290 = 81250.0*t38;
1253  t294 = 15625.0/6.0*t20;
1254  t296 = 1525.0/12.0+t288-105625.0/6.0*t33-t137-t290+t140-16250.0*t16+15625.0/3.0*t32-t251+6875.0/2.0*t1-t294-5125.0/4.0*eta;
1255  t300 = 243750.0*t14;
1256  t301 = 75000.0*t11;
1257  t308 = -325.0-105625.0/6.0*t47+128750.0/3.0*t44-t104-t300+t301+208125.0*t9+8500.0*xi-77500.0*t24+153750.0*t6+110000.0*t3-257500.0/3.0*t17;
1258  t310 = 50000.0*t21;
1259  t312 = 162500.0*t38;
1260  t316 = 15625.0/3.0*t20;
1261  t318 = -208125.0*t13-t310+105625.0/3.0*t33+t116+t312-t119+46250.0*t16-83125.0/6.0*t32+t231-7500.0*t1+t316+6125.0/2.0*eta;
1262  t328 = 2675.0/6.0+284375.0/6.0*t47-280000.0/3.0*t44+t81+t300-t301-270000.0*t9-25375.0/2.0*xi+98750.0*t24-236250.0*t6-130625.0*t3+560000.0/3.0*t17;
1263  t335 = 270000.0*t13+t310-284375.0/3.0*t33-t94-t312+t97-235000.0/3.0*t16+58750.0/3.0*t32-t211+8125.0*t1-t316-7375.0/2.0*eta;
1264  t346 = -641875.0/12.0*t47+270625.0/3.0*t44-t57-t278+t279+331875.0/2.0*t9+10000.0*xi-62500.0*t24+180000.0*t6+75625.0*t3-541250.0/3.0*t17-331875.0/2.0*t13;
1265  t352 = -1925.0/6.0-t288+641875.0/6.0*t33+t70+t290-t73+85000.0*t16-536875.0/12.0*t32+t188-4375.0*t1+t294+8875.0/4.0*eta;
1266  t365 = 96096.0*t47-168168.0*t44+480480.0*t22+432432.0*t18+576576.0*t14-576576.0*t11-504504.0*t9-25200.0*xi+277200.0*t24-493416.0*t6-421344.0*t3;
1267  t377 = 336336.0*t17+504504.0*t13+384384.0*t21-192192.0*t33-432432.0*t29-384384.0*t38-240240.0*t36+138600.0*t5-164472.0*t16+96096.0*t32-48048.0*t40;
1268  t382 = 741312.0*t18;
1269  t390 = -27456.0*t47+261888.0*t44-2196480.0*t22-t382-2635776.0*t14+2635776.0*t11+2021184.0*t9+115200.0*xi-1267200.0*t24+1970496.0*t6+1926144.0*t3;
1270  t395 = 741312.0*t29;
1271  t398 = 504000.0*t5;
1272  t402 = -523776.0*t17-2021184.0*t13-1757184.0*t21+54912.0*t33+t395+1757184.0*t38+1098240.0*t36-t398+403392.0*t16+67584.0*t32-109824.0*t40;
1273  t414 = -439296.0*t47+673728.0*t44+1098240.0*t22-t382-164736.0*t14-405504.0*t11+785664.0*t9+80640.0*xi-403200.0*t24+1210176.0*t6+405504.0*t3;
1274  t424 = -1347456.0*t17-785664.0*t13+270336.0*t21+878592.0*t33+t395+109824.0*t38-549120.0*t36-t398+656832.0*t16-439296.0*t32+219648.0*t40;
1275  t426 = 72072.0*t47;
1276  t429 = 272844.0*t18;
1277  t437 = t426-124146.0*t44+94380.0*t22+t429+329472.0*t14-157212.0*t11-320958.0*t9-12870.0*xi+110880.0*t24-284922.0*t6-155628.0*t3;
1278  t441 = 144144.0*t33;
1279  t442 = 272844.0*t29;
1280  t445 = 97020.0*t5;
1281  t449 = 248292.0*t17+320958.0*t13+104808.0*t21-t441-t442-219648.0*t38-47190.0*t36+t445-122034.0*t16+72072.0*t32-36036.0*t40;
1282  t454 = 432432.0*t14;
1283  t461 = t426-127116.0*t44+180180.0*t22+298584.0*t18+t454-260172.0*t11-381348.0*t9-15840.0*xi+154440.0*t24-337392.0*t6-224928.0*t3;
1284  t466 = 288288.0*t38;
1285  t472 = 254232.0*t17+381348.0*t13+173448.0*t21-t441-298584.0*t29-t466-90090.0*t36+102960.0*t5-112464.0*t16+43362.0*t32-18018.0*t40;
1286  t483 = 54912.0*t47-106986.0*t44+360360.0*t22+t429+t454-432432.0*t11-372438.0*t9-18810.0*xi+209880.0*t24-366102.0*t6-317988.0*t3;
1287  t492 = 213972.0*t17+372438.0*t13+288288.0*t21-109824.0*t33-t442-t466-180180.0*t36+t445-94974.0*t16+26202.0*t32-9438.0*t40;
1288  t497 = 5189184.0*t18;
1289  t510 = 5189184.0*t29;
1290  t513 = 1710720.0*t5;
1291  t516 = 2014848.0*t3-4701312.0*t17-5246208.0*t13-1106688.0*t21+2306304.0*t33+t510+2965248.0*t38+411840.0*t36-t513+1894464.0*t16-413952.0*t32;
1292  t518 = 832832.0*t47;
1293  t520 = 1281280.0*t22;
1294  t521 = 4365504.0*t18;
1295  t522 = 4996992.0*t14;
1296  t529 = -t518+1924736.0*t44-t520-t521-t522+2399232.0*t11+5489088.0*t9+232320.0*xi-1943040.0*t24+5028672.0*t6+2669568.0*t3;
1297  t533 = 1665664.0*t33;
1298  t534 = 4365504.0*t29;
1299  t535 = 3331328.0*t38;
1300  t536 = 640640.0*t36;
1301  t537 = 1668480.0*t5;
1302  t540 = 128128.0*t40;
1303  t541 = -3849472.0*t17-5489088.0*t13-1599488.0*t21+t533+t534+t535+t536-t537+1771264.0*t16-463232.0*t32+t540;
1304  t550 = -t518+1829696.0*t44-t520-t521-t522+2779392.0*t11+5774208.0*t9+264000.0*xi-2323200.0*t24+5313792.0*t6+3239808.0*t3;
1305  t556 = -3659392.0*t17-5774208.0*t13-1852928.0*t21+t533+t534+t535+t536-t537+1676224.0*t16-399872.0*t32+t540;
1306  t575 = -741312.0*t47+1748736.0*t44-t497-6918912.0*t14+2483712.0*t11+7051968.0*t9+295680.0*xi-2661120.0*t24+5683392.0*t6+3725568.0*t3-3497472.0*t17-7051968.0*t13-1655808.0*t21+1482624.0*t33+t510+4612608.0*t38-t513+1609344.0*t16-276672.0*t32+82368.0*t40;
1307  t576 = 2882880.0*t47;
1308  t579 = 5765760.0*t18;
1309  t587 = t576-3843840.0*t44+343200.0*t22+t579+2882880.0*t14-823680.0*t11-4324320.0*t9-171600.0*xi+1029600.0*t24-4942080.0*t6-1235520.0*t3;
1310  t591 = 5765760.0*t33;
1311  t592 = 5765760.0*t29;
1312  t595 = 2059200.0*t5;
1313  t599 = 7687680.0*t17+4324320.0*t13+549120.0*t21-t591-t592-1921920.0*t38-171600.0*t36+t595-3363360.0*t16+2882880.0*t32-1441440.0*t40;
1314  t603 = 9369360.0*t18;
1315  t611 = t576-4444440.0*t44+1201200.0*t22+t603+7001280.0*t14-2368080.0*t11-8185320.0*t9-257400.0*xi+2059200.0*t24-7258680.0*t6-2779920.0*t3;
1316  t615 = 9369360.0*t29;
1317  t618 = 2574000.0*t5;
1318  t622 = 8888880.0*t17+8185320.0*t13+1578720.0*t21-t591-t615-4667520.0*t38-600600.0*t36+t618-3277560.0*t16+1681680.0*t32-720720.0*t40;
1319  t635 = 2082080.0*t47-3901040.0*t44+3203200.0*t22+10913760.0*t18+12492480.0*t14-5285280.0*t11-11703120.0*t9-343200.0*xi+3432000.0*t24-8888880.0*t6-5010720.0*t3;
1320  t647 = 7802080.0*t17+11703120.0*t13+3523520.0*t21-4164160.0*t33-10913760.0*t29-8328320.0*t38-1601600.0*t36+2745600.0*t5-2962960.0*t16+880880.0*t32-320320.0*t40;
1321  t652 = 17297280.0*t14;
1322  t659 = 1166880.0*t47-2728440.0*t44+7207200.0*t22+t603+t652-10090080.0*t11-13333320.0*t9-429000.0*xi+5148000.0*t24-9832680.0*t6-7927920.0*t3;
1323  t664 = 11531520.0*t38;
1324  t669 = 5456880.0*t17+13333320.0*t13+6726720.0*t21-2333760.0*t33-t615-t664-3603600.0*t36+t618-2419560.0*t16+394680.0*t32-120120.0*t40;
1325  t680 = 480480.0*t47-1441440.0*t44+14414400.0*t22+t579+t652-17297280.0*t11-11531520.0*t9-514800.0*xi+7207200.0*t24-10090080.0*t6-11531520.0*t3;
1326  t689 = 2882880.0*t17+11531520.0*t13+11531520.0*t21-960960.0*t33-t592-t664-7207200.0*t36+t595-1647360.0*t16+137280.0*t32-34320.0*t40;
1327 
1328  values[0] = t31+t52;
1329  values[1] = -t54+110625.0/2.0*t44-t56-t57-641875.0/2.0*t14+536875.0/2.0*t11+270625.0*t9+59125.0/4.0*xi-151250.0*t24+255000.0*t6+t76;
1330  values[2] = t78-90000.0*t44+t80+t81+284375.0*t14-117500.0*t11-280000.0*t9-10625.0*xi+90000.0*t24-235000.0*t6+t100;
1331  values[3] = -t78+69375.0*t44-t103-t104-105625.0*t14+83125.0*t11+128750.0*t9+12875.0/2.0*xi-55000.0*t24+138750.0*t6+t121;
1332  values[4] = t54-24375.0*t44+t124+t125+105625.0/2.0*t14-31250.0*t11-43750.0*t9-4625.0/2.0*xi+20000.0*t24-48750.0*t6+t142;
1333  values[5] = t155+t166;
1334  values[6] = t178+t189;
1335  values[7] = t201+t212;
1336  values[8] = t223+t232;
1337  values[9] = t243+t252;
1338  values[10] = t264+t274;
1339  values[11] = t287+t296;
1340  values[12] = t308+t318;
1341  values[13] = t328+t335;
1342  values[14] = t346+t352;
1343  values[15] = t365+t377;
1344  values[16] = t390+t402;
1345  values[17] = t414+t424;
1346  values[18] = t437+t449;
1347  values[19] = t461+t472;
1348  values[20] = t483+t492;
1349  values[21] = -1153152.0*t47+2350656.0*t44-823680.0*t22-t497-4447872.0*t14+1660032.0*t11+5246208.0*t9+200640.0*xi-1520640.0*t24+4828032.0*t6+t516;
1350  values[22] = t529+t541;
1351  values[23] = t550+t556;
1352  values[24] = t575;
1353  values[25] = t587+t599;
1354  values[26] = t611+t622;
1355  values[27] = t635+t647;
1356  values[28] = t659+t669;
1357  values[29] = t680+t689;
1358 }
1359 
1360 // ***********************************************************************
1361 
1362 TBaseFunct2D *BF_C_T_UL5_2D_Obj = new TBaseFunct2D
1363  (30, BF_C_T_UL5_2D, BFUnitTriangle,
1364  C_T_UL5_2D_Funct, C_T_UL5_2D_DeriveXi,
1365  C_T_UL5_2D_DeriveEta, C_T_UL5_2D_DeriveXiXi,
1366  C_T_UL5_2D_DeriveXiEta, C_T_UL5_2D_DeriveEtaEta, 7, 5,
1367  0, NULL);
Definition: BaseFunct2D.h:27