ParMooN
 All Classes Functions Variables Friends Pages
BF_N_T_P5_2D.h
1 // ***********************************************************************
2 // P5 element, nonconforming, 2D
3 // ***********************************************************************
4 
5 // base function values
6 static void N_T_P5_2D_Funct(double xi, double eta, double *values)
7 {
8  double t1, t2, t3, t4, t5, t6, t7, t8, t9, t11, t12, t13, t14, t15;
9  double t16, t17, t18, t19, t21, t23, t25, t27, t29, t31, t33, t35, t37;
10  double t39, t41, t54, t55, t59, t68, t70, t73, t74, t77, t78, t80, t81;
11  double t82, t95, t106, t107, t108, t120, t129, t131, t135, t143, t147;
12  double t156, t170, t181, t182, t194, t206, t218, t219, t223, t229, t230;
13  double t231, t232, t250, t251, t263, t271, t273, t277, t285, t289, t298;
14  double t312, t323, t324, t336, t348, t362, t366, t374, t375, t377, t381;
15  double t389, t392, t393, t399, t400, t401, t410, t413, t422, t423, t424;
16  double t425, t427, t431, t435, t439, t443, t449, t450, t451, t452, t464;
17  double t466, t468, t474, t475, t476, t477, t479, t483, t487, t489, t499;
18  double t501, t505, t511, t514, t521, t532, t540, t550, t559, t572, t586;
19 
20  t1 = 30.0*eta;
21  t2 = eta*eta;
22  t3 = 210.0*t2;
23  t4 = t2*eta;
24  t5 = 560.0*t4;
25  t6 = t2*t2;
26  t7 = 630.0*t6;
27  t8 = t6*eta;
28  t9 = 252.0*t8;
29  t11 = 30.0*xi;
30  t12 = xi*xi;
31  t13 = 210.0*t12;
32  t14 = t12*xi;
33  t15 = 560.0*t14;
34  t16 = t12*t12;
35  t17 = 630.0*t16;
36  t18 = t16*xi;
37  t19 = 252.0*t18;
38  t21 = xi*eta;
39  t23 = t12*eta;
40  t25 = xi*t2;
41  t27 = t14*eta;
42  t29 = t12*t2;
43  t31 = xi*t4;
44  t33 = t16*eta;
45  t35 = t14*t2;
46  t37 = t12*t4;
47  t39 = xi*t6;
48  t41 = t9-420.0*t21+1680.0*t23+1680.0*t25-2520.0*t27-3780.0*t29-2520.0*t31+1260.0*t33+2520.0*t35+2520.0*t37+1260.0*t39;
49  t54 = 40.0*t21;
50  t55 = -20.0/21.0+4.0/7.0*xi+240.0/7.0*eta+10.0*t12-300.0*t2-80.0/3.0*t14+3200.0/3.0*t4+30.0*t16-1580.0*t6-12.0*t18+780.0*t8-t54;
51  t59 = 540.0*t29;
52  t68 = (1.0-2.0*xi-eta)*(xi-eta);
53  t70 = 2.0*eta-1.0+xi;
54  t73 = t68*t70*t12*eta;
55  t74 = 66.0*t73;
56  t77 = t68*t70*xi*t2;
57  t78 = 66.0*t77;
58  t80 = t68*t70*t4;
59  t81 = 220.0*t80;
60  t82 = -160.0*t23+400.0*t25+294.0*t27+t59-1474.0*t31-162.0*t33-522.0*t35-258.0*t37+1542.0*t39-t74-t78-t81;
61  t95 = 40.0/7.0*xi-40.0/7.0*eta-100.0*t12+100.0*t2+1600.0/3.0*t14-1600.0/3.0*t4-980.0*t16+980.0*t6+540.0*t18-540.0*t8+560.0*t23;
62  t106 = t68*t70*t14;
63  t107 = 220.0*t106;
64  t108 = -560.0*t25-2146.0*t27+2146.0*t31+1998.0*t33+1398.0*t35-1398.0*t37-1998.0*t39+594.0*t73+594.0*t77+t107+t81;
65  t120 = 20.0/21.0-240.0/7.0*xi-4.0/7.0*eta+300.0*t12-10.0*t2-3200.0/3.0*t14+80.0/3.0*t4+1580.0*t16-30.0*t6-780.0*t18+12.0*t8+t54;
66  t129 = -400.0*t23+160.0*t25+1474.0*t27-t59-294.0*t31-1542.0*t33+258.0*t35+522.0*t37+162.0*t39-t74-t78-t107;
67  t131 = 8046.0/77.0*t21;
68  t135 = 2088.0/77.0*t29;
69  t143 = t131+23616.0/77.0*t23-83760.0/77.0*t27-76176.0/77.0*t25-t135-180900.0/77.0*t39+1440.0/7.0*t18-30720.0/77.0*t16+19200.0/77.0*t14-4338.0/77.0*t12+18.0/77.0*xi-5184.0/7.0*t8;
70  t147 = 1500.0/7.0*t80;
71  t156 = 459.0/539.0+114810.0/77.0*t6-10960.0/11.0*t4+21270.0/77.0*t2+t147-16770.0/539.0*eta+81720.0/77.0*t35+62640.0/77.0*t33-74700.0/77.0*t37+216012.0/77.0*t31+3240.0/7.0*t77+720.0/7.0*t73+480.0/7.0*t106;
72  t170 = 3690.0/77.0*t21-19680.0/77.0*t23+85728.0/77.0*t27-19680.0/77.0*t25-22284.0/77.0*t29-99720.0/77.0*t39-3816.0/7.0*t18+77190.0/77.0*t16-43280.0/77.0*t14+8730.0/77.0*t12-4230.0/539.0*xi-3816.0/7.0*t8;
73  t181 = 1500.0/7.0*t106;
74  t182 = 41.0/539.0+77190.0/77.0*t6-43280.0/77.0*t4+8730.0/77.0*t2+t147-4230.0/539.0*eta+14400.0/77.0*t35-99720.0/77.0*t33+14400.0/77.0*t37+85728.0/77.0*t31+180.0*t77-180.0*t73-t181;
75  t194 = t131-76176.0/77.0*t23+216012.0/77.0*t27+23616.0/77.0*t25-t135+62640.0/77.0*t39-5184.0/7.0*t18+114810.0/77.0*t16-10960.0/11.0*t14+21270.0/77.0*t12-16770.0/539.0*xi+1440.0/7.0*t8;
76  t206 = 459.0/539.0-30720.0/77.0*t6+19200.0/77.0*t4-4338.0/77.0*t2-480.0/7.0*t80+18.0/77.0*eta-74700.0/77.0*t35-180900.0/77.0*t33+81720.0/77.0*t37-83760.0/77.0*t31-720.0/7.0*t77-3240.0/7.0*t73-t181;
77  t218 = 810.0/7.0*t21;
78  t219 = -30.0/49.0+18.0/49.0*xi+1080.0/49.0*eta+360.0/7.0*t12-1350.0/7.0*t2-1380.0/7.0*t14+4800.0/7.0*t4+1710.0/7.0*t16-7110.0/7.0*t6-684.0/7.0*t18+3510.0/7.0*t8-t218;
79  t223 = 3240.0/7.0*t29;
80  t229 = 3762.0/7.0*t73;
81  t230 = 3762.0/7.0*t77;
82  t231 = 990.0/7.0*t80;
83  t232 = -3240.0/7.0*t23+8100.0/7.0*t25+1854.0*t27-t223-20808.0/7.0*t31-9234.0/7.0*t33-10854.0/7.0*t35+13644.0/7.0*t37+15444.0/7.0*t39-t229-t230-t231;
84  t250 = 990.0/7.0*t106;
85  t251 = 180.0/49.0*xi-180.0/49.0*eta-450.0/7.0*t12+450.0/7.0*t2+2400.0/7.0*t14-2400.0/7.0*t4-t17+t7+2430.0/7.0*t18-2430.0/7.0*t8-1782.0/7.0*t27+1782.0/7.0*t31+3636.0/7.0*t33-2214.0/7.0*t35+2214.0/7.0*t37-3636.0/7.0*t39-792.0/7.0*t73-792.0/7.0*t77+t250+t231;
86  t263 = -1080.0/49.0*xi-18.0/49.0*eta+1350.0/7.0*t12-360.0/7.0*t2-4800.0/7.0*t14+1380.0/7.0*t4+7110.0/7.0*t16-1710.0/7.0*t6-3510.0/7.0*t18+684.0/7.0*t8+t218-8100.0/7.0*t23;
87  t271 = 30.0/49.0+3240.0/7.0*t25+20808.0/7.0*t27+t223-1854.0*t31-15444.0/7.0*t33-13644.0/7.0*t35+10854.0/7.0*t37+9234.0/7.0*t39-t229-t230-t250;
88  t273 = 5540.0/77.0*t21;
89  t277 = 84320.0/77.0*t29;
90  t285 = t273+26920.0/77.0*t23-174970.0/77.0*t27-57780.0/77.0*t25+t277-92000.0/77.0*t39+6900.0/7.0*t18-141810.0/77.0*t16+81220.0/77.0*t14-16070.0/77.0*t12+760.0/77.0*xi-1614.0/7.0*t8;
91  t289 = 450.0/7.0*t80;
92  t298 = 467.0/1617.0+36060.0/77.0*t6-10480.0/33.0*t4+6920.0/77.0*t2+t289-5570.0/539.0*eta+27750.0/77.0*t35+192350.0/77.0*t33-135600.0/77.0*t37+133580.0/77.0*t31+2050.0/7.0*t77+3450.0/7.0*t73+2300.0/7.0*t106;
93  t312 = -510.0/77.0*t21+2720.0/77.0*t23+6530.0/77.0*t27+2720.0/77.0*t25-24580.0/77.0*t29-16980.0/77.0*t39-1086.0/7.0*t18+21540.0/77.0*t16-34640.0/231.0*t14+2080.0/77.0*t12-730.0/539.0*xi-1086.0/7.0*t8;
94  t323 = 450.0/7.0*t106;
95  t324 = -17.0/1617.0+21540.0/77.0*t6-34640.0/231.0*t4+2080.0/77.0*t2+t289-730.0/539.0*eta+20490.0/77.0*t35-16980.0/77.0*t33+20490.0/77.0*t37+6530.0/77.0*t31-100.0*t77+100.0*t73-t323;
96  t336 = t273-57780.0/77.0*t23+133580.0/77.0*t27+26920.0/77.0*t25+t277+192350.0/77.0*t39-1614.0/7.0*t18+36060.0/77.0*t16-10480.0/33.0*t14+6920.0/77.0*t12-5570.0/539.0*xi+6900.0/7.0*t8;
97  t348 = 467.0/1617.0-141810.0/77.0*t6+81220.0/77.0*t4-16070.0/77.0*t2-2300.0/7.0*t80+760.0/77.0*eta-135600.0/77.0*t35-92000.0/77.0*t33+27750.0/77.0*t37-174970.0/77.0*t31-3450.0/7.0*t77-2050.0/7.0*t73-t323;
98  t362 = -19860.0/11.0*t21+71424.0/11.0*t23-163344.0/11.0*t27+71424.0/11.0*t25-35928.0/11.0*t29+141300.0/11.0*t39+5580.0*t18-125700.0/11.0*t16+86400.0/11.0*t14-24900.0/11.0*t12+20400.0/77.0*xi+5580.0*t8;
99  t366 = 1500.0*t80;
100  t374 = 1500.0*t106;
101  t375 = -580.0/77.0-125700.0/11.0*t6+86400.0/11.0*t4-24900.0/11.0*t2-t366+20400.0/77.0*eta+13320.0/11.0*t35+141300.0/11.0*t33+13320.0/11.0*t37-163344.0/11.0*t31-1260.0*t77+1260.0*t73+t374;
102  t377 = 3504.0/11.0*t21;
103  t381 = 51984.0/11.0*t29;
104  t389 = t377+15192.0/11.0*t23-110676.0/11.0*t27-34704.0/11.0*t25+t381-62640.0/11.0*t39+3420.0*t18-66300.0/11.0*t16+33600.0/11.0*t14-5100.0/11.0*t12+600.0/77.0*xi-1440.0*t8;
105  t392 = 4800.0/11.0*t2;
106  t393 = 480.0*t80;
107  t399 = 720.0*t77;
108  t400 = 3240.0*t73;
109  t401 = 80.0/77.0+30720.0/11.0*t6-19200.0/11.0*t4+t392+t393-480.0/11.0*eta+19260.0/11.0*t35+125460.0/11.0*t33-81720.0/11.0*t37+83760.0/11.0*t31+t399+t400+t374;
110  t410 = 4800.0/11.0*t12;
111  t413 = t377-34704.0/11.0*t23+83760.0/11.0*t27+15192.0/11.0*t25+t381+125460.0/11.0*t39-1440.0*t18+30720.0/11.0*t16-19200.0/11.0*t14+t410-480.0/11.0*xi+3420.0*t8;
112  t422 = 3240.0*t77;
113  t423 = 720.0*t73;
114  t424 = 480.0*t106;
115  t425 = 80.0/77.0-66300.0/11.0*t6+33600.0/11.0*t4-5100.0/11.0*t2-t366+600.0/77.0*eta-81720.0/11.0*t35-62640.0/11.0*t33+19260.0/11.0*t37-110676.0/11.0*t31-t422-t423-t424;
116  t427 = 26328.0/11.0*t21;
117  t431 = 102456.0/11.0*t29;
118  t435 = 68800.0/11.0*t14;
119  t439 = t427-19680.0/11.0*t23-44688.0/11.0*t27-174912.0/11.0*t25+t431-333360.0/11.0*t39-3960.0*t18+93360.0/11.0*t16-t435+21600.0/11.0*t12-19080.0/77.0*xi-7200.0*t8;
120  t443 = 2160.0*t80;
121  t449 = 7200.0*t77;
122  t450 = 4680.0*t73;
123  t451 = 840.0*t106;
124  t452 = 580.0/77.0+158040.0/11.0*t6-104000.0/11.0*t4+28200.0/11.0*t2+t443-21720.0/77.0*eta+145080.0/11.0*t35+23040.0/11.0*t33-254880.0/11.0*t37+437904.0/11.0*t31+t449+t450-t451;
125  t464 = t427-174912.0/11.0*t23+437904.0/11.0*t27-19680.0/11.0*t25+t431+23040.0/11.0*t39-7200.0*t18+158040.0/11.0*t16-104000.0/11.0*t14+28200.0/11.0*t12-21720.0/77.0*xi-3960.0*t8;
126  t466 = 68800.0/11.0*t4;
127  t468 = 840.0*t80;
128  t474 = 4680.0*t77;
129  t475 = 7200.0*t73;
130  t476 = 2160.0*t106;
131  t477 = 580.0/77.0+93360.0/11.0*t6-t466+21600.0/11.0*t2+t468-19080.0/77.0*eta-254880.0/11.0*t35-333360.0/11.0*t33+145080.0/11.0*t37-44688.0/11.0*t31-t474-t475-t476;
132  t479 = 14064.0/11.0*t21;
133  t483 = 125640.0/11.0*t29;
134  t487 = 1600.0/11.0*t14;
135  t489 = -300.0/77.0-t479-18624.0/11.0*t23+171528.0/11.0*t27+114432.0/11.0*t25-t483+169560.0/11.0*t39-1080.0*t18+14160.0/11.0*t16+t487-t410+7320.0/77.0*xi;
136  t499 = 2160.0*t8-50520.0/11.0*t6+36800.0/11.0*t4-11400.0/11.0*t2-t393+9960.0/77.0*eta-112320.0/11.0*t35-159120.0/11.0*t33+259920.0/11.0*t37-261168.0/11.0*t31-t474-t475-t451;
137  t501 = 8904.0/11.0*t21;
138  t505 = 38304.0/11.0*t29;
139  t511 = 20.0/11.0+t501-30240.0/11.0*t23+94176.0/11.0*t27-46872.0/11.0*t25+t505-44280.0/11.0*t39-5760.0*t18+118440.0/11.0*t16-t435+15000.0/11.0*t12-8520.0/77.0*xi;
140  t514 = 1600.0/11.0*t4;
141  t521 = 720.0*t8-10920.0/11.0*t6+t514+1800.0/11.0*t2-t393-3240.0/77.0*eta+46080.0/11.0*t35-119520.0/11.0*t33-96480.0/11.0*t37+93648.0/11.0*t31+t422+t423-t476;
142  t532 = -300.0/77.0-t479+114432.0/11.0*t23-261168.0/11.0*t27-18624.0/11.0*t25-t483-159120.0/11.0*t39+2160.0*t18-50520.0/11.0*t16+36800.0/11.0*t14-11400.0/11.0*t12+9960.0/77.0*xi;
143  t540 = -1080.0*t8+14160.0/11.0*t6+t514-t392+t468+7320.0/77.0*eta+259920.0/11.0*t35+169560.0/11.0*t33-112320.0/11.0*t37+171528.0/11.0*t31+t449+t450+t424;
144  t550 = 20.0/11.0+t501-46872.0/11.0*t23+93648.0/11.0*t27-30240.0/11.0*t25+t505-119520.0/11.0*t39+720.0*t18-10920.0/11.0*t16+t487+1800.0/11.0*t12-3240.0/77.0*xi;
145  t559 = -5760.0*t8+118440.0/11.0*t6-t466+15000.0/11.0*t2+t443-8520.0/77.0*eta-96480.0/11.0*t35-44280.0/11.0*t33+46080.0/11.0*t37+94176.0/11.0*t31-t399-t400+t424;
146  t572 = -60.0/11.0-24864.0/11.0*t21+105504.0/11.0*t23-273420.0/11.0*t27+105504.0/11.0*t25-56700.0/11.0*t29+245700.0/11.0*t39+7560.0*t18-161280.0/11.0*t16+100800.0/11.0*t14-25200.0/11.0*t12+2520.0/11.0*xi;
147  t586 = 7560.0*t8-161280.0/11.0*t6+100800.0/11.0*t4-25200.0/11.0*t2-2520.0*t80+2520.0/11.0*eta+34020.0/11.0*t35+245700.0/11.0*t33+34020.0/11.0*t37-273420.0/11.0*t31-3780.0*t77+3780.0*t73+2520.0*t106;
148 
149  values[0] = 1.0-t1+t3-t5+t7-t9;
150  values[1] = -1.0+t11+t1-t13-t3+t15+t5-t17-t7+t19+t41;
151  values[2] = 1.0-t11+t13-t15+t17-t19;
152  values[3] = t55+t82;
153  values[4] = t95+t108;
154  values[5] = t120+t129;
155  values[6] = t143+t156;
156  values[7] = t170+t182;
157  values[8] = t194+t206;
158  values[9] = t219+t232;
159  values[10] = t251;
160  values[11] = t263+t271;
161  values[12] = t285+t298;
162  values[13] = t312+t324;
163  values[14] = t336+t348;
164  values[15] = t362+t375;
165  values[16] = t389+t401;
166  values[17] = t413+t425;
167  values[18] = t439+t452;
168  values[19] = t464+t477;
169  values[20] = t489+t499;
170  values[21] = t511+t521;
171  values[22] = t532+t540;
172  values[23] = t550+t559;
173  values[24] = t572+t586;
174 }
175 
176 // values of the derivatives in xi direction
177 static void N_T_P5_2D_DeriveXi(double xi, double eta, double *values)
178 {
179  double t1, t2, t3, t4, t5, t6, t7, t9, t11, t13, t15, t17, t19, t21;
180  double t23, t25, t27, t31, t38, t41, t45, t46, t47, t48, t49, t51, t52;
181  double t53, t54, t56, t57, t58, t59, t60, t61, t62, t63, t64, t65, t66;
182  double t67, t68, t69, t70, t73, t74, t75, t76, t77, t92, t99, t100, t101;
183  double t102, t103, t104, t105, t106, t109, t121, t124, t128, t135, t140;
184  double t149, t150, t151, t156, t172, t179, t180, t181, t182, t185, t199;
185  double t211, t215, t222, t225, t228, t229, t230, t231, t232, t233, t234;
186  double t235, t236, t237, t238, t239, t252, t257, t258, t259, t260, t263;
187  double t275, t278, t282, t289, t294, t303, t304, t305, t310, t326, t333;
188  double t334, t335, t336, t339, t353, t365, t381, t388, t389, t390, t391;
189  double t392, t393, t394, t397, t401, t408, t413, t414, t415, t416, t417;
190  double t418, t419, t420, t421, t422, t423, t424, t425, t438, t439, t440;
191  double t441, t442, t443, t444, t445, t446, t447, t448, t449, t450, t451;
192  double t452, t456, t463, t464, t468, t469, t470, t471, t472, t473, t474;
193  double t475, t476, t477, t478, t479, t480, t481, t482, t483, t484, t498;
194  double t499, t500, t501, t502, t503, t504, t505, t506, t507, t508, t509;
195  double t510, t511, t512, t513, t514, t518, t525, t526, t529, t530, t534;
196  double t541, t545, t546, t560, t561, t574, t575, t591, t607;
197 
198  t1 = 420.0*xi;
199  t2 = xi*xi;
200  t3 = 1680.0*t2;
201  t4 = t2*xi;
202  t5 = 2520.0*t4;
203  t6 = t2*t2;
204  t7 = 1260.0*t6;
205  t9 = xi*eta;
206  t11 = eta*eta;
207  t13 = t2*eta;
208  t15 = xi*t11;
209  t17 = t11*eta;
210  t19 = t4*eta;
211  t21 = t2*t11;
212  t23 = xi*t17;
213  t25 = t11*t11;
214  t27 = 30.0-t1+t3-t5+t7-420.0*eta+3360.0*t9+1680.0*t11-7560.0*t13-7560.0*t15-2520.0*t17+5040.0*t19+7560.0*t21+5040.0*t23+1260.0*t25;
215  t31 = 1080.0*t15;
216  t38 = 40.0*eta;
217  t41 = 4.0/7.0-1474.0*t17+882.0*t13+t31-648.0*t19-1566.0*t21-320.0*t9-516.0*t23+400.0*t11+1542.0*t25-t38-80.0*t2+120.0*t4;
218  t45 = 1.0-2.0*xi-eta;
219  t46 = xi-eta;
220  t47 = t45*t46;
221  t48 = t47*t13;
222  t49 = 66.0*t48;
223  t51 = 2.0*eta-1.0+xi;
224  t52 = t46*t51;
225  t53 = t52*t15;
226  t54 = 132.0*t53;
227  t56 = t47*t51*t11;
228  t57 = 66.0*t56;
229  t58 = t45*t51;
230  t59 = t58*t15;
231  t60 = 66.0*t59;
232  t61 = t47*t15;
233  t62 = 66.0*t61;
234  t63 = t52*t13;
235  t64 = 132.0*t63;
236  t65 = t47*t17;
237  t66 = 220.0*t65;
238  t67 = t58*t17;
239  t68 = 220.0*t67;
240  t69 = t52*t17;
241  t70 = 440.0*t69;
242  t73 = t47*t51*xi*eta;
243  t74 = 132.0*t73;
244  t75 = t58*t13;
245  t76 = 66.0*t75;
246  t77 = -60.0*t6+20.0*xi-t49+t54-t57-t60-t62+t64-t66-t68+t70-t74-t76;
247  t92 = 40.0/7.0+2146.0*t17-6438.0*t13+7992.0*t19+4194.0*t21+1120.0*t9-2796.0*t23-560.0*t11-1998.0*t25+1600.0*t2-3920.0*t4+2700.0*t6-200.0*xi+594.0*t48;
248  t99 = t47*t51*t2;
249  t100 = 660.0*t99;
250  t101 = t47*t4;
251  t102 = 220.0*t101;
252  t103 = t58*t4;
253  t104 = 220.0*t103;
254  t105 = t52*t4;
255  t106 = 440.0*t105;
256  t109 = -1188.0*t53+594.0*t56+594.0*t59+594.0*t61-1188.0*t63+t100+t102+t66+t68-t70+t104-t106+1188.0*t73+594.0*t75;
257  t121 = -240.0/7.0-294.0*t17+4422.0*t13-t31-6168.0*t19+774.0*t21-800.0*t9+1044.0*t23+160.0*t11+162.0*t25+t38-3200.0*t2+6320.0*t4;
258  t124 = -3900.0*t6+600.0*xi-t49+t54-t57-t60-t62+t64-t100-t102-t104+t106-t74-t76;
259  t128 = 4176.0/77.0*t15;
260  t135 = 8046.0/77.0*eta;
261  t140 = 18.0/77.0+216012.0/77.0*t17-251280.0/77.0*t13-t128+250560.0/77.0*t19+245160.0/77.0*t21+47232.0/77.0*t9-149400.0/77.0*t23-76176.0/77.0*t11-180900.0/77.0*t25+t135+57600.0/77.0*t2-122880.0/77.0*t4+7200.0/7.0*t6-8676.0/77.0*xi;
262  t149 = 1500.0/7.0*t65;
263  t150 = 1500.0/7.0*t67;
264  t151 = 3000.0/7.0*t69;
265  t156 = 720.0/7.0*t48-6480.0/7.0*t53+3240.0/7.0*t56+3240.0/7.0*t59+3240.0/7.0*t61-1440.0/7.0*t63+1440.0/7.0*t99+480.0/7.0*t101+t149+t150-t151+480.0/7.0*t103-960.0/7.0*t105+1440.0/7.0*t73+720.0/7.0*t75;
266  t172 = -4230.0/539.0+85728.0/77.0*t17+257184.0/77.0*t13-44568.0/77.0*t15-398880.0/77.0*t19+43200.0/77.0*t21-39360.0/77.0*t9+28800.0/77.0*t23-19680.0/77.0*t11-99720.0/77.0*t25+3690.0/77.0*eta-129840.0/77.0*t2+308760.0/77.0*t4-19080.0/7.0*t6+17460.0/77.0*xi;
267  t179 = 4500.0/7.0*t99;
268  t180 = 1500.0/7.0*t101;
269  t181 = 1500.0/7.0*t103;
270  t182 = 3000.0/7.0*t105;
271  t185 = -180.0*t48-360.0*t53+180.0*t56+180.0*t59+180.0*t61+360.0*t63-t179-t180+t149+t150-t151-t181+t182-360.0*t73-180.0*t75;
272  t199 = -16770.0/539.0-83760.0/77.0*t17+648036.0/77.0*t13-t128-723600.0/77.0*t19-224100.0/77.0*t21-152352.0/77.0*t9+163440.0/77.0*t23+23616.0/77.0*t11+62640.0/77.0*t25+t135-32880.0/11.0*t2+459240.0/77.0*t4-25920.0/7.0*t6+42540.0/77.0*xi;
273  t211 = -3240.0/7.0*t48+1440.0/7.0*t53-720.0/7.0*t56-720.0/7.0*t59-720.0/7.0*t61+6480.0/7.0*t63-t179-t180-480.0/7.0*t65-480.0/7.0*t67+960.0/7.0*t69-t181+t182-6480.0/7.0*t73-3240.0/7.0*t75;
274  t215 = 6480.0/7.0*t15;
275  t222 = 810.0/7.0*eta;
276  t225 = 18.0/49.0-20808.0/7.0*t17+5562.0*t13-t215-36936.0/7.0*t19-32562.0/7.0*t21-6480.0/7.0*t9+27288.0/7.0*t23+8100.0/7.0*t11+15444.0/7.0*t25-t222-4140.0/7.0*t2+6840.0/7.0*t4;
277  t228 = 3762.0/7.0*t48;
278  t229 = 7524.0/7.0*t53;
279  t230 = 3762.0/7.0*t56;
280  t231 = 3762.0/7.0*t59;
281  t232 = 3762.0/7.0*t61;
282  t233 = 7524.0/7.0*t63;
283  t234 = 990.0/7.0*t65;
284  t235 = 990.0/7.0*t67;
285  t236 = 1980.0/7.0*t69;
286  t237 = 7524.0/7.0*t73;
287  t238 = 3762.0/7.0*t75;
288  t239 = -3420.0/7.0*t6+720.0/7.0*xi-t228+t229-t230-t231-t232+t233-t234-t235+t236-t237-t238;
289  t252 = 180.0/49.0+1782.0/7.0*t17-5346.0/7.0*t13+14544.0/7.0*t19-6642.0/7.0*t21+4428.0/7.0*t23-3636.0/7.0*t25+7200.0/7.0*t2-t5+12150.0/7.0*t6-900.0/7.0*xi-792.0/7.0*t48+1584.0/7.0*t53;
290  t257 = 2970.0/7.0*t99;
291  t258 = 990.0/7.0*t101;
292  t259 = 990.0/7.0*t103;
293  t260 = 1980.0/7.0*t105;
294  t263 = -792.0/7.0*t56-792.0/7.0*t59-792.0/7.0*t61+1584.0/7.0*t63+t257+t258+t234+t235-t236+t259-t260-1584.0/7.0*t73-792.0/7.0*t75;
295  t275 = -1080.0/49.0-1854.0*t17+62424.0/7.0*t13+t215-61776.0/7.0*t19-40932.0/7.0*t21-16200.0/7.0*t9+21708.0/7.0*t23+3240.0/7.0*t11+9234.0/7.0*t25+t222-14400.0/7.0*t2+28440.0/7.0*t4;
296  t278 = -17550.0/7.0*t6+2700.0/7.0*xi-t228+t229-t230-t231-t232+t233-t257-t258-t259+t260-t237-t238;
297  t282 = 168640.0/77.0*t15;
298  t289 = 5540.0/77.0*eta;
299  t294 = 760.0/77.0+133580.0/77.0*t17-524910.0/77.0*t13+t282+769400.0/77.0*t19+83250.0/77.0*t21+53840.0/77.0*t9-271200.0/77.0*t23-57780.0/77.0*t11-92000.0/77.0*t25+t289+243660.0/77.0*t2-567240.0/77.0*t4+34500.0/7.0*t6-32140.0/77.0*xi;
300  t303 = 450.0/7.0*t65;
301  t304 = 450.0/7.0*t67;
302  t305 = 900.0/7.0*t69;
303  t310 = 3450.0/7.0*t48-4100.0/7.0*t53+2050.0/7.0*t56+2050.0/7.0*t59+2050.0/7.0*t61-6900.0/7.0*t63+6900.0/7.0*t99+2300.0/7.0*t101+t303+t304-t305+2300.0/7.0*t103-4600.0/7.0*t105+6900.0/7.0*t73+3450.0/7.0*t75;
304  t326 = -730.0/539.0+6530.0/77.0*t17+19590.0/77.0*t13-49160.0/77.0*t15-67920.0/77.0*t19+61470.0/77.0*t21+5440.0/77.0*t9+40980.0/77.0*t23+2720.0/77.0*t11-16980.0/77.0*t25-510.0/77.0*eta-34640.0/77.0*t2+86160.0/77.0*t4-5430.0/7.0*t6+4160.0/77.0*xi;
305  t333 = 1350.0/7.0*t99;
306  t334 = 450.0/7.0*t101;
307  t335 = 450.0/7.0*t103;
308  t336 = 900.0/7.0*t105;
309  t339 = 100.0*t48+200.0*t53-100.0*t56-100.0*t59-100.0*t61-200.0*t63-t333-t334+t303+t304-t305-t335+t336+200.0*t73+100.0*t75;
310  t353 = -5570.0/539.0-174970.0/77.0*t17+400740.0/77.0*t13+t282-368000.0/77.0*t19-406800.0/77.0*t21-115560.0/77.0*t9+55500.0/77.0*t23+26920.0/77.0*t11+192350.0/77.0*t25+t289-10480.0/11.0*t2+144240.0/77.0*t4-8070.0/7.0*t6+13840.0/77.0*xi;
311  t365 = -2050.0/7.0*t48+6900.0/7.0*t53-3450.0/7.0*t56-3450.0/7.0*t59-3450.0/7.0*t61+4100.0/7.0*t63-t333-t334-2300.0/7.0*t65-2300.0/7.0*t67+4600.0/7.0*t69-t335+t336-4100.0/7.0*t73-2050.0/7.0*t75;
312  t381 = 20400.0/77.0-163344.0/11.0*t17-490032.0/11.0*t13-71856.0/11.0*t15+565200.0/11.0*t19+39960.0/11.0*t21+142848.0/11.0*t9+26640.0/11.0*t23+71424.0/11.0*t11+141300.0/11.0*t25-19860.0/11.0*eta+259200.0/11.0*t2-502800.0/11.0*t4+27900.0*t6-49800.0/11.0*xi;
313  t388 = 4500.0*t99;
314  t389 = 1500.0*t101;
315  t390 = 1500.0*t65;
316  t391 = 1500.0*t67;
317  t392 = 3000.0*t69;
318  t393 = 1500.0*t103;
319  t394 = 3000.0*t105;
320  t397 = 1260.0*t48+2520.0*t53-1260.0*t56-1260.0*t59-1260.0*t61-2520.0*t63+t388+t389-t390-t391+t392+t393-t394+2520.0*t73+1260.0*t75;
321  t401 = 103968.0/11.0*t15;
322  t408 = 3504.0/11.0*eta;
323  t413 = 600.0/77.0+83760.0/11.0*t17-332028.0/11.0*t13+t401+501840.0/11.0*t19+57780.0/11.0*t21+30384.0/11.0*t9-163440.0/11.0*t23-34704.0/11.0*t11-62640.0/11.0*t25+t408+100800.0/11.0*t2-265200.0/11.0*t4+17100.0*t6-10200.0/11.0*xi;
324  t414 = 3240.0*t48;
325  t415 = 1440.0*t53;
326  t416 = 720.0*t56;
327  t417 = 720.0*t59;
328  t418 = 720.0*t61;
329  t419 = 6480.0*t63;
330  t420 = 480.0*t65;
331  t421 = 480.0*t67;
332  t422 = 960.0*t69;
333  t423 = 6480.0*t73;
334  t424 = 3240.0*t75;
335  t425 = t414-t415+t416+t417+t418-t419+t388+t389+t420+t421-t422+t393-t394+t423+t424;
336  t438 = 9600.0/11.0*xi;
337  t439 = -480.0/11.0-110676.0/11.0*t17+251280.0/11.0*t13+t401-250560.0/11.0*t19-245160.0/11.0*t21-69408.0/11.0*t9+38520.0/11.0*t23+15192.0/11.0*t11+125460.0/11.0*t25+t408-57600.0/11.0*t2+122880.0/11.0*t4-7200.0*t6+t438;
338  t440 = 720.0*t48;
339  t441 = 6480.0*t53;
340  t442 = 3240.0*t56;
341  t443 = 3240.0*t59;
342  t444 = 3240.0*t61;
343  t445 = 1440.0*t63;
344  t446 = 1440.0*t99;
345  t447 = 480.0*t101;
346  t448 = 480.0*t103;
347  t449 = 960.0*t105;
348  t450 = 1440.0*t73;
349  t451 = 720.0*t75;
350  t452 = -t440+t441-t442-t443-t444+t445-t446-t447-t390-t391+t392-t448+t449-t450-t451;
351  t456 = 204912.0/11.0*t15;
352  t463 = 26328.0/11.0*eta;
353  t464 = 206400.0/11.0*t2;
354  t468 = -19080.0/77.0+437904.0/11.0*t17-134064.0/11.0*t13+t456+92160.0/11.0*t19+435240.0/11.0*t21-39360.0/11.0*t9-509760.0/11.0*t23-174912.0/11.0*t11-333360.0/11.0*t25+t463-t464+373440.0/11.0*t4-19800.0*t6+43200.0/11.0*xi;
355  t469 = 4680.0*t48;
356  t470 = 14400.0*t53;
357  t471 = 7200.0*t56;
358  t472 = 7200.0*t59;
359  t473 = 7200.0*t61;
360  t474 = 9360.0*t63;
361  t475 = 2520.0*t99;
362  t476 = 840.0*t101;
363  t477 = 2160.0*t65;
364  t478 = 2160.0*t67;
365  t479 = 4320.0*t69;
366  t480 = 840.0*t103;
367  t481 = 1680.0*t105;
368  t482 = 9360.0*t73;
369  t483 = 4680.0*t75;
370  t484 = t469-t470+t471+t472+t473-t474-t475-t476+t477+t478-t479-t480+t481+t482+t483;
371  t498 = -21720.0/77.0-44688.0/11.0*t17+1313712.0/11.0*t13+t456-1333440.0/11.0*t19-764640.0/11.0*t21-349824.0/11.0*t9+290160.0/11.0*t23-19680.0/11.0*t11+23040.0/11.0*t25+t463-312000.0/11.0*t2+632160.0/11.0*t4-36000.0*t6+56400.0/11.0*xi;
372  t499 = 7200.0*t48;
373  t500 = 9360.0*t53;
374  t501 = 4680.0*t56;
375  t502 = 4680.0*t59;
376  t503 = 4680.0*t61;
377  t504 = 14400.0*t63;
378  t505 = 6480.0*t99;
379  t506 = 2160.0*t101;
380  t507 = 840.0*t65;
381  t508 = 840.0*t67;
382  t509 = 1680.0*t69;
383  t510 = 2160.0*t103;
384  t511 = 4320.0*t105;
385  t512 = 14400.0*t73;
386  t513 = 7200.0*t75;
387  t514 = -t499+t500-t501-t502-t503+t504-t505-t506+t507+t508-t509-t510+t511-t512-t513;
388  t518 = 251280.0/11.0*t15;
389  t525 = 14064.0/11.0*eta;
390  t526 = 4800.0/11.0*t2;
391  t529 = 7320.0/77.0-261168.0/11.0*t17+514584.0/11.0*t13-t518-636480.0/11.0*t19-336960.0/11.0*t21-37248.0/11.0*t9+519840.0/11.0*t23+114432.0/11.0*t11+169560.0/11.0*t25-t525+t526+56640.0/11.0*t4-5400.0*t6-t438;
392  t530 = -t499+t500-t501-t502-t503+t504-t475-t476-t420-t421+t422-t480+t481-t512-t513;
393  t534 = 76608.0/11.0*t15;
394  t541 = 8904.0/11.0*eta;
395  t545 = -8520.0/77.0+93648.0/11.0*t17+282528.0/11.0*t13+t534-478080.0/11.0*t19+138240.0/11.0*t21-60480.0/11.0*t9-192960.0/11.0*t23-46872.0/11.0*t11-44280.0/11.0*t25+t541-t464+473760.0/11.0*t4-28800.0*t6+30000.0/11.0*xi;
396  t546 = t440-t441+t442+t443+t444-t445-t505-t506-t420-t421+t422-t510+t511+t450+t451;
397  t560 = 9960.0/77.0+171528.0/11.0*t17-783504.0/11.0*t13-t518+678240.0/11.0*t19+779760.0/11.0*t21+228864.0/11.0*t9-224640.0/11.0*t23-18624.0/11.0*t11-159120.0/11.0*t25-t525+110400.0/11.0*t2-202080.0/11.0*t4+10800.0*t6-22800.0/11.0*xi;
398  t561 = t469-t470+t471+t472+t473-t474+t446+t447+t507+t508-t509+t448-t449+t482+t483;
399  t574 = -3240.0/77.0+94176.0/11.0*t17+280944.0/11.0*t13+t534-177120.0/11.0*t19-289440.0/11.0*t21-93744.0/11.0*t9+92160.0/11.0*t23-30240.0/11.0*t11-119520.0/11.0*t25+t541+t526-43680.0/11.0*t4+3600.0*t6+3600.0/11.0*xi;
400  t575 = -t414+t415-t416-t417-t418+t419+t446+t447+t477+t478-t479+t448-t449-t423-t424;
401  t591 = 2520.0/11.0-273420.0/11.0*t17-820260.0/11.0*t13-113400.0/11.0*t15+982800.0/11.0*t19+102060.0/11.0*t21+211008.0/11.0*t9+68040.0/11.0*t23+105504.0/11.0*t11+245700.0/11.0*t25-24864.0/11.0*eta+302400.0/11.0*t2-645120.0/11.0*t4+37800.0*t6-50400.0/11.0*xi;
402  t607 = 3780.0*t48+7560.0*t53-3780.0*t56-3780.0*t59-3780.0*t61-7560.0*t63+7560.0*t99+2520.0*t101-2520.0*t65-2520.0*t67+5040.0*t69+2520.0*t103-5040.0*t105+7560.0*t73+3780.0*t75;
403 
404  values[0] = 0.0;
405  values[1] = t27;
406  values[2] = -30.0+t1-t3+t5-t7;
407  values[3] = t41+t77;
408  values[4] = t92+t109;
409  values[5] = t121+t124;
410  values[6] = t140+t156;
411  values[7] = t172+t185;
412  values[8] = t199+t211;
413  values[9] = t225+t239;
414  values[10] = t252+t263;
415  values[11] = t275+t278;
416  values[12] = t294+t310;
417  values[13] = t326+t339;
418  values[14] = t353+t365;
419  values[15] = t381+t397;
420  values[16] = t413+t425;
421  values[17] = t439+t452;
422  values[18] = t468+t484;
423  values[19] = t498+t514;
424  values[20] = t529+t530;
425  values[21] = t545+t546;
426  values[22] = t560+t561;
427  values[23] = t574+t575;
428  values[24] = t591+t607;
429 }
430 
431 // values of the derivatives in eta direction
432 static void N_T_P5_2D_DeriveEta(double xi, double eta, double *values)
433 {
434  double t1, t2, t3, t4, t5, t6, t7, t10, t12, t14, t16, t18, t20, t22;
435  double t24, t26, t28, t29, t38, t39, t40, t41, t42, t44, t45, t46, t47;
436  double t49, t50, t51, t52, t53, t54, t55, t56, t57, t58, t60, t61, t63;
437  double t64, t65, t66, t67, t68, t69, t76, t77, t78, t79, t80, t95, t96;
438  double t97, t104, t105, t106, t107, t109, t118, t124, t126, t136, t141;
439  double t144, t145, t146, t147, t156, t172, t173, t181, t182, t184, t200;
440  double t211, t213, t221, t222, t223, t224, t225, t226, t227, t228, t230;
441  double t231, t232, t233, t238, t239, t240, t253, t254, t260, t261, t263;
442  double t272, t278, t280, t290, t295, t298, t299, t300, t301, t310, t326;
443  double t327, t335, t336, t338, t354, t365, t377, t382, t383, t385, t386;
444  double t387, t394, t395, t397, t399, t407, t408, t409, t410, t411, t412;
445  double t413, t414, t415, t416, t417, t418, t419, t424, t425, t426, t435;
446  double t436, t437, t438, t439, t440, t441, t442, t448, t449, t450, t451;
447  double t452, t454, t462, t463, t464, t465, t466, t467, t468, t469, t470;
448  double t472, t473, t474, t475, t480, t481, t482, t483, t484, t491, t493;
449  double t494, t495, t496, t497, t498, t499, t500, t501, t503, t504, t505;
450  double t510, t511, t512, t513, t514, t516, t524, t526, t531, t533, t539;
451  double t541, t543, t548, t556, t561, t569, t575, t591, t607;
452 
453  t1 = 420.0*eta;
454  t2 = eta*eta;
455  t3 = 1680.0*t2;
456  t4 = t2*eta;
457  t5 = 2520.0*t4;
458  t6 = t2*t2;
459  t7 = 1260.0*t6;
460  t10 = xi*xi;
461  t12 = t10*xi;
462  t14 = t10*t10;
463  t16 = xi*eta;
464  t18 = t10*eta;
465  t20 = xi*t2;
466  t22 = t12*eta;
467  t24 = t10*t2;
468  t26 = xi*t4;
469  t28 = 30.0-420.0*xi+1680.0*t10-2520.0*t12+1260.0*t14-t1+3360.0*t16+t3-7560.0*t18-7560.0*t20-t5+5040.0*t22+7560.0*t24+5040.0*t26+t7;
470  t29 = 1080.0*t18;
471  t38 = 1.0-2.0*xi-eta;
472  t39 = xi-eta;
473  t40 = t38*t39;
474  t41 = t40*t18;
475  t42 = 132.0*t41;
476  t44 = 2.0*eta-1.0+xi;
477  t45 = t39*t44;
478  t46 = t45*t20;
479  t47 = 66.0*t46;
480  t49 = t40*t44*t2;
481  t50 = 660.0*t49;
482  t51 = t38*t44;
483  t52 = t51*t20;
484  t53 = 66.0*t52;
485  t54 = 240.0/7.0+t29-4422.0*t20-1044.0*t22-774.0*t24+6168.0*t26+800.0*t16+3200.0*t2-6320.0*t4-t42+t47-t50+t53;
486  t55 = t40*t20;
487  t56 = 132.0*t55;
488  t57 = t45*t18;
489  t58 = 66.0*t57;
490  t60 = t40*t44*t10;
491  t61 = 66.0*t60;
492  t63 = t40*t4;
493  t64 = 440.0*t63;
494  t65 = t51*t4;
495  t66 = 220.0*t65;
496  t67 = t45*t4;
497  t68 = 220.0*t67;
498  t69 = 40.0*xi;
499  t76 = t40*t44*xi*eta;
500  t77 = 132.0*t76;
501  t78 = t51*t18;
502  t79 = 66.0*t78;
503  t80 = -t56+t58-t61-600.0*eta-t64+t66+t68-t69+3900.0*t6-160.0*t10+294.0*t12-162.0*t14-t77+t79;
504  t95 = 6438.0*t20+2796.0*t22-4194.0*t24-7992.0*t26-1120.0*t16-1600.0*t2+3920.0*t4+1188.0*t41-594.0*t46+t50-594.0*t52+1188.0*t55-594.0*t57+594.0*t60;
505  t96 = t40*t12;
506  t97 = 440.0*t96;
507  t104 = t51*t12;
508  t105 = 220.0*t104;
509  t106 = t45*t12;
510  t107 = 220.0*t106;
511  t109 = -40.0/7.0+t97+200.0*eta+t64-t66-t68-2700.0*t6+560.0*t10-2146.0*t12+1998.0*t14+1188.0*t76-t105-t107-594.0*t78;
512  t118 = -t29-882.0*t20+516.0*t22+1566.0*t24+648.0*t26+320.0*t16+80.0*t2-120.0*t4-t42+t47+t53-t56+t58;
513  t124 = -4.0/7.0-t61-t97-20.0*eta+t69+60.0*t6-400.0*t10+1474.0*t12-1542.0*t14-t77+t105+t107+t79;
514  t126 = 4176.0/77.0*t18;
515  t136 = 4500.0/7.0*t49;
516  t141 = -t126+648036.0/77.0*t20+163440.0/77.0*t22-224100.0/77.0*t24-723600.0/77.0*t26-152352.0/77.0*t16-32880.0/11.0*t2+459240.0/77.0*t4+1440.0/7.0*t41-3240.0/7.0*t46+t136-3240.0/7.0*t52+6480.0/7.0*t55-720.0/7.0*t57+720.0/7.0*t60;
517  t144 = 3000.0/7.0*t63;
518  t145 = 1500.0/7.0*t65;
519  t146 = 1500.0/7.0*t67;
520  t147 = 8046.0/77.0*xi;
521  t156 = -16770.0/539.0+960.0/7.0*t96+42540.0/77.0*eta+t144-t145-t146+t147-25920.0/7.0*t6+23616.0/77.0*t10-83760.0/77.0*t12+62640.0/77.0*t14+6480.0/7.0*t76-480.0/7.0*t104-480.0/7.0*t106-720.0/7.0*t78;
522  t172 = -44568.0/77.0*t18+257184.0/77.0*t20+28800.0/77.0*t22+43200.0/77.0*t24-398880.0/77.0*t26-39360.0/77.0*t16-129840.0/77.0*t2+308760.0/77.0*t4-360.0*t41-180.0*t46+t136-180.0*t52+360.0*t55+180.0*t57-180.0*t60;
523  t173 = 3000.0/7.0*t96;
524  t181 = 1500.0/7.0*t104;
525  t182 = 1500.0/7.0*t106;
526  t184 = -4230.0/539.0-t173+17460.0/77.0*eta+t144-t145-t146+3690.0/77.0*xi-19080.0/7.0*t6-19680.0/77.0*t10+85728.0/77.0*t12-99720.0/77.0*t14+360.0*t76+t181+t182+180.0*t78;
527  t200 = -t126-251280.0/77.0*t20-149400.0/77.0*t22+245160.0/77.0*t24+250560.0/77.0*t26+47232.0/77.0*t16+57600.0/77.0*t2-122880.0/77.0*t4-6480.0/7.0*t41+720.0/7.0*t46-1440.0/7.0*t49+720.0/7.0*t52-1440.0/7.0*t55+3240.0/7.0*t57-3240.0/7.0*t60;
528  t211 = 18.0/77.0-t173-8676.0/77.0*eta-960.0/7.0*t63+480.0/7.0*t65+480.0/7.0*t67+t147+7200.0/7.0*t6-76176.0/77.0*t10+216012.0/77.0*t12-180900.0/77.0*t14-1440.0/7.0*t76+t181+t182+3240.0/7.0*t78;
529  t213 = 6480.0/7.0*t18;
530  t221 = 7524.0/7.0*t41;
531  t222 = 3762.0/7.0*t46;
532  t223 = 2970.0/7.0*t49;
533  t224 = 3762.0/7.0*t52;
534  t225 = 7524.0/7.0*t55;
535  t226 = -t213-62424.0/7.0*t20-21708.0/7.0*t22+40932.0/7.0*t24+61776.0/7.0*t26+16200.0/7.0*t16+14400.0/7.0*t2-28440.0/7.0*t4-t221+t222-t223+t224-t225;
536  t227 = 3762.0/7.0*t57;
537  t228 = 3762.0/7.0*t60;
538  t230 = 1980.0/7.0*t63;
539  t231 = 990.0/7.0*t65;
540  t232 = 990.0/7.0*t67;
541  t233 = 810.0/7.0*xi;
542  t238 = 7524.0/7.0*t76;
543  t239 = 3762.0/7.0*t78;
544  t240 = 1080.0/49.0+t227-t228-2700.0/7.0*eta-t230+t231+t232-t233+17550.0/7.0*t6-3240.0/7.0*t10+1854.0*t12-9234.0/7.0*t14-t238+t239;
545  t253 = 5346.0/7.0*t20-4428.0/7.0*t22+6642.0/7.0*t24-14544.0/7.0*t26-7200.0/7.0*t2+t5-1584.0/7.0*t41+792.0/7.0*t46+t223+792.0/7.0*t52-1584.0/7.0*t55+792.0/7.0*t57-792.0/7.0*t60;
546  t254 = 1980.0/7.0*t96;
547  t260 = 990.0/7.0*t104;
548  t261 = 990.0/7.0*t106;
549  t263 = -180.0/49.0+t254+900.0/7.0*eta+t230-t231-t232-12150.0/7.0*t6-1782.0/7.0*t12+3636.0/7.0*t14-1584.0/7.0*t76-t260-t261+792.0/7.0*t78;
550  t272 = t213-5562.0*t20-27288.0/7.0*t22+32562.0/7.0*t24+36936.0/7.0*t26+6480.0/7.0*t16+4140.0/7.0*t2-6840.0/7.0*t4-t221+t222+t224-t225+t227;
551  t278 = -18.0/49.0-t228-t254-720.0/7.0*eta+t233+3420.0/7.0*t6-8100.0/7.0*t10+20808.0/7.0*t12-15444.0/7.0*t14-t238+t260+t261+t239;
552  t280 = 168640.0/77.0*t18;
553  t290 = 1350.0/7.0*t49;
554  t295 = t280+400740.0/77.0*t20+55500.0/77.0*t22-406800.0/77.0*t24-368000.0/77.0*t26-115560.0/77.0*t16-10480.0/11.0*t2+144240.0/77.0*t4+6900.0/7.0*t41-2050.0/7.0*t46+t290-2050.0/7.0*t52+4100.0/7.0*t55-3450.0/7.0*t57+3450.0/7.0*t60;
555  t298 = 900.0/7.0*t63;
556  t299 = 450.0/7.0*t65;
557  t300 = 450.0/7.0*t67;
558  t301 = 5540.0/77.0*xi;
559  t310 = -5570.0/539.0+4600.0/7.0*t96+13840.0/77.0*eta+t298-t299-t300+t301-8070.0/7.0*t6+26920.0/77.0*t10-174970.0/77.0*t12+192350.0/77.0*t14+4100.0/7.0*t76-2300.0/7.0*t104-2300.0/7.0*t106-3450.0/7.0*t78;
560  t326 = -49160.0/77.0*t18+19590.0/77.0*t20+40980.0/77.0*t22+61470.0/77.0*t24-67920.0/77.0*t26+5440.0/77.0*t16-34640.0/77.0*t2+86160.0/77.0*t4+200.0*t41+100.0*t46+t290+100.0*t52-200.0*t55-100.0*t57+100.0*t60;
561  t327 = 900.0/7.0*t96;
562  t335 = 450.0/7.0*t104;
563  t336 = 450.0/7.0*t106;
564  t338 = -730.0/539.0-t327+4160.0/77.0*eta+t298-t299-t300-510.0/77.0*xi-5430.0/7.0*t6+2720.0/77.0*t10+6530.0/77.0*t12-16980.0/77.0*t14-200.0*t76+t335+t336-100.0*t78;
565  t354 = t280-524910.0/77.0*t20-271200.0/77.0*t22+83250.0/77.0*t24+769400.0/77.0*t26+53840.0/77.0*t16+243660.0/77.0*t2-567240.0/77.0*t4-4100.0/7.0*t41+3450.0/7.0*t46-6900.0/7.0*t49+3450.0/7.0*t52-6900.0/7.0*t55+2050.0/7.0*t57-2050.0/7.0*t60;
566  t365 = 760.0/77.0-t327-32140.0/77.0*eta-4600.0/7.0*t63+2300.0/7.0*t65+2300.0/7.0*t67+t301+34500.0/7.0*t6-57780.0/77.0*t10+133580.0/77.0*t12-92000.0/77.0*t14-6900.0/7.0*t76+t335+t336+2050.0/7.0*t78;
567  t377 = 4500.0*t49;
568  t382 = -71856.0/11.0*t18-490032.0/11.0*t20+26640.0/11.0*t22+39960.0/11.0*t24+565200.0/11.0*t26+142848.0/11.0*t16+259200.0/11.0*t2-502800.0/11.0*t4+2520.0*t41+1260.0*t46-t377+1260.0*t52-2520.0*t55-1260.0*t57+1260.0*t60;
569  t383 = 3000.0*t96;
570  t385 = 3000.0*t63;
571  t386 = 1500.0*t65;
572  t387 = 1500.0*t67;
573  t394 = 1500.0*t104;
574  t395 = 1500.0*t106;
575  t397 = 20400.0/77.0+t383-49800.0/11.0*eta-t385+t386+t387-19860.0/11.0*xi+27900.0*t6+71424.0/11.0*t10-163344.0/11.0*t12+141300.0/11.0*t14-2520.0*t76-t394-t395-1260.0*t78;
576  t399 = 103968.0/11.0*t18;
577  t407 = 6480.0*t41;
578  t408 = 720.0*t46;
579  t409 = 1440.0*t49;
580  t410 = 720.0*t52;
581  t411 = 1440.0*t55;
582  t412 = 3240.0*t57;
583  t413 = 3240.0*t60;
584  t414 = t399+251280.0/11.0*t20+38520.0/11.0*t22-245160.0/11.0*t24-250560.0/11.0*t26-69408.0/11.0*t16-57600.0/11.0*t2+122880.0/11.0*t4+t407-t408+t409-t410+t411-t412+t413;
585  t415 = 9600.0/11.0*eta;
586  t416 = 960.0*t63;
587  t417 = 480.0*t65;
588  t418 = 480.0*t67;
589  t419 = 3504.0/11.0*xi;
590  t424 = 1440.0*t76;
591  t425 = 3240.0*t78;
592  t426 = -480.0/11.0+t383+t415+t416-t417-t418+t419-7200.0*t6+15192.0/11.0*t10-110676.0/11.0*t12+125460.0/11.0*t14+t424-t394-t395-t425;
593  t435 = 1440.0*t41;
594  t436 = 3240.0*t46;
595  t437 = 3240.0*t52;
596  t438 = 6480.0*t55;
597  t439 = 720.0*t57;
598  t440 = 720.0*t60;
599  t441 = t399-332028.0/11.0*t20-163440.0/11.0*t22+57780.0/11.0*t24+501840.0/11.0*t26+30384.0/11.0*t16+100800.0/11.0*t2-265200.0/11.0*t4-t435+t436-t377+t437-t438+t439-t440;
600  t442 = 960.0*t96;
601  t448 = 6480.0*t76;
602  t449 = 480.0*t104;
603  t450 = 480.0*t106;
604  t451 = 720.0*t78;
605  t452 = 600.0/77.0-t442-10200.0/11.0*eta-t385+t386+t387+t419+17100.0*t6-34704.0/11.0*t10+83760.0/11.0*t12-62640.0/11.0*t14-t448+t449+t450+t451;
606  t454 = 204912.0/11.0*t18;
607  t462 = 9360.0*t41;
608  t463 = 7200.0*t46;
609  t464 = 6480.0*t49;
610  t465 = 7200.0*t52;
611  t466 = 14400.0*t55;
612  t467 = 4680.0*t57;
613  t468 = 4680.0*t60;
614  t469 = t454+1313712.0/11.0*t20+290160.0/11.0*t22-764640.0/11.0*t24-1333440.0/11.0*t26-349824.0/11.0*t16-312000.0/11.0*t2+632160.0/11.0*t4+t462-t463+t464-t465+t466-t467+t468;
615  t470 = 1680.0*t96;
616  t472 = 4320.0*t63;
617  t473 = 2160.0*t65;
618  t474 = 2160.0*t67;
619  t475 = 26328.0/11.0*xi;
620  t480 = 14400.0*t76;
621  t481 = 840.0*t104;
622  t482 = 840.0*t106;
623  t483 = 4680.0*t78;
624  t484 = -21720.0/77.0-t470+56400.0/11.0*eta+t472-t473-t474+t475-36000.0*t6-19680.0/11.0*t10-44688.0/11.0*t12+23040.0/11.0*t14+t480+t481+t482-t483;
625  t491 = 206400.0/11.0*t2;
626  t493 = 14400.0*t41;
627  t494 = 4680.0*t46;
628  t495 = 2520.0*t49;
629  t496 = 4680.0*t52;
630  t497 = 9360.0*t55;
631  t498 = 7200.0*t57;
632  t499 = 7200.0*t60;
633  t500 = t454-134064.0/11.0*t20-509760.0/11.0*t22+435240.0/11.0*t24+92160.0/11.0*t26-39360.0/11.0*t16-t491+373440.0/11.0*t4-t493+t494+t495+t496-t497+t498-t499;
634  t501 = 4320.0*t96;
635  t503 = 1680.0*t63;
636  t504 = 840.0*t65;
637  t505 = 840.0*t67;
638  t510 = 9360.0*t76;
639  t511 = 2160.0*t104;
640  t512 = 2160.0*t106;
641  t513 = 7200.0*t78;
642  t514 = -19080.0/77.0-t501+43200.0/11.0*eta+t503-t504-t505+t475-19800.0*t6-174912.0/11.0*t10+437904.0/11.0*t12-333360.0/11.0*t14-t510+t511+t512+t513;
643  t516 = 251280.0/11.0*t18;
644  t524 = 9960.0/77.0-t516-783504.0/11.0*t20-224640.0/11.0*t22+779760.0/11.0*t24+678240.0/11.0*t26+228864.0/11.0*t16+110400.0/11.0*t2-202080.0/11.0*t4-t493+t494-t409+t496-t497+t498;
645  t526 = 14064.0/11.0*xi;
646  t531 = -t499-t470-22800.0/11.0*eta-t416+t417+t418-t526+10800.0*t6-18624.0/11.0*t10+171528.0/11.0*t12-159120.0/11.0*t14-t510+t481+t482+t513;
647  t533 = 76608.0/11.0*t18;
648  t539 = 4800.0/11.0*t2;
649  t541 = -3240.0/77.0+t533+280944.0/11.0*t20+92160.0/11.0*t22-289440.0/11.0*t24-177120.0/11.0*t26-93744.0/11.0*t16+t539-43680.0/11.0*t4+t435-t436-t409-t437+t438-t439;
650  t543 = 8904.0/11.0*xi;
651  t548 = t440-t501+3600.0/11.0*eta-t416+t417+t418+t543+3600.0*t6-30240.0/11.0*t10+94176.0/11.0*t12-119520.0/11.0*t14+t448+t511+t512-t451;
652  t556 = 7320.0/77.0-t516+514584.0/11.0*t20+519840.0/11.0*t22-336960.0/11.0*t24-636480.0/11.0*t26-37248.0/11.0*t16+t539+56640.0/11.0*t4+t462-t463+t495-t465+t466-t467;
653  t561 = t468+t442-t415+t503-t504-t505-t526-5400.0*t6+114432.0/11.0*t10-261168.0/11.0*t12+169560.0/11.0*t14+t480-t449-t450-t483;
654  t569 = -8520.0/77.0+t533+282528.0/11.0*t20-192960.0/11.0*t22+138240.0/11.0*t24-478080.0/11.0*t26-60480.0/11.0*t16-t491+473760.0/11.0*t4-t407+t408+t464+t410-t411+t412;
655  t575 = -t413+t442+30000.0/11.0*eta+t472-t473-t474+t543-28800.0*t6-46872.0/11.0*t10+93648.0/11.0*t12-44280.0/11.0*t14-t424-t449-t450+t425;
656  t591 = 2520.0/11.0-113400.0/11.0*t18-820260.0/11.0*t20+68040.0/11.0*t22+102060.0/11.0*t24+982800.0/11.0*t26+211008.0/11.0*t16+302400.0/11.0*t2-645120.0/11.0*t4+7560.0*t41+3780.0*t46-7560.0*t49+3780.0*t52-7560.0*t55-3780.0*t57;
657  t607 = 3780.0*t60+5040.0*t96-50400.0/11.0*eta-5040.0*t63+2520.0*t65+2520.0*t67-24864.0/11.0*xi+37800.0*t6+105504.0/11.0*t10-273420.0/11.0*t12+245700.0/11.0*t14-7560.0*t76-2520.0*t104-2520.0*t106-3780.0*t78;
658 
659  values[0] = -30.0+t1-t3+t5-t7;
660  values[1] = t28;
661  values[2] = 0.0;
662  values[3] = t54+t80;
663  values[4] = t95+t109;
664  values[5] = t118+t124;
665  values[6] = t141+t156;
666  values[7] = t172+t184;
667  values[8] = t200+t211;
668  values[9] = t226+t240;
669  values[10] = t253+t263;
670  values[11] = t272+t278;
671  values[12] = t295+t310;
672  values[13] = t326+t338;
673  values[14] = t354+t365;
674  values[15] = t382+t397;
675  values[16] = t414+t426;
676  values[17] = t441+t452;
677  values[18] = t469+t484;
678  values[19] = t500+t514;
679  values[20] = t524+t531;
680  values[21] = t541+t548;
681  values[22] = t556+t561;
682  values[23] = t569+t575;
683  values[24] = t591+t607;
684 }
685 
686 // values of derivatives in xi-xi direction
687 static void N_T_P5_2D_DeriveXiXi(double xi, double eta, double *values)
688 {
689  double t1, t2, t3, t4, t5, t7, t9, t11, t13, t15, t26, t28, t29, t30;
690  double t31, t32, t33, t34, t35, t37, t38, t40, t41, t42, t43, t44, t45;
691  double t46, t47, t48, t49, t50, t51, t53, t54, t56, t57, t59, t60, t61;
692  double t62, t63, t66, t67, t68, t70, t71, t72, t73, t74, t75, t76, t92;
693  double t96, t97, t100, t101, t102, t103, t106, t107, t108, t109, t110;
694  double t111, t112, t113, t114, t122, t125, t139, t142, t149, t151, t155;
695  double t160, t176, t180, t184, t185, t188, t189, t190, t191, t192, t209;
696  double t219, t227, t228, t229, t230, t231, t232, t233, t234, t235, t236;
697  double t237, t238, t239, t242, t243, t244, t245, t246, t247, t261, t265;
698  double t268, t269, t272, t273, t274, t275, t276, t284, t287, t301, t304;
699  double t311, t313, t317, t322, t338, t342, t346, t347, t350, t351, t352;
700  double t353, t354, t371, t381, t395, t399, t402, t406, t407, t408, t411;
701  double t412, t413, t414, t415, t416, t424, t425, t426, t427, t428, t429;
702  double t430, t431, t432, t433, t434, t435, t436, t439, t440, t441, t442;
703  double t443, t444, t452, t453, t454, t455, t456, t457, t458, t459, t460;
704  double t461, t462, t463, t464, t467, t468, t469, t470, t471, t472, t473;
705  double t474, t475, t480, t483, t484, t485, t486, t487, t488, t489, t490;
706  double t491, t492, t493, t494, t495, t496, t499, t500, t501, t502, t503;
707  double t504, t505, t506, t507, t508, t509, t510, t518, t519, t520, t521;
708  double t522, t523, t524, t525, t526, t527, t528, t529, t530, t531, t534;
709  double t535, t536, t537, t538, t539, t540, t541, t542, t543, t544, t549;
710  double t552, t555, t556, t563, t566, t567, t575, t578, t585, t588, t605;
711  double t623;
712 
713  t1 = 3360.0*xi;
714  t2 = xi*xi;
715  t3 = 7560.0*t2;
716  t4 = xi*t2;
717  t5 = 5040.0*t4;
718  t7 = xi*eta;
719  t9 = eta*eta;
720  t11 = t2*eta;
721  t13 = xi*t9;
722  t15 = t9*eta;
723  t26 = 1.0-2.0*xi-eta;
724  t28 = 2.0*eta-1.0+xi;
725  t29 = t26*t28;
726  t30 = t29*t7;
727  t31 = 264.0*t30;
728  t32 = t28*xi;
729  t33 = t32*t9;
730  t34 = 264.0*t33;
731  t35 = xi-eta;
732  t37 = t35*xi*t9;
733  t38 = 264.0*t37;
734  t40 = t26*xi*t9;
735  t41 = 132.0*t40;
736  t42 = t26*t35;
737  t43 = t42*t9;
738  t44 = 132.0*t43;
739  t45 = t29*t9;
740  t46 = 132.0*t45;
741  t47 = 20.0-320.0*eta-240.0*t4+360.0*t2-160.0*xi+1764.0*t7-1944.0*t11-t31+t34+t38-t41-t44-t46;
742  t48 = t28*t15;
743  t49 = 880.0*t48;
744  t50 = t42*t7;
745  t51 = 264.0*t50;
746  t53 = t42*t28*eta;
747  t54 = 132.0*t53;
748  t56 = t28*t2*eta;
749  t57 = 264.0*t56;
750  t59 = t26*t2*eta;
751  t60 = 132.0*t59;
752  t61 = t35*t28;
753  t62 = t61*t9;
754  t63 = 264.0*t62;
755  t66 = 1080.0*t9;
756  t67 = t26*t15;
757  t68 = 440.0*t67;
758  t70 = t35*t2*eta;
759  t71 = 264.0*t70;
760  t72 = t61*t7;
761  t73 = 528.0*t72;
762  t74 = t35*t15;
763  t75 = 880.0*t74;
764  t76 = t49-t51-t54+t57-t60+t63-516.0*t15-3132.0*t13+t66-t68+t71+t73+t75;
765  t92 = -200.0+1120.0*eta+10800.0*t4-11760.0*t2+3200.0*xi-12876.0*t7+23976.0*t11+2376.0*t30-2376.0*t33-2376.0*t37+1188.0*t40+1188.0*t43+1188.0*t45-t49+2376.0*t50+1188.0*t53;
766  t96 = t29*t2;
767  t97 = 1320.0*t96;
768  t100 = t35*t4;
769  t101 = 880.0*t100;
770  t102 = t42*t2;
771  t103 = 1320.0*t102;
772  t106 = t42*t32;
773  t107 = 1320.0*t106;
774  t108 = t28*t4;
775  t109 = 880.0*t108;
776  t110 = t61*t2;
777  t111 = 2640.0*t110;
778  t112 = t26*t4;
779  t113 = 440.0*t112;
780  t114 = -2376.0*t56+1188.0*t59-2376.0*t62+t97-2796.0*t15+8388.0*t13-t101+t68+t103-2376.0*t70-4752.0*t72-t75+t107-t109-t111+t113;
781  t122 = 600.0-800.0*eta-15600.0*t4+18960.0*t2-6400.0*xi+8844.0*t7-18504.0*t11-t31+t34+t38-t41-t44-t46-t51-t54;
782  t125 = t57-t60+t63-t97+1044.0*t15+1548.0*t13-t66+t101-t103+t71+t73-t107+t109+t111-t113;
783  t139 = 6000.0/7.0*t48;
784  t142 = -8676.0/77.0+47232.0/77.0*eta+28800.0/7.0*t4-368640.0/77.0*t2+115200.0/77.0*xi-502560.0/77.0*t7+751680.0/77.0*t11+2880.0/7.0*t30-12960.0/7.0*t33-12960.0/7.0*t37+6480.0/7.0*t40+6480.0/7.0*t43+6480.0/7.0*t45-t139+2880.0/7.0*t50+1440.0/7.0*t53;
785  t149 = 4176.0/77.0*t9;
786  t151 = 3000.0/7.0*t67;
787  t155 = 6000.0/7.0*t74;
788  t160 = -2880.0/7.0*t56+1440.0/7.0*t59-12960.0/7.0*t62+2880.0/7.0*t96-149400.0/77.0*t15+490320.0/77.0*t13-t149-1920.0/7.0*t100+t151+2880.0/7.0*t102-2880.0/7.0*t70-5760.0/7.0*t72-t155+2880.0/7.0*t106-1920.0/7.0*t108-5760.0/7.0*t110+960.0/7.0*t112;
789  t176 = 17460.0/77.0-39360.0/77.0*eta-76320.0/7.0*t4+926280.0/77.0*t2-259680.0/77.0*xi+514368.0/77.0*t7-1196640.0/77.0*t11-720.0*t30-720.0*t33-720.0*t37+360.0*t40+360.0*t43+360.0*t45-t139-720.0*t50-360.0*t53;
790  t180 = 9000.0/7.0*t96;
791  t184 = 6000.0/7.0*t100;
792  t185 = 9000.0/7.0*t102;
793  t188 = 9000.0/7.0*t106;
794  t189 = 6000.0/7.0*t108;
795  t190 = 18000.0/7.0*t110;
796  t191 = 3000.0/7.0*t112;
797  t192 = 720.0*t56-360.0*t59-720.0*t62-t180+28800.0/77.0*t15+86400.0/77.0*t13-44568.0/77.0*t9+t184+t151-t185+720.0*t70+1440.0*t72-t155-t188+t189+t190-t191;
798  t209 = 42540.0/77.0-152352.0/77.0*eta-103680.0/7.0*t4+1377720.0/77.0*t2-65760.0/11.0*xi+1296072.0/77.0*t7-2170800.0/77.0*t11-12960.0/7.0*t30+2880.0/7.0*t33+2880.0/7.0*t37-1440.0/7.0*t40-1440.0/7.0*t43-1440.0/7.0*t45+1920.0/7.0*t48-12960.0/7.0*t50-6480.0/7.0*t53;
799  t219 = 12960.0/7.0*t56-6480.0/7.0*t59+2880.0/7.0*t62-t180+163440.0/77.0*t15-448200.0/77.0*t13-t149+t184-960.0/7.0*t67-t185+12960.0/7.0*t70+25920.0/7.0*t72+1920.0/7.0*t74-t188+t189+t190-t191;
800  t227 = 15048.0/7.0*t30;
801  t228 = 15048.0/7.0*t33;
802  t229 = 15048.0/7.0*t37;
803  t230 = 7524.0/7.0*t40;
804  t231 = 7524.0/7.0*t43;
805  t232 = 7524.0/7.0*t45;
806  t233 = 720.0/7.0-6480.0/7.0*eta-13680.0/7.0*t4+20520.0/7.0*t2-8280.0/7.0*xi+11124.0*t7-110808.0/7.0*t11-t227+t228+t229-t230-t231-t232;
807  t234 = 3960.0/7.0*t48;
808  t235 = 15048.0/7.0*t50;
809  t236 = 7524.0/7.0*t53;
810  t237 = 15048.0/7.0*t56;
811  t238 = 7524.0/7.0*t59;
812  t239 = 15048.0/7.0*t62;
813  t242 = 6480.0/7.0*t9;
814  t243 = 1980.0/7.0*t67;
815  t244 = 15048.0/7.0*t70;
816  t245 = 30096.0/7.0*t72;
817  t246 = 3960.0/7.0*t74;
818  t247 = t234-t235-t236+t237-t238+t239+27288.0/7.0*t15-65124.0/7.0*t13-t242-t243+t244+t245+t246;
819  t261 = -900.0/7.0+48600.0/7.0*t4-t3+14400.0/7.0*xi-10692.0/7.0*t7+43632.0/7.0*t11-3168.0/7.0*t30+3168.0/7.0*t33+3168.0/7.0*t37-1584.0/7.0*t40-1584.0/7.0*t43-1584.0/7.0*t45-t234-3168.0/7.0*t50-1584.0/7.0*t53;
820  t265 = 5940.0/7.0*t96;
821  t268 = 3960.0/7.0*t100;
822  t269 = 5940.0/7.0*t102;
823  t272 = 5940.0/7.0*t106;
824  t273 = 3960.0/7.0*t108;
825  t274 = 11880.0/7.0*t110;
826  t275 = 1980.0/7.0*t112;
827  t276 = 3168.0/7.0*t56-1584.0/7.0*t59+3168.0/7.0*t62+t265+4428.0/7.0*t15-13284.0/7.0*t13-t268+t243+t269+3168.0/7.0*t70+6336.0/7.0*t72-t246+t272-t273-t274+t275;
828  t284 = 2700.0/7.0-16200.0/7.0*eta-70200.0/7.0*t4+85320.0/7.0*t2-28800.0/7.0*xi+124848.0/7.0*t7-185328.0/7.0*t11-t227+t228+t229-t230-t231-t232-t235-t236;
829  t287 = t237-t238+t239-t265+21708.0/7.0*t15-81864.0/7.0*t13+t242+t268-t269+t244+t245-t272+t273+t274-t275;
830  t301 = 1800.0/7.0*t48;
831  t304 = -32140.0/77.0+53840.0/77.0*eta+138000.0/7.0*t4-1701720.0/77.0*t2+487320.0/77.0*xi-1049820.0/77.0*t7+2308200.0/77.0*t11+13800.0/7.0*t30-8200.0/7.0*t33-8200.0/7.0*t37+4100.0/7.0*t40+4100.0/7.0*t43+4100.0/7.0*t45-t301+13800.0/7.0*t50+6900.0/7.0*t53;
832  t311 = 168640.0/77.0*t9;
833  t313 = 900.0/7.0*t67;
834  t317 = 1800.0/7.0*t74;
835  t322 = -13800.0/7.0*t56+6900.0/7.0*t59-8200.0/7.0*t62+13800.0/7.0*t96-271200.0/77.0*t15+166500.0/77.0*t13+t311-9200.0/7.0*t100+t313+13800.0/7.0*t102-13800.0/7.0*t70-27600.0/7.0*t72-t317+13800.0/7.0*t106-9200.0/7.0*t108-27600.0/7.0*t110+4600.0/7.0*t112;
836  t338 = 4160.0/77.0+5440.0/77.0*eta-21720.0/7.0*t4+258480.0/77.0*t2-69280.0/77.0*xi+39180.0/77.0*t7-203760.0/77.0*t11+400.0*t30+400.0*t33+400.0*t37-200.0*t40-200.0*t43-200.0*t45-t301+400.0*t50+200.0*t53;
837  t342 = 2700.0/7.0*t96;
838  t346 = 1800.0/7.0*t100;
839  t347 = 2700.0/7.0*t102;
840  t350 = 2700.0/7.0*t106;
841  t351 = 1800.0/7.0*t108;
842  t352 = 5400.0/7.0*t110;
843  t353 = 900.0/7.0*t112;
844  t354 = -400.0*t56+200.0*t59+400.0*t62-t342+40980.0/77.0*t15+122940.0/77.0*t13-49160.0/77.0*t9+t346+t313-t347-400.0*t70-800.0*t72-t317-t350+t351+t352-t353;
845  t371 = 13840.0/77.0-115560.0/77.0*eta-32280.0/7.0*t4+432720.0/77.0*t2-20960.0/11.0*xi+801480.0/77.0*t7-1104000.0/77.0*t11-8200.0/7.0*t30+13800.0/7.0*t33+13800.0/7.0*t37-6900.0/7.0*t40-6900.0/7.0*t43-6900.0/7.0*t45+9200.0/7.0*t48-8200.0/7.0*t50-4100.0/7.0*t53;
846  t381 = 8200.0/7.0*t56-4100.0/7.0*t59+13800.0/7.0*t62-t342+55500.0/77.0*t15-813600.0/77.0*t13+t311+t346-4600.0/7.0*t67-t347+8200.0/7.0*t70+16400.0/7.0*t72+9200.0/7.0*t74-t350+t351+t352-t353;
847  t395 = 6000.0*t48;
848  t399 = 142848.0/11.0*eta+111600.0*t4-1508400.0/11.0*t2+518400.0/11.0*xi-980064.0/11.0*t7+1695600.0/11.0*t11+5040.0*t30+5040.0*t33+5040.0*t37-2520.0*t40-2520.0*t43-2520.0*t45+t395+5040.0*t50+2520.0*t53-5040.0*t56;
849  t402 = 9000.0*t96;
850  t406 = 6000.0*t100;
851  t407 = 3000.0*t67;
852  t408 = 9000.0*t102;
853  t411 = 6000.0*t74;
854  t412 = 9000.0*t106;
855  t413 = 6000.0*t108;
856  t414 = 18000.0*t110;
857  t415 = 3000.0*t112;
858  t416 = -49800.0/11.0+2520.0*t59+5040.0*t62+t402+26640.0/11.0*t15+79920.0/11.0*t13-71856.0/11.0*t9-t406-t407+t408-5040.0*t70-10080.0*t72+t411+t412-t413-t414+t415;
859  t424 = 12960.0*t30;
860  t425 = 2880.0*t33;
861  t426 = 2880.0*t37;
862  t427 = 1440.0*t40;
863  t428 = 1440.0*t43;
864  t429 = 1440.0*t45;
865  t430 = 1920.0*t48;
866  t431 = 12960.0*t50;
867  t432 = 6480.0*t53;
868  t433 = -10200.0/11.0+30384.0/11.0*eta+68400.0*t4-795600.0/11.0*t2+201600.0/11.0*xi-664056.0/11.0*t7+1505520.0/11.0*t11+t424-t425-t426+t427+t428+t429-t430+t431+t432;
869  t434 = 12960.0*t56;
870  t435 = 6480.0*t59;
871  t436 = 2880.0*t62;
872  t439 = 103968.0/11.0*t9;
873  t440 = 960.0*t67;
874  t441 = 12960.0*t70;
875  t442 = 25920.0*t72;
876  t443 = 1920.0*t74;
877  t444 = -t434+t435-t436+t402-163440.0/11.0*t15+115560.0/11.0*t13+t439-t406+t440+t408-t441-t442-t443+t412-t413-t414+t415;
878  t452 = 2880.0*t30;
879  t453 = 12960.0*t33;
880  t454 = 12960.0*t37;
881  t455 = 6480.0*t40;
882  t456 = 6480.0*t43;
883  t457 = 6480.0*t45;
884  t458 = 2880.0*t50;
885  t459 = 1440.0*t53;
886  t460 = 9600.0/11.0-69408.0/11.0*eta-28800.0*t4+368640.0/11.0*t2-115200.0/11.0*xi+502560.0/11.0*t7-751680.0/11.0*t11-t452+t453+t454-t455-t456-t457+t395-t458-t459;
887  t461 = 2880.0*t56;
888  t462 = 1440.0*t59;
889  t463 = 12960.0*t62;
890  t464 = 2880.0*t96;
891  t467 = 1920.0*t100;
892  t468 = 2880.0*t102;
893  t469 = 2880.0*t70;
894  t470 = 5760.0*t72;
895  t471 = 2880.0*t106;
896  t472 = 1920.0*t108;
897  t473 = 5760.0*t110;
898  t474 = 960.0*t112;
899  t475 = t461-t462+t463-t464+38520.0/11.0*t15-490320.0/11.0*t13+t439+t467-t407-t468+t469+t470+t411-t471+t472+t473-t474;
900  t480 = 412800.0/11.0*xi;
901  t483 = 18720.0*t30;
902  t484 = 28800.0*t33;
903  t485 = 28800.0*t37;
904  t486 = 14400.0*t40;
905  t487 = 14400.0*t43;
906  t488 = 14400.0*t45;
907  t489 = 8640.0*t48;
908  t490 = 18720.0*t50;
909  t491 = 9360.0*t53;
910  t492 = 43200.0/11.0-39360.0/11.0*eta-79200.0*t4+1120320.0/11.0*t2-t480-268128.0/11.0*t7+276480.0/11.0*t11+t483-t484-t485+t486+t487+t488-t489+t490+t491;
911  t493 = 18720.0*t56;
912  t494 = 9360.0*t59;
913  t495 = 28800.0*t62;
914  t496 = 5040.0*t96;
915  t499 = 204912.0/11.0*t9;
916  t500 = 3360.0*t100;
917  t501 = 4320.0*t67;
918  t502 = 5040.0*t102;
919  t503 = 18720.0*t70;
920  t504 = 37440.0*t72;
921  t505 = 8640.0*t74;
922  t506 = 5040.0*t106;
923  t507 = 3360.0*t108;
924  t508 = 10080.0*t110;
925  t509 = 1680.0*t112;
926  t510 = -t493+t494-t495-t496-509760.0/11.0*t15+870480.0/11.0*t13+t499+t500+t501-t502-t503-t504-t505-t506+t507+t508-t509;
927  t518 = 28800.0*t30;
928  t519 = 18720.0*t33;
929  t520 = 18720.0*t37;
930  t521 = 9360.0*t40;
931  t522 = 9360.0*t43;
932  t523 = 9360.0*t45;
933  t524 = 3360.0*t48;
934  t525 = 28800.0*t50;
935  t526 = 14400.0*t53;
936  t527 = 56400.0/11.0-349824.0/11.0*eta-144000.0*t4+1896480.0/11.0*t2-624000.0/11.0*xi+2627424.0/11.0*t7-4000320.0/11.0*t11-t518+t519+t520-t521-t522-t523-t524-t525-t526;
937  t528 = 28800.0*t56;
938  t529 = 14400.0*t59;
939  t530 = 18720.0*t62;
940  t531 = 12960.0*t96;
941  t534 = 8640.0*t100;
942  t535 = 1680.0*t67;
943  t536 = 12960.0*t102;
944  t537 = 28800.0*t70;
945  t538 = 57600.0*t72;
946  t539 = 3360.0*t74;
947  t540 = 12960.0*t106;
948  t541 = 8640.0*t108;
949  t542 = 25920.0*t110;
950  t543 = 4320.0*t112;
951  t544 = t528-t529+t530-t531+290160.0/11.0*t15-1529280.0/11.0*t13+t499+t534+t535-t536+t537+t538-t539-t540+t541+t542-t543;
952  t549 = 9600.0/11.0*xi;
953  t552 = -9600.0/11.0-37248.0/11.0*eta-21600.0*t4+169920.0/11.0*t2+t549+1029168.0/11.0*t7-1909440.0/11.0*t11-t518+t519+t520-t521-t522-t523+t430-t525-t526;
954  t555 = 251280.0/11.0*t9;
955  t556 = t528-t529+t530-t496+519840.0/11.0*t15-673920.0/11.0*t13-t555+t500-t440-t502+t537+t538+t443-t506+t507+t508-t509;
956  t563 = 30000.0/11.0-60480.0/11.0*eta-115200.0*t4+1421280.0/11.0*t2-t480+565056.0/11.0*t7-1434240.0/11.0*t11+t452-t453-t454+t455+t456+t457+t430+t458+t459;
957  t566 = 76608.0/11.0*t9;
958  t567 = -t461+t462-t463-t531-192960.0/11.0*t15+276480.0/11.0*t13+t566+t534-t440-t536-t469-t470+t443-t540+t541+t542-t543;
959  t575 = -22800.0/11.0+228864.0/11.0*eta+43200.0*t4-606240.0/11.0*t2+220800.0/11.0*xi-1567008.0/11.0*t7+2034720.0/11.0*t11+t483-t484-t485+t486+t487+t488-t524+t490+t491;
960  t578 = -t493+t494-t495+t464-224640.0/11.0*t15+1559520.0/11.0*t13-t555-t467+t535+t468-t503-t504-t539+t471-t472-t473+t474;
961  t585 = 3600.0/11.0-93744.0/11.0*eta+14400.0*t4-131040.0/11.0*t2+t549+561888.0/11.0*t7-531360.0/11.0*t11-t424+t425+t426-t427-t428-t429-t489-t431-t432;
962  t588 = t434-t435+t436+t464+92160.0/11.0*t15-578880.0/11.0*t13+t566-t467+t501+t468+t441+t442-t505+t471-t472-t473+t474;
963  t605 = -50400.0/11.0+211008.0/11.0*eta+151200.0*t4-1935360.0/11.0*t2+604800.0/11.0*xi-1640520.0/11.0*t7+2948400.0/11.0*t11+15120.0*t30+15120.0*t33+15120.0*t37-7560.0*t40-7560.0*t43-7560.0*t45+10080.0*t48+15120.0*t50+7560.0*t53;
964  t623 = -15120.0*t56+7560.0*t59+15120.0*t62+15120.0*t96+68040.0/11.0*t15+204120.0/11.0*t13-113400.0/11.0*t9-10080.0*t100-5040.0*t67+15120.0*t102-15120.0*t70-30240.0*t72+10080.0*t74+15120.0*t106-10080.0*t108-30240.0*t110+5040.0*t112;
965 
966  values[0] = 0.0;
967  values[1] = -420.0+t1-t3+t5+3360.0*eta-15120.0*t7-7560.0*t9+15120.0*t11+15120.0*t13+5040.0*t15;
968  values[2] = 420.0-t1+t3-t5;
969  values[3] = t47+t76;
970  values[4] = t92+t114;
971  values[5] = t122+t125;
972  values[6] = t142+t160;
973  values[7] = t176+t192;
974  values[8] = t209+t219;
975  values[9] = t233+t247;
976  values[10] = t261+t276;
977  values[11] = t284+t287;
978  values[12] = t304+t322;
979  values[13] = t338+t354;
980  values[14] = t371+t381;
981  values[15] = t399+t416;
982  values[16] = t433+t444;
983  values[17] = t460+t475;
984  values[18] = t492+t510;
985  values[19] = t527+t544;
986  values[20] = t552+t556;
987  values[21] = t563+t567;
988  values[22] = t575+t578;
989  values[23] = t585+t588;
990  values[24] = t605+t623;
991 }
992 
993 // values of derivatives in eta-eta direction
994 static void N_T_P5_2D_DeriveXiEta(double xi, double eta, double *values)
995 {
996  double t2, t4, t7, t9, t11, t13, t15, t18, t22, t23, t24, t25, t26, t27;
997  double t29, t30, t31, t32, t33, t34, t36, t37, t38, t39, t40, t41, t42;
998  double t43, t47, t52, t53, t55, t56, t57, t58, t59, t60, t62, t63, t65;
999  double t66, t67, t69, t70, t72, t74, t76, t77, t79, t80, t81, t85, t86;
1000  double t87, t88, t89, t90, t98, t109, t113, t120, t127, t129, t131, t134;
1001  double t135, t136, t144, t145, t147, t165, t170, t171, t172, t178, t196;
1002  double t208, t226, t228, t231, t232, t233, t234, t235, t236, t237, t241;
1003  double t245, t246, t247, t248, t249, t254, t255, t256, t260, t261, t262;
1004  double t269, t279, t283, t290, t298, t300, t303, t304, t305, t313, t315;
1005  double t333, t338, t339, t340, t346, t364, t376, t394, t399, t400, t401;
1006  double t402, t403, t404, t409, t411, t424, t427, t429, t431, t434, t435;
1007  double t436, t437, t438, t439, t440, t442, t448, t449, t451, t452, t453;
1008  double t455, t457, t459, t460, t464, t465, t466, t467, t468, t469, t470;
1009  double t472, t478, t479, t480, t481, t482, t483, t488, t489, t490, t492;
1010  double t495, t496, t497, t498, t499, t500, t501, t502, t503, t504, t506;
1011  double t512, t513, t515, t516, t517, t520, t521, t522, t523, t527, t528;
1012  double t529, t530, t531, t532, t533, t534, t535, t536, t538, t544, t545;
1013  double t546, t547, t548, t549, t551, t554, t555, t557, t561, t570, t571;
1014  double t573, t577, t583, t584, t585, t586, t591, t600, t605, t611, t628;
1015  double t646;
1016 
1017  t2 = xi*xi;
1018  t4 = t2*xi;
1019  t7 = xi*eta;
1020  t9 = eta*eta;
1021  t11 = t2*eta;
1022  t13 = xi*t9;
1023  t15 = t9*eta;
1024  t18 = 2160.0*t7;
1025  t22 = 2.0*eta-1.0+xi;
1026  t23 = t22*t15;
1027  t24 = 220.0*t23;
1028  t25 = xi-eta;
1029  t26 = t25*t15;
1030  t27 = 1100.0*t26;
1031  t29 = 1.0-2.0*xi-eta;
1032  t30 = t29*t15;
1033  t31 = 220.0*t30;
1034  t32 = t29*t25;
1035  t33 = t32*t7;
1036  t34 = 396.0*t33;
1037  t36 = t32*t22*eta;
1038  t37 = 132.0*t36;
1039  t38 = t25*t22;
1040  t39 = t38*t7;
1041  t40 = 396.0*t39;
1042  t41 = t22*xi;
1043  t42 = t32*t41;
1044  t43 = 132.0*t42;
1045  t47 = -40.0+t18-3132.0*t11-1548.0*t13-t24+t27-t31-t34-t37+t40-t43+882.0*t2-320.0*xi+800.0*eta;
1046  t52 = t22*t2*eta;
1047  t53 = 66.0*t52;
1048  t55 = t29*t2*eta;
1049  t56 = 66.0*t55;
1050  t57 = t38*t9;
1051  t58 = 1386.0*t57;
1052  t59 = t41*t9;
1053  t60 = 66.0*t59;
1054  t62 = t25*xi*t9;
1055  t63 = 330.0*t62;
1056  t65 = t29*xi*t9;
1057  t66 = 66.0*t65;
1058  t67 = t32*t9;
1059  t69 = t29*t22;
1060  t70 = t69*t9;
1061  t72 = t38*t2;
1062  t74 = t32*t2;
1063  t76 = t69*t2;
1064  t77 = 66.0*t76;
1065  t79 = t25*t2*eta;
1066  t80 = 330.0*t79;
1067  t81 = -648.0*t4-4422.0*t9+6168.0*t15-t53-t56+t58-t60+t63-t66-792.0*t67-594.0*t70+132.0*t72-66.0*t74-t77+t80;
1068  t85 = t25*t4;
1069  t86 = 1100.0*t85;
1070  t87 = t22*t4;
1071  t88 = 220.0*t87;
1072  t89 = t29*t4;
1073  t90 = 220.0*t89;
1074  t98 = 8388.0*t11-8388.0*t13+t24-t27+t31-t86+t88+t90+3564.0*t33+1188.0*t36-3564.0*t39+1188.0*t42-6438.0*t2+1120.0*xi-1120.0*eta;
1075  t109 = 66.0*t70;
1076  t113 = 7992.0*t4+6438.0*t9-7992.0*t15+594.0*t52+594.0*t55-1914.0*t57+594.0*t59-2970.0*t62+594.0*t65+1848.0*t67+t109-1848.0*t72+1914.0*t74-t77-2970.0*t79;
1077  t120 = 40.0-t18+1548.0*t11+3132.0*t13+t86-t88-t90-t34-t37+t40-t43+4422.0*t2-800.0*xi+320.0*eta;
1078  t127 = 1386.0*t74;
1079  t129 = -6168.0*t4-882.0*t9+648.0*t15-t53-t56+66.0*t57-t60+t63-t66-132.0*t67+t109+792.0*t72-t127+594.0*t76+t80;
1080  t131 = 8352.0/77.0*t7;
1081  t134 = 1500.0/7.0*t23;
1082  t135 = 7500.0/7.0*t26;
1083  t136 = 1500.0/7.0*t30;
1084  t144 = t69*t7;
1085  t145 = 720.0*t144;
1086  t147 = 8046.0/77.0-t131+490320.0/77.0*t11-448200.0/77.0*t13+t134-t135+t136-2400.0/7.0*t85+480.0/7.0*t87+480.0/7.0*t89+9360.0/7.0*t33+6480.0/7.0*t36-14400.0/7.0*t39+1440.0/7.0*t42+t145-251280.0/77.0*t2;
1087  t165 = 47232.0/77.0*xi-152352.0/77.0*eta+250560.0/77.0*t4+648036.0/77.0*t9-723600.0/77.0*t15+720.0/7.0*t52+720.0/7.0*t55-12240.0/7.0*t57+3240.0/7.0*t59-16200.0/7.0*t62+3240.0/7.0*t65+10980.0/7.0*t67+180.0*t70-2880.0/7.0*t72+3600.0/7.0*t74-720.0/7.0*t76-3600.0/7.0*t79;
1088  t170 = 7500.0/7.0*t85;
1089  t171 = 1500.0/7.0*t87;
1090  t172 = 1500.0/7.0*t89;
1091  t178 = 3690.0/77.0-89136.0/77.0*t7+86400.0/77.0*t11+86400.0/77.0*t13+t134-t135+t136+t170-t171-t172-360.0*t33+360.0*t36-360.0*t39-360.0*t42+t145+257184.0/77.0*t2;
1092  t196 = -39360.0/77.0*xi-39360.0/77.0*eta-398880.0/77.0*t4+257184.0/77.0*t9-398880.0/77.0*t15-180.0*t52-180.0*t55-10260.0/7.0*t57+180.0*t59-900.0*t62+180.0*t65+7020.0/7.0*t67+3240.0/7.0*t70+7020.0/7.0*t72-10260.0/7.0*t74+3240.0/7.0*t76+900.0*t79;
1093  t208 = 8046.0/77.0-t131-448200.0/77.0*t11+490320.0/77.0*t13-480.0/7.0*t23+2400.0/7.0*t26-480.0/7.0*t30+t170-t171-t172-14400.0/7.0*t33-1440.0/7.0*t36+9360.0/7.0*t39-6480.0/7.0*t42+t145+648036.0/77.0*t2;
1094  t226 = -152352.0/77.0*xi+47232.0/77.0*eta-723600.0/77.0*t4-251280.0/77.0*t9+250560.0/77.0*t15-3240.0/7.0*t52-3240.0/7.0*t55+3600.0/7.0*t57-720.0/7.0*t59+3600.0/7.0*t62-720.0/7.0*t65-2880.0/7.0*t67-720.0/7.0*t70+10980.0/7.0*t72-12240.0/7.0*t74+180.0*t76+16200.0/7.0*t79;
1095  t228 = 12960.0/7.0*t7;
1096  t231 = 990.0/7.0*t23;
1097  t232 = 4950.0/7.0*t26;
1098  t233 = 990.0/7.0*t30;
1099  t234 = 22572.0/7.0*t33;
1100  t235 = 7524.0/7.0*t36;
1101  t236 = 22572.0/7.0*t39;
1102  t237 = 7524.0/7.0*t42;
1103  t241 = -810.0/7.0-t228-65124.0/7.0*t11+81864.0/7.0*t13-t231+t232-t233-t234-t235+t236-t237+5562.0*t2-6480.0/7.0*xi+16200.0/7.0*eta;
1104  t245 = 3762.0/7.0*t52;
1105  t246 = 3762.0/7.0*t55;
1106  t247 = 3762.0/7.0*t59;
1107  t248 = 18810.0/7.0*t62;
1108  t249 = 3762.0/7.0*t65;
1109  t254 = 3762.0/7.0*t76;
1110  t255 = 18810.0/7.0*t79;
1111  t256 = -36936.0/7.0*t4-62424.0/7.0*t9+61776.0/7.0*t15-t245-t246+t58-t247+t248-t249-10494.0/7.0*t67+792.0/7.0*t70+7524.0/7.0*t72-3762.0/7.0*t74-t254+t255;
1112  t260 = 4950.0/7.0*t85;
1113  t261 = 990.0/7.0*t87;
1114  t262 = 990.0/7.0*t89;
1115  t269 = -13284.0/7.0*t11+13284.0/7.0*t13+t231-t232+t233-t260+t261+t262-4752.0/7.0*t33-1584.0/7.0*t36+4752.0/7.0*t39-1584.0/7.0*t42-5346.0/7.0*t2+14544.0/7.0*t4;
1116  t279 = 3762.0/7.0*t70;
1117  t283 = 5346.0/7.0*t9-14544.0/7.0*t15-792.0/7.0*t52-792.0/7.0*t55-5148.0/7.0*t57-792.0/7.0*t59+3960.0/7.0*t62-792.0/7.0*t65+198.0*t67+t279-198.0*t72+5148.0/7.0*t74-t254+3960.0/7.0*t79;
1118  t290 = 810.0/7.0+t228-81864.0/7.0*t11+65124.0/7.0*t13+t260-t261-t262-t234-t235+t236-t237+62424.0/7.0*t2-16200.0/7.0*xi+6480.0/7.0*eta;
1119  t298 = -61776.0/7.0*t4-5562.0*t9+36936.0/7.0*t15-t245-t246+3762.0/7.0*t57-t247+t248-t249-7524.0/7.0*t67+t279+10494.0/7.0*t72-t127-792.0/7.0*t76+t255;
1120  t300 = 337280.0/77.0*t7;
1121  t303 = 450.0/7.0*t23;
1122  t304 = 2250.0/7.0*t26;
1123  t305 = 450.0/7.0*t30;
1124  t313 = 400.0*t144;
1125  t315 = 5540.0/77.0+t300+166500.0/77.0*t11-813600.0/77.0*t13+t303-t304+t305-11500.0/7.0*t85+2300.0/7.0*t87+2300.0/7.0*t89+17900.0/7.0*t33+4100.0/7.0*t36-15100.0/7.0*t39+6900.0/7.0*t42-t313-524910.0/77.0*t2;
1126  t333 = 53840.0/77.0*xi-115560.0/77.0*eta+769400.0/77.0*t4+400740.0/77.0*t9-368000.0/77.0*t15+3450.0/7.0*t52+3450.0/7.0*t55-4750.0/7.0*t57+2050.0/7.0*t59-10250.0/7.0*t62+2050.0/7.0*t65+5450.0/7.0*t67-100.0*t70-13800.0/7.0*t72+17250.0/7.0*t74-3450.0/7.0*t76-17250.0/7.0*t79;
1127  t338 = 2250.0/7.0*t85;
1128  t339 = 450.0/7.0*t87;
1129  t340 = 450.0/7.0*t89;
1130  t346 = -510.0/77.0-98320.0/77.0*t7+122940.0/77.0*t11+122940.0/77.0*t13+t303-t304+t305+t338-t339-t340+200.0*t33-200.0*t36+200.0*t39+200.0*t42-t313+19590.0/77.0*t2;
1131  t364 = 5440.0/77.0*xi+5440.0/77.0*eta-67920.0/77.0*t4+19590.0/77.0*t9-67920.0/77.0*t15+100.0*t52+100.0*t55-2000.0/7.0*t57-100.0*t59+500.0*t62-100.0*t65-50.0/7.0*t67+2050.0/7.0*t70-50.0/7.0*t72-2000.0/7.0*t74+2050.0/7.0*t76-500.0*t79;
1132  t376 = 5540.0/77.0+t300-813600.0/77.0*t11+166500.0/77.0*t13-2300.0/7.0*t23+11500.0/7.0*t26-2300.0/7.0*t30+t338-t339-t340-15100.0/7.0*t33-6900.0/7.0*t36+17900.0/7.0*t39-4100.0/7.0*t42-t313+400740.0/77.0*t2;
1133  t394 = -115560.0/77.0*xi+53840.0/77.0*eta-368000.0/77.0*t4-524910.0/77.0*t9+769400.0/77.0*t15-2050.0/7.0*t52-2050.0/7.0*t55+17250.0/7.0*t57-3450.0/7.0*t59+17250.0/7.0*t62-3450.0/7.0*t65-13800.0/7.0*t67-3450.0/7.0*t70+5450.0/7.0*t72-4750.0/7.0*t74-100.0*t76+10250.0/7.0*t79;
1134  t399 = 1500.0*t23;
1135  t400 = 7500.0*t26;
1136  t401 = 1500.0*t30;
1137  t402 = 7500.0*t85;
1138  t403 = 1500.0*t87;
1139  t404 = 1500.0*t89;
1140  t409 = 5040.0*t144;
1141  t411 = -19860.0/11.0-143712.0/11.0*t7+79920.0/11.0*t11+79920.0/11.0*t13-t399+t400-t401-t402+t403+t404+2520.0*t33-2520.0*t36+2520.0*t39+2520.0*t42-t409-490032.0/11.0*t2;
1142  t424 = 3240.0*t70;
1143  t427 = 3240.0*t76;
1144  t429 = 142848.0/11.0*xi+142848.0/11.0*eta+565200.0/11.0*t4-490032.0/11.0*t9+565200.0/11.0*t15+1260.0*t52+1260.0*t55+10260.0*t57-1260.0*t59+6300.0*t62-1260.0*t65-7020.0*t67-t424-7020.0*t72+10260.0*t74-t427-6300.0*t79;
1145  t431 = 207936.0/11.0*t7;
1146  t434 = 480.0*t23;
1147  t435 = 2400.0*t26;
1148  t436 = 480.0*t30;
1149  t437 = 14400.0*t33;
1150  t438 = 1440.0*t36;
1151  t439 = 9360.0*t39;
1152  t440 = 6480.0*t42;
1153  t442 = 3504.0/11.0+t431+115560.0/11.0*t11-490320.0/11.0*t13+t434-t435+t436-t402+t403+t404+t437+t438-t439+t440-t409-332028.0/11.0*t2;
1154  t448 = 3240.0*t52;
1155  t449 = 3240.0*t55;
1156  t451 = 720.0*t59;
1157  t452 = 3600.0*t62;
1158  t453 = 720.0*t65;
1159  t455 = 720.0*t70;
1160  t457 = 12240.0*t74;
1161  t459 = 16200.0*t79;
1162  t460 = 30384.0/11.0*xi-69408.0/11.0*eta+501840.0/11.0*t4+251280.0/11.0*t9-250560.0/11.0*t15+t448+t449-3600.0*t57+t451-t452+t453+2880.0*t67+t455-10980.0*t72+t457-1260.0*t76-t459;
1163  t464 = 2400.0*t85;
1164  t465 = 480.0*t87;
1165  t466 = 480.0*t89;
1166  t467 = 9360.0*t33;
1167  t468 = 6480.0*t36;
1168  t469 = 14400.0*t39;
1169  t470 = 1440.0*t42;
1170  t472 = 3504.0/11.0+t431-490320.0/11.0*t11+115560.0/11.0*t13-t399+t400-t401+t464-t465-t466-t467-t468+t469-t470-t409+251280.0/11.0*t2;
1171  t478 = 720.0*t52;
1172  t479 = 720.0*t55;
1173  t480 = 12240.0*t57;
1174  t481 = 3240.0*t59;
1175  t482 = 16200.0*t62;
1176  t483 = 3240.0*t65;
1177  t488 = 720.0*t76;
1178  t489 = 3600.0*t79;
1179  t490 = -69408.0/11.0*xi+30384.0/11.0*eta-250560.0/11.0*t4-332028.0/11.0*t9+501840.0/11.0*t15-t478-t479+t480-t481+t482-t483-10980.0*t67-1260.0*t70+2880.0*t72-3600.0*t74+t488+t489;
1180  t492 = 409824.0/11.0*t7;
1181  t495 = 2160.0*t23;
1182  t496 = 10800.0*t26;
1183  t497 = 2160.0*t30;
1184  t498 = 4200.0*t85;
1185  t499 = 840.0*t87;
1186  t500 = 840.0*t89;
1187  t501 = 33120.0*t33;
1188  t502 = 14400.0*t36;
1189  t503 = 38160.0*t39;
1190  t504 = 9360.0*t42;
1191  t506 = 26328.0/11.0+t492+870480.0/11.0*t11-1529280.0/11.0*t13+t495-t496+t497+t498-t499-t500+t501+t502-t503+t504+t409-134064.0/11.0*t2;
1192  t512 = 4680.0*t52;
1193  t513 = 4680.0*t55;
1194  t515 = 7200.0*t59;
1195  t516 = 36000.0*t62;
1196  t517 = 7200.0*t65;
1197  t520 = 360.0*t74;
1198  t521 = 7200.0*t76;
1199  t522 = 23400.0*t79;
1200  t523 = -39360.0/11.0*xi-349824.0/11.0*eta+92160.0/11.0*t4+1313712.0/11.0*t9-1333440.0/11.0*t15+t512+t513-20160.0*t57+t515-t516+t517+20880.0*t67-t455-6840.0*t72-t520+t521-t522;
1201  t527 = 840.0*t23;
1202  t528 = 4200.0*t26;
1203  t529 = 840.0*t30;
1204  t530 = 10800.0*t85;
1205  t531 = 2160.0*t87;
1206  t532 = 2160.0*t89;
1207  t533 = 38160.0*t33;
1208  t534 = 9360.0*t36;
1209  t535 = 33120.0*t39;
1210  t536 = 14400.0*t42;
1211  t538 = 26328.0/11.0+t492-1529280.0/11.0*t11+870480.0/11.0*t13+t527-t528+t529+t530-t531-t532-t533-t534+t535-t536+t409+1313712.0/11.0*t2;
1212  t544 = 7200.0*t52;
1213  t545 = 7200.0*t55;
1214  t546 = 360.0*t57;
1215  t547 = 4680.0*t59;
1216  t548 = 23400.0*t62;
1217  t549 = 4680.0*t65;
1218  t551 = 7200.0*t70;
1219  t554 = 36000.0*t79;
1220  t555 = -349824.0/11.0*xi-39360.0/11.0*eta-1333440.0/11.0*t4-134064.0/11.0*t9+92160.0/11.0*t15-t544-t545-t546-t547+t548-t549-6840.0*t67+t551+20880.0*t72-20160.0*t74-t488+t554;
1221  t557 = 502560.0/11.0*t7;
1222  t561 = -14064.0/11.0-t557-673920.0/11.0*t11+1559520.0/11.0*t13-t434+t435-t436+t498-t499-t500-t533-t534+t535-t536+t409+514584.0/11.0*t2;
1223  t570 = 4680.0*t76;
1224  t571 = -37248.0/11.0*xi+228864.0/11.0*eta-636480.0/11.0*t4-783504.0/11.0*t9+678240.0/11.0*t15-t544-t545+7560.0*t57-t547+t548-t549-10800.0*t67+t424+16920.0*t72-t457-t570+t554;
1225  t573 = 153216.0/11.0*t7;
1226  t577 = 8904.0/11.0+t573+276480.0/11.0*t11-578880.0/11.0*t13-t434+t435-t436+t530-t531-t532+t467+t468-t469+t470+t409+282528.0/11.0*t2;
1227  t583 = 5040.0*t67;
1228  t584 = 4680.0*t70;
1229  t585 = 5040.0*t72;
1230  t586 = -60480.0/11.0*xi-93744.0/11.0*eta-478080.0/11.0*t4+280944.0/11.0*t9-177120.0/11.0*t15+t478+t479-t546+t481-t482+t483+t583-t584+t585-t457+t521-t489;
1231  t591 = -14064.0/11.0-t557+1559520.0/11.0*t11-673920.0/11.0*t13+t527-t528+t529-t464+t465+t466+t501+t502-t503+t504+t409-783504.0/11.0*t2;
1232  t600 = 228864.0/11.0*xi-37248.0/11.0*eta+678240.0/11.0*t4+514584.0/11.0*t9-636480.0/11.0*t15+t512+t513-t480+t515-t516+t517+16920.0*t67-t584-10800.0*t72+7560.0*t74+t427-t522;
1233  t605 = 8904.0/11.0+t573-578880.0/11.0*t11+276480.0/11.0*t13+t495-t496+t497-t464+t465+t466-t437-t438+t439-t440+t409+280944.0/11.0*t2;
1234  t611 = -93744.0/11.0*xi-60480.0/11.0*eta-177120.0/11.0*t4+282528.0/11.0*t9-478080.0/11.0*t15-t448-t449-t480-t451+t452-t453+t583+t551+t585-t520-t570+t459;
1235  t628 = -24864.0/11.0-226800.0/11.0*t7+204120.0/11.0*t11+204120.0/11.0*t13-2520.0*t23+12600.0*t26-2520.0*t30-12600.0*t85+2520.0*t87+2520.0*t89+7560.0*t33-7560.0*t36+7560.0*t39+7560.0*t42-15120.0*t144-820260.0/11.0*t2;
1236  t646 = 211008.0/11.0*xi+211008.0/11.0*eta+982800.0/11.0*t4-820260.0/11.0*t9+982800.0/11.0*t15+3780.0*t52+3780.0*t55+18900.0*t57-3780.0*t59+18900.0*t62-3780.0*t65-15120.0*t67-3780.0*t70-15120.0*t72+18900.0*t74-3780.0*t76-18900.0*t79;
1237 
1238  values[0] = 0.0;
1239  values[1] = -420.0+3360.0*xi-7560.0*t2+5040.0*t4+3360.0*eta-15120.0*t7-7560.0*t9+15120.0*t11+15120.0*t13+5040.0*t15;
1240  values[2] = 0.0;
1241  values[3] = t47+t81;
1242  values[4] = t98+t113;
1243  values[5] = t120+t129;
1244  values[6] = t147+t165;
1245  values[7] = t178+t196;
1246  values[8] = t208+t226;
1247  values[9] = t241+t256;
1248  values[10] = t269+t283;
1249  values[11] = t290+t298;
1250  values[12] = t315+t333;
1251  values[13] = t346+t364;
1252  values[14] = t376+t394;
1253  values[15] = t411+t429;
1254  values[16] = t442+t460;
1255  values[17] = t472+t490;
1256  values[18] = t506+t523;
1257  values[19] = t538+t555;
1258  values[20] = t561+t571;
1259  values[21] = t577+t586;
1260  values[22] = t591+t600;
1261  values[23] = t605+t611;
1262  values[24] = t628+t646;
1263 }
1264 
1265 // values of derivatives in xi-eta direction
1266 static void N_T_P5_2D_DeriveEtaEta(double xi, double eta, double *values)
1267 {
1268  double t1, t2, t3, t4, t5, t8, t10, t12, t14, t16, t20, t21, t22, t24;
1269  double t26, t27, t28, t29, t30, t31, t32, t33, t35, t36, t37, t38, t39;
1270  double t40, t41, t42, t44, t45, t46, t47, t48, t49, t50, t51, t52, t53;
1271  double t54, t56, t57, t59, t60, t65, t66, t68, t69, t71, t72, t73, t74;
1272  double t75, t76, t80, t81, t82, t83, t84, t96, t108, t109, t110, t111;
1273  double t112, t113, t114, t119, t125, t127, t131, t132, t135, t136, t141;
1274  double t142, t155, t156, t160, t172, t185, t186, t187, t188, t204, t219;
1275  double t221, t222, t223, t224, t225, t226, t227, t228, t229, t230, t231;
1276  double t232, t233, t235, t236, t241, t242, t243, t244, t245, t249, t250;
1277  double t251, t262, t273, t274, t275, t276, t281, t287, t289, t293, t294;
1278  double t297, t298, t303, t304, t317, t318, t322, t334, t347, t348, t349;
1279  double t350, t366, t381, t383, t387, t388, t391, t392, t399, t411, t412;
1280  double t413, t414, t415, t416, t418, t419, t420, t421, t422, t423, t424;
1281  double t425, t426, t427, t428, t429, t430, t432, t433, t438, t439, t440;
1282  double t441, t442, t446, t447, t448, t450, t451, t452, t453, t454, t455;
1283  double t456, t457, t459, t464, t465, t466, t467, t468, t472, t473, t474;
1284  double t475, t477, t478, t479, t480, t481, t482, t483, t484, t485, t486;
1285  double t487, t488, t489, t491, t492, t497, t498, t499, t500, t501, t505;
1286  double t506, t507, t508, t509, t510, t512, t513, t514, t515, t516, t517;
1287  double t518, t519, t520, t521, t522, t523, t524, t526, t529, t531, t532;
1288  double t533, t534, t535, t539, t540, t541, t542, t543, t544, t547, t548;
1289  double t556, t559, t560, t563, t568, t571, t578, t581, t588, t605, t623;
1290 
1291  t1 = 3360.0*eta;
1292  t2 = eta*eta;
1293  t3 = 7560.0*t2;
1294  t4 = t2*eta;
1295  t5 = 5040.0*t4;
1296  t8 = xi*xi;
1297  t10 = t8*xi;
1298  t12 = xi*eta;
1299  t14 = t8*eta;
1300  t16 = xi*t2;
1301  t20 = 1.0-2.0*xi-eta;
1302  t21 = xi-eta;
1303  t22 = t20*t21;
1304  t24 = 2.0*eta-1.0+xi;
1305  t26 = t24*eta*t22;
1306  t27 = 1320.0*t26;
1307  t28 = t21*t24;
1308  t29 = t28*t12;
1309  t30 = 264.0*t29;
1310  t31 = t24*xi;
1311  t32 = t22*t31;
1312  t33 = 132.0*t32;
1313  t35 = t21*xi*t2;
1314  t36 = 264.0*t35;
1315  t37 = t24*t4;
1316  t38 = 440.0*t37;
1317  t39 = t28*t2;
1318  t40 = 1320.0*t39;
1319  t41 = t31*t2;
1320  t42 = 132.0*t41;
1321  t44 = t20*xi*t2;
1322  t45 = 264.0*t44;
1323  t46 = t22*t2;
1324  t47 = 2640.0*t46;
1325  t48 = t20*t24;
1326  t49 = t48*t2;
1327  t50 = 1320.0*t49;
1328  t51 = t28*t8;
1329  t52 = 132.0*t51;
1330  t53 = t22*t8;
1331  t54 = 264.0*t53;
1332  t56 = t21*t8*eta;
1333  t57 = 264.0*t56;
1334  t59 = -600.0-t27+t30-t33+t36-t38+t40-t42+t45-t47+t50+t52-t54+t57-1044.0*t10;
1335  t60 = 1080.0*t8;
1336  t65 = t48*t8;
1337  t66 = 132.0*t65;
1338  t68 = t24*t8*eta;
1339  t69 = 132.0*t68;
1340  t71 = t20*t8*eta;
1341  t72 = 264.0*t71;
1342  t73 = t48*t12;
1343  t74 = 264.0*t73;
1344  t75 = t22*t12;
1345  t76 = 528.0*t75;
1346  t80 = t21*t4;
1347  t81 = 880.0*t80;
1348  t82 = t20*t4;
1349  t83 = 880.0*t82;
1350  t84 = t60-18960.0*t2+15600.0*t4+6400.0*eta+800.0*xi+t66-t69+t72+t74-t76-8844.0*t12-1548.0*t14+18504.0*t16+t81+t83;
1351  t96 = 200.0+t27-2376.0*t29+1188.0*t32-2376.0*t35+t38-t40+1188.0*t41-2376.0*t44+t47-t50-1188.0*t51+2376.0*t53-2376.0*t56+2796.0*t10+11760.0*t2;
1352  t108 = t21*t10;
1353  t109 = 880.0*t108;
1354  t110 = t24*t10;
1355  t111 = 440.0*t110;
1356  t112 = t20*t10;
1357  t113 = 880.0*t112;
1358  t114 = -10800.0*t4-3200.0*eta-1120.0*xi-1188.0*t65+1188.0*t68-2376.0*t71-2376.0*t73+4752.0*t75+12876.0*t12-8388.0*t14-23976.0*t16-t81-t83-t109+t111-t113;
1359  t119 = -20.0+t30-t33+t36-t42+t45+t52-t54+t57+516.0*t10-t60-360.0*t2+240.0*t4;
1360  t125 = 160.0*eta+320.0*xi+t66-t69+t72+t74-t76-1764.0*t12+3132.0*t14+1944.0*t16+t109-t111+t113;
1361  t127 = 9000.0/7.0*t26;
1362  t131 = 3000.0/7.0*t37;
1363  t132 = 9000.0/7.0*t39;
1364  t135 = 18000.0/7.0*t46;
1365  t136 = 9000.0/7.0*t49;
1366  t141 = 4176.0/77.0*t8;
1367  t142 = 42540.0/77.0+t127-12960.0/7.0*t29+6480.0/7.0*t32-12960.0/7.0*t35+t131-t132+6480.0/7.0*t41-12960.0/7.0*t44+t135-t136-1440.0/7.0*t51+2880.0/7.0*t53-2880.0/7.0*t56+163440.0/77.0*t10-t141;
1368  t155 = 6000.0/7.0*t80;
1369  t156 = 6000.0/7.0*t82;
1370  t160 = 1377720.0/77.0*t2-103680.0/7.0*t4-65760.0/11.0*eta-152352.0/77.0*xi-1440.0/7.0*t65+1440.0/7.0*t68-2880.0/7.0*t71-12960.0/7.0*t73+25920.0/7.0*t75+1296072.0/77.0*t12-448200.0/77.0*t14-2170800.0/77.0*t16-t155-t156-1920.0/7.0*t108+960.0/7.0*t110-1920.0/7.0*t112;
1371  t172 = 17460.0/77.0+t127-720.0*t29+360.0*t32-720.0*t35+t131-t132+360.0*t41-720.0*t44+t135-t136+360.0*t51-720.0*t53+720.0*t56+28800.0/77.0*t10-44568.0/77.0*t8;
1372  t185 = 6000.0/7.0*t108;
1373  t186 = 3000.0/7.0*t110;
1374  t187 = 6000.0/7.0*t112;
1375  t188 = 926280.0/77.0*t2-76320.0/7.0*t4-259680.0/77.0*eta-39360.0/77.0*xi+360.0*t65-360.0*t68+720.0*t71-720.0*t73+1440.0*t75+514368.0/77.0*t12+86400.0/77.0*t14-1196640.0/77.0*t16-t155-t156+t185-t186+t187;
1376  t204 = -8676.0/77.0-2880.0/7.0*t26+2880.0/7.0*t29-1440.0/7.0*t32+2880.0/7.0*t35-960.0/7.0*t37+2880.0/7.0*t39-1440.0/7.0*t41+2880.0/7.0*t44-5760.0/7.0*t46+2880.0/7.0*t49+6480.0/7.0*t51-12960.0/7.0*t53+12960.0/7.0*t56-149400.0/77.0*t10-t141;
1377  t219 = -368640.0/77.0*t2+28800.0/7.0*t4+115200.0/77.0*eta+47232.0/77.0*xi+6480.0/7.0*t65-6480.0/7.0*t68+12960.0/7.0*t71+2880.0/7.0*t73-5760.0/7.0*t75-502560.0/77.0*t12+490320.0/77.0*t14+751680.0/77.0*t16+1920.0/7.0*t80+1920.0/7.0*t82+t185-t186+t187;
1378  t221 = 5940.0/7.0*t26;
1379  t222 = 15048.0/7.0*t29;
1380  t223 = 7524.0/7.0*t32;
1381  t224 = 15048.0/7.0*t35;
1382  t225 = 1980.0/7.0*t37;
1383  t226 = 5940.0/7.0*t39;
1384  t227 = 7524.0/7.0*t41;
1385  t228 = 15048.0/7.0*t44;
1386  t229 = 11880.0/7.0*t46;
1387  t230 = 5940.0/7.0*t49;
1388  t231 = 7524.0/7.0*t51;
1389  t232 = 15048.0/7.0*t53;
1390  t233 = 15048.0/7.0*t56;
1391  t235 = -2700.0/7.0-t221+t222-t223+t224-t225+t226-t227+t228-t229+t230+t231-t232+t233-21708.0/7.0*t10;
1392  t236 = 6480.0/7.0*t8;
1393  t241 = 7524.0/7.0*t65;
1394  t242 = 7524.0/7.0*t68;
1395  t243 = 15048.0/7.0*t71;
1396  t244 = 15048.0/7.0*t73;
1397  t245 = 30096.0/7.0*t75;
1398  t249 = 3960.0/7.0*t80;
1399  t250 = 3960.0/7.0*t82;
1400  t251 = -t236-85320.0/7.0*t2+70200.0/7.0*t4+28800.0/7.0*eta+16200.0/7.0*xi+t241-t242+t243+t244-t245-124848.0/7.0*t12+81864.0/7.0*t14+185328.0/7.0*t16+t249+t250;
1401  t262 = 900.0/7.0+t221+3168.0/7.0*t29-1584.0/7.0*t32+3168.0/7.0*t35+t225-t226-1584.0/7.0*t41+3168.0/7.0*t44+t229-t230+1584.0/7.0*t51-3168.0/7.0*t53+3168.0/7.0*t56-4428.0/7.0*t10;
1402  t273 = 3960.0/7.0*t108;
1403  t274 = 1980.0/7.0*t110;
1404  t275 = 3960.0/7.0*t112;
1405  t276 = t3-48600.0/7.0*t4-14400.0/7.0*eta+1584.0/7.0*t65-1584.0/7.0*t68+3168.0/7.0*t71+3168.0/7.0*t73-6336.0/7.0*t75+10692.0/7.0*t12+13284.0/7.0*t14-43632.0/7.0*t16-t249-t250-t273+t274-t275;
1406  t281 = -720.0/7.0+t222-t223+t224-t227+t228+t231-t232+t233-27288.0/7.0*t10+t236-20520.0/7.0*t2+13680.0/7.0*t4;
1407  t287 = 8280.0/7.0*eta+6480.0/7.0*xi+t241-t242+t243+t244-t245-11124.0*t12+65124.0/7.0*t14+110808.0/7.0*t16+t273-t274+t275;
1408  t289 = 2700.0/7.0*t26;
1409  t293 = 900.0/7.0*t37;
1410  t294 = 2700.0/7.0*t39;
1411  t297 = 5400.0/7.0*t46;
1412  t298 = 2700.0/7.0*t49;
1413  t303 = 168640.0/77.0*t8;
1414  t304 = 13840.0/77.0+t289-8200.0/7.0*t29+4100.0/7.0*t32-8200.0/7.0*t35+t293-t294+4100.0/7.0*t41-8200.0/7.0*t44+t297-t298-6900.0/7.0*t51+13800.0/7.0*t53-13800.0/7.0*t56+55500.0/77.0*t10+t303;
1415  t317 = 1800.0/7.0*t80;
1416  t318 = 1800.0/7.0*t82;
1417  t322 = 432720.0/77.0*t2-32280.0/7.0*t4-20960.0/11.0*eta-115560.0/77.0*xi-6900.0/7.0*t65+6900.0/7.0*t68-13800.0/7.0*t71-8200.0/7.0*t73+16400.0/7.0*t75+801480.0/77.0*t12-813600.0/77.0*t14-1104000.0/77.0*t16-t317-t318-9200.0/7.0*t108+4600.0/7.0*t110-9200.0/7.0*t112;
1418  t334 = 4160.0/77.0+t289+400.0*t29-200.0*t32+400.0*t35+t293-t294-200.0*t41+400.0*t44+t297-t298-200.0*t51+400.0*t53-400.0*t56+40980.0/77.0*t10-49160.0/77.0*t8;
1419  t347 = 1800.0/7.0*t108;
1420  t348 = 900.0/7.0*t110;
1421  t349 = 1800.0/7.0*t112;
1422  t350 = 258480.0/77.0*t2-21720.0/7.0*t4-69280.0/77.0*eta+5440.0/77.0*xi-200.0*t65+200.0*t68-400.0*t71+400.0*t73-800.0*t75+39180.0/77.0*t12+122940.0/77.0*t14-203760.0/77.0*t16-t317-t318+t347-t348+t349;
1423  t366 = -32140.0/77.0-13800.0/7.0*t26+13800.0/7.0*t29-6900.0/7.0*t32+13800.0/7.0*t35-4600.0/7.0*t37+13800.0/7.0*t39-6900.0/7.0*t41+13800.0/7.0*t44-27600.0/7.0*t46+13800.0/7.0*t49+4100.0/7.0*t51-8200.0/7.0*t53+8200.0/7.0*t56-271200.0/77.0*t10+t303;
1424  t381 = -1701720.0/77.0*t2+138000.0/7.0*t4+487320.0/77.0*eta+53840.0/77.0*xi+4100.0/7.0*t65-4100.0/7.0*t68+8200.0/7.0*t71+13800.0/7.0*t73-27600.0/7.0*t75-1049820.0/77.0*t12+166500.0/77.0*t14+2308200.0/77.0*t16+9200.0/7.0*t80+9200.0/7.0*t82+t347-t348+t349;
1425  t383 = 9000.0*t26;
1426  t387 = 3000.0*t37;
1427  t388 = 9000.0*t39;
1428  t391 = 18000.0*t46;
1429  t392 = 9000.0*t49;
1430  t399 = -t383+5040.0*t29-2520.0*t32+5040.0*t35-t387+t388-2520.0*t41+5040.0*t44-t391+t392-2520.0*t51+5040.0*t53-5040.0*t56+26640.0/11.0*t10-71856.0/11.0*t8-1508400.0/11.0*t2;
1431  t411 = 6000.0*t80;
1432  t412 = 6000.0*t82;
1433  t413 = 6000.0*t108;
1434  t414 = 3000.0*t110;
1435  t415 = 6000.0*t112;
1436  t416 = -49800.0/11.0+111600.0*t4+518400.0/11.0*eta+142848.0/11.0*xi-2520.0*t65+2520.0*t68-5040.0*t71+5040.0*t73-10080.0*t75-980064.0/11.0*t12+79920.0/11.0*t14+1695600.0/11.0*t16+t411+t412-t413+t414-t415;
1437  t418 = 2880.0*t26;
1438  t419 = 2880.0*t29;
1439  t420 = 1440.0*t32;
1440  t421 = 2880.0*t35;
1441  t422 = 960.0*t37;
1442  t423 = 2880.0*t39;
1443  t424 = 1440.0*t41;
1444  t425 = 2880.0*t44;
1445  t426 = 5760.0*t46;
1446  t427 = 2880.0*t49;
1447  t428 = 6480.0*t51;
1448  t429 = 12960.0*t53;
1449  t430 = 12960.0*t56;
1450  t432 = 103968.0/11.0*t8;
1451  t433 = 9600.0/11.0+t418-t419+t420-t421+t422-t423+t424-t425+t426-t427-t428+t429-t430+38520.0/11.0*t10+t432;
1452  t438 = 6480.0*t65;
1453  t439 = 6480.0*t68;
1454  t440 = 12960.0*t71;
1455  t441 = 2880.0*t73;
1456  t442 = 5760.0*t75;
1457  t446 = 1920.0*t80;
1458  t447 = 1920.0*t82;
1459  t448 = 368640.0/11.0*t2-28800.0*t4-115200.0/11.0*eta-69408.0/11.0*xi-t438+t439-t440-t441+t442+502560.0/11.0*t12-490320.0/11.0*t14-751680.0/11.0*t16-t446-t447-t413+t414-t415;
1460  t450 = 12960.0*t29;
1461  t451 = 6480.0*t32;
1462  t452 = 12960.0*t35;
1463  t453 = 6480.0*t41;
1464  t454 = 12960.0*t44;
1465  t455 = 1440.0*t51;
1466  t456 = 2880.0*t53;
1467  t457 = 2880.0*t56;
1468  t459 = -10200.0/11.0-t383+t450-t451+t452-t387+t388-t453+t454-t391+t392+t455-t456+t457-163440.0/11.0*t10+t432;
1469  t464 = 1440.0*t65;
1470  t465 = 1440.0*t68;
1471  t466 = 2880.0*t71;
1472  t467 = 12960.0*t73;
1473  t468 = 25920.0*t75;
1474  t472 = 1920.0*t108;
1475  t473 = 960.0*t110;
1476  t474 = 1920.0*t112;
1477  t475 = -795600.0/11.0*t2+68400.0*t4+201600.0/11.0*eta+30384.0/11.0*xi+t464-t465+t466+t467-t468-664056.0/11.0*t12+115560.0/11.0*t14+1505520.0/11.0*t16+t411+t412+t472-t473+t474;
1478  t477 = 12960.0*t26;
1479  t478 = 28800.0*t29;
1480  t479 = 14400.0*t32;
1481  t480 = 28800.0*t35;
1482  t481 = 4320.0*t37;
1483  t482 = 12960.0*t39;
1484  t483 = 14400.0*t41;
1485  t484 = 28800.0*t44;
1486  t485 = 25920.0*t46;
1487  t486 = 12960.0*t49;
1488  t487 = 9360.0*t51;
1489  t488 = 18720.0*t53;
1490  t489 = 18720.0*t56;
1491  t491 = 204912.0/11.0*t8;
1492  t492 = 56400.0/11.0+t477-t478+t479-t480+t481-t482+t483-t484+t485-t486-t487+t488-t489+290160.0/11.0*t10+t491;
1493  t497 = 9360.0*t65;
1494  t498 = 9360.0*t68;
1495  t499 = 18720.0*t71;
1496  t500 = 28800.0*t73;
1497  t501 = 57600.0*t75;
1498  t505 = 8640.0*t80;
1499  t506 = 8640.0*t82;
1500  t507 = 3360.0*t108;
1501  t508 = 1680.0*t110;
1502  t509 = 3360.0*t112;
1503  t510 = 1896480.0/11.0*t2-144000.0*t4-624000.0/11.0*eta-349824.0/11.0*xi-t497+t498-t499-t500+t501+2627424.0/11.0*t12-1529280.0/11.0*t14-4000320.0/11.0*t16-t505-t506+t507-t508+t509;
1504  t512 = 5040.0*t26;
1505  t513 = 18720.0*t29;
1506  t514 = 9360.0*t32;
1507  t515 = 18720.0*t35;
1508  t516 = 1680.0*t37;
1509  t517 = 5040.0*t39;
1510  t518 = 9360.0*t41;
1511  t519 = 18720.0*t44;
1512  t520 = 10080.0*t46;
1513  t521 = 5040.0*t49;
1514  t522 = 14400.0*t51;
1515  t523 = 28800.0*t53;
1516  t524 = 28800.0*t56;
1517  t526 = 43200.0/11.0+t512+t513-t514+t515+t516-t517-t518+t519+t520-t521+t522-t523+t524-509760.0/11.0*t10+t491;
1518  t529 = 412800.0/11.0*eta;
1519  t531 = 14400.0*t65;
1520  t532 = 14400.0*t68;
1521  t533 = 28800.0*t71;
1522  t534 = 18720.0*t73;
1523  t535 = 37440.0*t75;
1524  t539 = 3360.0*t80;
1525  t540 = 3360.0*t82;
1526  t541 = 8640.0*t108;
1527  t542 = 4320.0*t110;
1528  t543 = 8640.0*t112;
1529  t544 = 1120320.0/11.0*t2-79200.0*t4-t529-39360.0/11.0*xi+t531-t532+t533+t534-t535-268128.0/11.0*t12+870480.0/11.0*t14+276480.0/11.0*t16-t539-t540+t541-t542+t543;
1530  t547 = 251280.0/11.0*t8;
1531  t548 = -22800.0/11.0-t418+t513-t514+t515-t422+t423-t518+t519-t426+t427+t522-t523+t524-224640.0/11.0*t10-t547;
1532  t556 = -606240.0/11.0*t2+43200.0*t4+220800.0/11.0*eta+228864.0/11.0*xi+t531-t532+t533+t534-t535-1567008.0/11.0*t12+1559520.0/11.0*t14+2034720.0/11.0*t16+t446+t447+t507-t508+t509;
1533  t559 = 76608.0/11.0*t8;
1534  t560 = 3600.0/11.0-t418-t450+t451-t452-t422+t423+t453-t454-t426+t427-t455+t456-t457+92160.0/11.0*t10+t559;
1535  t563 = 9600.0/11.0*eta;
1536  t568 = -131040.0/11.0*t2+14400.0*t4+t563-93744.0/11.0*xi-t464+t465-t466-t467+t468+561888.0/11.0*t12-578880.0/11.0*t14-531360.0/11.0*t16+t446+t447+t541-t542+t543;
1537  t571 = -9600.0/11.0+t512-t478+t479-t480+t516-t517+t483-t484+t520-t521-t487+t488-t489+519840.0/11.0*t10-t547;
1538  t578 = 169920.0/11.0*t2-21600.0*t4+t563-37248.0/11.0*xi-t497+t498-t499-t500+t501+1029168.0/11.0*t12-673920.0/11.0*t14-1909440.0/11.0*t16-t539-t540-t472+t473-t474;
1539  t581 = 30000.0/11.0+t477+t419-t420+t421+t481-t482-t424+t425+t485-t486+t428-t429+t430-192960.0/11.0*t10+t559;
1540  t588 = 1421280.0/11.0*t2-115200.0*t4-t529-60480.0/11.0*xi+t438-t439+t440+t441-t442+565056.0/11.0*t12+276480.0/11.0*t14-1434240.0/11.0*t16-t505-t506-t472+t473-t474;
1541  t605 = -50400.0/11.0-15120.0*t26+15120.0*t29-7560.0*t32+15120.0*t35-5040.0*t37+15120.0*t39-7560.0*t41+15120.0*t44-30240.0*t46+15120.0*t49-7560.0*t51+15120.0*t53-15120.0*t56+68040.0/11.0*t10-113400.0/11.0*t8;
1542  t623 = -1935360.0/11.0*t2+151200.0*t4+604800.0/11.0*eta+211008.0/11.0*xi-7560.0*t65+7560.0*t68-15120.0*t71+15120.0*t73-30240.0*t75-1640520.0/11.0*t12+204120.0/11.0*t14+2948400.0/11.0*t16+10080.0*t80+10080.0*t82-10080.0*t108+5040.0*t110-10080.0*t112;
1543 
1544  values[0] = 420.0-t1+t3-t5;
1545  values[1] = -420.0+3360.0*xi-7560.0*t8+5040.0*t10+t1-15120.0*t12-t3+15120.0*t14+15120.0*t16+t5;
1546  values[2] = 0.0;
1547  values[3] = t59+t84;
1548  values[4] = t96+t114;
1549  values[5] = t119+t125;
1550  values[6] = t142+t160;
1551  values[7] = t172+t188;
1552  values[8] = t204+t219;
1553  values[9] = t235+t251;
1554  values[10] = t262+t276;
1555  values[11] = t281+t287;
1556  values[12] = t304+t322;
1557  values[13] = t334+t350;
1558  values[14] = t366+t381;
1559  values[15] = t399+t416;
1560  values[16] = t433+t448;
1561  values[17] = t459+t475;
1562  values[18] = t492+t510;
1563  values[19] = t526+t544;
1564  values[20] = t548+t556;
1565  values[21] = t560+t568;
1566  values[22] = t571+t578;
1567  values[23] = t581+t588;
1568  values[24] = t605+t623;
1569 }
1570 
1571 static int N_T_P5_2D_ChangeJ0[2] = { 3, 9 };
1572 static int N_T_P5_2D_ChangeJ1[2] = { 4, 10 };
1573 static int N_T_P5_2D_ChangeJ2[2] = { 5, 11 };
1574 
1575 static int *N_T_P5_2D_Change[3] = { N_T_P5_2D_ChangeJ0, N_T_P5_2D_ChangeJ1,
1576  N_T_P5_2D_ChangeJ2 };
1577 
1578 // ***********************************************************************
1579 
1580 TBaseFunct2D *BF_N_T_P5_2D_Obj = new TBaseFunct2D
1581  (25, BF_N_T_P5_2D, BFUnitTriangle,
1582  N_T_P5_2D_Funct, N_T_P5_2D_DeriveXi,
1583  N_T_P5_2D_DeriveEta, N_T_P5_2D_DeriveXiXi,
1584  N_T_P5_2D_DeriveXiEta, N_T_P5_2D_DeriveEtaEta, 6, 5,
1585  2, N_T_P5_2D_Change);
Definition: BaseFunct2D.h:27