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