aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dx.compute.c
blob: 46af5f250fe6e8ccd8bebfad7a4ff9345e329988 (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
fp t779;
fp t782;
fp t742;
fp t740;
fp t851;
fp t781;
fp t741;
fp t739;
fp t850;
fp t849;
fp t774;
fp t778;
fp t769;
fp t848;
fp t758;
fp t752;
fp t847;
fp t759;
fp t800;
fp t743;
fp t846;
fp t760;
fp t720;
fp t750;
fp t802;
fp t845;
fp t821;
fp t761;
fp t748;
fp t717;
fp t780;
fp t715;
fp t812;
fp t807;
fp t844;
fp t745;
fp t714;
fp t770;
fp t843;
fp t811;
fp t754;
fp t804;
fp t842;
fp t819;
fp t798;
fp t771;
fp t735;
fp t824;
fp t841;
fp t815;
fp t840;
fp t809;
fp t839;
fp t731;
fp t801;
fp t838;
fp t777;
fp t810;
fp t822;
fp t823;
fp t837;
fp t818;
fp t836;
fp t757;
fp t776;
fp t820;
fp t772;
fp t803;
fp t835;
fp t834;
fp t833;
fp t744;
fp t832;
fp t831;
fp t830;
fp t766;
fp t799;
fp t716;
fp t829;
fp t765;
fp t828;
fp t827;
fp t718;
fp t826;
fp t825;
fp t816;
fp t817;
fp t813;
fp t773;
fp t806;
fp t805;
fp t751;
fp t797;
fp t764;
fp t749;
fp t796;
fp t746;
fp t795;
fp t767;
fp t768;
fp t794;
fp t753;
fp t793;
fp t792;
fp t791;
fp t790;
fp t789;
fp t788;
fp t787;
fp t786;
fp t763;
fp t762;
fp t747;
      t779 = y*y;
      t782 = t779*t779;
      t742 = RATIONAL(1.0,288.0);
      t740 = RATIONAL(-1.0,288.0);
      t851 = t740*t782+t742*t779;
      t781 = t779*y;
      t741 = RATIONAL(-1.0,144.0);
      t739 = RATIONAL(1.0,144.0);
      t850 = t741*t781+t739*y;
      t849 = t741*y+t739*t781;
      t774 = RATIONAL(1.0,3.0);
      t778 = x*x;
      t769 = RATIONAL(-4.0,9.0);
      t848 = t769+t774*t778;
      t758 = RATIONAL(-1.0,12.0);
      t752 = RATIONAL(-1.0,48.0);
      t847 = t752*t782+t758;
      t759 = RATIONAL(1.0,24.0);
      t800 = t759*t778;
      t743 = RATIONAL(1.0,18.0);
      t846 = t743*y+t781*t800;
      t760 = RATIONAL(-1.0,18.0);
      t720 = t760*y;
      t750 = RATIONAL(-1.0,24.0);
      t802 = t778*t750;
      t845 = t720+t781*t802;
      t821 = x*t743;
      t761 = RATIONAL(1.0,36.0);
      t748 = RATIONAL(-1.0,36.0);
      t717 = t748*t782;
      t780 = t778*x;
      t715 = t760*t779;
      t812 = t780*t717+x*t715;
      t807 = t779*t780;
      t844 = t761*t807+t812+t782*t821;
      t745 = RATIONAL(1.0,72.0);
      t714 = t745*t782;
      t770 = RATIONAL(1.0,9.0);
      t843 = x*t714+t812+t770*t807;
      t811 = t778*t782;
      t754 = RATIONAL(1.0,48.0);
      t804 = t778*t754;
      t842 = t761*t782+t752*t811+(t804+t748)*t779;
      t819 = y*t780;
      t798 = t780*t743;
      t771 = RATIONAL(-1.0,9.0);
      t735 = t771*t781;
      t824 = x*y;
      t841 = t760*t819+t770*t824+t781*t798+x*t735;
      t815 = t780*t782;
      t840 = t741*t807+t739*t815+t851*x;
      t809 = t778*t779;
      t839 = t717+t761*t779+t782*t804+t752*t809;
      t731 = t770*t781;
      t801 = t781*t780;
      t838 = y*t798+t760*t801+t771*t824+x*t731;
      t777 = RATIONAL(-2.0,9.0);
      t810 = t780*t770;
      t822 = x*RATIONAL(8.0,9.0);
      t823 = x*t782;
      t837 = t779*t822+t769*t807+t777*t823+t782*t810;
      t818 = x*t781;
      t836 = y*t822+t769*t819+t777*t818+t770*t801;
      t757 = RATIONAL(1.0,12.0);
      t776 = RATIONAL(-1.0,3.0);
      t820 = y*t778;
      t772 = RATIONAL(4.0,9.0);
      t803 = t778*t781;
      t835 = t757*t803+t776*t820+t735+t772*y;
      t834 = t731+t758*t803+t848*y;
      t833 = t772*t819+RATIONAL(2.0,9.0)*t818+t771*t801+RATIONAL(-8.0,9.0)*t824
;
      t744 = RATIONAL(-1.0,72.0);
      t832 = t744*t819+t745*t801+t850*x;
      t831 = t754*t803+t752*t820+t850;
      t830 = t744*t801+t745*t819+t849*x;
      t766 = RATIONAL(-1.0,6.0);
      t799 = t778*t766;
      t716 = t744*t781;
      t829 = y*t799+t716+t846;
      t765 = RATIONAL(1.0,6.0);
      t828 = t765*t809+t714+t715+t782*t802;
      t827 = (t800+t744)*t782+(t743+t799)*t779;
      t718 = t745*t781;
      t826 = t718+t765*t820+t845;
      t825 = t752*t803+y*t804+t849;
      t816 = RATIONAL(5.0,48.0)*t779;
      t817 = t759*t815+x*t816;
      t813 = t766*t782;
      t773 = RATIONAL(-5.0,3.0);
      t806 = t780*t813+t773*x*t779;
      t805 = t778*RATIONAL(-1.0,96.0);
      t751 = RATIONAL(-5.0,48.0);
      t797 = t751*t823+t750*t807+t817;
      t764 = RATIONAL(2.0,3.0);
      t749 = RATIONAL(5.0,12.0);
      t796 = t749*t823+t764*t807+t806;
      t746 = RATIONAL(-5.0,24.0);
      t795 = t817+t765*t780+t746*t807+t847*x;
      t767 = RATIONAL(5.0,6.0);
      t768 = RATIONAL(-2.0,3.0);
      t794 = t774*t823+RATIONAL(4.0,3.0)*x+t768*t780+t767*t807+t806;
      t753 = RATIONAL(1.0,96.0);
      t793 = t753*t811+t779*t805+t840+t851;
      t792 = t760*t781+y*t802+t844+t846;
      t791 = y*t810+t748*t801+t843+(t718+t720)*x;
      t790 = t776*t809+t771*t782+t772*t779+t757*t811+t837;
      t789 = t761*t801+y*t821+x*t716+t771*t819+t843;
      t788 = t753*t809+t740*t779+t840+(t742+t805)*t782;
      t787 = t837+t758*t811+t770*t782+t848*t779;
      t786 = t743*t781+y*t800+t844+t845;
      t763 = RATIONAL(1.0,4.0);
      t762 = RATIONAL(-5.0,8.0);
      t747 = RATIONAL(5.0,24.0);
      coeffs_dx->coeff_m2_m2 = t788+t830+t831;
      coeffs_dx->coeff_m1_m2 = t786+t839+t841;
      coeffs_dx->coeff_0_m2 = (t757*y+t758*t781)*t780+(t746*y+t747*t781)*x+t797
;
      coeffs_dx->coeff_p1_m2 = t792+t841+t842;
      coeffs_dx->coeff_p2_m2 = t793+t825+t830;
      coeffs_dx->coeff_m2_m1 = t789+t826+t827;
      coeffs_dx->coeff_m1_m1 = t787+t833+t835;
      coeffs_dx->coeff_0_m1 = (t765*t781+t768*y)*t780+(RATIONAL(5.0,3.0)*y+
RATIONAL(-5.0,12.0)*t781)*x+t796;
      coeffs_dx->coeff_p1_m1 = t790+t833+t834;
      coeffs_dx->coeff_p2_m1 = t789+t828+t829;
      coeffs_dx->coeff_m2_0 = t757+t754*t782+t751*t779+(RATIONAL(5.0,16.0)*t779
+RATIONAL(-1.0,4.0)+RATIONAL(-1.0,16.0)*t782)*t778+t795;
      coeffs_dx->coeff_m1_0 = t768+t813+t767*t779+(RATIONAL(1.0,2.0)+t762*t779+
RATIONAL(1.0,8.0)*t782)*t778+t794;
      coeffs_dx->coeff_0_0 = (1.0+t763*t782+RATIONAL(-5.0,4.0)*t779)*t780+(
RATIONAL(25.0,8.0)*t779+RATIONAL(-5.0,2.0)+t762*t782)*x;
      coeffs_dx->coeff_p1_0 = RATIONAL(-5.0,6.0)*t779+t764+t765*t782+(RATIONAL(
5.0,8.0)*t779+RATIONAL(-1.0,2.0)+RATIONAL(-1.0,8.0)*t782)*t778+t794;
      coeffs_dx->coeff_p2_0 = t816+(t763+RATIONAL(-5.0,16.0)*t779+RATIONAL(1.0,
16.0)*t782)*t778+t795+t847;
      coeffs_dx->coeff_m2_p1 = t791+t827+t829;
      coeffs_dx->coeff_m1_p1 = t787+t834+t836;
      coeffs_dx->coeff_0_p1 = (t766*t781+t764*y)*t780+(t773*y+t749*t781)*x+t796
;
      coeffs_dx->coeff_p1_p1 = t790+t835+t836;
      coeffs_dx->coeff_p2_p1 = t791+t826+t828;
      coeffs_dx->coeff_m2_p2 = t788+t825+t832;
      coeffs_dx->coeff_m1_p2 = t792+t838+t839;
      coeffs_dx->coeff_0_p2 = (t757*t781+t758*y)*t780+(t747*y+t746*t781)*x+t797
;
      coeffs_dx->coeff_p1_p2 = t786+t838+t842;
      coeffs_dx->coeff_p2_p2 = t793+t831+t832;