aboutsummaryrefslogtreecommitdiff
path: root/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxy.compute.c
blob: c0d06300786a806ff773e19e208846c52406ae43 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
fp t1387;
fp t1428;
fp t1430;
fp t1390;
fp t1485;
fp t1396;
fp t1484;
fp t1397;
fp t1388;
fp t1483;
fp t1426;
fp t1403;
fp t1482;
fp t1458;
fp t1427;
fp t1429;
fp t1389;
fp t1361;
fp t1425;
fp t1481;
fp t1343;
fp t1460;
fp t1480;
fp t1450;
fp t1418;
fp t1398;
fp t1395;
fp t1344;
fp t1479;
fp t1448;
fp t1414;
fp t1417;
fp t1423;
fp t1369;
fp t1401;
fp t1478;
fp t1459;
fp t1416;
fp t1422;
fp t1421;
fp t1402;
fp t1477;
fp t1366;
fp t1367;
fp t1413;
fp t1452;
fp t1476;
fp t1391;
fp t1475;
fp t1350;
fp t1474;
fp t1356;
fp t1400;
fp t1473;
fp t1373;
fp t1472;
fp t1362;
fp t1471;
fp t1351;
fp t1380;
fp t1470;
fp t1455;
fp t1469;
fp t1399;
fp t1355;
fp t1468;
fp t1386;
fp t1415;
fp t1467;
fp t1466;
fp t1449;
fp t1347;
fp t1465;
fp t1461;
fp t1464;
fp t1412;
fp t1463;
fp t1411;
fp t1447;
fp t1462;
fp t1457;
fp t1456;
fp t1454;
fp t1453;
fp t1393;
fp t1451;
fp t1446;
fp t1445;
fp t1409;
fp t1420;
fp t1444;
fp t1394;
fp t1443;
fp t1442;
fp t1441;
fp t1440;
fp t1439;
fp t1438;
fp t1437;
fp t1436;
fp t1435;
fp t1434;
fp t1433;
fp t1424;
fp t1419;
fp t1410;
fp t1408;
fp t1407;
fp t1406;
fp t1405;
fp t1404;
fp t1392;
      t1387 = RATIONAL(1.0,144.0);
      t1428 = y*y;
      t1430 = y*t1428;
      t1390 = RATIONAL(-1.0,72.0);
      t1485 = t1387*y+t1390*t1430;
      t1396 = RATIONAL(-1.0,48.0);
      t1484 = t1396*t1428+t1387;
      t1397 = RATIONAL(1.0,48.0);
      t1388 = RATIONAL(-1.0,144.0);
      t1483 = t1397*t1428+t1388;
      t1426 = RATIONAL(-8.0,9.0);
      t1403 = RATIONAL(4.0,9.0);
      t1482 = t1426*y+t1403*t1430;
      t1458 = x*t1430;
      t1427 = x*x;
      t1429 = t1427*x;
      t1389 = RATIONAL(1.0,18.0);
      t1361 = t1389*t1429;
      t1425 = RATIONAL(2.0,9.0);
      t1481 = y*t1361+t1425*t1458;
      t1343 = t1389*t1430;
      t1460 = y*t1429;
      t1480 = x*t1343+t1425*t1460;
      t1450 = t1428*t1429;
      t1418 = RATIONAL(-1.0,9.0);
      t1398 = RATIONAL(1.0,12.0);
      t1395 = RATIONAL(-1.0,24.0);
      t1344 = t1395*t1428;
      t1479 = t1398*t1450+t1418*t1429+(t1389+t1344)*x;
      t1448 = t1427*t1430;
      t1414 = RATIONAL(1.0,6.0);
      t1417 = RATIONAL(1.0,9.0);
      t1423 = RATIONAL(-1.0,3.0);
      t1369 = t1423*t1427;
      t1401 = RATIONAL(-1.0,18.0);
      t1478 = t1401*t1430+t1414*t1448+(t1369+t1417)*y;
      t1459 = x*t1428;
      t1416 = RATIONAL(-2.0,3.0);
      t1422 = RATIONAL(8.0,9.0);
      t1421 = RATIONAL(1.0,3.0);
      t1402 = RATIONAL(-4.0,9.0);
      t1477 = t1416*t1459+t1421*t1450+t1422*x+t1402*t1429;
      t1366 = t1421*t1427;
      t1367 = t1421*t1428;
      t1413 = RATIONAL(-1.0,4.0);
      t1452 = t1413*t1427;
      t1476 = t1366+t1367+t1428*t1452+t1402;
      t1391 = RATIONAL(1.0,72.0);
      t1475 = t1391*t1429+t1395*t1450+t1483*x;
      t1350 = t1397*t1427;
      t1474 = t1391*t1430+t1395*t1448+(t1350+t1388)*y;
      t1356 = t1396*t1427;
      t1400 = RATIONAL(1.0,24.0);
      t1473 = t1400*t1448+y*t1356+t1485;
      t1373 = t1423*t1428;
      t1472 = t1401*t1429+t1414*t1450+(t1373+t1417)*x;
      t1362 = t1390*t1429;
      t1471 = t1362+t1400*t1450+t1484*x;
      t1351 = t1395*t1427;
      t1380 = t1418*t1430;
      t1470 = t1398*t1448+t1380+(t1351+t1389)*y;
      t1455 = t1429*t1430;
      t1469 = y*t1362+RATIONAL(1.0,36.0)*t1455+t1485*x;
      t1399 = RATIONAL(-1.0,12.0);
      t1355 = t1400*t1427;
      t1468 = t1417*t1430+t1399*t1448+(t1355+t1401)*y;
      t1386 = t1418*x;
      t1415 = RATIONAL(-1.0,6.0);
      t1467 = t1386+t1361+t1415*t1450+x*t1367;
      t1466 = t1343+t1415*t1448+(t1366+t1418)*y;
      t1449 = t1399*t1429;
      t1347 = t1400*t1428;
      t1465 = t1417*t1429+t1428*t1449+(t1401+t1347)*x;
      t1461 = x*y;
      t1464 = RATIONAL(16.0,9.0)*t1461+t1426*t1458+t1482*t1429;
      t1412 = RATIONAL(2.0,3.0);
      t1463 = t1412*t1459+t1423*t1450+t1426*x+t1403*t1429;
      t1411 = RATIONAL(1.0,4.0);
      t1447 = t1411*t1428;
      t1462 = t1427*t1447+t1369+t1373+t1403;
      t1457 = y*t1427;
      t1456 = t1429*t1380+y*t1386;
      t1454 = t1428*t1427;
      t1453 = t1416*t1455+RATIONAL(-10.0,3.0)*t1461;
      t1393 = RATIONAL(5.0,24.0);
      t1451 = t1393*t1461+t1414*t1455;
      t1446 = RATIONAL(1.0,8.0)*t1454+t1389+t1456;
      t1445 = t1401+RATIONAL(-1.0,8.0)*t1454+t1456;
      t1409 = RATIONAL(4.0,3.0);
      t1420 = RATIONAL(5.0,3.0);
      t1444 = t1420*t1458+t1409*t1460+t1453;
      t1394 = RATIONAL(-5.0,12.0);
      t1443 = t1394*t1458+y*t1449+t1451;
      t1442 = t1420*t1460+t1409*t1458+t1453;
      t1441 = t1399*t1458+t1394*t1460+t1451;
      t1440 = t1350+RATIONAL(-1.0,16.0)*t1454+t1469+t1483;
      t1439 = RATIONAL(1.0,16.0)*t1454+t1356+t1469+t1484;
      t1438 = t1351+t1415*t1428+t1446+t1481;
      t1437 = t1344+t1415*t1427+t1446+t1480;
      t1436 = t1347+t1414*t1427+t1445+t1480;
      t1435 = t1402*t1430+t1422*y+t1416*t1457+t1421*t1448+t1464;
      t1434 = t1414*t1428+t1355+t1445+t1481;
      t1433 = t1423*t1448+t1412*t1457+t1464+t1482;
      t1424 = RATIONAL(5.0,4.0);
      t1419 = RATIONAL(-5.0,3.0);
      t1410 = RATIONAL(1.0,2.0);
      t1408 = RATIONAL(-5.0,8.0);
      t1407 = RATIONAL(-5.0,4.0);
      t1406 = RATIONAL(-5.0,2.0);
      t1405 = RATIONAL(-1.0,2.0);
      t1404 = RATIONAL(5.0,8.0);
      t1392 = RATIONAL(-5.0,24.0);
      coeffs_dxy->coeff_m2_m2 = t1439+t1474+t1475;
      coeffs_dxy->coeff_m1_m2 = t1434+t1470+t1472;
      coeffs_dxy->coeff_0_m2 = (t1413*t1428+t1398)*t1429+(t1404*t1428+t1392)*x+
