aboutsummaryrefslogtreecommitdiff
path: root/src/GeneralizedPolynomial-Uniform/Lagrange-maximum-degree/2d.coeffs/2d.cube.order4.smooth0/coeffs-I.compute.c
blob: 51bec8da9cae8c93113a996c96cbbf905035377f (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
fp t411;
fp t412;
fp t464;
fp t414;
fp t404;
fp t356;
fp t413;
fp t360;
fp t433;
fp t442;
fp t463;
fp t407;
fp t365;
fp t364;
fp t436;
fp t462;
fp t435;
fp t461;
fp t438;
fp t389;
fp t398;
fp t406;
fp t362;
fp t460;
fp t401;
fp t368;
fp t395;
fp t397;
fp t459;
fp t439;
fp t359;
fp t458;
fp t370;
fp t291;
fp t396;
fp t440;
fp t457;
fp t405;
fp t322;
fp t409;
fp t456;
fp t402;
fp t400;
fp t327;
fp t455;
fp t383;
fp t454;
fp t311;
fp t453;
fp t390;
fp t307;
fp t452;
fp t451;
fp t384;
fp t408;
fp t318;
fp t450;
fp t399;
fp t296;
fp t449;
fp t448;
fp t387;
fp t403;
fp t330;
fp t447;
fp t416;
fp t354;
fp t415;
fp t355;
fp t314;
fp t446;
fp t369;
fp t358;
fp t445;
fp t376;
fp t379;
fp t444;
fp t443;
fp t441;
fp t437;
fp t434;
fp t410;
fp t366;
fp t394;
fp t372;
fp t432;
fp t371;
fp t431;
fp t391;
fp t393;
fp t430;
fp t429;
fp t428;
fp t306;
fp t427;
fp t426;
fp t316;
fp t425;
fp t424;
fp t423;
fp t422;
fp t421;
fp t392;
fp t386;
fp t385;
fp t382;
      t411 = x*x;
      t412 = y*y;
      t464 = t411+t412;
      t414 = t411*x;
      t404 = RATIONAL(-1.0,60.0);
      t356 = t404*t414;
      t413 = y*t412;
      t360 = t404*t413;
      t433 = y*t356+x*t360;
      t442 = x*y;
      t463 = RATIONAL(11.0,150.0)*t442+t433;
      t407 = RATIONAL(1.0,60.0);
      t365 = t407*t414;
      t364 = t407*t413;
      t436 = x*t364+y*t365;
      t462 = RATIONAL(-11.0,150.0)*t442+t436;
      t435 = t411*t412;
      t461 = RATIONAL(-99.0,1225.0)+RATIONAL(-1.0,98.0)*t435;
      t438 = y*t411;
      t389 = RATIONAL(31.0,210.0);
      t398 = RATIONAL(-1.0,140.0);
      t406 = RATIONAL(-1.0,30.0);
      t362 = t406*t413;
      t460 = t389*y+t398*t438+t362;
      t401 = RATIONAL(1.0,30.0);
      t368 = t401*t413;
      t395 = RATIONAL(-31.0,210.0);
      t397 = RATIONAL(1.0,140.0);
      t459 = t368+t395*y+t397*t438;
      t439 = x*t412;
      t359 = t401*t414;
      t458 = t395*x+t397*t439+t359;
      t370 = t406*t414;
      t291 = RATIONAL(21.0,200.0)*t442;
      t396 = RATIONAL(1.0,120.0);
      t440 = x*t413;
      t457 = t291+t396*t440+y*t370;
      t405 = RATIONAL(-1.0,70.0);
      t322 = t405*t438;
      t409 = RATIONAL(1.0,84.0);
      t456 = t322+t409*y+t364;
      t402 = RATIONAL(-1.0,84.0);
      t400 = RATIONAL(1.0,70.0);
      t327 = t400*t439;
      t455 = t327+t356+t402*x;
      t383 = RATIONAL(-11.0,105.0);
      t454 = t368+t322+t383*y;
      t311 = t405*t439;
      t453 = t409*x+t311+t365;
      t390 = RATIONAL(11.0,105.0);
      t307 = t400*t438;
      t452 = t307+t390*y+t362;
      t451 = t307+t402*y+t360;
      t384 = RATIONAL(31.0,420.0);
      t408 = RATIONAL(-1.0,35.0);
      t318 = t408*t439;
      t450 = t384*x+t356+t318;
      t399 = RATIONAL(-1.0,120.0);
      t296 = RATIONAL(-21.0,200.0)*t442;
      t449 = y*t359+t399*t440+t296;
      t448 = t389*x+t370+t398*t439;
      t387 = RATIONAL(-31.0,420.0);
      t403 = RATIONAL(1.0,35.0);
      t330 = t403*t439;
      t447 = t330+t387*x+t365;
      t416 = t412*t412;
      t354 = t396*t416;
      t415 = t411*t411;
      t355 = t396*t415;
      t314 = RATIONAL(1.0,49.0)*t435;
      t446 = RATIONAL(51.0,1225.0)+t314+t354+t355+t464*RATIONAL(-289.0,5880.0);
      t369 = t406*t416;
      t358 = t406*t415;
      t445 = RATIONAL(1.0,196.0)*t435+t369+RATIONAL(-24.0,1225.0)+t358+t464*
RATIONAL(181.0,1470.0);
      t376 = RATIONAL(71.0,5880.0);
      t379 = RATIONAL(113.0,735.0);
      t444 = t376*t412+t379*t411+t354+t358+t461;
      t443 = t369+t376*t411+t379*t412+t355+t461;
      t441 = t414*y;
      t437 = RATIONAL(-1.0,49.0)*t435+RATIONAL(96.0,1225.0);
      t434 = RATIONAL(1.0,98.0)*t435+RATIONAL(246.0,1225.0);
      t410 = RATIONAL(1.0,20.0);
      t366 = t410*t416;
      t394 = RATIONAL(-41.0,196.0);
      t372 = RATIONAL(191.0,5880.0);
      t432 = t366+t372*t411+t355+t394*t412+t437;
      t371 = t410*t415;
      t431 = t371+t394*t411+t372*t412+t354+t437;
      t391 = RATIONAL(83.0,735.0);
      t393 = RATIONAL(-53.0,196.0);
      t430 = t369+t391*t412+t371+t393*t411+t434;
      t429 = t366+t393*t412+t391*t411+t358+t434;
      t428 = RATIONAL(-37.0,300.0)*t442+t436+t445;
      t306 = t403*t438;
      t427 = t306+t387*y+t364+t446;
      t426 = RATIONAL(37.0,300.0)*t442+t433+t445;
      t316 = t408*t438;
      t425 = t384*y+t360+t316+t446;
      t424 = t390*x+t370+t327+t444;
      t423 = t383*x+t311+t359+t444;
      t422 = x*t368+t296+t399*t441+t443;
      t421 = t396*t441+t291+x*t362+t443;
      t392 = RATIONAL(-17.0,105.0);
      t386 = RATIONAL(17.0,105.0);
      t385 = RATIONAL(17.0,420.0);
      t382 = RATIONAL(-17.0,420.0);
      coeffs_I->coeff_m2_m2 = t425+t450+t462;
      coeffs_I->coeff_m1_m2 = t423+t451+t457;
      coeffs_I->coeff_0_m2 = t306+t382*y+t360+t431;
      coeffs_I->coeff_p1_m2 = t424+t449+t451;
      coeffs_I->coeff_p2_m2 = t425+t447+t463;
      coeffs_I->coeff_m2_m1 = t421+t454+t455;
      coeffs_I->coeff_m1_m1 = t426+t458+t459;
      coeffs_I->coeff_0_m1 = t392*y+t368+t307+t430;
      coeffs_I->coeff_p1_m1 = t428+t448+t459;
      coeffs_I->coeff_p2_m1 = t422+t453+t454;
      coeffs_I->coeff_m2_0 = t356+t330+t382*x+t432;
      coeffs_I->coeff_m1_0 = t359+t392*x+t327+t429;
      coeffs_I->coeff_0_0 = RATIONAL(541.0,1225.0)+t366+t314+t371+t464*RATIONAL
(-57.0,196.0);
      coeffs_I->coeff_p1_0 = t311+t370+t386*x+t429;
      coeffs_I->coeff_p2_0 = t318+t365+t385*x+t432;
      coeffs_I->coeff_m2_p1 = t422+t452+t455;
      coeffs_I->coeff_m1_p1 = t428+t458+t460;
      coeffs_I->coeff_0_p1 = t386*y+t362+t322+t430;
      coeffs_I->coeff_p1_p1 = t426+t448+t460;
      coeffs_I->coeff_p2_p1 = t421+t452+t453;
      coeffs_I->coeff_m2_p2 = t427+t450+t463;
      coeffs_I->coeff_m1_p2 = t423+t449+t456;
      coeffs_I->coeff_0_p2 = t316+t364+t385*y+t431;
      coeffs_I->coeff_p1_p2 = t424+t456+t457;
      coeffs_I->coeff_p2_p2 = t427+t447+t462;