aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-I.compute.c
blob: 0eefdd1cb0ba801c9f9989f62bd6febeed744e75 (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
fp t652;
fp t654;
fp t674;
fp t672;
fp t771;
fp t666;
fp t667;
fp t770;
fp t675;
fp t677;
fp t655;
fp t676;
fp t708;
fp t747;
fp t769;
fp t699;
fp t651;
fp t742;
fp t768;
fp t636;
fp t707;
fp t653;
fp t720;
fp t638;
fp t732;
fp t673;
fp t719;
fp t656;
fp t735;
fp t767;
fp t724;
fp t700;
fp t766;
fp t662;
fp t713;
fp t726;
fp t744;
fp t663;
fp t746;
fp t765;
fp t722;
fp t743;
fp t739;
fp t764;
fp t741;
fp t763;
fp t737;
fp t635;
fp t640;
fp t745;
fp t701;
fp t702;
fp t762;
fp t730;
fp t761;
fp t740;
fp t733;
fp t721;
fp t728;
fp t760;
fp t705;
fp t717;
fp t759;
fp t664;
fp t738;
fp t706;
fp t665;
fp t758;
fp t757;
fp t756;
fp t714;
fp t755;
fp t754;
fp t753;
fp t729;
fp t752;
fp t703;
fp t712;
fp t710;
fp t751;
fp t750;
fp t749;
fp t748;
fp t671;
fp t711;
fp t646;
fp t715;
fp t736;
fp t734;
fp t731;
fp t727;
fp t723;
fp t658;
fp t718;
fp t709;
fp t668;
fp t632;
fp t641;
fp t629;
fp t633;
fp t627;
fp t698;
fp t697;
fp t644;
fp t643;
fp t696;
fp t695;
fp t694;
fp t642;
fp t647;
fp t648;
fp t693;
fp t645;
fp t649;
fp t650;
fp t692;
fp t691;
fp t690;
fp t689;
fp t688;
fp t687;
fp t686;
fp t685;
fp t684;
fp t683;
fp t682;
fp t670;
fp t669;
fp t661;
fp t660;
fp t659;
fp t657;
      t652 = RATIONAL(1.0,72.0);
      t654 = RATIONAL(-1.0,18.0);
      t674 = x*x;
      t672 = t674*t674;
      t771 = t654*t674+t652*t672;
      t666 = RATIONAL(1.0,6.0);
      t667 = RATIONAL(2.0,3.0);
      t770 = t666*t672+t667;
      t675 = y*y;
      t677 = t675*y;
      t655 = RATIONAL(-1.0,72.0);
      t676 = t674*x;
      t708 = t655*t676;
      t747 = x*y;
      t769 = t677*t708+t654*t747;
      t699 = t652*t676;
      t651 = RATIONAL(1.0,18.0);
      t742 = x*t651;
      t768 = t677*t699+y*t742;
      t636 = RATIONAL(1.0,144.0);
      t707 = t636*t674;
      t653 = RATIONAL(1.0,36.0);
      t720 = t653*t672;
      t638 = RATIONAL(-1.0,144.0);
      t732 = t674*t675;
      t673 = t675*t675;
      t719 = t673*t672;
      t656 = RATIONAL(-1.0,36.0);
      t735 = t638*t719+t656*t732;
      t767 = t673*t707+t735+t675*t720;
      t724 = t672*t675;
      t700 = t653*t674;
      t766 = t636*t724+t673*t700+t735;
      t662 = RATIONAL(1.0,9.0);
      t713 = t677*t662;
      t726 = t656*t672;
      t744 = y*t674;
      t663 = RATIONAL(-4.0,9.0);
      t746 = y*t672;
      t765 = t674*t713+t662*t746+t663*t744+t677*t726;
      t722 = t677*t676;
      t743 = y*t676;
      t739 = x*t663;
      t764 = x*t713+t662*t743+t656*t722+y*t739;
      t741 = t677*x;
      t763 = t651*t743+t652*t741+t769;
      t737 = x*t673;
      t635 = RATIONAL(-1.0,288.0);
      t640 = RATIONAL(1.0,288.0);
      t745 = x*t675;
      t701 = t635*t673;
      t702 = t640*t676;
      t762 = t675*t702+t635*t745+t676*t701+t640*t737;
      t730 = t673*t674;
      t761 = (t730+t724)*RATIONAL(-1.0,576.0)+(t719+t732)*RATIONAL(1.0,576.0);
      t740 = x*t636;
      t733 = t677*t638;
      t721 = t636*t676;
      t728 = t638*t676;
      t760 = x*t733+t677*t721+(t728+t740)*y;
      t705 = t655*t673;
      t717 = t675*t676;
      t759 = t673*t699+x*t705+t675*t742+t654*t717;
      t664 = RATIONAL(4.0,9.0);
      t738 = y*t664;
      t706 = t653*t676;
      t665 = RATIONAL(-1.0,9.0);
      t758 = x*t738+t677*t706+(t743+t741)*t665;
      t757 = t673*t702+t635*t717+x*t701+t640*t745;
      t756 = t677*t740+y*t721+(t722+t747)*t638;
      t714 = t665*t674;
      t755 = t665*t746+t674*t738+(t714+t720)*t677;
      t754 = t654*t745+t676*t705+t652*t737+t651*t717;
      t753 = t655*t741+t654*t743+t768;
      t729 = t656*t673;
      t752 = t673*t721+t653*t745+x*t729+t638*t717;
      t703 = t664*t675;
      t712 = t653*t673;
      t710 = t665*t675;
      t751 = t674*t703+t673*t714+(t712+t710)*t672;
      t750 = y*t699+t651*t741+t769;
      t749 = t636*t717+t673*t728+x*t712+t656*t745;
      t748 = y*t708+t654*t741+t768;
      t671 = RATIONAL(-5.0,6.0);
      t711 = t671*t674;
      t646 = RATIONAL(-1.0,24.0);
      t715 = t646*t672;
      t736 = t673*t715+t675*t711;
      t734 = t677*t672;
      t731 = t651*t674;
      t727 = t640*t672;
      t723 = RATIONAL(1.0,96.0)*t719+RATIONAL(5.0,96.0)*t732;
      t658 = RATIONAL(5.0,24.0);
      t718 = t658*t674;
      t709 = t640*t674;
      t668 = RATIONAL(-1.0,6.0);
      t632 = t668*t672;
      t641 = RATIONAL(1.0,24.0);
      t629 = t641*t672;
      t633 = t668*t673;
      t627 = t641*t673;
      t698 = t677*t635;
      t697 = t736+t673*t718+t633+t770*t675;
      t644 = RATIONAL(-1.0,96.0);
      t643 = RATIONAL(-5.0,96.0);
      t696 = t627+t646*t675+t643*t730+t644*t724+t723;
      t695 = t646*t674+t644*t730+t643*t724+t629+t723;
      t694 = t658*t724+t667*t674+t666*t730+t632+t736;
      t642 = RATIONAL(1.0,48.0);
      t647 = RATIONAL(-5.0,48.0);
      t648 = RATIONAL(1.0,12.0);
      t693 = t647*t675+t642*t673+t648;
      t645 = RATIONAL(5.0,48.0);
      t649 = RATIONAL(-1.0,48.0);
      t650 = RATIONAL(-1.0,12.0);
      t692 = t650+t649*t672+t645*t674;
      t691 = t642*t672+t648+t647*t674;
      t690 = t650+t649*t673+t645*t675;
      t689 = t677*t731+t655*t734+t766+t771*y;
      t688 = t767+t636*t734+t674*t733+(t700+t726)*y;
      t687 = t676*t710+x*t703+t673*t706+t665*t737+t751;
      t686 = t751+t675*t739+t676*t729+(t737+t717)*t662;
      t685 = t677*t727+t674*t698+y*t709+t635*t746+t761;
      t684 = t677*t707+t672*t733+t656*t744+y*t720+t767;
      t683 = t677*t709+y*t727+t672*t698+t635*t744+t761;
      t682 = t766+y*t731+t655*t746+t771*t677;
      t670 = RATIONAL(-2.0,3.0);
      t669 = RATIONAL(5.0,6.0);
      t661 = RATIONAL(-5.0,4.0);
      t660 = RATIONAL(1.0,4.0);
      t659 = RATIONAL(-5.0,24.0);
      t657 = RATIONAL(-5.0,16.0);
      coeffs_I->coeff_m2_m2 = t683+t760+t762;
      coeffs_I->coeff_m1_m2 = t682+t750+t752;
      coeffs_I->coeff_0_m2 = t691*y+t692*t677+t696;
      coeffs_I->coeff_p1_m2 = t682+t748+t749;
      coeffs_I->coeff_p2_m2 = t683+t756+t757;
      coeffs_I->coeff_m2_m1 = t688+t759+t763;
      coeffs_I->coeff_m1_m1 = t686+t758+t765;
      coeffs_I->coeff_0_m1 = (t669*t674+t632+t670)*y+(t659*t674+t666+t629)*t677
+t697;
      coeffs_I->coeff_p1_m1 = t687+t764+t765;
      coeffs_I->coeff_p2_m1 = t688+t753+t754;
      coeffs_I->coeff_m2_0 = t690*t676+t693*x+t695;
      coeffs_I->coeff_m1_0 = (t666+t659*t675+t627)*t676+(t669*t675+t670+t633)*x
+t694;
      coeffs_I->coeff_0_0 = RATIONAL(1.0,1.0)+t660*t672+t661*t674+(RATIONAL(
25.0,16.0)*t674+t661+t657*t672)*t675+(t660+t657*t674+RATIONAL(1.0,16.0)*t672)*
t673;
      coeffs_I->coeff_p1_0 = (t658*t675+t646*t673+t668)*t676+(t667+t671*t675+
t666*t673)*x+t694;
      coeffs_I->coeff_p2_0 = t693*t676+t690*x+t695;
      coeffs_I->coeff_m2_p1 = t684+t753+t759;
      coeffs_I->coeff_m1_p1 = t686+t755+t764;
      coeffs_I->coeff_0_p1 = (t711+t770)*y+(t715+t718+t668)*t677+t697;
      coeffs_I->coeff_p1_p1 = t687+t755+t758;
      coeffs_I->coeff_p2_p1 = t684+t754+t763;
      coeffs_I->coeff_m2_p2 = t685+t756+t762;
      coeffs_I->coeff_m1_p2 = t689+t748+t752;
      coeffs_I->coeff_0_p2 = t692*y+t691*t677+t696;
      coeffs_I->coeff_p1_p2 = t689+t749+t750;
      coeffs_I->coeff_p2_p2 = t685+t757+t760;