aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-I.compute.c
blob: 1c237bcc93b7237d9e6aefe58bf8c80eb31d73ae (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
fp t519;
fp t520;
fp t522;
fp t525;
fp t616;
fp t513;
fp t512;
fp t615;
fp t507;
fp t592;
fp t521;
fp t523;
fp t524;
fp t492;
fp t549;
fp t614;
fp t490;
fp t588;
fp t493;
fp t570;
fp t613;
fp t574;
fp t488;
fp t496;
fp t556;
fp t526;
fp t563;
fp t582;
fp t508;
fp t552;
fp t484;
fp t612;
fp t557;
fp t554;
fp t611;
fp t487;
fp t559;
fp t590;
fp t562;
fp t489;
fp t579;
fp t610;
fp t583;
fp t565;
fp t553;
fp t609;
fp t581;
fp t608;
fp t558;
fp t517;
fp t518;
fp t576;
fp t580;
fp t607;
fp t584;
fp t606;
fp t586;
fp t605;
fp t561;
fp t564;
fp t604;
fp t578;
fp t603;
fp t577;
fp t602;
fp t560;
fp t568;
fp t551;
fp t601;
fp t585;
fp t550;
fp t600;
fp t599;
fp t589;
fp t555;
fp t598;
fp t597;
fp t596;
fp t595;
fp t594;
fp t593;
fp t591;
fp t587;
fp t506;
fp t479;
fp t477;
fp t514;
fp t480;
fp t483;
fp t510;
fp t572;
fp t571;
fp t497;
fp t569;
fp t567;
fp t548;
fp t547;
fp t495;
fp t546;
fp t499;
fp t502;
fp t545;
fp t544;
fp t543;
fp t505;
fp t500;
fp t503;
fp t542;
fp t501;
fp t504;
fp t498;
fp t541;
fp t540;
fp t539;
fp t538;
fp t537;
fp t536;
fp t535;
fp t534;
fp t533;
fp t532;
fp t531;
fp t516;
fp t515;
fp t511;
fp t509;
fp t494;
fp t491;
      t519 = RATIONAL(-1.0,9.0);
      t520 = RATIONAL(4.0,9.0);
      t522 = y*y;
      t525 = t522*t522;
      t616 = t520*t522+t519*t525;
      t513 = RATIONAL(1.0,6.0);
      t512 = RATIONAL(2.0,3.0);
      t615 = t512+t513*t525;
      t507 = RATIONAL(-1.0,18.0);
      t592 = x*y;
      t521 = x*x;
      t523 = t521*x;
      t524 = y*t522;
      t492 = RATIONAL(-1.0,72.0);
      t549 = t524*t492;
      t614 = t507*t592+t523*t549;
      t490 = RATIONAL(1.0,18.0);
      t588 = x*t490;
      t493 = RATIONAL(1.0,72.0);
      t570 = t523*t493;
      t613 = t524*t570+y*t588;
      t574 = t521*t525;
      t488 = RATIONAL(-1.0,144.0);
      t496 = RATIONAL(-1.0,36.0);
      t556 = t521*t496;
      t526 = t521*t521;
      t563 = t526*t525;
      t582 = t488*t563+t522*t556;
      t508 = RATIONAL(1.0,36.0);
      t552 = t526*t508;
      t484 = RATIONAL(1.0,144.0);
      t612 = t484*t574+t582+t522*t552;
      t557 = t522*t526;
      t554 = t525*t508;
      t611 = t484*t557+t582+t521*t554;
      t487 = RATIONAL(-1.0,288.0);
      t559 = t487*t521;
      t590 = y*t526;
      t562 = t526*t524;
      t489 = RATIONAL(1.0,288.0);
      t579 = t521*t489;
      t610 = t487*t562+y*t559+t524*t579+t489*t590;
      t583 = t523*t484;
      t565 = t524*t484;
      t553 = t523*t488;
      t609 = t488*t592+x*t565+y*t583+t524*t553;
      t581 = t522*t523;
      t608 = t519*t581+t523*t554+t616*x;
      t558 = t524*t496;
      t517 = RATIONAL(-4.0,9.0);
      t518 = RATIONAL(1.0,9.0);
      t576 = t518*t523;
      t580 = t518*t524;
      t607 = y*t576+t523*t558+x*t580+t517*t592;
      t584 = t521*t522;
      t606 = (t574+t557)*RATIONAL(-1.0,576.0)+(t584+t563)*RATIONAL(1.0,576.0);
      t586 = x*t524;
      t605 = t490*t586+y*t570+t614;
      t561 = t507*t521;
      t564 = t521*t490;
      t604 = t526*t549+t524*t564+y*t561+t493*t590;
      t578 = t524*t488;
      t603 = t484*t592+x*t578+y*t553+t523*t565;
      t577 = t523*t490;
      t602 = t493*t586+y*t577+t614;
      t560 = t525*t493;
      t568 = t507*t523;
      t551 = t525*t492;
      t601 = x*t551+t523*t560+(t568+t588)*t522;
      t585 = y*t523;
      t550 = t519*t524;
      t600 = x*t550+t520*t592+t519*t585+t508*t523*t524;
      t599 = t487*t590+t489*t562+y*t579+t524*t559;
      t589 = x*t522;
      t555 = t525*t496;
      t598 = t522*t576+t517*t589+t523*t555+t518*x*t525;
      t597 = x*t560+t523*t551+t522*t577+t507*t589;
      t596 = x*t549+y*t568+t613;
      t595 = t492*t590+t524*t561+y*t564+t493*t562;
      t594 = t492*t585+t507*t586+t613;
      t593 = t525*t552+t519*t557+t616*t521;
      t591 = y*t521;
      t587 = x*t487;
      t506 = RATIONAL(1.0,24.0);
      t479 = t506*t525;
      t477 = t506*t526;
      t514 = RATIONAL(-1.0,6.0);
      t480 = t514*t526;
      t483 = t514*t525;
      t510 = RATIONAL(-5.0,6.0);
      t572 = t510*t522;
      t571 = RATIONAL(1.0,96.0)*t563+RATIONAL(5.0,96.0)*t584;
      t497 = RATIONAL(-1.0,24.0);
      t569 = t497*t563+t521*t572;
      t567 = t525*t489;
      t548 = t487*t523;
      t547 = t522*t489;
      t495 = RATIONAL(5.0,24.0);
      t546 = t495*t557+t480+t569+t615*t521;
      t499 = RATIONAL(-5.0,96.0);
      t502 = RATIONAL(-1.0,96.0);
      t545 = t499*t557+t497*t521+t477+t502*t574+t571;
      t544 = t499*t574+t497*t522+t479+t502*t557+t571;
      t543 = t513*t557+t512*t522+t495*t574+t483+t569;
      t505 = RATIONAL(-1.0,12.0);
      t500 = RATIONAL(-1.0,48.0);
      t503 = RATIONAL(5.0,48.0);
      t542 = t505+t503*t521+t500*t526;
      t501 = RATIONAL(1.0,48.0);
      t504 = RATIONAL(1.0,12.0);
      t498 = RATIONAL(-5.0,48.0);
      t541 = t504+t501*t526+t498*t521;
      t540 = t500*t525+t503*t522+t505;
      t539 = t501*t525+t498*t522+t504;
      t538 = t612+t521*t565+t488*t562+(t552+t556)*y;
      t537 = t496*t590+t508*t591+t484*t562+t521*t578+t612;
      t536 = t525*t583+t508*t589+t522*t553+x*t555+t611;
      t535 = t522*t587+t523*t547+x*t567+t525*t548+t606;
      t534 = t525*t587+t522*t548+t523*t567+x*t547+t606;
      t533 = t496*t589+x*t554+t525*t553+t484*t581+t611;
      t532 = t517*t591+t526*t558+t518*t590+t521*t580+t593;
      t531 = t521*t550+t520*t591+t519*t590+t524*t552+t593;
      t516 = RATIONAL(-2.0,3.0);
      t515 = RATIONAL(5.0,6.0);
      t511 = RATIONAL(1.0,4.0);
      t509 = RATIONAL(-5.0,4.0);
      t494 = RATIONAL(-5.0,24.0);
      t491 = RATIONAL(-5.0,16.0);
      coeffs_I->coeff_m2_m2 = t535+t603+t610;
      coeffs_I->coeff_m1_m2 = t536+t595+t605;
      coeffs_I->coeff_0_m2 = t542*t524+t541*y+t544;
      coeffs_I->coeff_p1_m2 = t533+t594+t595;
      coeffs_I->coeff_p2_m2 = t534+t609+t610;
      coeffs_I->coeff_m2_m1 = t537+t601+t602;
      coeffs_I->coeff_m1_m1 = t532+t598+t600;
      coeffs_I->coeff_0_m1 = (t513+t477+t494*t521)*t524+(t516+t515*t521+t480)*y
+t543;
      coeffs_I->coeff_p1_m1 = t532+t607+t608;
      coeffs_I->coeff_p2_m1 = t537+t596+t597;
      coeffs_I->coeff_m2_0 = t540*t523+t539*x+t545;
      coeffs_I->coeff_m1_0 = (t479+t494*t522+t513)*t523+(t483+t516+t515*t522)*x
+t546;
      coeffs_I->coeff_0_0 = t511*t525+RATIONAL(1.0,1.0)+t509*t522+(t491*t522+
t511+RATIONAL(1.0,16.0)*t525)*t526+(t491*t525+RATIONAL(25.0,16.0)*t522+t509)*
t521;
      coeffs_I->coeff_p1_0 = (t497*t525+t495*t522+t514)*t523+(t572+t615)*x+t546
;
      coeffs_I->coeff_p2_0 = t539*t523+t540*x+t545;
      coeffs_I->coeff_m2_p1 = t538+t596+t601;
      coeffs_I->coeff_m1_p1 = t531+t598+t607;
      coeffs_I->coeff_0_p1 = (t497*t526+t495*t521+t514)*t524+(t512+t510*t521+
t513*t526)*y+t543;
      coeffs_I->coeff_p1_p1 = t531+t600+t608;
      coeffs_I->coeff_p2_p1 = t538+t597+t602;
      coeffs_I->coeff_m2_p2 = t535+t599+t609;
      coeffs_I->coeff_m1_p2 = t536+t594+t604;
      coeffs_I->coeff_0_p2 = t541*t524+t542*y+t544;
      coeffs_I->coeff_p1_p2 = t533+t604+t605;
      coeffs_I->coeff_p2_p2 = t534+t599+t603;