aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxx.compute.c
blob: 3f214f70f8781dbcd0525889132b0a2c5d334e38 (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
fp t1370;
fp t1373;
fp t1440;
fp t1340;
fp t1439;
fp t1338;
fp t1438;
fp t1374;
fp t1372;
fp t1412;
fp t1410;
fp t1437;
fp t1364;
fp t1415;
fp t1375;
fp t1341;
fp t1397;
fp t1261;
fp t1436;
fp t1435;
fp t1409;
fp t1414;
fp t1343;
fp t1265;
fp t1434;
fp t1365;
fp t1394;
fp t1433;
fp t1339;
fp t1393;
fp t1342;
fp t1432;
fp t1413;
fp t1417;
fp t1431;
fp t1311;
fp t1430;
fp t1429;
fp t1428;
fp t1371;
fp t1368;
fp t1427;
fp t1400;
fp t1426;
fp t1408;
fp t1369;
fp t1425;
fp t1424;
fp t1411;
fp t1423;
fp t1422;
fp t1392;
fp t1346;
fp t1404;
fp t1345;
fp t1421;
fp t1403;
fp t1344;
fp t1420;
fp t1419;
fp t1405;
fp t1362;
fp t1363;
fp t1366;
fp t1418;
fp t1416;
fp t1350;
fp t1399;
fp t1361;
fp t1407;
fp t1354;
fp t1322;
fp t1324;
fp t1353;
fp t1406;
fp t1402;
fp t1401;
fp t1356;
fp t1391;
fp t1357;
fp t1390;
fp t1351;
fp t1389;
fp t1347;
fp t1388;
fp t1352;
fp t1387;
fp t1355;
fp t1386;
fp t1385;
fp t1384;
fp t1383;
fp t1382;
fp t1367;
fp t1381;
fp t1380;
fp t1379;
fp t1378;
      t1370 = RATIONAL(-1.0,6.0);
      t1373 = x*x;
      t1440 = t1370*t1373+RATIONAL(1.0,9.0);
      t1340 = RATIONAL(-1.0,24.0);
      t1439 = t1373*t1340+RATIONAL(1.0,144.0);
      t1338 = RATIONAL(1.0,24.0);
      t1438 = t1373*t1338+RATIONAL(-1.0,144.0);
      t1374 = y*y;
      t1372 = t1374*t1374;
      t1412 = x*t1372;
      t1410 = x*t1374;
      t1437 = t1338*t1410+t1340*t1412;
      t1364 = RATIONAL(-1.0,3.0);
      t1415 = y*t1364;
      t1375 = t1374*y;
      t1341 = RATIONAL(1.0,12.0);
      t1397 = t1375*t1341;
      t1261 = x*t1397;
      t1436 = x*t1415+t1261;
      t1435 = t1340*t1410+t1338*t1412;
      t1409 = x*t1341;
      t1414 = x*t1375;
      t1343 = RATIONAL(-1.0,12.0);
      t1265 = t1343*t1414;
      t1434 = t1265+y*t1409;
      t1365 = RATIONAL(1.0,3.0);
      t1394 = t1374*t1365;
      t1433 = x*t1394+t1343*t1412;
      t1339 = RATIONAL(1.0,48.0);
      t1393 = t1372*t1339;
      t1342 = RATIONAL(-1.0,48.0);
      t1432 = t1342*t1410+x*t1393;
      t1413 = x*t1338;
      t1417 = x*y;
      t1431 = t1375*t1413+t1340*t1417;
      t1311 = t1342*t1372;
      t1430 = t1339*t1410+x*t1311;
      t1429 = t1340*t1414+y*t1413;
      t1428 = t1372*t1409+t1364*t1410;
      t1371 = RATIONAL(-2.0,3.0);
      t1368 = RATIONAL(1.0,6.0);
      t1427 = t1368*t1412+t1371*t1410;
      t1400 = t1375*t1368;
      t1426 = x*t1400+t1371*t1417;
      t1408 = t1370*x;
      t1369 = RATIONAL(2.0,3.0);
      t1425 = t1369*t1417+t1375*t1408;
      t1424 = t1372*t1408+t1369*t1410;
      t1411 = y*t1365;
      t1423 = t1265+x*t1411;
      t1422 = t1343*t1417+t1261;
      t1392 = t1373*t1343;
      t1346 = RATIONAL(-1.0,18.0);
      t1404 = t1346*t1374+t1372*t1392;
      t1345 = RATIONAL(1.0,72.0);
      t1421 = t1404+t1345*t1372+t1373*t1394;
      t1403 = t1373*t1374;
      t1344 = RATIONAL(1.0,18.0);
      t1420 = t1344*t1372+t1404+t1341*t1403;
      t1419 = t1342*t1403+t1373*t1393+RATIONAL(-1.0,288.0)*t1372+RATIONAL(1.0,
288.0)*t1374;
      t1405 = t1372*t1373;
      t1362 = RATIONAL(-2.0,9.0);
      t1363 = RATIONAL(-4.0,3.0);
      t1366 = RATIONAL(8.0,9.0);
      t1418 = t1362*t1372+t1366*t1374+t1365*t1405+t1363*t1403;
      t1416 = y*t1373;
      t1350 = RATIONAL(-1.0,2.0);
      t1399 = t1350*t1373;
      t1361 = RATIONAL(-5.0,3.0);
      t1407 = t1361*t1374+t1372*t1399;
      t1354 = RATIONAL(1.0,2.0);
      t1322 = t1354*t1373;
      t1324 = RATIONAL(-2.0,1.0)*t1373;
      t1353 = RATIONAL(1.0,8.0);
      t1406 = t1353*t1405+RATIONAL(5.0,48.0)*t1374;
      t1402 = t1373*t1375;
      t1401 = RATIONAL(2.0,1.0)*t1373;
      t1356 = RATIONAL(-1.0,4.0);
      t1391 = t1356*t1373+RATIONAL(5.0,24.0);
      t1357 = RATIONAL(1.0,4.0);
      t1390 = t1357*t1373+RATIONAL(-5.0,24.0);
      t1351 = RATIONAL(-1.0,8.0);
      t1389 = RATIONAL(-5.0,48.0)*t1372+t1351*t1403+t1406;
      t1347 = RATIONAL(5.0,12.0);
      t1388 = t1347*t1372+t1374*t1401+t1407;
      t1352 = RATIONAL(4.0,3.0);
      t1387 = t1352+RATIONAL(5.0,2.0)*t1403+t1365*t1372+t1324+t1407;
      t1355 = RATIONAL(-5.0,8.0);
      t1386 = t1355*t1403+t1322+t1311+t1343+t1406;
      t1385 = t1419+t1439*t1375+t1438*y;
      t1384 = t1419+t1438*t1375+t1439*y;
      t1383 = t1421+t1373*t1411+t1346*y+(t1392+t1345)*t1375;
      t1382 = t1344*y+RATIONAL(-1.0,72.0)*t1375+t1421+(t1415+t1397)*t1373;
      t1367 = RATIONAL(-1.0,9.0);
      t1381 = t1367*y+t1368*t1416+t1420+t1440*t1375;
      t1380 = RATIONAL(-8.0,9.0)*y+RATIONAL(2.0,9.0)*t1375+t1364*t1402+t1352*
t1416+t1418;
      t1379 = t1420+t1373*t1400+t1367*t1375+t1440*y;
      t1378 = t1363*t1416+t1362*t1375+t1365*t1402+y*t1366+t1418;
      coeffs_dxx->coeff_m2_m2 = t1385+t1430+t1431;
      coeffs_dxx->coeff_m1_m2 = t1379+t1434+t1435;
      coeffs_dxx->coeff_0_m2 = t1391*t1375+t1390*y+t1389;
      coeffs_dxx->coeff_p1_m2 = t1379+t1422+t1437;
      coeffs_dxx->coeff_p2_m2 = t1385+t1429+t1432;
      coeffs_dxx->coeff_m2_m1 = t1382+t1423+t1428;
      coeffs_dxx->coeff_m1_m1 = t1380+t1424+t1426;
      coeffs_dxx->coeff_0_m1 = (RATIONAL(-5.0,12.0)+t1322)*t1375+(RATIONAL(5.0,
3.0)+t1324)*y+t1388;
      coeffs_dxx->coeff_p1_m1 = t1380+t1425+t1427;
      coeffs_dxx->coeff_p2_m1 = t1382+t1433+t1436;
      coeffs_dxx->coeff_m2_0 = (t1351*t1372+RATIONAL(5.0,8.0)*t1374+t1350)*x+
t1386;
      coeffs_dxx->coeff_m1_0 = (t1357*t1372+1.0+RATIONAL(-5.0,4.0)*t1374)*x+
t1387;
      coeffs_dxx->coeff_0_0 = RATIONAL(25.0,8.0)*t1374+RATIONAL(-5.0,2.0)+t1355
*t1372+(RATIONAL(3.0,4.0)*t1372+RATIONAL(-15.0,4.0)*t1374+RATIONAL(3.0,1.0))*
t1373;
      coeffs_dxx->coeff_p1_0 = (RATIONAL(5.0,4.0)*t1374+t1356*t1372-1.0)*x+
t1387;
      coeffs_dxx->coeff_p2_0 = (t1353*t1372+t1355*t1374+t1354)*x+t1386;
      coeffs_dxx->coeff_m2_p1 = t1383+t1428+t1436;
      coeffs_dxx->coeff_m1_p1 = t1378+t1424+t1425;
      coeffs_dxx->coeff_0_p1 = (t1347+t1399)*t1375+(t1401+t1361)*y+t1388;
      coeffs_dxx->coeff_p1_p1 = t1378+t1426+t1427;
      coeffs_dxx->coeff_p2_p1 = t1383+t1423+t1433;
      coeffs_dxx->coeff_m2_p2 = t1384+t1429+t1430;
      coeffs_dxx->coeff_m1_p2 = t1381+t1422+t1435;
      coeffs_dxx->coeff_0_p2 = t1390*t1375+t1391*y+t1389;
      coeffs_dxx->coeff_p1_p2 = t1381+t1434+t1437;
      coeffs_dxx->coeff_p2_p2 = t1384+t1431+t1432;