aboutsummaryrefslogtreecommitdiff
path: root/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-I.compute.c
blob: cfe391a8a07061677eaf9fa713a1385025471347 (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
fp t417,
   t418,
   t468,
   t414,
   t411,
   t361,
   t416,
   t373,
   t435,
   t446,
   t467,
   t410,
   t367,
   t368,
   t437,
   t466,
   t465,
   t441,
   t464,
   t463,
   t395,
   t445,
   t400,
   t403,
   t366,
   t462,
   t397,
   t404,
   t362,
   t401,
   t461,
   t389,
   t409,
   t316,
   t460,
   t402,
   t312,
   t384,
   t459,
   t412,
   t443,
   t407,
   t309,
   t458,
   t295,
   t399,
   t442,
   t457,
   t406,
   t405,
   t319,
   t456,
   t444,
   t370,
   t455,
   t315,
   t454,
   t325,
   t453,
   t358,
   t398,
   t294,
   t452,
   t451,
   t413,
   t364,
   t450,
   t357,
   t380,
   t381,
   t449,
   t332,
   t415,
   t356,
   t448,
   t447,
   t440,
   t360,
   t439,
   t438,
   t436,
   t390,
   t408,
   t363,
   t392,
   t434,
   t386,
   t374,
   t433,
   t372,
   t432,
   t431,
   t317,
   t430,
   t429,
   t326,
   t428,
   t427,
   t385,
   t426,
   t425,
   t388,
   t424,
   t423,
   t394,
   t393,
   t391,
   t387;
      t417 = x*x;
      t418 = y*y;
      t468 = t418+t417;
      t414 = t417*x;
      t411 = RATIONAL(1.0,60.0);
      t361 = t411*t414;
      t416 = y*t418;
      t373 = t411*t416;
      t435 = x*t373+y*t361;
      t446 = x*y;
      t467 = RATIONAL(-11.0,150.0)*t446+t435;
      t410 = RATIONAL(-1.0,60.0);
      t367 = t410*t416;
      t368 = t410*t414;
      t437 = x*t367+y*t368;
      t466 = RATIONAL(37.0,300.0)*t446+t437;
      t465 = t437+RATIONAL(11.0,150.0)*t446;
      t441 = t417*t418;
      t464 = RATIONAL(-1.0,98.0)*t441+RATIONAL(-99.0,1225.0);
      t463 = RATIONAL(-37.0,300.0)*t446+t435;
      t395 = RATIONAL(-31.0,210.0);
      t445 = x*t418;
      t400 = RATIONAL(1.0,140.0);
      t403 = RATIONAL(1.0,30.0);
      t366 = t403*t414;
      t462 = t400*t445+t395*x+t366;
      t397 = RATIONAL(31.0,210.0);
      t404 = RATIONAL(-1.0,30.0);
      t362 = t404*t414;
      t401 = RATIONAL(-1.0,140.0);
      t461 = t362+t401*t445+t397*x;
      t389 = RATIONAL(31.0,420.0);
      t409 = RATIONAL(-1.0,35.0);
      t316 = t409*t445;
      t460 = t389*x+t316+t368;
      t402 = RATIONAL(1.0,35.0);
      t312 = t402*t445;
      t384 = RATIONAL(-31.0,420.0);
      t459 = t312+t384*x+t361;
      t412 = RATIONAL(1.0,84.0);
      t443 = y*t417;
      t407 = RATIONAL(-1.0,70.0);
      t309 = t407*t443;
      t458 = t309+t412*y+t373;
      t295 = RATIONAL(21.0,200.0)*t446;
      t399 = RATIONAL(1.0,120.0);
      t442 = x*t416;
      t457 = t295+t399*t442+y*t362;
      t406 = RATIONAL(-1.0,84.0);
      t405 = RATIONAL(1.0,70.0);
      t319 = t405*t443;
      t456 = t406*y+t319+t367;
      t444 = t414*y;
      t370 = t404*t416;
      t455 = t295+t399*t444+x*t370;
      t315 = t405*t445;
      t454 = t406*x+t315+t368;
      t325 = t407*t445;
      t453 = t361+t412*x+t325;
      t358 = t403*t416;
      t398 = RATIONAL(-1.0,120.0);
      t294 = RATIONAL(-21.0,200.0)*t446;
      t452 = t398*t444+t294+x*t358;
      t451 = y*t366+t398*t442+t294;
      t413 = t417*t417;
      t364 = t404*t413;
      t450 = RATIONAL(-24.0,1225.0)+t364+RATIONAL(1.0,196.0)*t441+t468*RATIONAL
(181.0,1470.0);
      t357 = t399*t413;
      t380 = RATIONAL(71.0,5880.0);
      t381 = RATIONAL(113.0,735.0);
      t449 = t357+t380*t417+t381*t418+t464;
      t332 = RATIONAL(1.0,49.0)*t441;
      t415 = t418*t418;
      t356 = t399*t415;
      t448 = RATIONAL(51.0,1225.0)+t357+t332+t356+t468*RATIONAL(-289.0,5880.0);
      t447 = t380*t418+t364+t381*t417+t356+t464;
      t440 = RATIONAL(1.0,98.0)*t441+RATIONAL(246.0,1225.0);
      t360 = t404*t415;
      t439 = t360+t370;
      t438 = t358+t360;
      t436 = RATIONAL(-1.0,49.0)*t441+RATIONAL(96.0,1225.0);
      t390 = RATIONAL(83.0,735.0);
      t408 = RATIONAL(1.0,20.0);
      t363 = t408*t413;
      t392 = RATIONAL(-53.0,196.0);
      t434 = t360+t390*t418+t363+t392*t417+t440;
      t386 = RATIONAL(-41.0,196.0);
      t374 = RATIONAL(191.0,5880.0);
      t433 = t386*t417+t374*t418+t363+t356+t436;
      t372 = t408*t415;
      t432 = t357+t374*t417+t386*t418+t372+t436;
      t431 = t364+t392*t418+t390*t417+t372+t440;
      t317 = t402*t443;
      t430 = t384*y+t317+t373+t448;
      t429 = t395*y+t400*t443+t438+t450;
      t326 = t409*t443;
      t428 = t389*y+t367+t326+t448;
      t427 = t397*y+t401*t443+t439+t450;
      t385 = RATIONAL(11.0,105.0);
      t426 = t385*x+t315+t362+t447;
      t425 = t385*y+t319+t439+t449;
      t388 = RATIONAL(-11.0,105.0);
      t424 = t388*y+t309+t438+t449;
      t423 = t388*x+t366+t325+t447;
      t394 = RATIONAL(17.0,420.0);
      t393 = RATIONAL(17.0,105.0);
      t391 = RATIONAL(-17.0,105.0);
      t387 = RATIONAL(-17.0,420.0);
      coeff_I_m2_m2 = t428+t460+t467;
      coeff_I_m1_m2 = t423+t456+t457;
      coeff_I_0_m2 = t317+t367+t387*y+t433;
      coeff_I_p1_m2 = t426+t451+t456;
      coeff_I_p2_m2 = t428+t459+t465;
      coeff_I_m2_m1 = t424+t454+t455;
      coeff_I_m1_m1 = t429+t462+t466;
      coeff_I_0_m1 = t319+t391*y+t358+t434;
      coeff_I_p1_m1 = t429+t461+t463;
      coeff_I_p2_m1 = t424+t452+t453;
      coeff_I_m2_0 = t368+t312+t387*x+t432;
      coeff_I_m1_0 = t366+t315+t391*x+t431;
      coeff_I_0_0 = RATIONAL(541.0,1225.0)+t363+t372+t332+t468*RATIONAL(-57.0,
196.0);
      coeff_I_p1_0 = t325+t362+t393*x+t431;
      coeff_I_p2_0 = t394*x+t316+t361+t432;
      coeff_I_m2_p1 = t425+t452+t454;
      coeff_I_m1_p1 = t427+t462+t463;
      coeff_I_0_p1 = t309+t370+t393*y+t434;
      coeff_I_p1_p1 = t427+t461+t466;
      coeff_I_p2_p1 = t425+t453+t455;
      coeff_I_m2_p2 = t430+t460+t465;
      coeff_I_m1_p2 = t423+t451+t458;
      coeff_I_0_p2 = t394*y+t326+t373+t433;
      coeff_I_p1_p2 = t426+t457+t458;
      coeff_I_p2_p2 = t430+t459+t467;