t1443;
      coeffs_dxy->coeff_p1_m2 = t1438+t1468+t1472;
      coeffs_dxy->coeff_p2_m2 = t1440+t1473+t1475;
      coeffs_dxy->coeff_m2_m1 = t1436+t1478+t1479;
      coeffs_dxy->coeff_m1_m1 = t1433+t1462+t1463;
      coeffs_dxy->coeff_0_m1 = (t1410*t1428+t1416)*t1429+(t1420+t1407*t1428)*x+
t1444;
      coeffs_dxy->coeff_p1_m1 = t1435+t1463+t1476;
      coeffs_dxy->coeff_p2_m1 = t1437+t1466+t1479;
      coeffs_dxy->coeff_m2_0 = (t1452+t1398)*t1430+(t1392+t1404*t1427)*y+t1441;
      coeffs_dxy->coeff_m1_0 = (t1410*t1427+t1416)*t1430+(t1407*t1427+t1420)*y+
t1442;
      coeffs_dxy->coeff_0_0 = (t1406*y+t1430)*t1429+(RATIONAL(25.0,4.0)*y+t1406
*t1430)*x;
      coeffs_dxy->coeff_p1_0 = (t1405*t1427+t1412)*t1430+(t1424*t1427+t1419)*y+
t1442;
      coeffs_dxy->coeff_p2_0 = (t1399+t1411*t1427)*t1430+(t1393+t1408*t1427)*y+
t1441;
      coeffs_dxy->coeff_m2_p1 = t1437+t1465+t1478;
      coeffs_dxy->coeff_m1_p1 = t1433+t1476+t1477;
      coeffs_dxy->coeff_0_p1 = (t1405*t1428+t1412)*t1429+(t1419+t1424*t1428)*x+
t1444;
      coeffs_dxy->coeff_p1_p1 = t1435+t1462+t1477;
      coeffs_dxy->coeff_p2_p1 = t1436+t1465+t1466;
      coeffs_dxy->coeff_m2_p2 = t1440+t1471+t1474;
      coeffs_dxy->coeff_m1_p2 = t1438+t1467+t1470;
      coeffs_dxy->coeff_0_p2 = (t1447+t1399)*t1429+(t1393+t1408*t1428)*x+t1443;
      coeffs_dxy->coeff_p1_p2 = t1434+t1467+t1468;
      coeffs_dxy->coeff_p2_p2 = t1439+t1471+t1473;