ParMooN
 All Classes Functions Variables Friends Pages
BF_C_Q_M7_2D.h
1 // ***********************************************************************
2 // M7 element with bubbles, conforming, 2D
3 // ***********************************************************************
4 
5 // base function values
6 static void C_Q_M7_2D_Funct(double xi, double eta, double *values)
7 {
8  double t1 = xi*xi;
9  double t2 = eta*eta;
10  double t3 = t2*eta;
11  double t4 = t1*t3;
12  double t6 = t1*xi;
13  double t7 = t6*t2;
14  double t9 = t1*t1;
15  double t10 = t9*eta;
16  double t12 = xi*t3;
17  double t14 = t1*t2;
18  double t16 = t6*eta;
19  double t18 = xi*t2;
20  double t20 = t1*eta;
21  double t22 = xi*eta;
22  double t24 = t2*t2;
23  double t25 = xi*t24;
24  double t27 = t9*xi;
25  double t28 = t27*eta;
26  double t29 = 184877.0/207360.0*t28;
27  double t30 = t9*t2;
28  double t32 = t6*t3;
29  double t34 = t9*t1;
30  double t35 = t34*eta;
31  double t36 = 117649.0/184320.0*t35;
32  double t37 = t27*t2;
33  double t38 = 85829.0/76800.0*t37;
34  double t42 = 12113881.0/4147200.0*t27;
35  double t43 = 967443.0/320000.0*t4+967443.0/320000.0*t7+2529463.0/921600.0*t10+971869.0/5184000.0*t12+225377.0/300000.0*t14+971869.0/5184000.0*t16-45493571.0/23040000.0*t18-45493571.0/23040000.0*t20-802049.0/864000.0*t22+2529463.0/921600.0*t25+t29-98177.0/384000.0*t30-94061.0/96000.0*t32-t36-t38+1478113.0/69120000.0*eta+1478113.0/69120000.0*xi-49159793.0/115200000.0*t1+t42;
36  double t53 = (1.0-t1)*(63.0/8.0*t27-35.0/4.0*t6+15.0/8.0*xi)*(1.0+eta);
37  double t54 = 16807.0/207360.0*t53;
38  double t56 = t24*eta;
39  double t63 = (1.0-t2)*(63.0/8.0*t56-35.0/4.0*t3+15.0/8.0*eta)*(1.0+xi);
40  double t64 = 16807.0/207360.0*t63;
41  double t65 = 117649.0/184320.0*t34;
42  double t66 = t9*t6;
43  double t67 = 117649.0/92160.0*t66;
44  double t69 = 12113881.0/4147200.0*t56;
45  double t72 = t24*t2;
46  double t73 = 117649.0/184320.0*t72;
47  double t74 = t24*t3;
48  double t75 = 117649.0/92160.0*t74;
49  double t76 = t6*t24;
50  double t78 = t1*t56;
51  double t79 = 85829.0/76800.0*t78;
52  double t80 = xi*t72;
53  double t81 = 117649.0/184320.0*t80;
54  double t82 = t9*t3;
55  double t84 = t1*t24;
56  double t86 = xi*t56;
57  double t87 = 184877.0/207360.0*t86;
58  double t88 = -80149.0/6400000.0-172582507.0/103680000.0*t6-922751.0/4608000.0*t9-t54-t64+t65-t67-49159793.0/115200000.0*t2+t69-172582507.0/103680000.0*t3-922751.0/4608000.0*t24+t73-t75-87887.0/38400.0*t76-t79-t81-87887.0/38400.0*t82-98177.0/384000.0*t84+t87;
59  double t94 = 141169.0/100000.0*t14;
60  double t98 = 1526693.0/1728000.0*t22;
61  double t100 = 84035.0/20736.0*t28;
62  double t102 = 556003.0/288000.0*t32;
63  double t103 = 117649.0/36864.0*t35;
64  double t104 = 16807.0/7680.0*t37;
65  double t106 = -633619.0/160000.0*t4-3142223.0/720000.0*t7-1370971.0/184320.0*t10-556003.0/288000.0*t12-t94+16428671.0/5184000.0*t16+6266267.0/2880000.0*t18+93337111.0/23040000.0*t20+t98-136171.0/57600.0*t25-t100+45619.0/76800.0*t30+t102+t103+t104+1501213.0/7680000.0*eta;
66  double t109 = 14101073.0/829440.0*t27;
67  double t112 = 117649.0/207360.0*t53;
68  double t113 = 117649.0/36864.0*t34;
69  double t114 = 823543.0/92160.0*t66;
70  double t116 = 14651.0/25600.0*t56;
71  double t120 = 14651.0/25600.0*t78;
72  double t123 = 517293.0/12800000.0-30836729.0/34560000.0*xi-11258387.0/115200000.0*t1-t109+928658437.0/103680000.0*t6+2994047.0/921600.0*t9+t112-t113+t114+7849849.0/9600000.0*t2-t116+700553.0/1920000.0*t3-324919.0/384000.0*t24+136171.0/57600.0*t76+t120+55223.0/15360.0*t82+324919.0/384000.0*t84;
73  double t129 = 2277569.0/800000.0*t14;
74  double t133 = 3412507.0/576000.0*t22;
75  double t135 = 16807.0/2160.0*t28;
76  double t137 = 22981.0/8000.0*t32;
77  double t138 = 117649.0/20480.0*t35;
78  double t139 = 16807.0/3840.0*t37;
79  double t141 = -1021993.0/160000.0*t4-1934177.0/240000.0*t7+64827.0/20480.0*t10-22981.0/8000.0*t12-t129-23683121.0/1728000.0*t16+1767479.0/480000.0*t18+5465313.0/2560000.0*t20+t133-6517.0/1600.0*t25+t135+2401.0/1600.0*t30+t137-t138+t139+1137437.0/2560000.0*eta;
80  double t144 = 12285917.0/276480.0*t27;
81  double t147 = 117649.0/69120.0*t53;
82  double t148 = 117649.0/20480.0*t34;
83  double t149 = 823543.0/30720.0*t66;
84  double t151 = 14651.0/12800.0*t56;
85  double t155 = 14651.0/12800.0*t78;
86  double t158 = -1368129.0/12800000.0+42876029.0/11520000.0*xi+59592379.0/12800000.0*t1+t144-737881837.0/34560000.0*t6-1054039.0/102400.0*t9-t147+t148-t149+1077069.0/800000.0*t2-t151+60809.0/80000.0*t3-83153.0/64000.0*t24+6517.0/1600.0*t76+t155+7203.0/1280.0*t82+83153.0/64000.0*t84;
87  double t164 = 12372059.0/1280000.0*t14;
88  double t168 = 139601.0/10800.0*t22;
89  double t170 = 386561.0/41472.0*t28;
90  double t172 = 62083.0/6400.0*t32;
91  double t173 = 16807.0/1024.0*t37;
92  double t174 = 45521.0/320000.0+3889473.0/256000.0*t4+46991.0/2000.0*t7+362551.0/36864.0*t10+62083.0/6400.0*t12+t164+23065721.0/1036800.0*t16-906549.0/128000.0*t18-3445729.0/288000.0*t20-t168+31899.0/5120.0*t25-t170-69629.0/10240.0*t30-t172+t103-t173;
93  double t178 = 11378339.0/165888.0*t27;
94  double t181 = 117649.0/41472.0*t53;
95  double t182 = 823543.0/18432.0*t66;
96  double t184 = 14651.0/20480.0*t56;
97  double t188 = 14651.0/20480.0*t78;
98  double t191 = -67963.0/64000.0*eta-45345629.0/6912000.0*xi-52741003.0/5760000.0*t1-t178+631843387.0/20736000.0*t6+2249737.0/184320.0*t9+t181-t113+t182-1834217.0/640000.0*t2+t184+1519.0/32000.0*t3+309533.0/102400.0*t24-31899.0/5120.0*t76-t188-31213.0/2048.0*t82-309533.0/102400.0*t84;
99  double t197 = 2978171.0/320000.0*t14;
100  double t201 = 1693391.0/172800.0*t22;
101  double t204 = 7889.0/1200.0*t32;
102  double t205 = -211043.0/2560000.0-287581.0/32000.0*t4-3850861.0/192000.0*t7-545027.0/36864.0*t10-7889.0/1200.0*t12-t197-19824371.0/1036800.0*t16+1399097.0/384000.0*t18+50965439.0/4608000.0*t20+t201-4459.0/1920.0*t25+t170+40817.0/5120.0*t30+t204+t103+t173;
103  double t212 = 1127.0/2560.0*t56;
104  double t216 = 1127.0/2560.0*t78;
105  double t219 = 272979.0/512000.0*eta+42104279.0/6912000.0*xi+134554637.0/23040000.0*t1+t178-622119337.0/20736000.0*t6-472997.0/184320.0*t9-t181-t113-t182+854217.0/640000.0*t2-t212-49.0/125.0*t3-12201.0/12800.0*t24+4459.0/1920.0*t76+t216+2401.0/256.0*t82+12201.0/12800.0*t84;
106  double t225 = 13696039.0/4800000.0*t14;
107  double t229 = 1971907.0/576000.0*t22;
108  double t232 = 4459.0/12000.0*t32;
109  double t234 = 1083733.0/960000.0*t4+212317.0/40000.0*t7+199283.0/20480.0*t10+4459.0/12000.0*t12+t225+19361321.0/1728000.0*t16-446929.0/480000.0*t18-10357963.0/2560000.0*t20-t229+18179.0/19200.0*t25-t135-31213.0/12800.0*t30-t232-t138-t139+153713.0/2560000.0*eta;
110  double t240 = 14651.0/76800.0*t56;
111  double t244 = 14651.0/76800.0*t78;
112  double t245 = 2401.0/2560.0*t82;
113  double t247 = -285621.0/12800000.0-38554229.0/11520000.0*xi+8088871.0/12800000.0*t1-t144+724916437.0/34560000.0*t6-650671.0/102400.0*t9+t147+t148+t149-497791.0/1200000.0*t2+t240-91679.0/480000.0*t3+144893.0/384000.0*t24-18179.0/19200.0*t76-t244-t245-144893.0/384000.0*t84;
114  double t253 = 3234931.0/2400000.0*t14;
115  double t257 = 2074807.0/1728000.0*t22;
116  double t260 = 14749.0/96000.0*t32;
117  double t262 = -169393.0/480000.0*t4+994357.0/480000.0*t7-765919.0/184320.0*t10-14749.0/96000.0*t12+t253-27233171.0/5184000.0*t16+37387.0/320000.0*t18+22266139.0/23040000.0*t20+t257-3087.0/12800.0*t25+t100-26411.0/19200.0*t30+t260+t103-t104-6321.0/2560000.0*eta;
118  double t268 = 637.0/12800.0*t56;
119  double t272 = 637.0/12800.0*t78;
120  double t275 = 1013957.0/12800000.0+41641229.0/34560000.0*xi-242622863.0/115200000.0*t1+t109-961071937.0/103680000.0*t6+4809203.0/921600.0*t9-t112-t113-t114+5537.0/200000.0*t2-t268+12887.0/320000.0*t3-1519.0/16000.0*t24+3087.0/12800.0*t76+t272+2401.0/7680.0*t82+1519.0/16000.0*t84;
121  double t280 = 7123181.0/5184000.0*t12;
122  double t281 = 20141003.0/19200000.0*t14;
123  double t282 = 7123181.0/5184000.0*t16;
124  double t285 = 273563.0/432000.0*t22;
125  double t288 = 83771.0/144000.0*t32;
126  double t292 = 2126359.0/3840000.0*t4-940003.0/720000.0*t7+947497.0/921600.0*t10+t280-t281+t282+5917771.0/23040000.0*t18-3953677.0/11520000.0*t20-t285-243521.0/307200.0*t25-t29+594629.0/768000.0*t30-t288-t36+t38-6820919.0/34560000.0*eta-17668213.0/69120000.0*xi+36656459.0/57600000.0*t1;
127  double t296 = 5077697.0/8294400.0*t56;
128  double t300 = 85829.0/307200.0*t78;
129  double t303 = -585327.0/12800000.0-t42+196867657.0/103680000.0*t6-5668649.0/4608000.0*t9+t54+t64+t65+t67+39583543.0/115200000.0*t2-t296+21015133.0/25920000.0*t3-2159737.0/2304000.0*t24+t73+77597.0/230400.0*t76-t300+t81-87887.0/153600.0*t82+246511.0/512000.0*t84-t87;
130  double t309 = 4969237.0/4800000.0*t14;
131  double t313 = 326879.0/216000.0*t22;
132  double t316 = 29841.0/64000.0*t32;
133  double t317 = 14651.0/25600.0*t37;
134  double t321 = -6341041.0/5760000.0*t4+7623959.0/5760000.0*t7-136171.0/230400.0*t10-14426923.0/2592000.0*t12+t309-29841.0/64000.0*t16-10863937.0/7680000.0*t18+6379457.0/11520000.0*t20+t313+818741.0/184320.0*t25-41797.0/192000.0*t30+t316-t317+37473779.0/34560000.0*eta+419979.0/2560000.0*xi-1367737.0/4800000.0*t1;
135  double t322 = 14651.0/25600.0*t27;
136  double t325 = 117649.0/207360.0*t63;
137  double t327 = 2907611.0/829440.0*t56;
138  double t330 = 117649.0/36864.0*t72;
139  double t332 = 16807.0/30720.0*t78;
140  double t333 = 117649.0/36864.0*t80;
141  double t336 = 84035.0/20736.0*t86;
142  double t337 = 1013957.0/12800000.0+t322-927227.0/1280000.0*t6+41797.0/192000.0*t9-t325-170592863.0/115200000.0*t2+t327-59484089.0/12960000.0*t3+4232963.0/921600.0*t24-t330-55223.0/92160.0*t76+t332-t333+136171.0/230400.0*t82-2401.0/3200.0*t84+t336;
143  double t343 = 2564513.0/1600000.0*t14;
144  double t347 = 1251607.0/576000.0*t22;
145  double t350 = 14063.0/16000.0*t32;
146  double t351 = 14651.0/12800.0*t37;
147  double t355 = -1841567.0/960000.0*t4+180859.0/80000.0*t7-6517.0/6400.0*t10+17200421.0/1728000.0*t12+t343-14063.0/16000.0*t16+5099087.0/2560000.0*t18+1582259.0/1920000.0*t20-t347-160867.0/20480.0*t25-34153.0/64000.0*t30+t350-t351-16460227.0/5760000.0*eta+303163.0/2560000.0*xi-1827651.0/3200000.0*t1;
148  double t356 = 14651.0/12800.0*t27;
149  double t359 = 117649.0/69120.0*t63;
150  double t361 = 1226911.0/138240.0*t56;
151  double t364 = 117649.0/20480.0*t72;
152  double t365 = 2401.0/2560.0*t76;
153  double t366 = 16807.0/15360.0*t78;
154  double t367 = 117649.0/20480.0*t80;
155  double t370 = 16807.0/2160.0*t86;
156  double t371 = -285621.0/12800000.0+t356-423311.0/320000.0*t6+34153.0/64000.0*t9+t359+26096371.0/12800000.0*t2-t361+50686139.0/4320000.0*t3-794731.0/102400.0*t24+t364-t365+t366+t367+6517.0/6400.0*t82-26411.0/25600.0*t84-t370;
157  double t377 = 1454467.0/640000.0*t14;
158  double t381 = 1911539.0/691200.0*t22;
159  double t384 = 35329.0/76800.0*t32;
160  double t385 = 14651.0/20480.0*t37;
161  double t388 = -211043.0/2560000.0+112847.0/16000.0*t4-771211.0/192000.0*t7+31899.0/20480.0*t10-25062667.0/2073600.0*t12-t377+35329.0/76800.0*t16+1815107.0/2304000.0*t18-1510229.0/512000.0*t20+t381+26411.0/9216.0*t25-11221.0/51200.0*t30-t384+t385+31917179.0/6912000.0*eta-236621.0/512000.0*xi;
162  double t390 = 14651.0/20480.0*t27;
163  double t393 = 117649.0/41472.0*t63;
164  double t395 = 4453855.0/331776.0*t56;
165  double t399 = 16807.0/4096.0*t78;
166  double t402 = 386561.0/41472.0*t86;
167  double t403 = 207809.0/1280000.0*t1-t390+756021.0/512000.0*t6+11221.0/51200.0*t9-t393-501613.0/23040000.0*t2+t395-748234949.0/41472000.0*t3+607453.0/184320.0*t24-t330+31213.0/12288.0*t76-t399-t333-31899.0/20480.0*t82+21609.0/10240.0*t84+t402;
168  double t409 = 367059.0/1280000.0*t14;
169  double t413 = 612941.0/172800.0*t22;
170  double t416 = 1029.0/3200.0*t32;
171  double t417 = 1127.0/2560.0*t37;
172  double t420 = 45521.0/320000.0-4756381.0/768000.0*t4+50617.0/24000.0*t7-4459.0/7680.0*t10+13341671.0/1036800.0*t12+t409+1029.0/3200.0*t16-19272239.0/4608000.0*t18+3210137.0/1536000.0*t20-t413+256907.0/36864.0*t25-147.0/1600.0*t30-t416-t417-65455033.0/13824000.0*eta+207221.0/512000.0*xi;
173  double t422 = 1127.0/2560.0*t27;
174  double t431 = 166257.0/2560000.0*t1+t422-34937.0/64000.0*t6+147.0/1600.0*t9+t393-62402137.0/23040000.0*t2-t395+376548487.0/20736000.0*t3+530621.0/92160.0*t24-t330-2401.0/1536.0*t76+t399-t333+4459.0/7680.0*t82-7203.0/20480.0*t84-t402;
175  double t437 = 335209.0/1200000.0*t14;
176  double t441 = 1611757.0/576000.0*t22;
177  double t444 = 24353.0/96000.0*t32;
178  double t445 = 14651.0/76800.0*t37;
179  double t449 = 98441.0/80000.0*t4-211729.0/480000.0*t7+18179.0/76800.0*t10-18280871.0/1728000.0*t12+t437+24353.0/96000.0*t16+25791689.0/7680000.0*t18-261709.0/1920000.0*t20+t441-549829.0/61440.0*t25+9359.0/192000.0*t30-t444+t445+4250113.0/1440000.0*eta-393813.0/2560000.0*xi+229957.0/2400000.0*t1;
180  double t450 = 14651.0/76800.0*t27;
181  double t459 = -1368129.0/12800000.0-t450+546791.0/1920000.0*t6-9359.0/192000.0*t9-t359+35582379.0/12800000.0*t2+t361-204365231.0/17280000.0*t3-861959.0/102400.0*t24+t364+2401.0/15360.0*t76-t366+t367-18179.0/76800.0*t82-2401.0/6400.0*t84+t370;
182  double t465 = 371273.0/800000.0*t14;
183  double t469 = 1714657.0/1728000.0*t22;
184  double t472 = 7889.0/144000.0*t32;
185  double t473 = 637.0/12800.0*t37;
186  double t477 = 506611.0/960000.0*t4+356083.0/2880000.0*t7-3087.0/51200.0*t10+26152721.0/5184000.0*t12+t465-7889.0/144000.0*t16-5661313.0/7680000.0*t18+74431.0/3840000.0*t20-t469+717899.0/184320.0*t25-3479.0/256000.0*t30+t472-t473-17386327.0/17280000.0*eta+259063.0/7680000.0*xi-269059.0/6400000.0*t1;
187  double t478 = 637.0/12800.0*t27;
188  double t487 = 517293.0/12800000.0+t478-137347.0/1920000.0*t6+3479.0/256000.0*t9+t325-128307137.0/115200000.0*t2-t327+467769337.0/103680000.0*t3+3930437.0/921600.0*t24-t330-2401.0/46080.0*t76-t332-t333+3087.0/51200.0*t82-21609.0/51200.0*t84-t336;
189  double t501 = 85829.0/307200.0*t37;
190  double t505 = -716051.0/5760000.0*t4-716051.0/5760000.0*t7-13827.0/25600.0*t10-10198837.0/10368000.0*t12-1343687.0/3200000.0*t14-10198837.0/10368000.0*t16-496649.0/5760000.0*t18-496649.0/5760000.0*t20+839329.0/3456000.0*t22-13827.0/25600.0*t25+t29+84457.0/256000.0*t30+73481.0/384000.0*t32+t36+t501+9594313.0/69120000.0*eta+9594313.0/69120000.0*xi+18298957.0/115200000.0*t1;
191  double t506 = 5077697.0/8294400.0*t27;
192  double t515 = -80149.0/6400000.0+t506-155978489.0/207360000.0*t6-3621101.0/4608000.0*t9-t54-t64+t65+18298957.0/115200000.0*t2+t296-155978489.0/207360000.0*t3-3621101.0/4608000.0*t24+t73+77597.0/921600.0*t76+t300+t81+77597.0/921600.0*t82+84457.0/256000.0*t84+t87;
193  double t526 = 16807.0/30720.0*t37;
194  double t528 = 356083.0/2880000.0*t4+506611.0/960000.0*t7+717899.0/184320.0*t10-7889.0/144000.0*t12+t465+26152721.0/5184000.0*t16+74431.0/3840000.0*t18-5661313.0/7680000.0*t20-t469-3087.0/51200.0*t25-t100-21609.0/51200.0*t30+t472-t103-t526+259063.0/7680000.0*eta;
195  double t531 = 2907611.0/829440.0*t27;
196  double t540 = 517293.0/12800000.0-17386327.0/17280000.0*xi-128307137.0/115200000.0*t1-t531+467769337.0/103680000.0*t6+3930437.0/921600.0*t9+t112-t113-269059.0/6400000.0*t2+t268-137347.0/1920000.0*t3+3479.0/256000.0*t24+3087.0/51200.0*t76-t272-2401.0/46080.0*t82-3479.0/256000.0*t84;
197  double t551 = 16807.0/15360.0*t37;
198  double t553 = -211729.0/480000.0*t4+98441.0/80000.0*t7-549829.0/61440.0*t10+24353.0/96000.0*t12+t437-18280871.0/1728000.0*t16-261709.0/1920000.0*t18+25791689.0/7680000.0*t20+t441+18179.0/76800.0*t25+t135-2401.0/6400.0*t30-t444+t138-t551-393813.0/2560000.0*eta;
199  double t556 = 1226911.0/138240.0*t27;
200  double t565 = -1368129.0/12800000.0+4250113.0/1440000.0*xi+35582379.0/12800000.0*t1+t556-204365231.0/17280000.0*t6-861959.0/102400.0*t9-t147+t148+229957.0/2400000.0*t2-t240+546791.0/1920000.0*t3-9359.0/192000.0*t24-18179.0/76800.0*t76+t244+2401.0/15360.0*t82+9359.0/192000.0*t84;
201  double t576 = 16807.0/4096.0*t37;
202  double t577 = 45521.0/320000.0+50617.0/24000.0*t4-4756381.0/768000.0*t7+256907.0/36864.0*t10+1029.0/3200.0*t12+t409+13341671.0/1036800.0*t16+3210137.0/1536000.0*t18-19272239.0/4608000.0*t20-t413-4459.0/7680.0*t25-t170-7203.0/20480.0*t30-t416-t103+t576;
203  double t581 = 4453855.0/331776.0*t27;
204  double t590 = 207221.0/512000.0*eta-65455033.0/13824000.0*xi-62402137.0/23040000.0*t1-t581+376548487.0/20736000.0*t6+530621.0/92160.0*t9+t181-t113+166257.0/2560000.0*t2+t212-34937.0/64000.0*t3+147.0/1600.0*t24+4459.0/7680.0*t76-t216-2401.0/1536.0*t82-147.0/1600.0*t84;
205  double t601 = -211043.0/2560000.0-771211.0/192000.0*t4+112847.0/16000.0*t7+26411.0/9216.0*t10+35329.0/76800.0*t12-t377-25062667.0/2073600.0*t16-1510229.0/512000.0*t18+1815107.0/2304000.0*t20+t381+31899.0/20480.0*t25+t170+21609.0/10240.0*t30-t384-t103-t576;
206  double t613 = -236621.0/512000.0*eta+31917179.0/6912000.0*xi-501613.0/23040000.0*t1+t581-748234949.0/41472000.0*t6+607453.0/184320.0*t9-t181-t113+207809.0/1280000.0*t2-t184+756021.0/512000.0*t3+11221.0/51200.0*t24-31899.0/20480.0*t76+t188+31213.0/12288.0*t82-11221.0/51200.0*t84;
207  double t625 = 180859.0/80000.0*t4-1841567.0/960000.0*t7-160867.0/20480.0*t10-14063.0/16000.0*t12+t343+17200421.0/1728000.0*t16+1582259.0/1920000.0*t18+5099087.0/2560000.0*t20-t347-6517.0/6400.0*t25-t135-26411.0/25600.0*t30+t350+t138+t551+303163.0/2560000.0*eta;
208  double t635 = -285621.0/12800000.0-16460227.0/5760000.0*xi+26096371.0/12800000.0*t1-t556+50686139.0/4320000.0*t6-794731.0/102400.0*t9+t147+t148-1827651.0/3200000.0*t2+t151-423311.0/320000.0*t3+34153.0/64000.0*t24+6517.0/6400.0*t76-t155-t245-34153.0/64000.0*t84;
209  double t647 = 7623959.0/5760000.0*t4-6341041.0/5760000.0*t7+818741.0/184320.0*t10-29841.0/64000.0*t12+t309-14426923.0/2592000.0*t16+6379457.0/11520000.0*t18-10863937.0/7680000.0*t20+t313-136171.0/230400.0*t25+t100-2401.0/3200.0*t30+t316-t103+t526+419979.0/2560000.0*eta;
210  double t658 = 1013957.0/12800000.0+37473779.0/34560000.0*xi-170592863.0/115200000.0*t1+t531-59484089.0/12960000.0*t6+4232963.0/921600.0*t9-t112-t113-1367737.0/4800000.0*t2+t116-927227.0/1280000.0*t3+41797.0/192000.0*t24+136171.0/230400.0*t76-t120-55223.0/92160.0*t82-41797.0/192000.0*t84;
211  double t670 = -940003.0/720000.0*t4+2126359.0/3840000.0*t7-243521.0/307200.0*t10+t280-t281+t282-3953677.0/11520000.0*t18+5917771.0/23040000.0*t20-t285+947497.0/921600.0*t25-t29+246511.0/512000.0*t30-t288+t36-t501-17668213.0/69120000.0*eta-6820919.0/34560000.0*xi+39583543.0/115200000.0*t1;
212  double t679 = -585327.0/12800000.0-t506+21015133.0/25920000.0*t6-2159737.0/2304000.0*t9+t54+t64+t65+36656459.0/57600000.0*t2-t69+196867657.0/103680000.0*t3-5668649.0/4608000.0*t24+t73+t75-87887.0/153600.0*t76+t79-t81+77597.0/230400.0*t82+594629.0/768000.0*t84-t87;
213  double t693 = 994357.0/480000.0*t4-169393.0/480000.0*t7-3087.0/12800.0*t10-27233171.0/5184000.0*t12+t253-14749.0/96000.0*t16+22266139.0/23040000.0*t18+37387.0/320000.0*t20+t257-765919.0/184320.0*t25+1519.0/16000.0*t30+t260+t473+41641229.0/34560000.0*eta-6321.0/2560000.0*xi+5537.0/200000.0*t1;
214  double t697 = 14101073.0/829440.0*t56;
215  double t700 = 823543.0/92160.0*t74;
216  double t702 = 16807.0/7680.0*t78;
217  double t705 = 1013957.0/12800000.0-t478+12887.0/320000.0*t6-1519.0/16000.0*t9-t325-242622863.0/115200000.0*t2+t697-961071937.0/103680000.0*t3+4809203.0/921600.0*t24-t330-t700+2401.0/7680.0*t76-t702+t333+3087.0/12800.0*t82-26411.0/19200.0*t84+t336;
218  double t719 = 212317.0/40000.0*t4+1083733.0/960000.0*t7+18179.0/19200.0*t10+19361321.0/1728000.0*t12+t225+4459.0/12000.0*t16-10357963.0/2560000.0*t18-446929.0/480000.0*t20-t229+199283.0/20480.0*t25-144893.0/384000.0*t30-t232-t445-38554229.0/11520000.0*eta+153713.0/2560000.0*xi-497791.0/1200000.0*t1;
219  double t723 = 12285917.0/276480.0*t56;
220  double t726 = 823543.0/30720.0*t74;
221  double t727 = 16807.0/3840.0*t78;
222  double t730 = -285621.0/12800000.0+t450-91679.0/480000.0*t6+144893.0/384000.0*t9+t359+8088871.0/12800000.0*t2-t723+724916437.0/34560000.0*t3-650671.0/102400.0*t24+t364+t726-t365-t727-t367-18179.0/19200.0*t82-31213.0/12800.0*t84-t370;
223  double t743 = -211043.0/2560000.0-3850861.0/192000.0*t4-287581.0/32000.0*t7-4459.0/1920.0*t10-19824371.0/1036800.0*t12-t197-7889.0/1200.0*t16+50965439.0/4608000.0*t18+1399097.0/384000.0*t20+t201-545027.0/36864.0*t25+12201.0/12800.0*t30+t204+t417+42104279.0/6912000.0*eta+272979.0/512000.0*xi;
224  double t748 = 11378339.0/165888.0*t56;
225  double t751 = 823543.0/18432.0*t74;
226  double t753 = 16807.0/1024.0*t78;
227  double t756 = 854217.0/640000.0*t1-t422-49.0/125.0*t6-12201.0/12800.0*t9-t393+134554637.0/23040000.0*t2+t748-622119337.0/20736000.0*t3-472997.0/184320.0*t24-t330-t751+2401.0/256.0*t76+t753+t333+4459.0/1920.0*t82+40817.0/5120.0*t84+t402;
228  double t769 = 45521.0/320000.0+46991.0/2000.0*t4+3889473.0/256000.0*t7+31899.0/5120.0*t10+23065721.0/1036800.0*t12+t164+62083.0/6400.0*t16-3445729.0/288000.0*t18-906549.0/128000.0*t20-t168+362551.0/36864.0*t25-309533.0/102400.0*t30-t172-t385-45345629.0/6912000.0*eta-67963.0/64000.0*xi;
229  double t779 = -1834217.0/640000.0*t1+t390+1519.0/32000.0*t6+309533.0/102400.0*t9+t393-52741003.0/5760000.0*t2-t748+631843387.0/20736000.0*t3+2249737.0/184320.0*t24-t330+t751-31213.0/2048.0*t76-t753+t333-31899.0/5120.0*t82-69629.0/10240.0*t84-t402;
230  double t793 = -1934177.0/240000.0*t4-1021993.0/160000.0*t7-6517.0/1600.0*t10-23683121.0/1728000.0*t12-t129-22981.0/8000.0*t16+5465313.0/2560000.0*t18+1767479.0/480000.0*t20+t133+64827.0/20480.0*t25+83153.0/64000.0*t30+t137+t351+42876029.0/11520000.0*eta+1137437.0/2560000.0*xi+1077069.0/800000.0*t1;
231  double t802 = -1368129.0/12800000.0-t356+60809.0/80000.0*t6-83153.0/64000.0*t9-t359+59592379.0/12800000.0*t2+t723-737881837.0/34560000.0*t3-1054039.0/102400.0*t24+t364-t726+7203.0/1280.0*t76+t727-t367+6517.0/1600.0*t82+2401.0/1600.0*t84+t370;
232  double t816 = -3142223.0/720000.0*t4-633619.0/160000.0*t7-136171.0/57600.0*t10+16428671.0/5184000.0*t12-t94-556003.0/288000.0*t16+93337111.0/23040000.0*t18+6266267.0/2880000.0*t20+t98-1370971.0/184320.0*t25+324919.0/384000.0*t30+t102+t317-30836729.0/34560000.0*eta+1501213.0/7680000.0*xi+7849849.0/9600000.0*t1;
233  double t825 = 517293.0/12800000.0-t322+700553.0/1920000.0*t6-324919.0/384000.0*t9+t325-11258387.0/115200000.0*t2-t697+928658437.0/103680000.0*t3+2994047.0/921600.0*t24-t330+t700+55223.0/15360.0*t76+t702+t333+136171.0/57600.0*t82+45619.0/76800.0*t84-t336;
234  double t837 = 15625.0/2048.0*t30;
235  double t839 = 78125.0/12288.0*t37;
236  double t840 = 71875.0/3072.0*t4+71875.0/3072.0*t7+78125.0/4096.0*t10+15625.0/1024.0*t12+15625.0/1024.0*t14+15625.0/1024.0*t16-209375.0/12288.0*t18-209375.0/12288.0*t20-15625.0/1024.0*t22+78125.0/4096.0*t25-t837-15625.0/1024.0*t32-t839;
237  double t843 = 15625.0/2048.0*t1;
238  double t844 = 78125.0/12288.0*t27;
239  double t846 = 15625.0/2048.0*t9;
240  double t847 = 15625.0/2048.0*t2;
241  double t848 = 78125.0/12288.0*t56;
242  double t850 = 15625.0/2048.0*t24;
243  double t852 = 78125.0/12288.0*t78;
244  double t854 = 15625.0/2048.0*t84;
245  double t855 = -3125.0/1536.0*eta-3125.0/1536.0*xi-t843+t844-53125.0/12288.0*t6+t846-t847+t848-53125.0/12288.0*t3+t850-78125.0/4096.0*t76-t852-78125.0/4096.0*t82-t854;
246  double t857 = 34375.0/1536.0*t4;
247  double t858 = 34375.0/1536.0*t7;
248  double t860 = 15625.0/768.0*t12;
249  double t861 = 15625.0/768.0*t14;
250  double t862 = 15625.0/768.0*t16;
251  double t865 = 15625.0/768.0*t22;
252  double t866 = 78125.0/6144.0*t25;
253  double t868 = 15625.0/768.0*t32;
254  double t869 = -t857-t858-78125.0/3072.0*t10-t860-t861-t862+59375.0/6144.0*t18+34375.0/1536.0*t20+t865-t866+78125.0/6144.0*t30+t868;
255  double t870 = 78125.0/6144.0*t37;
256  double t871 = 3125.0/1024.0*eta;
257  double t872 = 3125.0/1024.0*xi;
258  double t874 = 78125.0/6144.0*t27;
259  double t878 = 78125.0/6144.0*t76;
260  double t880 = t870+t871+t872+78125.0/6144.0*t1-t874+59375.0/6144.0*t6-78125.0/6144.0*t9+t847-3125.0/1024.0*t3-t850+t878+78125.0/3072.0*t82+t854;
261  double t882 = 78125.0/6144.0*t10;
262  double t886 = -t857-t858-t882-t860-t861-t862+34375.0/1536.0*t18+59375.0/6144.0*t20+t865-78125.0/3072.0*t25+t837+t868;
263  double t889 = 78125.0/6144.0*t56;
264  double t893 = 78125.0/6144.0*t78;
265  double t894 = 78125.0/6144.0*t82;
266  double t896 = t871+t872+t843-3125.0/1024.0*t6-t846+78125.0/6144.0*t2-t889+59375.0/6144.0*t3-78125.0/6144.0*t24+78125.0/3072.0*t76+t893+t894+78125.0/6144.0*t84;
267  double t898 = 3125.0/2048.0*xi;
268  double t899 = 3125.0/2048.0*eta;
269  double t900 = 15625.0/1536.0*t22;
270  double t903 = 15625.0/1536.0*t16;
271  double t904 = 15625.0/1536.0*t14;
272  double t905 = 15625.0/1536.0*t12;
273  double t906 = 34375.0/3072.0*t7;
274  double t907 = 34375.0/3072.0*t4;
275  double t908 = 15625.0/1536.0*t30;
276  double t909 = 15625.0/1536.0*t1;
277  double t912 = 15625.0/1536.0*t9;
278  double t913 = 15625.0/1536.0*t32;
279  double t914 = -t898-t899-t900-34375.0/3072.0*t20+3125.0/2048.0*t18+t903+t904+t905+t882+t906+t907-t908-t909-34375.0/3072.0*t6+3125.0/2048.0*t3+t912+t874-t913-t870-t894;
280  double t922 = -3125.0/512.0*xi-3125.0/512.0*eta-t865-3125.0/288.0*t20-3125.0/288.0*t18+t862+t861+t860+78125.0/4608.0*t10+3125.0/288.0*t7+3125.0/288.0*t4;
281  double t924 = 15625.0/1536.0*t2;
282  double t927 = 15625.0/1536.0*t24;
283  double t928 = 15625.0/1536.0*t84;
284  double t931 = 78125.0/4608.0*t25-t908-t909-t924+3125.0/512.0*t6+3125.0/512.0*t3+t912+t927-t868-t928-78125.0/4608.0*t82-78125.0/4608.0*t76;
285  double t937 = -t898-t899-t900+3125.0/2048.0*t20-34375.0/3072.0*t18+t903+t904+t905+t906+t907+t866-t924+3125.0/2048.0*t6-34375.0/3072.0*t3+t927+t889-t913-t928-t878-t893;
286  double t940 = 8125.0/2048.0*t14;
287  double t947 = -3125.0/12288.0*xi+3125.0/12288.0*t18-t940-40625.0/6144.0*t7+15625.0/4096.0*t30+8125.0/2048.0*t1+625.0/4096.0*t2+40625.0/6144.0*t6-15625.0/4096.0*t9-t844+t839-625.0/4096.0;
288  double t950 = 15625.0/2304.0*t22;
289  double t953 = 15625.0/2304.0*t16;
290  double t954 = 6875.0/1536.0*t14;
291  double t955 = 15625.0/2304.0*t12;
292  double t965 = 15625.0/2304.0*t32;
293  double t967 = 3125.0/768.0*xi+3125.0/3072.0*eta+t950+34375.0/4608.0*t20-3125.0/768.0*t18-t953-t954-t955-78125.0/9216.0*t10+3125.0/768.0*t7-34375.0/4608.0*t4+15625.0/3072.0*t30+6875.0/1536.0*t1-625.0/1024.0*t2-3125.0/768.0*t6-3125.0/3072.0*t3-15625.0/3072.0*t9+t965+78125.0/9216.0*t82+625.0/1024.0;
294  double t982 = 3125.0/3072.0*xi+3125.0/768.0*eta+t950-3125.0/768.0*t20+34375.0/4608.0*t18-t953-t954-t955-34375.0/4608.0*t7+3125.0/768.0*t4-78125.0/9216.0*t25-625.0/1024.0*t1+6875.0/1536.0*t2-3125.0/3072.0*t6-3125.0/768.0*t3-15625.0/3072.0*t24+t965+15625.0/3072.0*t84+78125.0/9216.0*t76+625.0/1024.0;
295  double t991 = -3125.0/12288.0*eta+3125.0/12288.0*t20-t940-40625.0/6144.0*t4+625.0/4096.0*t1+8125.0/2048.0*t2+40625.0/6144.0*t3-15625.0/4096.0*t24-t848+15625.0/4096.0*t84+t852-625.0/4096.0;
296 
297  values[0] = t43+t88;
298  values[1] = t106+t123;
299  values[2] = t141+t158;
300  values[3] = t174+t191;
301  values[4] = t205+t219;
302  values[5] = t234+t247;
303  values[6] = t262+t275;
304  values[7] = t292+t303;
305  values[8] = t321+t337;
306  values[9] = t355+t371;
307  values[10] = t388+t403;
308  values[11] = t420+t431;
309  values[12] = t449+t459;
310  values[13] = t477+t487;
311  values[14] = t505+t515;
312  values[15] = t528+t540;
313  values[16] = t553+t565;
314  values[17] = t577+t590;
315  values[18] = t601+t613;
316  values[19] = t625+t635;
317  values[20] = t647+t658;
318  values[21] = t670+t679;
319  values[22] = t693+t705;
320  values[23] = t719+t730;
321  values[24] = t743+t756;
322  values[25] = t769+t779;
323  values[26] = t793+t802;
324  values[27] = t816+t825;
325  values[28] = t840+t855;
326  values[29] = t869+t880;
327  values[30] = t886+t896;
328  values[31] = t914;
329  values[32] = t922+t931;
330  values[33] = t937;
331  values[34] = t947;
332  values[35] = t967;
333  values[36] = t982;
334  values[37] = t991;
335 }
336 
337 // values of the derivatives in xi direction
338 static void C_Q_M7_2D_DeriveXi(double xi, double eta, double *values)
339 {
340  double t1 = eta*eta;
341  double t4 = t1*t1;
342  double t5 = t4*eta;
343  double t6 = 184877.0/207360.0*t5;
344  double t7 = t1*eta;
345  double t9 = xi*xi;
346  double t11 = t9*t9;
347  double t12 = 12113881.0/829440.0*t11;
348  double t13 = t9*xi;
349  double t15 = t11*xi;
350  double t16 = 117649.0/30720.0*t15;
351  double t17 = t11*t9;
352  double t18 = 823543.0/92160.0*t17;
353  double t20 = t4*t1;
354  double t21 = 117649.0/184320.0*t20;
355  double t27 = 1.0+eta;
356  double t28 = (1.0-t9)*(315.0/8.0*t11-105.0/4.0*t9+15.0/8.0)*t27;
357  double t29 = 16807.0/207360.0*t28;
358  double t35 = (1.0-t1)*(63.0/8.0*t5-35.0/4.0*t7+15.0/8.0*eta);
359  double t36 = 16807.0/207360.0*t35;
360  double t37 = xi*eta;
361  double t39 = 1478113.0/69120000.0-45493571.0/23040000.0*t1-802049.0/864000.0*eta+t6+971869.0/5184000.0*t7-172582507.0/34560000.0*t9+t12-922751.0/1152000.0*t13+t16-t18+2529463.0/921600.0*t4-t21-t29-t36-45493571.0/11520000.0*t37;
362  double t40 = t9*eta;
363  double t42 = xi*t1;
364  double t44 = t13*eta;
365  double t46 = t9*t1;
366  double t48 = xi*t7;
367  double t50 = t11*eta;
368  double t51 = 184877.0/41472.0*t50;
369  double t52 = t13*t1;
370  double t54 = t9*t7;
371  double t56 = xi*t4;
372  double t58 = t15*eta;
373  double t59 = 117649.0/30720.0*t58;
374  double t60 = t11*t1;
375  double t61 = 85829.0/15360.0*t60;
376  double t62 = t13*t7;
377  double t64 = t9*t4;
378  double t66 = xi*t5;
379  double t67 = 85829.0/38400.0*t66;
380  double t74 = xi*(63.0/8.0*t15-35.0/4.0*t13+15.0/8.0*xi)*t27;
381  double t75 = 16807.0/103680.0*t74;
382  double t76 = 971869.0/1728000.0*t40+225377.0/150000.0*t42+2529463.0/230400.0*t44+2902329.0/320000.0*t46+967443.0/160000.0*t48+t51-98177.0/96000.0*t52-94061.0/32000.0*t54-98177.0/192000.0*t56-t59-t61-87887.0/9600.0*t62-87887.0/12800.0*t64-t67-49159793.0/57600000.0*xi+t75;
383  double t79 = 1526693.0/1728000.0*eta;
384  double t82 = 14101073.0/165888.0*t11;
385  double t84 = 117649.0/6144.0*t15;
386  double t85 = 5764801.0/92160.0*t17;
387  double t87 = 117649.0/207360.0*t28;
388  double t90 = 141169.0/50000.0*t42;
389  double t92 = 6266267.0/2880000.0*t1+t79-556003.0/288000.0*t7+928658437.0/34560000.0*t9-t82+2994047.0/230400.0*t13-t84+t85-136171.0/57600.0*t4+t87+93337111.0/11520000.0*t37+16428671.0/1728000.0*t40-t90-1370971.0/46080.0*t44;
390  double t95 = 420175.0/20736.0*t50;
391  double t97 = 556003.0/96000.0*t54;
392  double t99 = 117649.0/6144.0*t58;
393  double t100 = 16807.0/1536.0*t60;
394  double t103 = 14651.0/12800.0*t66;
395  double t105 = 117649.0/103680.0*t74;
396  double t106 = -30836729.0/34560000.0-3142223.0/240000.0*t46-633619.0/80000.0*t48-t95+45619.0/19200.0*t52+t97+324919.0/192000.0*t56+t99+t100+55223.0/3840.0*t62+136171.0/19200.0*t64+t103-11258387.0/57600000.0*xi-t105;
397  double t109 = 3412507.0/576000.0*eta;
398  double t112 = 12285917.0/55296.0*t11;
399  double t114 = 352947.0/10240.0*t15;
400  double t115 = 5764801.0/30720.0*t17;
401  double t117 = 117649.0/69120.0*t28;
402  double t120 = 2277569.0/400000.0*t42;
403  double t122 = 1767479.0/480000.0*t1+t109-22981.0/8000.0*t7-737881837.0/11520000.0*t9+t112-1054039.0/25600.0*t13+t114-t115-6517.0/1600.0*t4-t117+5465313.0/1280000.0*t37-23683121.0/576000.0*t40-t120+64827.0/5120.0*t44;
404  double t125 = 16807.0/432.0*t50;
405  double t127 = 68943.0/8000.0*t54;
406  double t129 = 352947.0/10240.0*t58;
407  double t130 = 16807.0/768.0*t60;
408  double t133 = 14651.0/6400.0*t66;
409  double t135 = 117649.0/34560.0*t74;
410  double t136 = 42876029.0/11520000.0-1934177.0/80000.0*t46-1021993.0/80000.0*t48+t125+2401.0/400.0*t52+t127+83153.0/32000.0*t56-t129+t130+7203.0/320.0*t62+19551.0/1600.0*t64+t133+59592379.0/6400000.0*xi+t135;
411  double t139 = 139601.0/10800.0*eta;
412  double t142 = 56891695.0/165888.0*t11;
413  double t144 = 5764801.0/18432.0*t17;
414  double t146 = 117649.0/41472.0*t28;
415  double t149 = 12372059.0/640000.0*t42;
416  double t150 = -45345629.0/6912000.0-906549.0/128000.0*t1-t139+62083.0/6400.0*t7+631843387.0/6912000.0*t9-t142+2249737.0/46080.0*t13-t84+t144+31899.0/5120.0*t4+t146-3445729.0/144000.0*t37+23065721.0/345600.0*t40+t149;
417  double t154 = 1932805.0/41472.0*t50;
418  double t156 = 186249.0/6400.0*t54;
419  double t158 = 84035.0/1024.0*t60;
420  double t161 = 14651.0/10240.0*t66;
421  double t163 = 117649.0/20736.0*t74;
422  double t164 = 362551.0/9216.0*t44+140973.0/2000.0*t46+3889473.0/128000.0*t48-t154-69629.0/2560.0*t52-t156-309533.0/51200.0*t56+t99-t158-31213.0/512.0*t62-95697.0/5120.0*t64-t161-52741003.0/2880000.0*xi-t163;
423  double t167 = 1693391.0/172800.0*eta;
424  double t174 = 2978171.0/160000.0*t42;
425  double t175 = 42104279.0/6912000.0+1399097.0/384000.0*t1+t167-7889.0/1200.0*t7-622119337.0/6912000.0*t9+t142-472997.0/46080.0*t13-t84-t144-4459.0/1920.0*t4-t146+50965439.0/2304000.0*t37-19824371.0/345600.0*t40-t174;
426  double t180 = 7889.0/400.0*t54;
427  double t184 = 1127.0/1280.0*t66;
428  double t186 = -545027.0/9216.0*t44-3850861.0/64000.0*t46-287581.0/16000.0*t48+t154+40817.0/1280.0*t52+t180+12201.0/6400.0*t56+t99+t158+2401.0/64.0*t62+4459.0/640.0*t64+t184+134554637.0/11520000.0*xi+t163;
429  double t189 = 1971907.0/576000.0*eta;
430  double t196 = 13696039.0/2400000.0*t42;
431  double t198 = -446929.0/480000.0*t1-t189+4459.0/12000.0*t7+724916437.0/11520000.0*t9-t112-650671.0/25600.0*t13+t114+t115+18179.0/19200.0*t4+t117-10357963.0/1280000.0*t37+19361321.0/576000.0*t40+t196+199283.0/5120.0*t44;
432  double t202 = 4459.0/4000.0*t54;
433  double t204 = 2401.0/640.0*t62;
434  double t206 = 14651.0/38400.0*t66;
435  double t208 = -38554229.0/11520000.0+636951.0/40000.0*t46+1083733.0/480000.0*t48-t125-31213.0/3200.0*t52-t202-144893.0/192000.0*t56-t129-t130-t204-18179.0/6400.0*t64-t206+8088871.0/6400000.0*xi-t135;
436  double t211 = 2074807.0/1728000.0*eta;
437  double t218 = 3234931.0/1200000.0*t42;
438  double t220 = 37387.0/320000.0*t1+t211-14749.0/96000.0*t7-961071937.0/34560000.0*t9+t82+4809203.0/230400.0*t13-t84-t85-3087.0/12800.0*t4-t87+22266139.0/11520000.0*t37-27233171.0/1728000.0*t40+t218-765919.0/46080.0*t44;
439  double t224 = 14749.0/32000.0*t54;
440  double t228 = 637.0/6400.0*t66;
441  double t230 = 41641229.0/34560000.0+994357.0/160000.0*t46-169393.0/240000.0*t48+t95-26411.0/4800.0*t52+t224+1519.0/8000.0*t56+t99-t100+2401.0/1920.0*t62+9261.0/12800.0*t64+t228-242622863.0/57600000.0*xi+t105;
442  double t233 = 273563.0/432000.0*eta;
443  double t234 = 7123181.0/5184000.0*t7;
444  double t239 = 7123181.0/1728000.0*t40;
445  double t240 = 5917771.0/23040000.0*t1-t233-t6+t234+196867657.0/34560000.0*t9-t12-5668649.0/1152000.0*t13+t16+t18-243521.0/307200.0*t4+t21+t29+t36-3953677.0/5760000.0*t37+t239;
446  double t241 = 20141003.0/9600000.0*t42;
447  double t246 = 83771.0/48000.0*t54;
448  double t250 = 85829.0/153600.0*t66;
449  double t252 = -17668213.0/69120000.0-t241+947497.0/230400.0*t44-940003.0/240000.0*t46+2126359.0/1920000.0*t48-t51+594629.0/192000.0*t52-t246+246511.0/256000.0*t56-t59+t61-87887.0/38400.0*t62+77597.0/76800.0*t64-t250+36656459.0/28800000.0*xi-t75;
450  double t255 = 326879.0/216000.0*eta;
451  double t256 = 84035.0/20736.0*t5;
452  double t259 = 14651.0/5120.0*t11;
453  double t262 = 117649.0/36864.0*t20;
454  double t263 = 117649.0/207360.0*t35;
455  double t265 = 419979.0/2560000.0-10863937.0/7680000.0*t1+t255+t256-14426923.0/2592000.0*t7-2781681.0/1280000.0*t9+t259+41797.0/48000.0*t13+818741.0/184320.0*t4-t262-t263+6379457.0/5760000.0*t37;
456  double t267 = 4969237.0/2400000.0*t42;
457  double t272 = 89523.0/64000.0*t54;
458  double t274 = 14651.0/5120.0*t60;
459  double t277 = 16807.0/15360.0*t66;
460  double t279 = -89523.0/64000.0*t40+t267-136171.0/57600.0*t44+7623959.0/1920000.0*t46-6341041.0/2880000.0*t48-41797.0/48000.0*t52+t272-2401.0/1600.0*t56-t274+136171.0/57600.0*t62-55223.0/30720.0*t64+t277-1367737.0/2400000.0*xi;
461  double t282 = 1251607.0/576000.0*eta;
462  double t283 = 16807.0/2160.0*t5;
463  double t286 = 14651.0/2560.0*t11;
464  double t289 = 117649.0/20480.0*t20;
465  double t290 = 117649.0/69120.0*t35;
466  double t293 = 5099087.0/2560000.0*t1-t282-t283+17200421.0/1728000.0*t7-1269933.0/320000.0*t9+t286+34153.0/16000.0*t13-160867.0/20480.0*t4+t289+t290+1582259.0/960000.0*t37-42189.0/16000.0*t40;
467  double t294 = 2564513.0/800000.0*t42;
468  double t299 = 42189.0/16000.0*t54;
469  double t301 = 14651.0/2560.0*t60;
470  double t303 = 7203.0/2560.0*t64;
471  double t304 = 16807.0/7680.0*t66;
472  double t306 = 303163.0/2560000.0+t294-6517.0/1600.0*t44+542577.0/80000.0*t46-1841567.0/480000.0*t48-34153.0/16000.0*t52+t299-26411.0/12800.0*t56-t301+6517.0/1600.0*t62-t303+t304-1827651.0/1600000.0*xi;
473  double t309 = 1911539.0/691200.0*eta;
474  double t310 = 386561.0/41472.0*t5;
475  double t313 = 14651.0/4096.0*t11;
476  double t316 = 117649.0/41472.0*t35;
477  double t319 = 1815107.0/2304000.0*t1+t309+t310-25062667.0/2073600.0*t7+2268063.0/512000.0*t9-t313+11221.0/12800.0*t13+26411.0/9216.0*t4-t262-t316-1510229.0/256000.0*t37+35329.0/25600.0*t40;
478  double t320 = 1454467.0/320000.0*t42;
479  double t325 = 35329.0/25600.0*t54;
480  double t327 = 14651.0/4096.0*t60;
481  double t330 = 16807.0/2048.0*t66;
482  double t332 = -236621.0/512000.0-t320+31899.0/5120.0*t44-771211.0/64000.0*t46+112847.0/8000.0*t48-11221.0/12800.0*t52-t325+21609.0/5120.0*t56+t327-31899.0/5120.0*t62+31213.0/4096.0*t64-t330+207809.0/640000.0*xi;
483  double t335 = 612941.0/172800.0*eta;
484  double t338 = 1127.0/512.0*t11;
485  double t342 = 207221.0/512000.0-19272239.0/4608000.0*t1-t335-t310+13341671.0/1036800.0*t7-104811.0/64000.0*t9+t338+147.0/400.0*t13+256907.0/36864.0*t4-t262+t316+3210137.0/768000.0*t37;
486  double t344 = 367059.0/640000.0*t42;
487  double t349 = 3087.0/3200.0*t54;
488  double t351 = 1127.0/512.0*t60;
489  double t355 = 3087.0/3200.0*t40+t344-4459.0/1920.0*t44+50617.0/8000.0*t46-4756381.0/384000.0*t48-147.0/400.0*t52-t349-7203.0/10240.0*t56-t351+4459.0/1920.0*t62-2401.0/512.0*t64+t330+166257.0/1280000.0*xi;
490  double t358 = 1611757.0/576000.0*eta;
491  double t361 = 14651.0/15360.0*t11;
492  double t365 = -393813.0/2560000.0+25791689.0/7680000.0*t1+t358+t283-18280871.0/1728000.0*t7+546791.0/640000.0*t9-t361-9359.0/48000.0*t13-549829.0/61440.0*t4+t289-t290-261709.0/960000.0*t37;
493  double t367 = 335209.0/600000.0*t42;
494  double t372 = 24353.0/32000.0*t54;
495  double t374 = 14651.0/15360.0*t60;
496  double t378 = 24353.0/32000.0*t40+t367+18179.0/19200.0*t44-211729.0/160000.0*t46+98441.0/40000.0*t48+9359.0/48000.0*t52-t372-2401.0/3200.0*t56+t374-18179.0/19200.0*t62+2401.0/5120.0*t64-t304+229957.0/1200000.0*xi;
497  double t381 = 1714657.0/1728000.0*eta;
498  double t384 = 637.0/2560.0*t11;
499  double t389 = -5661313.0/7680000.0*t1-t381-t256+26152721.0/5184000.0*t7-137347.0/640000.0*t9+t384+3479.0/64000.0*t13+717899.0/184320.0*t4-t262+t263+74431.0/1920000.0*t37-7889.0/48000.0*t40;
500  double t390 = 371273.0/400000.0*t42;
501  double t395 = 7889.0/48000.0*t54;
502  double t397 = 637.0/2560.0*t60;
503  double t401 = 259063.0/7680000.0+t390-3087.0/12800.0*t44+356083.0/960000.0*t46+506611.0/480000.0*t48-3479.0/64000.0*t52+t395-21609.0/25600.0*t56-t397+3087.0/12800.0*t62-2401.0/15360.0*t64-t277-269059.0/3200000.0*xi;
504  double t407 = 5077697.0/1658880.0*t11;
505  double t413 = -496649.0/5760000.0*t1+839329.0/3456000.0*eta+t6-10198837.0/10368000.0*t7-155978489.0/69120000.0*t9+t407-3621101.0/1152000.0*t13+t16-13827.0/25600.0*t4+t21-t29-t36-496649.0/2880000.0*t37-10198837.0/3456000.0*t40-1343687.0/1600000.0*t42;
506  double t420 = 85829.0/61440.0*t60;
507  double t424 = 9594313.0/69120000.0-13827.0/6400.0*t44-716051.0/1920000.0*t46-716051.0/2880000.0*t48+t51+84457.0/64000.0*t52+73481.0/128000.0*t54+84457.0/128000.0*t56+t59+t420+77597.0/230400.0*t62+77597.0/307200.0*t64+t250+18298957.0/57600000.0*xi+t75;
508  double t429 = 2907611.0/165888.0*t11;
509  double t435 = 74431.0/3840000.0*t1-t381-7889.0/144000.0*t7+467769337.0/34560000.0*t9-t429+3930437.0/230400.0*t13-t84-3087.0/51200.0*t4+t87-5661313.0/3840000.0*t37+26152721.0/1728000.0*t40+t390+717899.0/46080.0*t44;
510  double t440 = 16807.0/6144.0*t60;
511  double t444 = -17386327.0/17280000.0+506611.0/320000.0*t46+356083.0/1440000.0*t48-t95-21609.0/12800.0*t52+t395-3479.0/128000.0*t56-t99-t440-2401.0/11520.0*t62+9261.0/51200.0*t64-t228-128307137.0/57600000.0*xi-t105;
512  double t449 = 1226911.0/27648.0*t11;
513  double t454 = 4250113.0/1440000.0-261709.0/1920000.0*t1+t358+24353.0/96000.0*t7-204365231.0/5760000.0*t9+t449-861959.0/25600.0*t13+t114+18179.0/76800.0*t4-t117+25791689.0/3840000.0*t37-18280871.0/576000.0*t40+t367;
514  double t460 = 16807.0/3072.0*t60;
515  double t464 = -549829.0/15360.0*t44+295323.0/80000.0*t46-211729.0/240000.0*t48+t125-2401.0/1600.0*t52-t372+9359.0/96000.0*t56+t129-t460+2401.0/3840.0*t62-18179.0/25600.0*t64+t206+35582379.0/6400000.0*xi+t135;
516  double t469 = 22269275.0/331776.0*t11;
517  double t474 = -65455033.0/13824000.0+3210137.0/1536000.0*t1-t335+1029.0/3200.0*t7+376548487.0/6912000.0*t9-t469+530621.0/23040.0*t13-t84-4459.0/7680.0*t4+t146-19272239.0/2304000.0*t37+13341671.0/345600.0*t40+t344;
518  double t480 = 84035.0/4096.0*t60;
519  double t484 = 256907.0/9216.0*t44-4756381.0/256000.0*t46+50617.0/12000.0*t48-t154-7203.0/5120.0*t52-t349-147.0/800.0*t56-t99+t480-2401.0/384.0*t62+4459.0/2560.0*t64-t184-62402137.0/11520000.0*xi-t163;
520  double t494 = -1510229.0/512000.0*t1+t309+35329.0/76800.0*t7-748234949.0/13824000.0*t9+t469+607453.0/46080.0*t13-t84+31899.0/20480.0*t4-t146+1815107.0/1152000.0*t37-25062667.0/691200.0*t40-t320+26411.0/2304.0*t44;
521  double t502 = 31917179.0/6912000.0+338541.0/16000.0*t46-771211.0/96000.0*t48+t154+21609.0/2560.0*t52-t325-11221.0/25600.0*t56-t99-t480+31213.0/3072.0*t62-95697.0/20480.0*t64+t161-501613.0/11520000.0*xi+t163;
522  double t512 = 1582259.0/1920000.0*t1-t282-14063.0/16000.0*t7+50686139.0/1440000.0*t9-t449-794731.0/25600.0*t13+t114-6517.0/6400.0*t4+t117+5099087.0/1280000.0*t37+17200421.0/576000.0*t40+t294-160867.0/5120.0*t44;
523  double t519 = -16460227.0/5760000.0-1841567.0/320000.0*t46+180859.0/40000.0*t48-t125-26411.0/6400.0*t52+t299-34153.0/32000.0*t56+t129+t460-t204+19551.0/6400.0*t64-t133+26096371.0/6400000.0*xi-t135;
524  double t528 = 37473779.0/34560000.0+6379457.0/11520000.0*t1+t255-29841.0/64000.0*t7-59484089.0/4320000.0*t9+t429+4232963.0/230400.0*t13-t84-136171.0/230400.0*t4-t87-10863937.0/3840000.0*t37-14426923.0/864000.0*t40+t267;
525  double t537 = 818741.0/46080.0*t44-6341041.0/1920000.0*t46+7623959.0/2880000.0*t48+t95-2401.0/800.0*t52+t272-41797.0/96000.0*t56-t99+t440-55223.0/23040.0*t62+136171.0/76800.0*t64-t103-170592863.0/57600000.0*xi+t105;
526  double t544 = -3953677.0/11520000.0*t1-t233-t6+t234+21015133.0/8640000.0*t9-t407-2159737.0/576000.0*t13+t16+947497.0/921600.0*t4-t21+t29+t36+5917771.0/11520000.0*t37+t239-t241;
527  double t553 = -6820919.0/34560000.0-243521.0/76800.0*t44+2126359.0/1280000.0*t46-940003.0/360000.0*t48-t51+246511.0/128000.0*t52-t246+594629.0/384000.0*t56+t59-t420+77597.0/57600.0*t62-87887.0/51200.0*t64+t67+39583543.0/57600000.0*xi-t75;
528  double t562 = 22266139.0/23040000.0*t1+t211+t256-27233171.0/5184000.0*t7+38661.0/320000.0*t9-t384-1519.0/4000.0*t13-765919.0/184320.0*t4+t262-t263+37387.0/160000.0*t37-14749.0/32000.0*t40;
529  double t570 = 16807.0/3840.0*t66;
530  double t572 = -6321.0/2560000.0+t218-3087.0/3200.0*t44-169393.0/160000.0*t46+994357.0/240000.0*t48+1519.0/4000.0*t52+t224-26411.0/9600.0*t56+t397+3087.0/3200.0*t62+2401.0/2560.0*t64-t570+5537.0/100000.0*xi;
531  double t581 = -10357963.0/2560000.0*t1-t189-t283+19361321.0/1728000.0*t7-91679.0/160000.0*t9+t361+144893.0/96000.0*t13+199283.0/20480.0*t4-t289+t290-446929.0/240000.0*t37+4459.0/4000.0*t40;
532  double t588 = 16807.0/1920.0*t66;
533  double t590 = 153713.0/2560000.0+t196+18179.0/4800.0*t44+1083733.0/320000.0*t46+212317.0/20000.0*t48-144893.0/96000.0*t52-t202-31213.0/6400.0*t56-t374-18179.0/4800.0*t62-t303-t588-497791.0/600000.0*xi;
534  double t598 = 272979.0/512000.0+50965439.0/4608000.0*t1+t167+t310-19824371.0/1036800.0*t7-147.0/125.0*t9-t338-12201.0/3200.0*t13-545027.0/36864.0*t4+t262-t316+1399097.0/192000.0*t37;
535  double t607 = 16807.0/512.0*t66;
536  double t609 = -7889.0/400.0*t40-t174-4459.0/480.0*t44-862743.0/32000.0*t46-3850861.0/96000.0*t48+12201.0/3200.0*t52+t180+40817.0/2560.0*t56+t351+4459.0/480.0*t62+7203.0/256.0*t64+t607+854217.0/320000.0*xi;
537  double t617 = -67963.0/64000.0-3445729.0/288000.0*t1-t139-t310+23065721.0/1036800.0*t7+4557.0/32000.0*t9+t313+309533.0/25600.0*t13+362551.0/36864.0*t4+t262+t316-906549.0/64000.0*t37;
538  double t627 = 186249.0/6400.0*t40+t149+31899.0/1280.0*t44+11668419.0/256000.0*t46+46991.0/1000.0*t48-309533.0/25600.0*t52-t156-69629.0/5120.0*t56-t327-31899.0/1280.0*t62-93639.0/2048.0*t64-t607-1834217.0/320000.0*xi;
539  double t636 = 5465313.0/2560000.0*t1+t109+t283-23683121.0/1728000.0*t7+182427.0/80000.0*t9-t286-83153.0/16000.0*t13+64827.0/20480.0*t4-t289-t290+1767479.0/240000.0*t37-68943.0/8000.0*t40;
540  double t645 = 1137437.0/2560000.0-t120-6517.0/400.0*t44-3065979.0/160000.0*t46-1934177.0/120000.0*t48+83153.0/16000.0*t52+t127+2401.0/800.0*t56+t301+6517.0/400.0*t62+21609.0/1280.0*t64+t588+1077069.0/400000.0*xi;
541  double t654 = 93337111.0/23040000.0*t1+t79-t256+16428671.0/5184000.0*t7+700553.0/640000.0*t9-t259-324919.0/96000.0*t13-1370971.0/184320.0*t4+t262+t263+6266267.0/1440000.0*t37-556003.0/96000.0*t40;
542  double t663 = 1501213.0/7680000.0-t90-136171.0/14400.0*t44-1900857.0/160000.0*t46-3142223.0/360000.0*t48+324919.0/96000.0*t52+t97+45619.0/38400.0*t56+t274+136171.0/14400.0*t62+55223.0/5120.0*t64+t570+7849849.0/4800000.0*xi;
543  double t665 = 15625.0/1024.0*xi;
544  double t667 = 15625.0/512.0*t13;
545  double t668 = 390625.0/12288.0*t11;
546  double t675 = -3125.0/1536.0-t665-53125.0/4096.0*t9+t667+t668-78125.0/1024.0*t62-15625.0/1024.0*eta-209375.0/6144.0*t37-209375.0/12288.0*t1+46875.0/1024.0*t40+15625.0/512.0*t42;
547  double t681 = 15625.0/512.0*t52;
548  double t683 = 15625.0/1024.0*t56;
549  double t684 = 390625.0/12288.0*t60;
550  double t686 = 78125.0/6144.0*t66;
551  double t687 = 15625.0/1024.0*t7+78125.0/1024.0*t44+71875.0/1024.0*t46+71875.0/1536.0*t48+78125.0/4096.0*t4-t681-46875.0/1024.0*t54-t683-t684-234375.0/4096.0*t64-t686;
552  double t692 = 390625.0/6144.0*t11;
553  double t694 = 15625.0/768.0*eta;
554  double t697 = 15625.0/256.0*t40;
555  double t699 = 15625.0/384.0*t42;
556  double t700 = 15625.0/768.0*t7;
557  double t702 = 34375.0/512.0*t46;
558  double t703 = 34375.0/768.0*t48;
559  double t704 = 78125.0/6144.0*t4;
560  double t706 = 15625.0/256.0*t54;
561  double t707 = 390625.0/6144.0*t60;
562  double t708 = 78125.0/2048.0*t64;
563  double t709 = -t699-t700-78125.0/768.0*t44-t702-t703-t704+78125.0/1536.0*t52+t706+t683+t707+t708;
564  double t712 = 78125.0/1536.0*t62;
565  double t715 = 78125.0/1536.0*t44;
566  double t719 = 78125.0/3072.0*t66;
567  double t720 = 3125.0/1024.0+t665-9375.0/1024.0*t9-t667+t712+t694+59375.0/3072.0*t37+34375.0/1536.0*t1-t697-t699-t700-t715-t702-t703-78125.0/3072.0*t4+t681+t706+78125.0/3072.0*t56+78125.0/1024.0*t64+t719;
568  double t721 = 15625.0/1536.0*eta;
569  double t724 = 15625.0/512.0*t40;
570  double t725 = 15625.0/768.0*t42;
571  double t726 = 15625.0/1536.0*t7;
572  double t727 = 34375.0/1024.0*t46;
573  double t728 = 34375.0/1536.0*t48;
574  double t729 = 15625.0/384.0*t52;
575  double t730 = 15625.0/768.0*xi;
576  double t732 = 15625.0/384.0*t13;
577  double t733 = 15625.0/512.0*t54;
578  double t734 = -3125.0/2048.0-t721-34375.0/1536.0*t37+3125.0/2048.0*t1+t724+t725+t726+t715+t727+t728-t729-t730-34375.0/1024.0*t9+t732+t692-t733-t707-t712;
579  double t742 = 15625.0/768.0*t56;
580  double t745 = -3125.0/512.0-t694-3125.0/144.0*t37-3125.0/288.0*t1+t697+t699+t700+78125.0/1152.0*t44+3125.0/96.0*t46+3125.0/144.0*t48+78125.0/4608.0*t4-t729-t730+9375.0/512.0*t9+t732-t706-t742-78125.0/1152.0*t62-78125.0/1536.0*t64;
581  double t749 = -3125.0/2048.0-t721+3125.0/1024.0*t37-34375.0/3072.0*t1+t724+t725+t726+t727+t728+t704+9375.0/2048.0*t9-t733-t742-t708-t719;
582  double t751 = 8125.0/1024.0*t42;
583  double t758 = 15625.0/2304.0*eta;
584  double t761 = 15625.0/768.0*t40;
585  double t762 = 6875.0/768.0*t42;
586  double t763 = 15625.0/2304.0*t7;
587  double t771 = 15625.0/768.0*t54;
588  double t773 = 3125.0/768.0+t758+34375.0/2304.0*t37-3125.0/768.0*t1-t761-t762-t763-78125.0/2304.0*t44+3125.0/256.0*t46-34375.0/2304.0*t48+15625.0/768.0*t52+6875.0/768.0*xi-3125.0/256.0*t9-15625.0/768.0*t13+t771+78125.0/2304.0*t62;
589  double t783 = 3125.0/3072.0+t758-3125.0/384.0*t37+34375.0/4608.0*t1-t761-t762-t763-34375.0/1536.0*t46+3125.0/384.0*t48-78125.0/9216.0*t4-625.0/512.0*xi-3125.0/1024.0*t9+t771+15625.0/1536.0*t56+78125.0/3072.0*t64;
590 
591  values[0] = t39+t76;
592  values[1] = t92+t106;
593  values[2] = t122+t136;
594  values[3] = t150+t164;
595  values[4] = t175+t186;
596  values[5] = t198+t208;
597  values[6] = t220+t230;
598  values[7] = t240+t252;
599  values[8] = t265+t279;
600  values[9] = t293+t306;
601  values[10] = t319+t332;
602  values[11] = t342+t355;
603  values[12] = t365+t378;
604  values[13] = t389+t401;
605  values[14] = t413+t424;
606  values[15] = t435+t444;
607  values[16] = t454+t464;
608  values[17] = t474+t484;
609  values[18] = t494+t502;
610  values[19] = t512+t519;
611  values[20] = t528+t537;
612  values[21] = t544+t553;
613  values[22] = t562+t572;
614  values[23] = t581+t590;
615  values[24] = t598+t609;
616  values[25] = t617+t627;
617  values[26] = t636+t645;
618  values[27] = t654+t663;
619  values[28] = t675+t687;
620  values[29] = 3125.0/1024.0+78125.0/3072.0*xi+59375.0/2048.0*t9-78125.0/1536.0*t13-t692+78125.0/768.0*t62+t694+34375.0/768.0*t37+59375.0/6144.0*t1-t697+t709;
621  values[30] = t720;
622  values[31] = t734;
623  values[32] = t745;
624  values[33] = t749;
625  values[34] = -3125.0/12288.0+3125.0/12288.0*t1-t751-40625.0/2048.0*t46+15625.0/1024.0*t52+8125.0/1024.0*xi+40625.0/2048.0*t9-15625.0/1024.0*t13-t668+t684;
626  values[35] = t773;
627  values[36] = t783;
628  values[37] = 3125.0/6144.0*t37-t751-40625.0/3072.0*t48+625.0/2048.0*xi+15625.0/2048.0*t56+t686;
629 }
630 
631 // values of the derivatives in eta direction
632 static void C_Q_M7_2D_DeriveEta(double xi, double eta, double *values)
633 {
634  double t1 = eta*eta;
635  double t2 = t1*t1;
636  double t3 = t2*eta;
637  double t4 = 117649.0/30720.0*t3;
638  double t5 = xi*xi;
639  double t7 = t5*t5;
640  double t8 = xi*t7;
641  double t10 = xi*t5;
642  double t14 = (1.0-t5)*(63.0/8.0*t8-35.0/4.0*t10+15.0/8.0*xi);
643  double t15 = 16807.0/207360.0*t14;
644  double t16 = xi*t3;
645  double t17 = 117649.0/30720.0*t16;
646  double t18 = t5*t2;
647  double t19 = 85829.0/15360.0*t18;
648  double t20 = t1*eta;
649  double t21 = t10*t20;
650  double t28 = 1.0+xi;
651  double t29 = eta*(63.0/8.0*t3-35.0/4.0*t20+15.0/8.0*eta)*t28;
652  double t30 = 16807.0/103680.0*t29;
653  double t36 = (1.0-t1)*(315.0/8.0*t2-105.0/4.0*t1+15.0/8.0)*t28;
654  double t37 = 16807.0/207360.0*t36;
655  double t38 = t7*t1;
656  double t40 = xi*t2;
657  double t41 = 184877.0/41472.0*t40;
658  double t42 = t8*eta;
659  double t43 = 85829.0/38400.0*t42;
660  double t44 = t5*t20;
661  double t46 = t7*eta;
662  double t48 = t10*t1;
663  double t50 = xi*t20;
664  double t52 = 1478113.0/69120000.0+t4-t15-t17-t19-87887.0/9600.0*t21+t30-t37-87887.0/12800.0*t38+t41-t43-98177.0/96000.0*t44-98177.0/192000.0*t46-94061.0/32000.0*t48+2529463.0/230400.0*t50;
665  double t53 = t5*t1;
666  double t55 = t5*eta;
667  double t57 = xi*t1;
668  double t59 = xi*eta;
669  double t64 = t10*eta;
670  double t68 = 12113881.0/829440.0*t2;
671  double t70 = 184877.0/207360.0*t8;
672  double t71 = t2*t1;
673  double t72 = 823543.0/92160.0*t71;
674  double t73 = t7*t5;
675  double t74 = 117649.0/184320.0*t73;
676  double t76 = 2902329.0/320000.0*t53+225377.0/150000.0*t55+971869.0/1728000.0*t57-45493571.0/11520000.0*t59-802049.0/864000.0*xi-922751.0/1152000.0*t20-49159793.0/57600000.0*eta+967443.0/160000.0*t64-45493571.0/23040000.0*t5+971869.0/5184000.0*t10+t68+2529463.0/921600.0*t7+t70-t72-t74-172582507.0/34560000.0*t1;
677  double t78 = 117649.0/207360.0*t14;
678  double t79 = 14651.0/5120.0*t18;
679  double t82 = 16807.0/3840.0*t42;
680  double t85 = 556003.0/96000.0*t48;
681  double t88 = 141169.0/50000.0*t55;
682  double t90 = t78+t79+136171.0/14400.0*t21+55223.0/5120.0*t38+t82+324919.0/96000.0*t44+45619.0/38400.0*t46+t85-136171.0/14400.0*t50-1900857.0/160000.0*t53-t88-556003.0/96000.0*t57;
683  double t92 = 1526693.0/1728000.0*xi;
684  double t98 = 14651.0/5120.0*t2;
685  double t100 = 84035.0/20736.0*t8;
686  double t101 = 117649.0/36864.0*t73;
687  double t103 = 1501213.0/7680000.0+6266267.0/1440000.0*t59+t92-324919.0/96000.0*t20+7849849.0/4800000.0*eta-3142223.0/360000.0*t64+93337111.0/23040000.0*t5+16428671.0/5184000.0*t10-t98-1370971.0/184320.0*t7-t100+t101+700553.0/640000.0*t1;
688  double t105 = 117649.0/69120.0*t14;
689  double t106 = 14651.0/2560.0*t18;
690  double t109 = 16807.0/1920.0*t42;
691  double t112 = 68943.0/8000.0*t48;
692  double t115 = 2277569.0/400000.0*t55;
693  double t117 = -t105+t106+6517.0/400.0*t21+21609.0/1280.0*t38+t109+83153.0/16000.0*t44+2401.0/800.0*t46+t112-6517.0/400.0*t50-3065979.0/160000.0*t53-t115-68943.0/8000.0*t57;
694  double t119 = 3412507.0/576000.0*xi;
695  double t125 = 14651.0/2560.0*t2;
696  double t127 = 16807.0/2160.0*t8;
697  double t128 = 117649.0/20480.0*t73;
698  double t130 = 1137437.0/2560000.0+1767479.0/240000.0*t59+t119-83153.0/16000.0*t20+1077069.0/400000.0*eta-1934177.0/120000.0*t64+5465313.0/2560000.0*t5-23683121.0/1728000.0*t10-t125+64827.0/20480.0*t7+t127-t128+182427.0/80000.0*t1;
699  double t132 = 117649.0/41472.0*t14;
700  double t133 = 14651.0/4096.0*t18;
701  double t136 = 16807.0/512.0*t42;
702  double t139 = 186249.0/6400.0*t48;
703  double t142 = 12372059.0/640000.0*t55;
704  double t143 = -67963.0/64000.0+t132-t133-31899.0/1280.0*t21-93639.0/2048.0*t38-t136-309533.0/25600.0*t44-69629.0/5120.0*t46-t139+31899.0/1280.0*t50+11668419.0/256000.0*t53+t142;
705  double t146 = 139601.0/10800.0*xi;
706  double t152 = 14651.0/4096.0*t2;
707  double t154 = 386561.0/41472.0*t8;
708  double t156 = 186249.0/6400.0*t57-906549.0/64000.0*t59-t146+309533.0/25600.0*t20-1834217.0/320000.0*eta+46991.0/1000.0*t64-3445729.0/288000.0*t5+23065721.0/1036800.0*t10+t152+362551.0/36864.0*t7-t154+t101+4557.0/32000.0*t1;
709  double t158 = 1127.0/512.0*t18;
710  double t163 = 7889.0/400.0*t48;
711  double t166 = 2978171.0/160000.0*t55;
712  double t167 = 272979.0/512000.0-t132+t158+4459.0/480.0*t21+7203.0/256.0*t38+t136+12201.0/3200.0*t44+40817.0/2560.0*t46+t163-4459.0/480.0*t50-862743.0/32000.0*t53-t166;
713  double t170 = 1693391.0/172800.0*xi;
714  double t176 = 1127.0/512.0*t2;
715  double t179 = -7889.0/400.0*t57+1399097.0/192000.0*t59+t170-12201.0/3200.0*t20+854217.0/320000.0*eta-3850861.0/96000.0*t64+50965439.0/4608000.0*t5-19824371.0/1036800.0*t10-t176-545027.0/36864.0*t7+t154+t101-147.0/125.0*t1;
716  double t181 = 14651.0/15360.0*t18;
717  double t183 = 7203.0/2560.0*t38;
718  double t186 = 4459.0/4000.0*t48;
719  double t189 = 13696039.0/2400000.0*t55;
720  double t191 = t105-t181-18179.0/4800.0*t21-t183-t109-144893.0/96000.0*t44-31213.0/6400.0*t46-t186+18179.0/4800.0*t50+1083733.0/320000.0*t53+t189+4459.0/4000.0*t57;
721  double t193 = 1971907.0/576000.0*xi;
722  double t199 = 14651.0/15360.0*t2;
723  double t202 = 153713.0/2560000.0-446929.0/240000.0*t59-t193+144893.0/96000.0*t20-497791.0/600000.0*eta+212317.0/20000.0*t64-10357963.0/2560000.0*t5+19361321.0/1728000.0*t10+t199+199283.0/20480.0*t7-t127-t128-91679.0/160000.0*t1;
724  double t204 = 637.0/2560.0*t18;
725  double t209 = 14749.0/32000.0*t48;
726  double t212 = 3234931.0/1200000.0*t55;
727  double t214 = -t78+t204+3087.0/3200.0*t21+2401.0/2560.0*t38-t82+1519.0/4000.0*t44-26411.0/9600.0*t46+t209-3087.0/3200.0*t50-169393.0/160000.0*t53+t212-14749.0/32000.0*t57;
728  double t216 = 2074807.0/1728000.0*xi;
729  double t222 = 637.0/2560.0*t2;
730  double t225 = -6321.0/2560000.0+37387.0/160000.0*t59+t216-1519.0/4000.0*t20+5537.0/100000.0*eta+994357.0/240000.0*t64+22266139.0/23040000.0*t5-27233171.0/5184000.0*t10-t222-765919.0/184320.0*t7+t100+t101+38661.0/320000.0*t1;
731  double t227 = 85829.0/61440.0*t18;
732  double t232 = 83771.0/48000.0*t48;
733  double t235 = t4+t15+t17-t227+77597.0/57600.0*t21-t30+t37-87887.0/51200.0*t38-t41+t43+246511.0/128000.0*t44+594629.0/384000.0*t46-t232-243521.0/76800.0*t50+2126359.0/1280000.0*t53;
734  double t236 = 20141003.0/9600000.0*t55;
735  double t237 = 7123181.0/1728000.0*t57;
736  double t239 = 273563.0/432000.0*xi;
737  double t244 = 7123181.0/5184000.0*t10;
738  double t245 = 5077697.0/1658880.0*t2;
739  double t248 = -6820919.0/34560000.0-t236+t237+5917771.0/11520000.0*t59-t239-2159737.0/576000.0*t20+39583543.0/57600000.0*eta-940003.0/360000.0*t64-3953677.0/11520000.0*t5+t244-t245+947497.0/921600.0*t7-t70-t74+21015133.0/8640000.0*t1;
740  double t250 = 117649.0/6144.0*t3;
741  double t251 = 117649.0/6144.0*t16;
742  double t252 = 16807.0/6144.0*t18;
743  double t254 = 117649.0/103680.0*t29;
744  double t255 = 117649.0/207360.0*t36;
745  double t257 = 420175.0/20736.0*t40;
746  double t258 = 14651.0/12800.0*t42;
747  double t261 = 89523.0/64000.0*t48;
748  double t262 = 37473779.0/34560000.0-t250-t251+t252-55223.0/23040.0*t21+t254-t255+136171.0/76800.0*t38+t257-t258-2401.0/800.0*t44-41797.0/96000.0*t46+t261;
749  double t265 = 4969237.0/2400000.0*t55;
750  double t268 = 326879.0/216000.0*xi;
751  double t274 = 2907611.0/165888.0*t2;
752  double t277 = 818741.0/46080.0*t50-6341041.0/1920000.0*t53+t265-14426923.0/864000.0*t57-10863937.0/3840000.0*t59+t268+4232963.0/230400.0*t20-170592863.0/57600000.0*eta+7623959.0/2880000.0*t64+6379457.0/11520000.0*t5-29841.0/64000.0*t10+t274-136171.0/230400.0*t7-59484089.0/4320000.0*t1;
753  double t279 = 352947.0/10240.0*t3;
754  double t280 = 352947.0/10240.0*t16;
755  double t281 = 16807.0/3072.0*t18;
756  double t282 = 2401.0/640.0*t21;
757  double t283 = 117649.0/34560.0*t29;
758  double t284 = 117649.0/69120.0*t36;
759  double t286 = 16807.0/432.0*t40;
760  double t287 = 14651.0/6400.0*t42;
761  double t290 = 42189.0/16000.0*t48;
762  double t292 = t279+t280+t281-t282-t283+t284+19551.0/6400.0*t38-t286-t287-26411.0/6400.0*t44-34153.0/32000.0*t46+t290-160867.0/5120.0*t50;
763  double t294 = 2564513.0/800000.0*t55;
764  double t297 = 1251607.0/576000.0*xi;
765  double t303 = 1226911.0/27648.0*t2;
766  double t306 = -16460227.0/5760000.0-1841567.0/320000.0*t53+t294+17200421.0/576000.0*t57+5099087.0/1280000.0*t59-t297-794731.0/25600.0*t20+26096371.0/6400000.0*eta+180859.0/40000.0*t64+1582259.0/1920000.0*t5-14063.0/16000.0*t10-t303-6517.0/6400.0*t7+50686139.0/1440000.0*t1;
767  double t308 = 84035.0/4096.0*t18;
768  double t310 = 117649.0/20736.0*t29;
769  double t311 = 117649.0/41472.0*t36;
770  double t313 = 1932805.0/41472.0*t40;
771  double t314 = 14651.0/10240.0*t42;
772  double t317 = 35329.0/25600.0*t48;
773  double t319 = -t250-t251-t308+31213.0/3072.0*t21+t310-t311-95697.0/20480.0*t38+t313+t314+21609.0/2560.0*t44-11221.0/25600.0*t46-t317+26411.0/2304.0*t50;
774  double t321 = 1454467.0/320000.0*t55;
775  double t324 = 1911539.0/691200.0*xi;
776  double t330 = 22269275.0/331776.0*t2;
777  double t333 = 31917179.0/6912000.0+338541.0/16000.0*t53-t321-25062667.0/691200.0*t57+1815107.0/1152000.0*t59+t324+607453.0/46080.0*t20-501613.0/11520000.0*eta-771211.0/96000.0*t64-1510229.0/512000.0*t5+35329.0/76800.0*t10+t330+31899.0/20480.0*t7-748234949.0/13824000.0*t1;
778  double t337 = 1127.0/1280.0*t42;
779  double t340 = 3087.0/3200.0*t48;
780  double t341 = -65455033.0/13824000.0-t250-t251+t308-2401.0/384.0*t21-t310+t311+4459.0/2560.0*t38-t313-t337-7203.0/5120.0*t44-147.0/800.0*t46-t340;
781  double t344 = 367059.0/640000.0*t55;
782  double t347 = 612941.0/172800.0*xi;
783  double t355 = 256907.0/9216.0*t50-4756381.0/256000.0*t53+t344+13341671.0/345600.0*t57-19272239.0/2304000.0*t59-t347+530621.0/23040.0*t20-62402137.0/11520000.0*eta+50617.0/12000.0*t64+3210137.0/1536000.0*t5+1029.0/3200.0*t10-t330-4459.0/7680.0*t7+376548487.0/6912000.0*t1;
784  double t359 = 14651.0/38400.0*t42;
785  double t362 = 24353.0/32000.0*t48;
786  double t363 = 4250113.0/1440000.0+t279+t280-t281+2401.0/3840.0*t21+t283-t284-18179.0/25600.0*t38+t286+t359-2401.0/1600.0*t44+9359.0/96000.0*t46-t362;
787  double t366 = 335209.0/600000.0*t55;
788  double t369 = 1611757.0/576000.0*xi;
789  double t377 = -549829.0/15360.0*t50+295323.0/80000.0*t53+t366-18280871.0/576000.0*t57+25791689.0/3840000.0*t59+t369-861959.0/25600.0*t20+35582379.0/6400000.0*eta-211729.0/240000.0*t64-261709.0/1920000.0*t5+24353.0/96000.0*t10+t303+18179.0/76800.0*t7-204365231.0/5760000.0*t1;
790  double t381 = 637.0/6400.0*t42;
791  double t384 = 7889.0/48000.0*t48;
792  double t386 = -t250-t251-t252-2401.0/11520.0*t21-t254+t255+9261.0/51200.0*t38-t257-t381-21609.0/12800.0*t44-3479.0/128000.0*t46+t384+717899.0/46080.0*t50;
793  double t388 = 371273.0/400000.0*t55;
794  double t391 = 1714657.0/1728000.0*xi;
795  double t399 = -17386327.0/17280000.0+506611.0/320000.0*t53+t388+26152721.0/1728000.0*t57-5661313.0/3840000.0*t59-t391+3930437.0/230400.0*t20-128307137.0/57600000.0*eta+356083.0/1440000.0*t64+74431.0/3840000.0*t5-7889.0/144000.0*t10-t274-3087.0/51200.0*t7+467769337.0/34560000.0*t1;
796  double t403 = 85829.0/153600.0*t42;
797  double t409 = t4-t15+t17+t227+77597.0/230400.0*t21+t30-t37+77597.0/307200.0*t38+t41+t403+84457.0/64000.0*t44+84457.0/128000.0*t46+73481.0/128000.0*t48-13827.0/6400.0*t50-716051.0/1920000.0*t53;
798  double t421 = 9594313.0/69120000.0-1343687.0/1600000.0*t55-10198837.0/3456000.0*t57-496649.0/2880000.0*t59+839329.0/3456000.0*xi-3621101.0/1152000.0*t20+18298957.0/57600000.0*eta-716051.0/2880000.0*t64-496649.0/5760000.0*t5-10198837.0/10368000.0*t10+t245-13827.0/25600.0*t7+t70+t74-155978489.0/69120000.0*t1;
799  double t425 = 16807.0/15360.0*t42;
800  double t431 = t78-t204+3087.0/12800.0*t21-2401.0/15360.0*t38-t425-3479.0/64000.0*t44-21609.0/25600.0*t46+t384-3087.0/12800.0*t50+356083.0/960000.0*t53+t388-7889.0/48000.0*t57;
801  double t440 = 259063.0/7680000.0+74431.0/1920000.0*t59-t391+3479.0/64000.0*t20-269059.0/3200000.0*eta+506611.0/480000.0*t64-5661313.0/7680000.0*t5+26152721.0/5184000.0*t10+t222+717899.0/184320.0*t7-t100-t101-137347.0/640000.0*t1;
802  double t444 = 16807.0/7680.0*t42;
803  double t449 = -393813.0/2560000.0-t105+t181-18179.0/19200.0*t21+2401.0/5120.0*t38-t444+9359.0/48000.0*t44-2401.0/3200.0*t46-t362+18179.0/19200.0*t50-211729.0/160000.0*t53+t366;
804  double t459 = 24353.0/32000.0*t57-261709.0/960000.0*t59+t369-9359.0/48000.0*t20+229957.0/1200000.0*eta+98441.0/40000.0*t64+25791689.0/7680000.0*t5-18280871.0/1728000.0*t10-t199-549829.0/61440.0*t7+t127+t128+546791.0/640000.0*t1;
805  double t463 = 16807.0/2048.0*t42;
806  double t468 = 207221.0/512000.0+t132-t158+4459.0/1920.0*t21-2401.0/512.0*t38+t463-147.0/400.0*t44-7203.0/10240.0*t46-t340-4459.0/1920.0*t50+50617.0/8000.0*t53+t344;
807  double t478 = 3087.0/3200.0*t57+3210137.0/768000.0*t59-t347+147.0/400.0*t20+166257.0/1280000.0*eta-4756381.0/384000.0*t64-19272239.0/4608000.0*t5+13341671.0/1036800.0*t10+t176+256907.0/36864.0*t7-t154-t101-104811.0/64000.0*t1;
808  double t487 = -t132+t133-31899.0/5120.0*t21+31213.0/4096.0*t38-t463-11221.0/12800.0*t44+21609.0/5120.0*t46-t317+31899.0/5120.0*t50-771211.0/64000.0*t53-t321+35329.0/25600.0*t57;
809  double t496 = -236621.0/512000.0-1510229.0/256000.0*t59+t324+11221.0/12800.0*t20+207809.0/640000.0*eta+112847.0/8000.0*t64+1815107.0/2304000.0*t5-25062667.0/2073600.0*t10-t152+26411.0/9216.0*t7+t154-t101+2268063.0/512000.0*t1;
810  double t504 = t105-t106+6517.0/1600.0*t21-t183+t444-34153.0/16000.0*t44-26411.0/12800.0*t46+t290-6517.0/1600.0*t50+542577.0/80000.0*t53+t294-42189.0/16000.0*t57;
811  double t513 = 303163.0/2560000.0+1582259.0/960000.0*t59-t297+34153.0/16000.0*t20-1827651.0/1600000.0*eta-1841567.0/480000.0*t64+5099087.0/2560000.0*t5+17200421.0/1728000.0*t10+t125-160867.0/20480.0*t7-t127+t128-1269933.0/320000.0*t1;
812  double t521 = 419979.0/2560000.0-t78-t79+136171.0/57600.0*t21-55223.0/30720.0*t38+t425-41797.0/48000.0*t44-2401.0/1600.0*t46+t261-136171.0/57600.0*t50+7623959.0/1920000.0*t53+t265;
813  double t531 = -89523.0/64000.0*t57+6379457.0/5760000.0*t59+t268+41797.0/48000.0*t20-1367737.0/2400000.0*eta-6341041.0/2880000.0*t64-10863937.0/7680000.0*t5-14426923.0/2592000.0*t10+t98+818741.0/184320.0*t7+t100-t101-2781681.0/1280000.0*t1;
814  double t539 = t4+t15-t17+t19-87887.0/38400.0*t21-t30+t37+77597.0/76800.0*t38-t41-t403+594629.0/192000.0*t44+246511.0/256000.0*t46-t232+947497.0/230400.0*t50-940003.0/240000.0*t53;
815  double t547 = -17668213.0/69120000.0-t236+t237-3953677.0/5760000.0*t59-t239-5668649.0/1152000.0*t20+36656459.0/28800000.0*eta+2126359.0/1920000.0*t64+5917771.0/23040000.0*t5+t244-t68-243521.0/307200.0*t7-t70+t72+t74+196867657.0/34560000.0*t1;
816  double t549 = 16807.0/1536.0*t18;
817  double t556 = -t250+t251-t549+2401.0/1920.0*t21+t254-t255+9261.0/12800.0*t38+t257+t381-26411.0/4800.0*t44+1519.0/8000.0*t46+t209-765919.0/46080.0*t50+994357.0/160000.0*t53;
818  double t564 = 14101073.0/165888.0*t2;
819  double t566 = 5764801.0/92160.0*t71;
820  double t568 = 41641229.0/34560000.0+t212-27233171.0/1728000.0*t57+22266139.0/11520000.0*t59+t216+4809203.0/230400.0*t20-242622863.0/57600000.0*eta-169393.0/240000.0*t64+37387.0/320000.0*t5-14749.0/96000.0*t10+t564-3087.0/12800.0*t7-t566-961071937.0/34560000.0*t1;
821  double t570 = 16807.0/768.0*t18;
822  double t576 = t279-t280-t570-t282-t283+t284-18179.0/6400.0*t38-t286-t359-31213.0/3200.0*t44-144893.0/192000.0*t46-t186+199283.0/5120.0*t50+636951.0/40000.0*t53;
823  double t584 = 12285917.0/55296.0*t2;
824  double t586 = 5764801.0/30720.0*t71;
825  double t588 = -38554229.0/11520000.0+t189+19361321.0/576000.0*t57-10357963.0/1280000.0*t59-t193-650671.0/25600.0*t20+8088871.0/6400000.0*eta+1083733.0/480000.0*t64-446929.0/480000.0*t5+4459.0/12000.0*t10-t584+18179.0/19200.0*t7+t586+724916437.0/11520000.0*t1;
826  double t590 = 84035.0/1024.0*t18;
827  double t596 = 42104279.0/6912000.0-t250+t251+t590+2401.0/64.0*t21+t310-t311+4459.0/640.0*t38+t313+t337+40817.0/1280.0*t44+12201.0/6400.0*t46+t163-545027.0/9216.0*t50;
828  double t605 = 56891695.0/165888.0*t2;
829  double t607 = 5764801.0/18432.0*t71;
830  double t609 = -3850861.0/64000.0*t53-t166-19824371.0/345600.0*t57+50965439.0/2304000.0*t59+t170-472997.0/46080.0*t20+134554637.0/11520000.0*eta-287581.0/16000.0*t64+1399097.0/384000.0*t5-7889.0/1200.0*t10+t605-4459.0/1920.0*t7-t607-622119337.0/6912000.0*t1;
831  double t616 = -45345629.0/6912000.0-t250+t251-t590-31213.0/512.0*t21-t310+t311-95697.0/5120.0*t38-t313-t314-69629.0/2560.0*t44-309533.0/51200.0*t46-t139+362551.0/9216.0*t50;
832  double t627 = 140973.0/2000.0*t53+t142+23065721.0/345600.0*t57-3445729.0/144000.0*t59-t146+2249737.0/46080.0*t20-52741003.0/2880000.0*eta+3889473.0/128000.0*t64-906549.0/128000.0*t5+62083.0/6400.0*t10-t605+31899.0/5120.0*t7+t607+631843387.0/6912000.0*t1;
833  double t635 = t279-t280+t570+7203.0/320.0*t21+t283-t284+19551.0/1600.0*t38+t286+t287+2401.0/400.0*t44+83153.0/32000.0*t46+t112+64827.0/5120.0*t50-1934177.0/80000.0*t53;
834  double t645 = 42876029.0/11520000.0-t115-23683121.0/576000.0*t57+5465313.0/1280000.0*t59+t119-1054039.0/25600.0*t20+59592379.0/6400000.0*eta-1021993.0/80000.0*t64+1767479.0/480000.0*t5-22981.0/8000.0*t10+t584-6517.0/1600.0*t7-t586-737881837.0/11520000.0*t1;
835  double t653 = -t250+t251+t549+55223.0/3840.0*t21-t254+t255+136171.0/19200.0*t38-t257+t258+45619.0/19200.0*t44+324919.0/192000.0*t46+t85-1370971.0/46080.0*t50-3142223.0/240000.0*t53;
836  double t663 = -30836729.0/34560000.0-t88+16428671.0/1728000.0*t57+93337111.0/11520000.0*t59+t92+2994047.0/230400.0*t20-11258387.0/57600000.0*eta-633619.0/80000.0*t64+6266267.0/2880000.0*t5-556003.0/288000.0*t10-t564-136171.0/57600.0*t7+t566+928658437.0/34560000.0*t1;
837  double t665 = 15625.0/1024.0*eta;
838  double t666 = 390625.0/12288.0*t2;
839  double t668 = 15625.0/512.0*t20;
840  double t675 = -3125.0/1536.0-t665+t666-53125.0/4096.0*t1+t668-234375.0/4096.0*t38-15625.0/1024.0*xi-209375.0/12288.0*t5-209375.0/6144.0*t59+15625.0/1024.0*t10+15625.0/512.0*t55;
841  double t681 = 15625.0/1024.0*t46;
842  double t683 = 15625.0/512.0*t44;
843  double t684 = 78125.0/6144.0*t42;
844  double t686 = 390625.0/12288.0*t18;
845  double t687 = 46875.0/1024.0*t57+78125.0/4096.0*t7+71875.0/1536.0*t64+71875.0/1024.0*t53+78125.0/1024.0*t50-t681-46875.0/1024.0*t48-t683-t684-78125.0/1024.0*t21-t686;
846  double t691 = 15625.0/768.0*xi;
847  double t694 = 15625.0/768.0*t10;
848  double t695 = 15625.0/384.0*t55;
849  double t696 = 15625.0/256.0*t57;
850  double t698 = 34375.0/768.0*t64;
851  double t699 = 34375.0/512.0*t53;
852  double t700 = 78125.0/1536.0*t50;
853  double t702 = 15625.0/256.0*t48;
854  double t703 = 78125.0/3072.0*t42;
855  double t704 = 78125.0/1536.0*t21;
856  double t705 = 3125.0/1024.0+t665-9375.0/1024.0*t1-t668+78125.0/1024.0*t38+t691+34375.0/1536.0*t5+59375.0/3072.0*t59-t694-t695-t696-78125.0/3072.0*t7-t698-t699-t700+78125.0/3072.0*t46+t702+t683+t703+t704;
857  double t707 = 390625.0/6144.0*t2;
858  double t710 = 78125.0/2048.0*t38;
859  double t714 = 78125.0/6144.0*t7;
860  double t718 = 390625.0/6144.0*t18;
861  double t719 = -t695-t696-t714-t698-t699-78125.0/768.0*t50+t681+t702+78125.0/1536.0*t44+78125.0/768.0*t21+t718;
862  double t721 = 15625.0/1536.0*xi;
863  double t724 = 15625.0/1536.0*t10;
864  double t725 = 15625.0/768.0*t55;
865  double t726 = 15625.0/512.0*t57;
866  double t727 = 34375.0/1536.0*t64;
867  double t728 = 34375.0/1024.0*t53;
868  double t729 = 15625.0/768.0*t46;
869  double t731 = 15625.0/512.0*t48;
870  double t732 = -3125.0/2048.0-t721-34375.0/3072.0*t5+3125.0/1024.0*t59+t724+t725+t726+t714+t727+t728-t729+9375.0/2048.0*t1-t731-t703-t710;
871  double t739 = 15625.0/768.0*eta;
872  double t741 = 15625.0/384.0*t20;
873  double t742 = 15625.0/384.0*t44;
874  double t745 = -3125.0/512.0-t691-3125.0/288.0*t5-3125.0/144.0*t59+t694+t695+t696+78125.0/4608.0*t7+3125.0/144.0*t64+3125.0/96.0*t53+78125.0/1152.0*t50-t729-t739+9375.0/512.0*t1+t741-t702-t742-78125.0/1536.0*t38-78125.0/1152.0*t21;
875  double t749 = -3125.0/2048.0-t721+3125.0/2048.0*t5-34375.0/1536.0*t59+t724+t725+t726+t727+t728+t700-t739-34375.0/1024.0*t1+t741+t707-t731-t742-t704-t718;
876  double t751 = 8125.0/1024.0*t55;
877  double t756 = 15625.0/2304.0*xi;
878  double t759 = 15625.0/2304.0*t10;
879  double t760 = 6875.0/768.0*t55;
880  double t761 = 15625.0/768.0*t57;
881  double t768 = 15625.0/768.0*t48;
882  double t770 = 3125.0/3072.0+t756+34375.0/4608.0*t5-3125.0/384.0*t59-t759-t760-t761-78125.0/9216.0*t7+3125.0/384.0*t64-34375.0/1536.0*t53+15625.0/1536.0*t46-625.0/512.0*eta-3125.0/1024.0*t1+t768+78125.0/3072.0*t38;
883  double t781 = 3125.0/768.0+t756-3125.0/768.0*t5+34375.0/2304.0*t59-t759-t760-t761-34375.0/2304.0*t64+3125.0/256.0*t53-78125.0/2304.0*t50+6875.0/768.0*eta-3125.0/256.0*t1-15625.0/768.0*t20+t768+15625.0/768.0*t44+78125.0/2304.0*t21;
884 
885  values[0] = t52+t76;
886  values[1] = t90+t103;
887  values[2] = t117+t130;
888  values[3] = t143+t156;
889  values[4] = t167+t179;
890  values[5] = t191+t202;
891  values[6] = t214+t225;
892  values[7] = t235+t248;
893  values[8] = t262+t277;
894  values[9] = t292+t306;
895  values[10] = t319+t333;
896  values[11] = t341+t355;
897  values[12] = t363+t377;
898  values[13] = t386+t399;
899  values[14] = t409+t421;
900  values[15] = t431+t440;
901  values[16] = t449+t459;
902  values[17] = t468+t478;
903  values[18] = t487+t496;
904  values[19] = t504+t513;
905  values[20] = t521+t531;
906  values[21] = t539+t547;
907  values[22] = t556+t568;
908  values[23] = t576+t588;
909  values[24] = t596+t609;
910  values[25] = t616+t627;
911  values[26] = t635+t645;
912  values[27] = t653+t663;
913  values[28] = t675+t687;
914  values[29] = t705;
915  values[30] = 3125.0/1024.0+78125.0/3072.0*eta-t707+59375.0/2048.0*t1-78125.0/1536.0*t20+t710+t691+59375.0/6144.0*t5+34375.0/768.0*t59-t694+t719;
916  values[31] = t732;
917  values[32] = t745;
918  values[33] = t749;
919  values[34] = 3125.0/6144.0*t59-t751-40625.0/3072.0*t64+15625.0/2048.0*t46+625.0/2048.0*eta+t684;
920  values[35] = t770;
921  values[36] = t781;
922  values[37] = -3125.0/12288.0+3125.0/12288.0*t5-t751-40625.0/2048.0*t53+8125.0/1024.0*eta+40625.0/2048.0*t1-15625.0/1024.0*t20-t666+15625.0/1024.0*t44+t686;
923 }
924 
925 // values of the derivatives in xi-xi direction
926 static void C_Q_M7_2D_DeriveXiXi(double xi, double eta, double *values)
927 {
928  double t3 = xi*xi;
929  double t4 = t3*t3;
930  double t5 = t4*xi;
931  double t7 = xi*t3;
932  double t11 = 1.0+eta;
933  double t12 = (63.0/8.0*t5-35.0/4.0*t7+15.0/8.0*xi)*t11;
934  double t13 = 16807.0/103680.0*t12;
935  double t15 = 12113881.0/207360.0*t7;
936  double t16 = 117649.0/6144.0*t4;
937  double t17 = 823543.0/15360.0*t5;
938  double t18 = eta*eta;
939  double t20 = t18*t18;
940  double t21 = t20*eta;
941  double t22 = 85829.0/38400.0*t21;
942  double t23 = t18*eta;
943  double t26 = -49159793.0/57600000.0-172582507.0/17280000.0*xi-45493571.0/11520000.0*eta+t13-922751.0/384000.0*t3+t15+t16-t17+225377.0/150000.0*t18-t22+967443.0/160000.0*t23-98177.0/192000.0*t20;
944  double t27 = xi*eta;
945  double t29 = t3*eta;
946  double t31 = xi*t18;
947  double t33 = t7*eta;
948  double t34 = 184877.0/10368.0*t33;
949  double t35 = t3*t18;
950  double t37 = xi*t23;
951  double t39 = t4*eta;
952  double t40 = 117649.0/6144.0*t39;
953  double t41 = t7*t18;
954  double t42 = 85829.0/3840.0*t41;
955  double t43 = t3*t23;
956  double t45 = xi*t20;
957  double t51 = xi*(315.0/8.0*t4-105.0/4.0*t3+15.0/8.0)*t11;
958  double t52 = 16807.0/51840.0*t51;
959  double t58 = (1.0-t3)*(315.0/2.0*t7-105.0/2.0*xi)*t11;
960  double t59 = 16807.0/207360.0*t58;
961  double t60 = 971869.0/864000.0*t27+2529463.0/76800.0*t29+2902329.0/160000.0*t31+t34-98177.0/32000.0*t35-94061.0/16000.0*t37-t40-t42-87887.0/3200.0*t43-87887.0/6400.0*t45+t52-t59;
962  double t64 = 117649.0/103680.0*t12;
963  double t66 = 14101073.0/41472.0*t7;
964  double t67 = 588245.0/6144.0*t4;
965  double t68 = 5764801.0/15360.0*t5;
966  double t69 = 141169.0/50000.0*t18;
967  double t70 = 14651.0/12800.0*t21;
968  double t73 = -11258387.0/57600000.0+928658437.0/17280000.0*xi+93337111.0/11520000.0*eta-t64+2994047.0/76800.0*t3-t66-t67+t68-t69+t70-633619.0/80000.0*t23+324919.0/192000.0*t20;
969  double t77 = 420175.0/5184.0*t33;
970  double t79 = 556003.0/48000.0*t37;
971  double t80 = 588245.0/6144.0*t39;
972  double t81 = 16807.0/384.0*t41;
973  double t84 = 117649.0/51840.0*t51;
974  double t85 = 117649.0/207360.0*t58;
975  double t86 = 16428671.0/864000.0*t27-1370971.0/15360.0*t29-3142223.0/120000.0*t31-t77+45619.0/6400.0*t35+t79+t80+t81+55223.0/1280.0*t43+136171.0/9600.0*t45-t84+t85;
976  double t90 = 117649.0/34560.0*t12;
977  double t92 = 12285917.0/13824.0*t7;
978  double t93 = 352947.0/2048.0*t4;
979  double t94 = 5764801.0/5120.0*t5;
980  double t95 = 2277569.0/400000.0*t18;
981  double t96 = 14651.0/6400.0*t21;
982  double t99 = 59592379.0/6400000.0-737881837.0/5760000.0*xi+5465313.0/1280000.0*eta+t90-3162117.0/25600.0*t3+t92+t93-t94-t95+t96-1021993.0/80000.0*t23+83153.0/32000.0*t20;
983  double t103 = 16807.0/108.0*t33;
984  double t105 = 68943.0/4000.0*t37;
985  double t106 = 352947.0/2048.0*t39;
986  double t107 = 16807.0/192.0*t41;
987  double t110 = 117649.0/17280.0*t51;
988  double t111 = 117649.0/69120.0*t58;
989  double t112 = -23683121.0/288000.0*t27+194481.0/5120.0*t29-1934177.0/40000.0*t31+t103+7203.0/400.0*t35+t105-t106+t107+21609.0/320.0*t43+19551.0/800.0*t45+t110-t111;
990  double t116 = 117649.0/20736.0*t12;
991  double t118 = 56891695.0/41472.0*t7;
992  double t119 = 5764801.0/3072.0*t5;
993  double t120 = 12372059.0/640000.0*t18;
994  double t121 = 14651.0/10240.0*t21;
995  double t124 = -52741003.0/2880000.0+631843387.0/3456000.0*xi-3445729.0/144000.0*eta-t116+2249737.0/15360.0*t3-t118-t67+t119+t120-t121+3889473.0/128000.0*t23-309533.0/51200.0*t20;
996  double t128 = 1932805.0/10368.0*t33;
997  double t130 = 186249.0/3200.0*t37;
998  double t131 = 84035.0/256.0*t41;
999  double t134 = 117649.0/10368.0*t51;
1000  double t135 = 117649.0/41472.0*t58;
1001  double t136 = 23065721.0/172800.0*t27+362551.0/3072.0*t29+140973.0/1000.0*t31-t128-208887.0/2560.0*t35-t130+t80-t131-93639.0/512.0*t43-95697.0/2560.0*t45-t134+t135;
1002  double t141 = 2978171.0/160000.0*t18;
1003  double t142 = 1127.0/1280.0*t21;
1004  double t145 = 134554637.0/11520000.0-622119337.0/3456000.0*xi+50965439.0/2304000.0*eta+t116-472997.0/15360.0*t3+t118-t67-t119-t141+t142-287581.0/16000.0*t23+12201.0/6400.0*t20;
1005  double t150 = 7889.0/200.0*t37;
1006  double t153 = -19824371.0/172800.0*t27-545027.0/3072.0*t29-3850861.0/32000.0*t31+t128+122451.0/1280.0*t35+t150+t80+t131+7203.0/64.0*t43+4459.0/320.0*t45+t134-t135;
1007  double t158 = 13696039.0/2400000.0*t18;
1008  double t159 = 14651.0/38400.0*t21;
1009  double t162 = 8088871.0/6400000.0+724916437.0/5760000.0*xi-10357963.0/1280000.0*eta-t90-1952013.0/25600.0*t3-t92+t93+t94+t158-t159+1083733.0/480000.0*t23-144893.0/192000.0*t20;
1010  double t167 = 4459.0/2000.0*t37;
1011  double t168 = 7203.0/640.0*t43;
1012  double t170 = 19361321.0/288000.0*t27+597849.0/5120.0*t29+636951.0/20000.0*t31-t103-93639.0/3200.0*t35-t167-t106-t107-t168-18179.0/3200.0*t45-t110+t111;
1013  double t175 = 3234931.0/1200000.0*t18;
1014  double t176 = 637.0/6400.0*t21;
1015  double t179 = -242622863.0/57600000.0-961071937.0/17280000.0*xi+22266139.0/11520000.0*eta+t64+4809203.0/76800.0*t3+t66-t67-t68+t175+t176-169393.0/240000.0*t23+1519.0/8000.0*t20;
1016  double t184 = 14749.0/16000.0*t37;
1017  double t187 = -27233171.0/864000.0*t27-765919.0/15360.0*t29+994357.0/80000.0*t31+t77-26411.0/1600.0*t35+t184+t80-t81+2401.0/640.0*t43+9261.0/6400.0*t45+t84-t85;
1018  double t192 = 20141003.0/9600000.0*t18;
1019  double t193 = 85829.0/153600.0*t21;
1020  double t196 = 7123181.0/864000.0*t27;
1021  double t197 = 196867657.0/17280000.0*xi-3953677.0/5760000.0*eta-t13-5668649.0/384000.0*t3-t15+t16+t17-t192-t193+2126359.0/1920000.0*t23+246511.0/256000.0*t20+t196;
1022  double t201 = 83771.0/24000.0*t37;
1023  double t204 = 36656459.0/28800000.0+947497.0/76800.0*t29-940003.0/120000.0*t31-t34+594629.0/64000.0*t35-t201-t40+t42-87887.0/12800.0*t43+77597.0/38400.0*t45-t52+t59;
1024  double t208 = 14651.0/1280.0*t7;
1025  double t211 = 4969237.0/2400000.0*t18;
1026  double t216 = 89523.0/32000.0*t37;
1027  double t218 = 14651.0/1280.0*t41;
1028  double t221 = 16807.0/15360.0*t21;
1029  double t222 = -1367737.0/2400000.0-2781681.0/640000.0*xi+41797.0/16000.0*t3+t208+6379457.0/5760000.0*eta-89523.0/32000.0*t27+t211-136171.0/19200.0*t29+7623959.0/960000.0*t31-6341041.0/2880000.0*t23-41797.0/16000.0*t35+t216-2401.0/1600.0*t20-t218+136171.0/19200.0*t43-55223.0/15360.0*t45+t221;
1030  double t225 = 14651.0/640.0*t7;
1031  double t228 = 2564513.0/800000.0*t18;
1032  double t233 = 42189.0/8000.0*t37;
1033  double t235 = 14651.0/640.0*t41;
1034  double t237 = 7203.0/1280.0*t45;
1035  double t238 = 16807.0/7680.0*t21;
1036  double t239 = -1827651.0/1600000.0-1269933.0/160000.0*xi+102459.0/16000.0*t3+t225+1582259.0/960000.0*eta-42189.0/8000.0*t27+t228-19551.0/1600.0*t29+542577.0/40000.0*t31-1841567.0/480000.0*t23-102459.0/16000.0*t35+t233-26411.0/12800.0*t20-t235+19551.0/1600.0*t43-t237+t238;
1037  double t242 = 14651.0/1024.0*t7;
1038  double t245 = 1454467.0/320000.0*t18;
1039  double t250 = 35329.0/12800.0*t37;
1040  double t252 = 14651.0/1024.0*t41;
1041  double t255 = 16807.0/2048.0*t21;
1042  double t256 = 207809.0/640000.0+2268063.0/256000.0*xi+33663.0/12800.0*t3-t242-1510229.0/256000.0*eta+35329.0/12800.0*t27-t245+95697.0/5120.0*t29-771211.0/32000.0*t31+112847.0/8000.0*t23-33663.0/12800.0*t35-t250+21609.0/5120.0*t20+t252-95697.0/5120.0*t43+31213.0/2048.0*t45-t255;
1043  double t259 = 1127.0/128.0*t7;
1044  double t262 = 367059.0/640000.0*t18;
1045  double t267 = 3087.0/1600.0*t37;
1046  double t269 = 1127.0/128.0*t41;
1047  double t272 = 166257.0/1280000.0-104811.0/32000.0*xi+441.0/400.0*t3+t259+3210137.0/768000.0*eta+3087.0/1600.0*t27+t262-4459.0/640.0*t29+50617.0/4000.0*t31-4756381.0/384000.0*t23-441.0/400.0*t35-t267-7203.0/10240.0*t20-t269+4459.0/640.0*t43-2401.0/256.0*t45+t255;
1048  double t275 = 14651.0/3840.0*t7;
1049  double t278 = 335209.0/600000.0*t18;
1050  double t283 = 24353.0/16000.0*t37;
1051  double t285 = 14651.0/3840.0*t41;
1052  double t288 = 229957.0/1200000.0+546791.0/320000.0*xi-9359.0/16000.0*t3-t275-261709.0/960000.0*eta+24353.0/16000.0*t27+t278+18179.0/6400.0*t29-211729.0/80000.0*t31+98441.0/40000.0*t23+9359.0/16000.0*t35-t283-2401.0/3200.0*t20+t285-18179.0/6400.0*t43+2401.0/2560.0*t45-t238;
1053  double t291 = 637.0/640.0*t7;
1054  double t294 = 371273.0/400000.0*t18;
1055  double t299 = 7889.0/24000.0*t37;
1056  double t301 = 637.0/640.0*t41;
1057  double t304 = -269059.0/3200000.0-137347.0/320000.0*xi+10437.0/64000.0*t3+t291+74431.0/1920000.0*eta-7889.0/24000.0*t27+t294-9261.0/12800.0*t29+356083.0/480000.0*t31+506611.0/480000.0*t23-10437.0/64000.0*t35+t299-21609.0/25600.0*t20-t301+9261.0/12800.0*t43-2401.0/7680.0*t45-t221;
1058  double t307 = 5077697.0/414720.0*t7;
1059  double t312 = 18298957.0/57600000.0+t52-t59-155978489.0/34560000.0*xi-3621101.0/384000.0*t3+t307+t16-496649.0/2880000.0*eta-10198837.0/1728000.0*t27-1343687.0/1600000.0*t18-41481.0/6400.0*t29;
1060  double t318 = 85829.0/15360.0*t41;
1061  double t321 = -716051.0/960000.0*t31-716051.0/2880000.0*t23+t34+253371.0/64000.0*t35+73481.0/64000.0*t37+84457.0/128000.0*t20+t40+t318+77597.0/76800.0*t43+77597.0/153600.0*t45+t193+t13;
1062  double t325 = 2907611.0/41472.0*t7;
1063  double t329 = -128307137.0/57600000.0-t84+t85+467769337.0/17280000.0*xi+3930437.0/76800.0*t3-t325-t67-5661313.0/3840000.0*eta+26152721.0/864000.0*t27+t294+717899.0/15360.0*t29;
1064  double t334 = 16807.0/1536.0*t41;
1065  double t337 = 506611.0/160000.0*t31+356083.0/1440000.0*t23-t77-64827.0/12800.0*t35+t299-3479.0/128000.0*t20-t80-t334-2401.0/3840.0*t43+9261.0/25600.0*t45-t176-t64;
1066  double t341 = 1226911.0/6912.0*t7;
1067  double t345 = 35582379.0/6400000.0+t110-t111-204365231.0/2880000.0*xi-2585877.0/25600.0*t3+t341+t93+25791689.0/3840000.0*eta-18280871.0/288000.0*t27+t278-549829.0/5120.0*t29;
1068  double t350 = 16807.0/768.0*t41;
1069  double t353 = 295323.0/40000.0*t31-211729.0/240000.0*t23+t103-7203.0/1600.0*t35-t283+9359.0/96000.0*t20+t106-t350+2401.0/1280.0*t43-18179.0/12800.0*t45+t159+t90;
1070  double t357 = 22269275.0/82944.0*t7;
1071  double t361 = -62402137.0/11520000.0-t134+t135+376548487.0/3456000.0*xi+530621.0/7680.0*t3-t357-t67-19272239.0/2304000.0*eta+13341671.0/172800.0*t27+t262+256907.0/3072.0*t29;
1072  double t366 = 84035.0/1024.0*t41;
1073  double t369 = -4756381.0/128000.0*t31+50617.0/12000.0*t23-t128-21609.0/5120.0*t35-t267-147.0/800.0*t20-t80+t366-2401.0/128.0*t43+4459.0/1280.0*t45-t142-t116;
1074  double t376 = -501613.0/11520000.0+t134-t135-748234949.0/6912000.0*xi+607453.0/15360.0*t3+t357-t67+1815107.0/1152000.0*eta-25062667.0/345600.0*t27-t245+26411.0/768.0*t29;
1075  double t383 = 338541.0/8000.0*t31-771211.0/96000.0*t23+t128+64827.0/2560.0*t35-t250-11221.0/25600.0*t20-t80-t366+31213.0/1024.0*t43-95697.0/10240.0*t45+t121+t116;
1076  double t390 = 26096371.0/6400000.0-t110+t111+50686139.0/720000.0*xi-2384193.0/25600.0*t3-t341+t93+5099087.0/1280000.0*eta+17200421.0/288000.0*t27+t228-482601.0/5120.0*t29;
1077  double t396 = -1841567.0/160000.0*t31+180859.0/40000.0*t23-t103-79233.0/6400.0*t35+t233-34153.0/32000.0*t20+t106+t350-t168+19551.0/3200.0*t45-t96-t90;
1078  double t403 = -170592863.0/57600000.0+t84-t85-59484089.0/2160000.0*xi+4232963.0/76800.0*t3+t325-t67-10863937.0/3840000.0*eta-14426923.0/432000.0*t27+t211+818741.0/15360.0*t29;
1079  double t410 = -6341041.0/960000.0*t31+7623959.0/2880000.0*t23+t77-7203.0/800.0*t35+t216-41797.0/96000.0*t20-t80+t334-55223.0/7680.0*t43+136171.0/38400.0*t45-t70+t64;
1080  double t416 = 39583543.0/57600000.0-t52+t59+21015133.0/4320000.0*xi-2159737.0/192000.0*t3-t307+t16+5917771.0/11520000.0*eta+t196-t192-243521.0/25600.0*t29;
1081  double t423 = 2126359.0/640000.0*t31-940003.0/360000.0*t23-t34+739533.0/128000.0*t35-t201+594629.0/384000.0*t20+t40-t318+77597.0/19200.0*t43-87887.0/25600.0*t45+t22-t13;
1082  double t436 = 16807.0/3840.0*t21;
1083  double t437 = 5537.0/100000.0+38661.0/160000.0*xi-4557.0/4000.0*t3-t291+37387.0/160000.0*eta-14749.0/16000.0*t27+t175-9261.0/3200.0*t29-169393.0/80000.0*t31+994357.0/240000.0*t23+4557.0/4000.0*t35+t184-26411.0/9600.0*t20+t301+9261.0/3200.0*t43+2401.0/1280.0*t45-t436;
1084  double t448 = 16807.0/1920.0*t21;
1085  double t449 = -497791.0/600000.0-91679.0/80000.0*xi+144893.0/32000.0*t3+t275-446929.0/240000.0*eta+4459.0/2000.0*t27+t158+18179.0/1600.0*t29+1083733.0/160000.0*t31+212317.0/20000.0*t23-144893.0/32000.0*t35-t167-31213.0/6400.0*t20-t285-18179.0/1600.0*t43-t237-t448;
1086  double t461 = 16807.0/512.0*t21;
1087  double t462 = 854217.0/320000.0-294.0/125.0*xi-36603.0/3200.0*t3-t259+1399097.0/192000.0*eta-7889.0/200.0*t27-t141-4459.0/160.0*t29-862743.0/16000.0*t31-3850861.0/96000.0*t23+36603.0/3200.0*t35+t150+40817.0/2560.0*t20+t269+4459.0/160.0*t43+7203.0/128.0*t45+t461;
1088  double t474 = -1834217.0/320000.0+4557.0/16000.0*xi+928599.0/25600.0*t3+t242-906549.0/64000.0*eta+186249.0/3200.0*t27+t120+95697.0/1280.0*t29+11668419.0/128000.0*t31+46991.0/1000.0*t23-928599.0/25600.0*t35-t130-69629.0/5120.0*t20-t252-95697.0/1280.0*t43-93639.0/1024.0*t45-t461;
1089  double t486 = 1077069.0/400000.0+182427.0/40000.0*xi-249459.0/16000.0*t3-t225+1767479.0/240000.0*eta-68943.0/4000.0*t27-t95-19551.0/400.0*t29-3065979.0/80000.0*t31-1934177.0/120000.0*t23+249459.0/16000.0*t35+t105+2401.0/800.0*t20+t235+19551.0/400.0*t43+21609.0/640.0*t45+t448;
1090  double t498 = 7849849.0/4800000.0+700553.0/320000.0*xi-324919.0/32000.0*t3-t208+6266267.0/1440000.0*eta-556003.0/48000.0*t27-t69-136171.0/4800.0*t29-1900857.0/80000.0*t31-3142223.0/360000.0*t23+324919.0/32000.0*t35+t79+45619.0/38400.0*t20+t218+136171.0/4800.0*t43+55223.0/2560.0*t45+t436;
1091  double t500 = 46875.0/512.0*t3;
1092  double t501 = 390625.0/3072.0*t7;
1093  double t509 = 46875.0/512.0*t35;
1094  double t511 = 15625.0/1024.0*t20;
1095  double t512 = 390625.0/3072.0*t41;
1096  double t514 = 78125.0/6144.0*t21;
1097  double t515 = -15625.0/1024.0-53125.0/2048.0*xi+t500+t501-234375.0/1024.0*t43-209375.0/6144.0*eta+46875.0/512.0*t27+15625.0/512.0*t18+234375.0/1024.0*t29+71875.0/512.0*t31+71875.0/1536.0*t23-t509-46875.0/512.0*t37-t511-t512-234375.0/2048.0*t45-t514;
1098  double t518 = 390625.0/1536.0*t7;
1099  double t521 = 15625.0/128.0*t27;
1100  double t522 = 15625.0/384.0*t18;
1101  double t524 = 34375.0/256.0*t31;
1102  double t525 = 34375.0/768.0*t23;
1103  double t527 = 15625.0/128.0*t37;
1104  double t528 = 390625.0/1536.0*t41;
1105  double t529 = 78125.0/1024.0*t45;
1106  double t530 = 78125.0/3072.0+59375.0/1024.0*xi-78125.0/512.0*t3-t518+78125.0/256.0*t43+34375.0/768.0*eta-t521-t522-78125.0/256.0*t29-t524-t525+78125.0/512.0*t35+t527+t511+t528+t529;
1107  double t532 = 78125.0/512.0*t43;
1108  double t534 = 78125.0/512.0*t29;
1109  double t537 = 78125.0/3072.0*t21;
1110  double t538 = 15625.0/1024.0-9375.0/512.0*xi-t500+t532+59375.0/3072.0*eta-t521-t522-t534-t524-t525+t509+t527+78125.0/3072.0*t20+78125.0/512.0*t45+t537;
1111  double t540 = 15625.0/256.0*t27;
1112  double t541 = 15625.0/768.0*t18;
1113  double t542 = 34375.0/512.0*t31;
1114  double t543 = 34375.0/1536.0*t23;
1115  double t544 = 15625.0/128.0*t35;
1116  double t546 = 15625.0/128.0*t3;
1117  double t547 = 15625.0/256.0*t37;
1118  double t548 = -34375.0/1536.0*eta+t540+t541+t534+t542+t543-t544-15625.0/768.0-34375.0/512.0*xi+t546+t518-t547-t528-t532;
1119  double t554 = 15625.0/768.0*t20;
1120  double t557 = -3125.0/144.0*eta+t521+t522+78125.0/384.0*t29+3125.0/48.0*t31+3125.0/144.0*t23-t544-15625.0/768.0+9375.0/256.0*xi+t546-t527-t554-78125.0/384.0*t43-78125.0/768.0*t45;
1121  double t561 = 8125.0/1024.0*t18;
1122  double t568 = 15625.0/384.0*t27;
1123  double t569 = 6875.0/768.0*t18;
1124  double t576 = 15625.0/384.0*t37;
1125  double t578 = 34375.0/2304.0*eta-t568-t569-78125.0/768.0*t29+3125.0/128.0*t31-34375.0/2304.0*t23+15625.0/256.0*t35+6875.0/768.0-3125.0/128.0*xi-15625.0/256.0*t3+t576+78125.0/768.0*t43;
1126 
1127  values[0] = t26+t60;
1128  values[1] = t73+t86;
1129  values[2] = t99+t112;
1130  values[3] = t124+t136;
1131  values[4] = t145+t153;
1132  values[5] = t162+t170;
1133  values[6] = t179+t187;
1134  values[7] = t197+t204;
1135  values[8] = t222;
1136  values[9] = t239;
1137  values[10] = t256;
1138  values[11] = t272;
1139  values[12] = t288;
1140  values[13] = t304;
1141  values[14] = t312+t321;
1142  values[15] = t329+t337;
1143  values[16] = t345+t353;
1144  values[17] = t361+t369;
1145  values[18] = t376+t383;
1146  values[19] = t390+t396;
1147  values[20] = t403+t410;
1148  values[21] = t416+t423;
1149  values[22] = t437;
1150  values[23] = t449;
1151  values[24] = t462;
1152  values[25] = t474;
1153  values[26] = t486;
1154  values[27] = t498;
1155  values[28] = t515;
1156  values[29] = t530;
1157  values[30] = t538;
1158  values[31] = t548;
1159  values[32] = t557;
1160  values[33] = 3125.0/1024.0*eta+t540+t541+t542+t543+9375.0/1024.0*xi-t547-t554-t529-t537;
1161  values[34] = -t561-40625.0/1024.0*t31+46875.0/1024.0*t35+8125.0/1024.0+40625.0/1024.0*xi-46875.0/1024.0*t3-t501+t512;
1162  values[35] = t578;
1163  values[36] = -3125.0/384.0*eta-t568-t569-34375.0/768.0*t31+3125.0/384.0*t23-625.0/512.0-3125.0/512.0*xi+t576+15625.0/1536.0*t20+78125.0/1536.0*t45;
1164  values[37] = 3125.0/6144.0*eta-t561-40625.0/3072.0*t23+625.0/2048.0+15625.0/2048.0*t20+t514;
1165 }
1166 
1167 // values of the derivatives in xi-eta direction
1168 static void C_Q_M7_2D_DeriveXiEta(double xi, double eta, double *values)
1169 {
1170  double t1 = eta*eta;
1171  double t3 = t1*t1;
1172  double t7 = (1.0-t1)*(315.0/8.0*t3-105.0/4.0*t1+15.0/8.0);
1173  double t8 = 16807.0/207360.0*t7;
1174  double t9 = t3*eta;
1175  double t11 = t1*eta;
1176  double t15 = eta*(63.0/8.0*t9-35.0/4.0*t11+15.0/8.0*eta);
1177  double t16 = 16807.0/103680.0*t15;
1178  double t17 = xi*eta;
1179  double t19 = xi*xi;
1180  double t20 = t19*eta;
1181  double t23 = t19*t1;
1182  double t25 = xi*t11;
1183  double t27 = t19*t19;
1184  double t28 = t27*eta;
1185  double t29 = 85829.0/7680.0*t28;
1186  double t30 = t19*xi;
1187  double t31 = t30*t1;
1188  double t33 = t19*t11;
1189  double t35 = xi*t3;
1190  double t36 = 85829.0/7680.0*t35;
1191  double t37 = -802049.0/864000.0-t8+t16+225377.0/75000.0*t17+2902329.0/160000.0*t20-45493571.0/11520000.0*eta-282183.0/32000.0*t23-98177.0/48000.0*t25-t29-87887.0/3200.0*t31-87887.0/3200.0*t33-t36;
1192  double t39 = 184877.0/41472.0*t3;
1193  double t40 = 117649.0/30720.0*t9;
1194  double t41 = t30*eta;
1195  double t43 = t27*xi;
1196  double t48 = xi*(63.0/8.0*t43-35.0/4.0*t30+15.0/8.0*xi);
1197  double t49 = 16807.0/103680.0*t48;
1198  double t54 = (1.0-t19)*(315.0/8.0*t27-105.0/4.0*t19+15.0/8.0);
1199  double t55 = 16807.0/207360.0*t54;
1200  double t59 = 184877.0/41472.0*t27;
1201  double t60 = xi*t1;
1202  double t63 = 117649.0/30720.0*t43;
1203  double t64 = 971869.0/1728000.0*t1+t39-t40-98177.0/48000.0*t41+t49-t55+2529463.0/230400.0*t11+971869.0/1728000.0*t19+2529463.0/230400.0*t30+t59+2902329.0/160000.0*t60-45493571.0/11520000.0*xi-t63;
1204  double t66 = 117649.0/207360.0*t54;
1205  double t72 = 141169.0/25000.0*t17;
1206  double t77 = 420175.0/20736.0*t27;
1207  double t79 = 556003.0/32000.0*t23;
1208  double t81 = 117649.0/6144.0*t43;
1209  double t82 = 16807.0/768.0*t28;
1210  double t85 = 14651.0/2560.0*t35;
1211  double t86 = 117649.0/103680.0*t48;
1212  double t87 = -1900857.0/80000.0*t60-t77+45619.0/9600.0*t41+t79+324919.0/48000.0*t25+t81+t82+55223.0/1280.0*t31+136171.0/4800.0*t33+t85-t86;
1213  double t89 = 117649.0/69120.0*t54;
1214  double t95 = 2277569.0/200000.0*t17;
1215  double t100 = 16807.0/432.0*t27;
1216  double t102 = 206829.0/8000.0*t23;
1217  double t104 = 352947.0/10240.0*t43;
1218  double t105 = 16807.0/384.0*t28;
1219  double t108 = 14651.0/1280.0*t35;
1220  double t109 = 117649.0/34560.0*t48;
1221  double t110 = -3065979.0/80000.0*t60+t100+2401.0/200.0*t41+t102+83153.0/8000.0*t25-t104+t105+21609.0/320.0*t31+19551.0/400.0*t33+t108+t109;
1222  double t112 = 117649.0/41472.0*t54;
1223  double t118 = 12372059.0/320000.0*t17;
1224  double t123 = 1932805.0/41472.0*t27;
1225  double t125 = 558747.0/6400.0*t23;
1226  double t127 = 84035.0/512.0*t28;
1227  double t130 = 14651.0/2048.0*t35;
1228  double t131 = 117649.0/20736.0*t48;
1229  double t132 = 11668419.0/128000.0*t60-t123-69629.0/1280.0*t41-t125-309533.0/12800.0*t25+t81-t127-93639.0/512.0*t31-95697.0/1280.0*t33-t130-t131;
1230  double t139 = 2978171.0/80000.0*t17;
1231  double t145 = 23667.0/400.0*t23;
1232  double t149 = 1127.0/256.0*t35;
1233  double t150 = -862743.0/16000.0*t60+t123+40817.0/640.0*t41+t145+12201.0/1600.0*t25+t81+t127+7203.0/64.0*t31+4459.0/160.0*t33+t149+t131;
1234  double t157 = 13696039.0/1200000.0*t17;
1235  double t163 = 13377.0/4000.0*t23;
1236  double t165 = 7203.0/640.0*t31;
1237  double t167 = 14651.0/7680.0*t35;
1238  double t168 = 1083733.0/160000.0*t60-t100-31213.0/1600.0*t41-t163-144893.0/48000.0*t25-t104-t105-t165-18179.0/1600.0*t33-t167-t109;
1239  double t175 = 3234931.0/600000.0*t17;
1240  double t181 = 44247.0/32000.0*t23;
1241  double t185 = 637.0/1280.0*t35;
1242  double t186 = -169393.0/80000.0*t60+t77-26411.0/2400.0*t41+t181+1519.0/2000.0*t25+t81-t82+2401.0/640.0*t31+9261.0/3200.0*t33+t185+t86;
1243  double t188 = 20141003.0/4800000.0*t17;
1244  double t191 = 83771.0/16000.0*t23;
1245  double t195 = 85829.0/30720.0*t35;
1246  double t196 = 7123181.0/1728000.0*t1;
1247  double t197 = t8-t16-t188-940003.0/120000.0*t20+5917771.0/11520000.0*eta-t191+246511.0/64000.0*t25+t29-87887.0/12800.0*t31+77597.0/19200.0*t33-t195+t196;
1248  double t200 = 7123181.0/1728000.0*t19;
1249  double t204 = -273563.0/432000.0-t39+t40+594629.0/96000.0*t41-t49+t55-243521.0/76800.0*t11+t200+947497.0/230400.0*t30-t59+2126359.0/640000.0*t60-3953677.0/5760000.0*xi-t63;
1250  double t206 = 117649.0/103680.0*t15;
1251  double t207 = 117649.0/207360.0*t7;
1252  double t209 = 420175.0/20736.0*t3;
1253  double t212 = 117649.0/6144.0*t9;
1254  double t216 = 4969237.0/1200000.0*t17;
1255  double t221 = 268569.0/64000.0*t23;
1256  double t223 = 14651.0/2560.0*t28;
1257  double t226 = 16807.0/3072.0*t35;
1258  double t227 = t216-136171.0/57600.0*t30+7623959.0/960000.0*t20-6341041.0/960000.0*t60-41797.0/24000.0*t41+t221-2401.0/400.0*t25-t223+136171.0/19200.0*t31-55223.0/7680.0*t33+t226;
1259  double t229 = 117649.0/34560.0*t15;
1260  double t230 = 117649.0/69120.0*t7;
1261  double t232 = 16807.0/432.0*t3;
1262  double t235 = 352947.0/10240.0*t9;
1263  double t239 = 2564513.0/400000.0*t17;
1264  double t244 = 126567.0/16000.0*t23;
1265  double t246 = 14651.0/1280.0*t28;
1266  double t248 = 7203.0/640.0*t33;
1267  double t249 = 16807.0/1536.0*t35;
1268  double t250 = t239-6517.0/1600.0*t30+542577.0/40000.0*t20-1841567.0/160000.0*t60-34153.0/8000.0*t41+t244-26411.0/3200.0*t25-t246+19551.0/1600.0*t31-t248+t249;
1269  double t252 = 117649.0/20736.0*t15;
1270  double t253 = 117649.0/41472.0*t7;
1271  double t255 = 1932805.0/41472.0*t3;
1272  double t261 = 1454467.0/160000.0*t17;
1273  double t266 = 105987.0/25600.0*t23;
1274  double t268 = 14651.0/2048.0*t28;
1275  double t271 = 84035.0/2048.0*t35;
1276  double t272 = -t261+31899.0/5120.0*t30-771211.0/32000.0*t20+338541.0/8000.0*t60-11221.0/6400.0*t41-t266+21609.0/1280.0*t25+t268-95697.0/5120.0*t31+31213.0/1024.0*t33-t271;
1277  double t280 = 367059.0/320000.0*t17;
1278  double t285 = 9261.0/3200.0*t23;
1279  double t287 = 1127.0/256.0*t28;
1280  double t290 = t280-4459.0/1920.0*t30+50617.0/4000.0*t20-4756381.0/128000.0*t60-147.0/200.0*t41-t285-7203.0/2560.0*t25-t287+4459.0/640.0*t31-2401.0/128.0*t33+t271;
1281  double t298 = 335209.0/300000.0*t17;
1282  double t303 = 73059.0/32000.0*t23;
1283  double t305 = 14651.0/7680.0*t28;
1284  double t308 = t298+18179.0/19200.0*t30-211729.0/80000.0*t20+295323.0/40000.0*t60+9359.0/24000.0*t41-t303-2401.0/800.0*t25+t305-18179.0/6400.0*t31+2401.0/1280.0*t33-t249;
1285  double t316 = 371273.0/200000.0*t17;
1286  double t321 = 7889.0/16000.0*t23;
1287  double t323 = 637.0/1280.0*t28;
1288  double t326 = t316-3087.0/12800.0*t30+356083.0/480000.0*t20+506611.0/160000.0*t60-3479.0/32000.0*t41+t321-21609.0/6400.0*t25-t323+9261.0/12800.0*t31-2401.0/3840.0*t33-t226;
1289  double t333 = 85829.0/30720.0*t28;
1290  double t337 = -t8+t16-1343687.0/800000.0*t17-716051.0/960000.0*t20-496649.0/2880000.0*eta+220443.0/128000.0*t23+84457.0/32000.0*t25+t333+77597.0/76800.0*t31+77597.0/76800.0*t33+t195-10198837.0/3456000.0*t1;
1291  double t344 = 839329.0/3456000.0+t39+t40+84457.0/32000.0*t41+t49-t55-13827.0/6400.0*t11-10198837.0/3456000.0*t19-13827.0/6400.0*t30+t59-716051.0/960000.0*t60-496649.0/2880000.0*xi+t63;
1292  double t357 = 16807.0/3072.0*t28;
1293  double t360 = 356083.0/480000.0*t60-t77-21609.0/6400.0*t41+t321-3479.0/32000.0*t25-t81-t357-2401.0/3840.0*t31+9261.0/12800.0*t33-t185-t86;
1294  double t373 = 16807.0/1536.0*t28;
1295  double t376 = -211729.0/80000.0*t60+t100-2401.0/800.0*t41-t303+9359.0/24000.0*t25+t104-t373+2401.0/1280.0*t31-18179.0/6400.0*t33+t167+t109;
1296  double t389 = 84035.0/2048.0*t28;
1297  double t392 = 50617.0/4000.0*t60-t123-7203.0/2560.0*t41-t285-147.0/200.0*t25-t81+t389-2401.0/128.0*t31+4459.0/640.0*t33-t149-t131;
1298  double t407 = -771211.0/32000.0*t60+t123+21609.0/1280.0*t41-t266-11221.0/6400.0*t25-t81-t389+31213.0/1024.0*t31-95697.0/5120.0*t33+t130+t131;
1299  double t421 = 542577.0/40000.0*t60-t100-26411.0/3200.0*t41+t244-34153.0/8000.0*t25+t104+t373-t165+19551.0/1600.0*t33-t108-t109;
1300  double t436 = 7623959.0/960000.0*t60+t77-2401.0/400.0*t41+t221-41797.0/24000.0*t25-t81+t357-55223.0/7680.0*t31+136171.0/19200.0*t33-t85+t86;
1301  double t443 = t8-t16-t188+2126359.0/640000.0*t20-3953677.0/5760000.0*eta-t191+594629.0/96000.0*t25-t333+77597.0/19200.0*t31-87887.0/12800.0*t33+t36+t196;
1302  double t449 = -273563.0/432000.0-t39-t40+246511.0/64000.0*t41-t49+t55+947497.0/230400.0*t11+t200-243521.0/76800.0*t30-t59-940003.0/120000.0*t60+5917771.0/11520000.0*xi+t63;
1303  double t464 = 16807.0/768.0*t35;
1304  double t465 = t175-3087.0/3200.0*t30-169393.0/80000.0*t20+994357.0/80000.0*t60+1519.0/2000.0*t41+t181-26411.0/2400.0*t25+t323+9261.0/3200.0*t31+2401.0/640.0*t33-t464;
1305  double t479 = 16807.0/384.0*t35;
1306  double t480 = t157+18179.0/4800.0*t30+1083733.0/160000.0*t20+636951.0/20000.0*t60-144893.0/48000.0*t41-t163-31213.0/1600.0*t25-t305-18179.0/1600.0*t31-t248-t479;
1307  double t495 = 84035.0/512.0*t35;
1308  double t496 = -t139-4459.0/480.0*t30-862743.0/16000.0*t20-3850861.0/32000.0*t60+12201.0/1600.0*t41+t145+40817.0/640.0*t25+t287+4459.0/160.0*t31+7203.0/64.0*t33+t495;
1309  double t511 = t118+31899.0/1280.0*t30+11668419.0/128000.0*t20+140973.0/1000.0*t60-309533.0/12800.0*t41-t125-69629.0/1280.0*t25-t268-95697.0/1280.0*t31-93639.0/512.0*t33-t495;
1310  double t526 = -t95-6517.0/400.0*t30-3065979.0/80000.0*t20-1934177.0/40000.0*t60+83153.0/8000.0*t41+t102+2401.0/200.0*t25+t246+19551.0/400.0*t31+21609.0/320.0*t33+t479;
1311  double t541 = -t72-136171.0/14400.0*t30-1900857.0/80000.0*t20-3142223.0/120000.0*t60+324919.0/48000.0*t41+t79+45619.0/9600.0*t25+t223+136171.0/4800.0*t31+55223.0/1280.0*t33+t464;
1312  double t553 = 15625.0/256.0*t41;
1313  double t555 = 15625.0/256.0*t25;
1314  double t556 = 390625.0/6144.0*t28;
1315  double t558 = 390625.0/6144.0*t35;
1316  double t559 = -234375.0/1024.0*t31-15625.0/1024.0-209375.0/6144.0*xi-209375.0/6144.0*eta+46875.0/1024.0*t19+15625.0/256.0*t17+46875.0/1024.0*t1+78125.0/1024.0*t30+71875.0/512.0*t20+71875.0/512.0*t60+78125.0/1024.0*t11-t553-140625.0/1024.0*t23-t555-t556-234375.0/1024.0*t33-t558;
1317  double t563 = 15625.0/256.0*t19;
1318  double t564 = 15625.0/192.0*t17;
1319  double t565 = 15625.0/256.0*t1;
1320  double t567 = 34375.0/256.0*t20;
1321  double t568 = 34375.0/256.0*t60;
1322  double t569 = 78125.0/1536.0*t11;
1323  double t571 = 46875.0/256.0*t23;
1324  double t572 = 390625.0/3072.0*t28;
1325  double t573 = 78125.0/512.0*t33;
1326  double t574 = 78125.0/256.0*t31+15625.0/768.0+34375.0/768.0*xi+59375.0/3072.0*eta-t563-t564-t565-78125.0/768.0*t30-t567-t568-t569+78125.0/768.0*t41+t571+t555+t572+t573;
1327  double t575 = 78125.0/512.0*t31;
1328  double t578 = 78125.0/1536.0*t30;
1329  double t582 = 390625.0/3072.0*t35;
1330  double t583 = t575+15625.0/768.0+59375.0/3072.0*xi+34375.0/768.0*eta-t563-t564-t565-t578-t567-t568-78125.0/768.0*t11+t553+t571+78125.0/768.0*t25+78125.0/256.0*t33+t582;
1331  double t586 = 15625.0/512.0*t19;
1332  double t587 = 15625.0/384.0*t17;
1333  double t588 = 15625.0/512.0*t1;
1334  double t589 = 34375.0/512.0*t20;
1335  double t590 = 34375.0/512.0*t60;
1336  double t591 = 15625.0/192.0*t41;
1337  double t592 = 46875.0/512.0*t23;
1338  double t593 = -15625.0/1536.0-34375.0/1536.0*xi+3125.0/1024.0*eta+t586+t587+t588+t578+t589+t590-t591-t592-t572-t575;
1339  double t600 = 15625.0/192.0*t25;
1340  double t603 = -15625.0/768.0-3125.0/144.0*xi-3125.0/144.0*eta+t563+t564+t565+78125.0/1152.0*t30+3125.0/48.0*t20+3125.0/48.0*t60+78125.0/1152.0*t11-t591-t571-t600-78125.0/384.0*t31-78125.0/384.0*t33;
1341  double t606 = -15625.0/1536.0+3125.0/1024.0*xi-34375.0/1536.0*eta+t586+t587+t588+t589+t590+t569-t592-t600-t573-t582;
1342  double t608 = 8125.0/512.0*t17;
1343  double t614 = 15625.0/768.0*t19;
1344  double t615 = 6875.0/384.0*t17;
1345  double t616 = 15625.0/768.0*t1;
1346  double t621 = 15625.0/256.0*t23;
1347  double t623 = 15625.0/2304.0+34375.0/2304.0*xi-3125.0/384.0*eta-t614-t615-t616-78125.0/2304.0*t30+3125.0/128.0*t20-34375.0/768.0*t60+15625.0/384.0*t41+t621+78125.0/768.0*t31;
1348  double t631 = 15625.0/2304.0-3125.0/384.0*xi+34375.0/2304.0*eta-t614-t615-t616-34375.0/768.0*t20+3125.0/128.0*t60-78125.0/2304.0*t11+t621+15625.0/384.0*t25+78125.0/768.0*t33;
1349 
1350  values[0] = t37+t64;
1351  values[1] = 1526693.0/1728000.0+t66+6266267.0/1440000.0*eta-556003.0/96000.0*t1-136171.0/14400.0*t11+93337111.0/11520000.0*xi+16428671.0/1728000.0*t19-t72-1370971.0/46080.0*t30-3142223.0/120000.0*t20+t87;
1352  values[2] = 3412507.0/576000.0-t89+1767479.0/240000.0*eta-68943.0/8000.0*t1-6517.0/400.0*t11+5465313.0/1280000.0*xi-23683121.0/576000.0*t19-t95+64827.0/5120.0*t30-1934177.0/40000.0*t20+t110;
1353  values[3] = -139601.0/10800.0+t112-906549.0/64000.0*eta+186249.0/6400.0*t1+31899.0/1280.0*t11-3445729.0/144000.0*xi+23065721.0/345600.0*t19+t118+362551.0/9216.0*t30+140973.0/1000.0*t20+t132;
1354  values[4] = 1693391.0/172800.0-t112+1399097.0/192000.0*eta-7889.0/400.0*t1-4459.0/480.0*t11+50965439.0/2304000.0*xi-19824371.0/345600.0*t19-t139-545027.0/9216.0*t30-3850861.0/32000.0*t20+t150;
1355  values[5] = -1971907.0/576000.0+t89-446929.0/240000.0*eta+4459.0/4000.0*t1+18179.0/4800.0*t11-10357963.0/1280000.0*xi+19361321.0/576000.0*t19+t157+199283.0/5120.0*t30+636951.0/20000.0*t20+t168;
1356  values[6] = 2074807.0/1728000.0-t66+37387.0/160000.0*eta-14749.0/32000.0*t1-3087.0/3200.0*t11+22266139.0/11520000.0*xi-27233171.0/1728000.0*t19+t175-765919.0/46080.0*t30+994357.0/80000.0*t20+t186;
1357  values[7] = t197+t204;
1358  values[8] = 326879.0/216000.0+t206-t207-10863937.0/3840000.0*eta+t209-14426923.0/864000.0*t1+818741.0/46080.0*t11-t212+6379457.0/5760000.0*xi-89523.0/64000.0*t19+t227;
1359  values[9] = -1251607.0/576000.0-t229+t230+5099087.0/1280000.0*eta-t232+17200421.0/576000.0*t1-160867.0/5120.0*t11+t235+1582259.0/960000.0*xi-42189.0/16000.0*t19+t250;
1360  values[10] = 1911539.0/691200.0+t252-t253+1815107.0/1152000.0*eta+t255-25062667.0/691200.0*t1+26411.0/2304.0*t11-t212-1510229.0/256000.0*xi+35329.0/25600.0*t19+t272;
1361  values[11] = -612941.0/172800.0-t252+t253-19272239.0/2304000.0*eta-t255+13341671.0/345600.0*t1+256907.0/9216.0*t11-t212+3210137.0/768000.0*xi+3087.0/3200.0*t19+t290;
1362  values[12] = 1611757.0/576000.0+t229-t230+25791689.0/3840000.0*eta+t232-18280871.0/576000.0*t1-549829.0/15360.0*t11+t235-261709.0/960000.0*xi+24353.0/32000.0*t19+t308;
1363  values[13] = -1714657.0/1728000.0-t206+t207-5661313.0/3840000.0*eta-t209+26152721.0/1728000.0*t1+717899.0/46080.0*t11-t212+74431.0/1920000.0*xi-7889.0/48000.0*t19+t326;
1364  values[14] = t337+t344;
1365  values[15] = -1714657.0/1728000.0+t66+74431.0/1920000.0*eta-7889.0/48000.0*t1-3087.0/12800.0*t11-5661313.0/3840000.0*xi+26152721.0/1728000.0*t19+t316+717899.0/46080.0*t30+506611.0/160000.0*t20+t360;
1366  values[16] = 1611757.0/576000.0-t89-261709.0/960000.0*eta+24353.0/32000.0*t1+18179.0/19200.0*t11+25791689.0/3840000.0*xi-18280871.0/576000.0*t19+t298-549829.0/15360.0*t30+295323.0/40000.0*t20+t376;
1367  values[17] = -612941.0/172800.0+t112+3210137.0/768000.0*eta+3087.0/3200.0*t1-4459.0/1920.0*t11-19272239.0/2304000.0*xi+13341671.0/345600.0*t19+t280+256907.0/9216.0*t30-4756381.0/128000.0*t20+t392;
1368  values[18] = 1911539.0/691200.0-t112-1510229.0/256000.0*eta+35329.0/25600.0*t1+31899.0/5120.0*t11+1815107.0/1152000.0*xi-25062667.0/691200.0*t19-t261+26411.0/2304.0*t30+338541.0/8000.0*t20+t407;
1369  values[19] = -1251607.0/576000.0+t89+1582259.0/960000.0*eta-42189.0/16000.0*t1-6517.0/1600.0*t11+5099087.0/1280000.0*xi+17200421.0/576000.0*t19+t239-160867.0/5120.0*t30-1841567.0/160000.0*t20+t421;
1370  values[20] = 326879.0/216000.0-t66+6379457.0/5760000.0*eta-89523.0/64000.0*t1-136171.0/57600.0*t11-10863937.0/3840000.0*xi-14426923.0/864000.0*t19+t216+818741.0/46080.0*t30-6341041.0/960000.0*t20+t436;
1371  values[21] = t443+t449;
1372  values[22] = 2074807.0/1728000.0+t206-t207+22266139.0/11520000.0*eta+t209-27233171.0/1728000.0*t1-765919.0/46080.0*t11+t212+37387.0/160000.0*xi-14749.0/32000.0*t19+t465;
1373  values[23] = -1971907.0/576000.0-t229+t230-10357963.0/1280000.0*eta-t232+19361321.0/576000.0*t1+199283.0/5120.0*t11-t235-446929.0/240000.0*xi+4459.0/4000.0*t19+t480;
1374  values[24] = 1693391.0/172800.0+t252-t253+50965439.0/2304000.0*eta+t255-19824371.0/345600.0*t1-545027.0/9216.0*t11+t212+1399097.0/192000.0*xi-7889.0/400.0*t19+t496;
1375  values[25] = -139601.0/10800.0-t252+t253-3445729.0/144000.0*eta-t255+23065721.0/345600.0*t1+362551.0/9216.0*t11+t212-906549.0/64000.0*xi+186249.0/6400.0*t19+t511;
1376  values[26] = 3412507.0/576000.0+t229-t230+5465313.0/1280000.0*eta+t232-23683121.0/576000.0*t1+64827.0/5120.0*t11-t235+1767479.0/240000.0*xi-68943.0/8000.0*t19+t526;
1377  values[27] = 1526693.0/1728000.0-t206+t207+93337111.0/11520000.0*eta-t209+16428671.0/1728000.0*t1-1370971.0/46080.0*t11+t212+6266267.0/1440000.0*xi-556003.0/96000.0*t19+t541;
1378  values[28] = t559;
1379  values[29] = t574;
1380  values[30] = t583;
1381  values[31] = t593;
1382  values[32] = t603;
1383  values[33] = t606;
1384  values[34] = 3125.0/6144.0*eta-t608-40625.0/1024.0*t20+15625.0/512.0*t41+t556;
1385  values[35] = t623;
1386  values[36] = t631;
1387  values[37] = 3125.0/6144.0*xi-t608-40625.0/1024.0*t60+15625.0/512.0*t25+t558;
1388 }
1389 
1390 // values of the derivatives in eta-eta direction
1391 static void C_Q_M7_2D_DeriveEtaEta(double xi, double eta, double *values)
1392 {
1393  double t1 = eta*eta;
1394  double t3 = t1*eta;
1395  double t8 = 1.0+xi;
1396  double t9 = (1.0-t1)*(315.0/2.0*t3-105.0/2.0*eta)*t8;
1397  double t10 = 16807.0/207360.0*t9;
1398  double t13 = t1*t1;
1399  double t18 = eta*(315.0/8.0*t13-105.0/4.0*t1+15.0/8.0)*t8;
1400  double t19 = 16807.0/51840.0*t18;
1401  double t20 = xi*xi;
1402  double t22 = t20*xi;
1403  double t24 = t20*t20;
1404  double t26 = t24*xi;
1405  double t27 = 85829.0/38400.0*t26;
1406  double t29 = t13*eta;
1407  double t30 = 823543.0/15360.0*t29;
1408  double t31 = 12113881.0/207360.0*t3;
1409  double t32 = 117649.0/6144.0*t13;
1410  double t33 = -t10-45493571.0/11520000.0*xi-172582507.0/17280000.0*eta+t19+225377.0/150000.0*t20+967443.0/160000.0*t22-98177.0/192000.0*t24-t27-922751.0/384000.0*t1-t30+t31+t32;
1411  double t34 = xi*eta;
1412  double t36 = t20*eta;
1413  double t38 = xi*t1;
1414  double t40 = t22*eta;
1415  double t42 = t20*t1;
1416  double t44 = xi*t3;
1417  double t45 = 184877.0/10368.0*t44;
1418  double t46 = t24*eta;
1419  double t48 = t22*t1;
1420  double t50 = t20*t3;
1421  double t51 = 85829.0/3840.0*t50;
1422  double t52 = xi*t13;
1423  double t53 = 117649.0/6144.0*t52;
1424  double t58 = (63.0/8.0*t29-35.0/4.0*t3+15.0/8.0*eta)*t8;
1425  double t59 = 16807.0/103680.0*t58;
1426  double t60 = -49159793.0/57600000.0+971869.0/864000.0*t34+2902329.0/160000.0*t36+2529463.0/76800.0*t38-94061.0/16000.0*t40-98177.0/32000.0*t42+t45-87887.0/6400.0*t46-87887.0/3200.0*t48-t51-t53+t59;
1427  double t64 = 14651.0/1280.0*t3;
1428  double t66 = 141169.0/50000.0*t20;
1429  double t72 = 556003.0/48000.0*t40;
1430  double t74 = 16807.0/3840.0*t26;
1431  double t77 = 14651.0/1280.0*t50;
1432  double t78 = 7849849.0/4800000.0+700553.0/320000.0*eta-324919.0/32000.0*t1-t64+6266267.0/1440000.0*xi-t66-556003.0/48000.0*t34-3142223.0/360000.0*t22-1900857.0/80000.0*t36-136171.0/4800.0*t38+45619.0/38400.0*t24+t72+324919.0/32000.0*t42+t74+55223.0/2560.0*t46+136171.0/4800.0*t48+t77;
1433  double t81 = 14651.0/640.0*t3;
1434  double t83 = 2277569.0/400000.0*t20;
1435  double t89 = 68943.0/4000.0*t40;
1436  double t91 = 16807.0/1920.0*t26;
1437  double t94 = 14651.0/640.0*t50;
1438  double t95 = 1077069.0/400000.0+182427.0/40000.0*eta-249459.0/16000.0*t1-t81+1767479.0/240000.0*xi-t83-68943.0/4000.0*t34-1934177.0/120000.0*t22-3065979.0/80000.0*t36-19551.0/400.0*t38+2401.0/800.0*t24+t89+249459.0/16000.0*t42+t91+21609.0/640.0*t46+19551.0/400.0*t48+t94;
1439  double t98 = 14651.0/1024.0*t3;
1440  double t100 = 12372059.0/640000.0*t20;
1441  double t106 = 186249.0/3200.0*t40;
1442  double t108 = 16807.0/512.0*t26;
1443  double t111 = 14651.0/1024.0*t50;
1444  double t112 = -1834217.0/320000.0+4557.0/16000.0*eta+928599.0/25600.0*t1+t98-906549.0/64000.0*xi+t100+186249.0/3200.0*t34+46991.0/1000.0*t22+11668419.0/128000.0*t36+95697.0/1280.0*t38-69629.0/5120.0*t24-t106-928599.0/25600.0*t42-t108-93639.0/1024.0*t46-95697.0/1280.0*t48-t111;
1445  double t115 = 1127.0/128.0*t3;
1446  double t117 = 2978171.0/160000.0*t20;
1447  double t123 = 7889.0/200.0*t40;
1448  double t127 = 1127.0/128.0*t50;
1449  double t128 = 854217.0/320000.0-294.0/125.0*eta-36603.0/3200.0*t1-t115+1399097.0/192000.0*xi-t117-7889.0/200.0*t34-3850861.0/96000.0*t22-862743.0/16000.0*t36-4459.0/160.0*t38+40817.0/2560.0*t24+t123+36603.0/3200.0*t42+t108+7203.0/128.0*t46+4459.0/160.0*t48+t127;
1450  double t131 = 14651.0/3840.0*t3;
1451  double t133 = 13696039.0/2400000.0*t20;
1452  double t139 = 4459.0/2000.0*t40;
1453  double t141 = 7203.0/1280.0*t46;
1454  double t143 = 14651.0/3840.0*t50;
1455  double t144 = -497791.0/600000.0-91679.0/80000.0*eta+144893.0/32000.0*t1+t131-446929.0/240000.0*xi+t133+4459.0/2000.0*t34+212317.0/20000.0*t22+1083733.0/160000.0*t36+18179.0/1600.0*t38-31213.0/6400.0*t24-t139-144893.0/32000.0*t42-t91-t141-18179.0/1600.0*t48-t143;
1456  double t147 = 637.0/640.0*t3;
1457  double t149 = 3234931.0/1200000.0*t20;
1458  double t155 = 14749.0/16000.0*t40;
1459  double t159 = 637.0/640.0*t50;
1460  double t160 = 5537.0/100000.0+38661.0/160000.0*eta-4557.0/4000.0*t1-t147+37387.0/160000.0*xi+t149-14749.0/16000.0*t34+994357.0/240000.0*t22-169393.0/80000.0*t36-9261.0/3200.0*t38-26411.0/9600.0*t24+t155+4557.0/4000.0*t42-t74+2401.0/1280.0*t46+9261.0/3200.0*t48+t159;
1461  double t163 = 5077697.0/414720.0*t3;
1462  double t165 = 20141003.0/9600000.0*t20;
1463  double t166 = 7123181.0/864000.0*t34;
1464  double t170 = 39583543.0/57600000.0+21015133.0/4320000.0*eta+t32-2159737.0/192000.0*t1-t163+5917771.0/11520000.0*xi-t165+t166-940003.0/360000.0*t22+2126359.0/640000.0*t36-243521.0/25600.0*t38;
1465  double t172 = 83771.0/24000.0*t40;
1466  double t176 = 85829.0/15360.0*t50;
1467  double t177 = 594629.0/384000.0*t24-t172+739533.0/128000.0*t42-t45+t27-87887.0/25600.0*t46+77597.0/19200.0*t48-t176+t53-t59-t19+t10;
1468  double t180 = 588245.0/6144.0*t13;
1469  double t182 = 2907611.0/41472.0*t3;
1470  double t184 = 4969237.0/2400000.0*t20;
1471  double t189 = -170592863.0/57600000.0-59484089.0/2160000.0*eta-t180+4232963.0/76800.0*t1+t182-10863937.0/3840000.0*xi+t184-14426923.0/432000.0*t34+7623959.0/2880000.0*t22-6341041.0/960000.0*t36+818741.0/15360.0*t38;
1472  double t191 = 89523.0/32000.0*t40;
1473  double t193 = 420175.0/5184.0*t44;
1474  double t194 = 14651.0/12800.0*t26;
1475  double t197 = 16807.0/1536.0*t50;
1476  double t198 = 588245.0/6144.0*t52;
1477  double t199 = 117649.0/103680.0*t58;
1478  double t200 = 117649.0/51840.0*t18;
1479  double t201 = 117649.0/207360.0*t9;
1480  double t202 = -41797.0/96000.0*t24+t191-7203.0/800.0*t42+t193-t194+136171.0/38400.0*t46-55223.0/7680.0*t48+t197-t198+t199+t200-t201;
1481  double t205 = 352947.0/2048.0*t13;
1482  double t207 = 1226911.0/6912.0*t3;
1483  double t209 = 2564513.0/800000.0*t20;
1484  double t214 = 26096371.0/6400000.0+50686139.0/720000.0*eta+t205-2384193.0/25600.0*t1-t207+5099087.0/1280000.0*xi+t209+17200421.0/288000.0*t34+180859.0/40000.0*t22-1841567.0/160000.0*t36-482601.0/5120.0*t38;
1485  double t216 = 42189.0/8000.0*t40;
1486  double t218 = 16807.0/108.0*t44;
1487  double t219 = 14651.0/6400.0*t26;
1488  double t221 = 7203.0/640.0*t48;
1489  double t222 = 16807.0/768.0*t50;
1490  double t223 = 352947.0/2048.0*t52;
1491  double t224 = 117649.0/34560.0*t58;
1492  double t225 = 117649.0/17280.0*t18;
1493  double t226 = 117649.0/69120.0*t9;
1494  double t227 = -34153.0/32000.0*t24+t216-79233.0/6400.0*t42-t218-t219+19551.0/3200.0*t46-t221+t222+t223-t224-t225+t226;
1495  double t231 = 22269275.0/82944.0*t3;
1496  double t233 = 1454467.0/320000.0*t20;
1497  double t238 = -501613.0/11520000.0-748234949.0/6912000.0*eta-t180+607453.0/15360.0*t1+t231+1815107.0/1152000.0*xi-t233-25062667.0/345600.0*t34-771211.0/96000.0*t22+338541.0/8000.0*t36+26411.0/768.0*t38;
1498  double t240 = 35329.0/12800.0*t40;
1499  double t242 = 1932805.0/10368.0*t44;
1500  double t243 = 14651.0/10240.0*t26;
1501  double t246 = 84035.0/1024.0*t50;
1502  double t247 = 117649.0/20736.0*t58;
1503  double t248 = 117649.0/10368.0*t18;
1504  double t249 = 117649.0/41472.0*t9;
1505  double t250 = -11221.0/25600.0*t24-t240+64827.0/2560.0*t42+t242+t243-95697.0/10240.0*t46+31213.0/1024.0*t48-t246-t198+t247+t248-t249;
1506  double t255 = 367059.0/640000.0*t20;
1507  double t260 = -62402137.0/11520000.0+376548487.0/3456000.0*eta-t180+530621.0/7680.0*t1-t231-19272239.0/2304000.0*xi+t255+13341671.0/172800.0*t34+50617.0/12000.0*t22-4756381.0/128000.0*t36+256907.0/3072.0*t38;
1508  double t262 = 3087.0/1600.0*t40;
1509  double t264 = 1127.0/1280.0*t26;
1510  double t267 = -147.0/800.0*t24-t262-21609.0/5120.0*t42-t242-t264+4459.0/1280.0*t46-2401.0/128.0*t48+t246-t198-t247-t248+t249;
1511  double t272 = 335209.0/600000.0*t20;
1512  double t277 = 35582379.0/6400000.0-204365231.0/2880000.0*eta+t205-2585877.0/25600.0*t1+t207+25791689.0/3840000.0*xi+t272-18280871.0/288000.0*t34-211729.0/240000.0*t22+295323.0/40000.0*t36-549829.0/5120.0*t38;
1513  double t279 = 24353.0/16000.0*t40;
1514  double t281 = 14651.0/38400.0*t26;
1515  double t284 = 9359.0/96000.0*t24-t279-7203.0/1600.0*t42+t218+t281-18179.0/12800.0*t46+2401.0/1280.0*t48-t222+t223+t224+t225-t226;
1516  double t289 = 371273.0/400000.0*t20;
1517  double t294 = -128307137.0/57600000.0+467769337.0/17280000.0*eta-t180+3930437.0/76800.0*t1-t182-5661313.0/3840000.0*xi+t289+26152721.0/864000.0*t34+356083.0/1440000.0*t22+506611.0/160000.0*t36+717899.0/15360.0*t38;
1518  double t296 = 7889.0/24000.0*t40;
1519  double t298 = 637.0/6400.0*t26;
1520  double t301 = -3479.0/128000.0*t24+t296-64827.0/12800.0*t42-t193-t298+9261.0/25600.0*t46-2401.0/3840.0*t48-t197-t198-t199-t200+t201;
1521  double t311 = 18298957.0/57600000.0-155978489.0/34560000.0*eta+t32-3621101.0/384000.0*t1+t163-496649.0/2880000.0*xi-1343687.0/1600000.0*t20-10198837.0/1728000.0*t34-716051.0/2880000.0*t22-716051.0/960000.0*t36-41481.0/6400.0*t38;
1522  double t315 = 85829.0/153600.0*t26;
1523  double t318 = 84457.0/128000.0*t24+73481.0/64000.0*t40+253371.0/64000.0*t42+t45+t315+77597.0/153600.0*t46+77597.0/76800.0*t48+t176+t53+t59+t19-t10;
1524  double t329 = 16807.0/15360.0*t26;
1525  double t332 = -269059.0/3200000.0-137347.0/320000.0*eta+10437.0/64000.0*t1+t147+74431.0/1920000.0*xi+t289-7889.0/24000.0*t34+506611.0/480000.0*t22+356083.0/480000.0*t36-9261.0/12800.0*t38-21609.0/25600.0*t24+t296-10437.0/64000.0*t42-t329-2401.0/7680.0*t46+9261.0/12800.0*t48-t159;
1526  double t342 = 16807.0/7680.0*t26;
1527  double t345 = 229957.0/1200000.0+546791.0/320000.0*eta-9359.0/16000.0*t1-t131-261709.0/960000.0*xi+t272+24353.0/16000.0*t34+98441.0/40000.0*t22-211729.0/80000.0*t36+18179.0/6400.0*t38-2401.0/3200.0*t24-t279+9359.0/16000.0*t42-t342+2401.0/2560.0*t46-18179.0/6400.0*t48+t143;
1528  double t355 = 16807.0/2048.0*t26;
1529  double t358 = 166257.0/1280000.0-104811.0/32000.0*eta+441.0/400.0*t1+t115+3210137.0/768000.0*xi+t255+3087.0/1600.0*t34-4756381.0/384000.0*t22+50617.0/4000.0*t36-4459.0/640.0*t38-7203.0/10240.0*t24-t262-441.0/400.0*t42+t355-2401.0/256.0*t46+4459.0/640.0*t48-t127;
1530  double t370 = 207809.0/640000.0+2268063.0/256000.0*eta+33663.0/12800.0*t1-t98-1510229.0/256000.0*xi-t233+35329.0/12800.0*t34+112847.0/8000.0*t22-771211.0/32000.0*t36+95697.0/5120.0*t38+21609.0/5120.0*t24-t240-33663.0/12800.0*t42-t355+31213.0/2048.0*t46-95697.0/5120.0*t48+t111;
1531  double t381 = -1827651.0/1600000.0-1269933.0/160000.0*eta+102459.0/16000.0*t1+t81+1582259.0/960000.0*xi+t209-42189.0/8000.0*t34-1841567.0/480000.0*t22+542577.0/40000.0*t36-19551.0/1600.0*t38-26411.0/12800.0*t24+t216-102459.0/16000.0*t42+t342-t141+19551.0/1600.0*t48-t94;
1532  double t393 = -1367737.0/2400000.0-2781681.0/640000.0*eta+41797.0/16000.0*t1+t64+6379457.0/5760000.0*xi+t184-89523.0/32000.0*t34-6341041.0/2880000.0*t22+7623959.0/960000.0*t36-136171.0/19200.0*t38-2401.0/1600.0*t24+t191-41797.0/16000.0*t42+t329-55223.0/15360.0*t46+136171.0/19200.0*t48-t77;
1533  double t399 = t10-3953677.0/5760000.0*xi+196867657.0/17280000.0*eta-t19-t165+2126359.0/1920000.0*t22+246511.0/256000.0*t24-t315-5668649.0/384000.0*t1+t30-t31+t32;
1534  double t405 = 36656459.0/28800000.0+t166-940003.0/120000.0*t36+947497.0/76800.0*t38-t172+594629.0/64000.0*t42-t45+77597.0/38400.0*t46-87887.0/12800.0*t48+t51-t53-t59;
1535  double t412 = 5764801.0/15360.0*t29;
1536  double t413 = 14101073.0/41472.0*t3;
1537  double t414 = -242622863.0/57600000.0-t201+22266139.0/11520000.0*xi-961071937.0/17280000.0*eta+t200+t149-169393.0/240000.0*t22+1519.0/8000.0*t24+t298+4809203.0/76800.0*t1-t412+t413;
1538  double t421 = 16807.0/384.0*t50;
1539  double t422 = -t180-27233171.0/864000.0*t34+994357.0/80000.0*t36-765919.0/15360.0*t38+t155-26411.0/1600.0*t42+t193+9261.0/6400.0*t46+2401.0/640.0*t48-t421+t198+t199;
1540  double t429 = 5764801.0/5120.0*t29;
1541  double t430 = 12285917.0/13824.0*t3;
1542  double t431 = 8088871.0/6400000.0+t226-10357963.0/1280000.0*xi+724916437.0/5760000.0*eta-t225+t133+1083733.0/480000.0*t22-144893.0/192000.0*t24-t281-1952013.0/25600.0*t1+t429-t430;
1543  double t437 = 16807.0/192.0*t50;
1544  double t438 = t205+19361321.0/288000.0*t34+636951.0/20000.0*t36+597849.0/5120.0*t38-t139-93639.0/3200.0*t42-t218-18179.0/3200.0*t46-t221-t437-t223-t224;
1545  double t445 = 5764801.0/3072.0*t29;
1546  double t446 = 56891695.0/41472.0*t3;
1547  double t447 = 134554637.0/11520000.0-t249+50965439.0/2304000.0*xi-622119337.0/3456000.0*eta+t248-t117-287581.0/16000.0*t22+12201.0/6400.0*t24+t264-472997.0/15360.0*t1-t445+t446;
1548  double t454 = 84035.0/256.0*t50;
1549  double t455 = -t180-19824371.0/172800.0*t34-3850861.0/32000.0*t36-545027.0/3072.0*t38+t123+122451.0/1280.0*t42+t242+4459.0/320.0*t46+7203.0/64.0*t48+t454+t198+t247;
1550  double t462 = -52741003.0/2880000.0+t249-3445729.0/144000.0*xi+631843387.0/3456000.0*eta-t248+t100+3889473.0/128000.0*t22-309533.0/51200.0*t24-t243+2249737.0/15360.0*t1+t445-t446;
1551  double t469 = -t180+23065721.0/172800.0*t34+140973.0/1000.0*t36+362551.0/3072.0*t38-t106-208887.0/2560.0*t42-t242-95697.0/2560.0*t46-93639.0/512.0*t48-t454+t198-t247;
1552  double t476 = 59592379.0/6400000.0-t226+5465313.0/1280000.0*xi-737881837.0/5760000.0*eta+t225-t83-1021993.0/80000.0*t22+83153.0/32000.0*t24+t219-3162117.0/25600.0*t1-t429+t430;
1553  double t483 = t205-23683121.0/288000.0*t34-1934177.0/40000.0*t36+194481.0/5120.0*t38+t89+7203.0/400.0*t42+t218+19551.0/800.0*t46+21609.0/320.0*t48+t437-t223+t224;
1554  double t490 = -11258387.0/57600000.0+t201+93337111.0/11520000.0*xi+928658437.0/17280000.0*eta-t200-t66-633619.0/80000.0*t22+324919.0/192000.0*t24+t194+2994047.0/76800.0*t1+t412-t413;
1555  double t497 = -t180+16428671.0/864000.0*t34-3142223.0/120000.0*t36-1370971.0/15360.0*t38+t72+45619.0/6400.0*t42-t193+136171.0/9600.0*t46+55223.0/1280.0*t48+t421+t198-t199;
1556  double t499 = 390625.0/3072.0*t3;
1557  double t501 = 46875.0/512.0*t1;
1558  double t509 = 15625.0/1024.0*t24;
1559  double t511 = 46875.0/512.0*t42;
1560  double t512 = 78125.0/6144.0*t26;
1561  double t514 = 390625.0/3072.0*t50;
1562  double t515 = -15625.0/1024.0+t499-53125.0/2048.0*eta+t501-234375.0/2048.0*t46-209375.0/6144.0*xi+15625.0/512.0*t20+46875.0/512.0*t34+71875.0/1536.0*t22+71875.0/512.0*t36+234375.0/1024.0*t38-t509-46875.0/512.0*t40-t511-t512-234375.0/1024.0*t48-t514;
1563  double t519 = 15625.0/384.0*t20;
1564  double t520 = 15625.0/128.0*t34;
1565  double t521 = 34375.0/768.0*t22;
1566  double t522 = 34375.0/256.0*t36;
1567  double t523 = 78125.0/512.0*t38;
1568  double t525 = 15625.0/128.0*t40;
1569  double t526 = 78125.0/3072.0*t26;
1570  double t527 = 78125.0/512.0*t48;
1571  double t528 = 15625.0/1024.0-9375.0/512.0*eta-t501+78125.0/512.0*t46+59375.0/3072.0*xi-t519-t520-t521-t522-t523+78125.0/3072.0*t24+t525+t511+t526+t527;
1572  double t529 = 390625.0/1536.0*t3;
1573  double t532 = 78125.0/1024.0*t46;
1574  double t537 = 390625.0/1536.0*t50;
1575  double t538 = 78125.0/3072.0-t529+59375.0/1024.0*eta-78125.0/512.0*t1+t532+34375.0/768.0*xi-t519-t520-t521-t522-78125.0/256.0*t38+t509+t525+78125.0/512.0*t42+78125.0/256.0*t48+t537;
1576  double t540 = 15625.0/768.0*t20;
1577  double t541 = 15625.0/256.0*t34;
1578  double t542 = 34375.0/1536.0*t22;
1579  double t543 = 34375.0/512.0*t36;
1580  double t544 = 15625.0/768.0*t24;
1581  double t546 = 15625.0/256.0*t40;
1582  double t553 = 15625.0/128.0*t1;
1583  double t554 = 15625.0/128.0*t42;
1584  double t557 = -3125.0/144.0*xi+t519+t520+3125.0/144.0*t22+3125.0/48.0*t36+78125.0/384.0*t38-t544-15625.0/768.0+9375.0/256.0*eta+t553-t525-t554-78125.0/768.0*t46-78125.0/384.0*t48;
1585  double t560 = -34375.0/1536.0*xi+t540+t541+t542+t543+t523-15625.0/768.0-34375.0/512.0*eta+t553+t529-t546-t554-t527-t537;
1586  double t562 = 8125.0/1024.0*t20;
1587  double t567 = 6875.0/768.0*t20;
1588  double t568 = 15625.0/384.0*t34;
1589  double t573 = 15625.0/384.0*t40;
1590  double t584 = 34375.0/2304.0*xi-t567-t568-34375.0/2304.0*t22+3125.0/128.0*t36-78125.0/768.0*t38+6875.0/768.0-3125.0/128.0*eta-15625.0/256.0*t1+t573+15625.0/256.0*t42+78125.0/768.0*t48;
1591 
1592  values[0] = t33+t60;
1593  values[1] = t78;
1594  values[2] = t95;
1595  values[3] = t112;
1596  values[4] = t128;
1597  values[5] = t144;
1598  values[6] = t160;
1599  values[7] = t170+t177;
1600  values[8] = t189+t202;
1601  values[9] = t214+t227;
1602  values[10] = t238+t250;
1603  values[11] = t260+t267;
1604  values[12] = t277+t284;
1605  values[13] = t294+t301;
1606  values[14] = t311+t318;
1607  values[15] = t332;
1608  values[16] = t345;
1609  values[17] = t358;
1610  values[18] = t370;
1611  values[19] = t381;
1612  values[20] = t393;
1613  values[21] = t399+t405;
1614  values[22] = t414+t422;
1615  values[23] = t431+t438;
1616  values[24] = t447+t455;
1617  values[25] = t462+t469;
1618  values[26] = t476+t483;
1619  values[27] = t490+t497;
1620  values[28] = t515;
1621  values[29] = t528;
1622  values[30] = t538;
1623  values[31] = 3125.0/1024.0*xi+t540+t541+t542+t543-t544+9375.0/1024.0*eta-t546-t526-t532;
1624  values[32] = t557;
1625  values[33] = t560;
1626  values[34] = 3125.0/6144.0*xi-t562-40625.0/3072.0*t22+15625.0/2048.0*t24+625.0/2048.0+t512;
1627  values[35] = -3125.0/384.0*xi-t567-t568+3125.0/384.0*t22-34375.0/768.0*t36+15625.0/1536.0*t24-625.0/512.0-3125.0/512.0*eta+t573+78125.0/1536.0*t46;
1628  values[36] = t584;
1629  values[37] = -t562-40625.0/1024.0*t36+8125.0/1024.0+40625.0/1024.0*eta-46875.0/1024.0*t1-t499+46875.0/1024.0*t42+t514;
1630 }
1631 
1632 // ***********************************************************************
1633 
1634 TBaseFunct2D *BF_C_Q_M7_2D_Obj = new TBaseFunct2D
1635  (38, BF_C_Q_M7_2D, BFUnitSquare,
1636  C_Q_M7_2D_Funct, C_Q_M7_2D_DeriveXi,
1637  C_Q_M7_2D_DeriveEta, C_Q_M7_2D_DeriveXiXi,
1638  C_Q_M7_2D_DeriveXiEta, C_Q_M7_2D_DeriveEtaEta, 7, 7,
1639  0, NULL);
Definition: BaseFunct2D.h:27