aboutsummaryrefslogtreecommitdiff
path: root/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeffs-I.compute.c
blob: a5b63608c7164fb26400d1929ee60d121ab35422 (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
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
fp t469;
fp t472;
fp t517;
fp t471;
fp t470;
fp t516;
fp t515;
fp t514;
fp t409;
fp t513;
fp t416;
fp t512;
fp t511;
fp t435;
fp t510;
fp t509;
fp t423;
fp t508;
fp t507;
fp t506;
fp t413;
fp t487;
fp t368;
fp t489;
fp t445;
fp t505;
fp t410;
fp t504;
fp t499;
fp t503;
fp t502;
fp t449;
fp t501;
fp t453;
fp t491;
fp t401;
fp t441;
fp t500;
fp t442;
fp t490;
fp t408;
fp t498;
fp t497;
fp t421;
fp t366;
fp t496;
fp t495;
fp t403;
fp t494;
fp t404;
fp t493;
fp t424;
fp t492;
fp t434;
fp t488;
fp t450;
fp t371;
fp t486;
fp t485;
fp t484;
fp t389;
fp t483;
fp t444;
fp t482;
fp t481;
fp t480;
fp t479;
fp t426;
fp t428;
fp t478;
fp t475;
fp t468;
fp t467;
fp t466;
fp t465;
fp t464;
fp t463;
fp t462;
fp t461;
fp t460;
fp t459;
fp t458;
fp t457;
fp t456;
fp t454;
fp t452;
fp t451;
fp t448;
fp t447;
fp t446;
fp t443;
fp t440;
fp t439;
fp t438;
fp t437;
fp t436;
fp t433;
fp t432;
fp t431;
fp t430;
fp t429;
fp t427;
fp t425;
fp t422;
fp t420;
fp t419;
fp t418;
fp t417;
fp t415;
fp t414;
fp t412;
fp t411;
fp t407;
fp t406;
fp t405;
fp t402;
fp t400;
fp t396;
fp t395;
fp t394;
fp t393;
fp t391;
fp t386;
fp t383;
fp t382;
fp t380;
fp t379;
fp t374;
fp t373;
fp t372;
fp t369;
fp t364;
fp t362;
fp t361;
fp t359;
fp t358;
fp t357;
fp t356;
fp t355;
fp t354;
fp t352;
      t469 = y*y;
      t472 = t469*y;
      t517 = t472+y;
      t471 = x*x;
      t470 = x*t471;
      t516 = t470+x;
      t515 = x*y;
      t514 = x*t472;
      t409 = RATIONAL(1.0,144.0);
      t513 = x*t409;
      t416 = RATIONAL(-1.0,144.0);
      t512 = t416*x;
      t511 = y*t470;
      t435 = RATIONAL(1.0,18.0);
      t510 = t435*x;
      t509 = t435*y;
      t423 = RATIONAL(-1.0,18.0);
      t508 = t423*x;
      t507 = t470*t469;
      t506 = t471*t472;
      t413 = RATIONAL(-7.0,144.0);
      t487 = t413*t470;
      t368 = t472*t487;
      t489 = t469*t471;
      t445 = RATIONAL(-1.0,24.0);
      t505 = t445*t489+t368;
      t410 = RATIONAL(7.0,144.0);
      t504 = t410*t472;
      t499 = t470*t472;
      t503 = RATIONAL(-35.0,9.0)*t489+RATIONAL(-16.0,9.0)*t499;
      t502 = t410*t470;
      t449 = RATIONAL(1.0,9.0);
      t501 = t449*t472;
      t453 = RATIONAL(-1.0,9.0);
      t491 = t453*t472;
      t401 = t470*t491;
      t441 = RATIONAL(-7.0,36.0);
      t500 = t441*t489+t401;
      t442 = RATIONAL(7.0,18.0);
      t490 = t449*t470;
      t408 = t472*t490;
      t498 = t442*t489+t408;
      t497 = t413*t472;
      t421 = RATIONAL(5.0,48.0);
      t366 = t410*t499;
      t496 = t421*t489+t366;
      t495 = t453*t470;
      t403 = RATIONAL(7.0,9.0)*t499;
      t494 = RATIONAL(7.0,6.0)*t489+t403;
      t404 = RATIONAL(-7.0,9.0)*t499;
      t493 = t404+RATIONAL(-35.0,12.0)*t489;
      t424 = RATIONAL(-1.0,72.0);
      t492 = t424*t489+t416*t499;
      t434 = RATIONAL(5.0,36.0);
      t488 = t408+t434*t489;
      t450 = RATIONAL(-5.0,6.0);
      t371 = y*t508;
      t486 = t450*t489+t404+t371;
      t485 = RATIONAL(-4.0,9.0)*t515+t403+RATIONAL(25.0,12.0)*t489;
      t484 = t368+RATIONAL(-5.0,24.0)*t489+t371;
      t389 = x*t509;
      t483 = RATIONAL(-5.0,18.0)*t489+t389+t401;
      t444 = RATIONAL(1.0,12.0);
      t482 = y*t512+t444*t489+t366;
      t481 = t389+RATIONAL(-49.0,144.0)*t499+RATIONAL(-5.0,8.0)*t489;
      t480 = t516*t469;
      t479 = t517*t471;
      t426 = RATIONAL(7.0,36.0);
      t428 = RATIONAL(-1.0,12.0);
      t478 = t428+t426*t471;
      t475 = t426*t469+t428;
      t468 = RATIONAL(5.0,3.0);
      t467 = RATIONAL(-8.0,9.0);
      t466 = RATIONAL(5.0,4.0);
      t465 = RATIONAL(-2.0,3.0);
      t464 = RATIONAL(8.0,9.0);
      t463 = RATIONAL(-5.0,3.0);
      t462 = RATIONAL(5.0,6.0);
      t461 = RATIONAL(1.0,3.0);
      t460 = RATIONAL(-1.0,6.0);
      t459 = RATIONAL(-1.0,2.0);
      t458 = RATIONAL(-2.0,9.0);
      t457 = RATIONAL(-1.0,3.0);
      t456 = RATIONAL(2.0,9.0);
      t454 = RATIONAL(-7.0,3.0);
      t452 = RATIONAL(-4.0,3.0);
      t451 = RATIONAL(2.0,3.0);
      t448 = RATIONAL(1.0,72.0);
      t447 = RATIONAL(14.0,9.0);
      t446 = RATIONAL(-7.0,12.0);
      t443 = RATIONAL(-10.0,9.0);
      t440 = RATIONAL(1.0,24.0);
      t439 = RATIONAL(-7.0,24.0);
      t438 = RATIONAL(-7.0,18.0);
      t437 = RATIONAL(-5.0,48.0);
      t436 = RATIONAL(7.0,24.0);
      t433 = RATIONAL(-5.0,36.0);
      t432 = RATIONAL(-28.0,9.0);
      t431 = RATIONAL(28.0,9.0);
      t430 = RATIONAL(20.0,9.0);
      t429 = RATIONAL(-14.0,9.0);
      t427 = RATIONAL(7.0,12.0);
      t425 = RATIONAL(-7.0,72.0);
      t422 = RATIONAL(10.0,9.0);
      t420 = RATIONAL(7.0,72.0);
      t419 = RATIONAL(-20.0,9.0);
      t418 = RATIONAL(-35.0,48.0);
      t417 = RATIONAL(49.0,36.0);
      t415 = RATIONAL(-35.0,36.0);
      t414 = RATIONAL(35.0,48.0);
      t412 = RATIONAL(35.0,36.0);
      t411 = RATIONAL(-49.0,36.0);
      t407 = RATIONAL(4.0,9.0)*t515;
      t406 = y*t490;
      t405 = t467*t511;
      t402 = y*t495;
      t400 = t464*t511;
      t396 = x*t501;
      t395 = t467*t514;
      t394 = t464*t514;
      t393 = x*t491;
      t391 = t442*t511;
      t386 = t472*t510;
      t383 = t438*t514;
      t382 = t438*t511;
      t380 = t470*t509;
      t379 = RATIONAL(16.0,9.0)*t499;
      t374 = t423*t511;
      t373 = t472*t508;
      t372 = t442*t514;
      t369 = y*t487;
      t364 = y*t502;
      t362 = t409*t511;
      t361 = t472*t512;
      t359 = x*t504;
      t358 = t416*t511;
      t357 = t472*t513;
      t356 = x*t497;
      t355 = t409*t499;
      t354 = y*t513;
      t352 = RATIONAL(49.0,144.0)*t499;
      coeffs_I->coeff_m2_m2 = t354+t357+t362+t355+RATIONAL(1.0,36.0)*t489+(t479
+t480)*t424;
      coeffs_I->coeff_m1_m2 = t369+t373+t421*t479+(t420*t470+t449*x)*t469+t484;
      coeffs_I->coeff_0_m2 = t406+(t458*t470+t460)*t469+t498+t517*(t444+t441*
t471);
      coeffs_I->coeff_p1_m2 = t386+t402+t434*t479+(t456*t470+t453*x)*t469+t483;
      coeffs_I->coeff_p2_m2 = t361+t364+t445*t479+(t425*t470+t448*x)*t469+t482;
      coeffs_I->coeff_p3_m2 = t358+t448*t507+t409*t479+t492;
      coeffs_I->coeff_m2_m1 = t356+t374+(t420*t472+t449*y)*t471+t421*t480+t484;
      coeffs_I->coeff_m1_m1 = t372+t407+t352+t391+(t450*y+t418*t472)*t471+(t450
*x+RATIONAL(25.0,16.0)*t471+t418*t470)*t469;
      coeffs_I->coeff_0_m1 = t405+(t417*t471+t446)*t472+(t468*t470+t466)*t469+(
t447*t471+t465)*y+t493;
      coeffs_I->coeff_p1_m1 = t400+t383+(t415*t472+t443*y)*t471+(t462*x+t463*
t470)*t469+t485;
      coeffs_I->coeff_p2_m1 = t359+t382+(t436*t472+t461*y)*t471+(t414*t470+t437
*x)*t469+t481;
      coeffs_I->coeff_p3_m1 = t437*t507+t380+(t423*y+t497)*t471+t496;
      coeffs_I->coeff_m2_0 = t396+(t460+t458*t472)*t471+t498+t516*(t444+t441*
t469);
      coeffs_I->coeff_m1_0 = t395+(t468*t472+t466)*t471+(t417*t469+t446)*t470+(
t465+t447*t469)*x+t493;
      coeffs_I->coeff_0_0 = t379+RATIONAL(1.0,1.0)+(t470+t472)*RATIONAL(4.0,3.0
)+(t454+t432*t472)*t471+(RATIONAL(49.0,9.0)*t471+t454+t432*t470)*t469;
      coeffs_I->coeff_p1_0 = t394+(t430*t472+t468)*t471+(t452+t431*t469)*t470+(
t429*t469+t451)*x+t503;
      coeffs_I->coeff_p2_0 = t393+(t465*t472+t459)*t471+(t411*t469+t427)*t470+
t475*x+t494;
      coeffs_I->coeff_p3_0 = (t444+t501)*t471+t475*t470+t500;
      coeffs_I->coeff_m2_p1 = t393+t380+(t456*t472+t453*y)*t471+t434*t480+t483;
      coeffs_I->coeff_m1_p1 = t382+t394+(t462*y+t463*t472)*t471+(t443*x+t415*
t470)*t469+t485;
      coeffs_I->coeff_0_p1 = t400+(t452+t431*t471)*t472+(t468+t430*t470)*t469+(
t451+t429*t471)*y+t503;
      coeffs_I->coeff_p1_p1 = t379+t405+t407+t395+(t422*y+t419*t472)*t471+(t422
*x+t419*t470+RATIONAL(25.0,9.0)*t471)*t469;
      coeffs_I->coeff_p2_p1 = t391+t396+(t457*y+t451*t472)*t471+(t412*t470+t433
*x)*t469+t486;
      coeffs_I->coeff_p3_p1 = t433*t507+t374+(t509+t491)*t471+t488;
      coeffs_I->coeff_m2_p2 = t359+t358+(t425*t472+t448*y)*t471+t445*t480+t482;
      coeffs_I->coeff_m1_p2 = t383+t364+(t414*t472+t437*y)*t471+(t461*x+t436*
t470)*t469+t481;
      coeffs_I->coeff_0_p2 = t402+(t411*t471+t427)*t472+(t459+t465*t470)*t469+
t478*y+t494;
      coeffs_I->coeff_p1_p2 = t372+t406+(t412*t472+t433*y)*t471+(t457*x+t451*
t470)*t469+t486;
      coeffs_I->coeff_p2_p2 = t369+t354+t352+t356+(t439*t472+t440*y)*t471+(t440
*x+t439*t470+RATIONAL(1.0,4.0)*t471)*t469;
      coeffs_I->coeff_p3_p2 = t440*t507+t362+(t416*y+t504)*t471+t505;
      coeffs_I->coeff_m2_p3 = t361+t448*t506+t409*t480+t492;
      coeffs_I->coeff_m1_p3 = t386+t437*t506+(t508+t487)*t469+t496;
      coeffs_I->coeff_0_p3 = t478*t472+(t490+t444)*t469+t500;
      coeffs_I->coeff_p1_p3 = t373+t433*t506+(t510+t495)*t469+t488;
      coeffs_I->coeff_p2_p3 = t440*t506+t357+(t512+t502)*t469+t505;
      coeffs_I->coeff_p3_p3 = t355+t409*t489+(t507+t506)*t416;