aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dx.compute.c
blob: 74e206a16755309be773a89d334990bc10a365e2 (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
fp t935;
fp t933;
fp t894;
fp t897;
fp t1007;
fp t920;
fp t934;
fp t936;
fp t1006;
fp t937;
fp t895;
fp t896;
fp t1005;
fp t1004;
fp t932;
fp t1003;
fp t906;
fp t907;
fp t1002;
fp t899;
fp t956;
fp t913;
fp t1001;
fp t911;
fp t873;
fp t955;
fp t903;
fp t1000;
fp t908;
fp t879;
fp t898;
fp t999;
fp t919;
fp t964;
fp t909;
fp t870;
fp t965;
fp t871;
fp t971;
fp t998;
fp t910;
fp t874;
fp t972;
fp t997;
fp t927;
fp t974;
fp t954;
fp t889;
fp t996;
fp t961;
fp t995;
fp t922;
fp t994;
fp t923;
fp t973;
fp t924;
fp t892;
fp t905;
fp t993;
fp t975;
fp t992;
fp t976;
fp t978;
fp t991;
fp t925;
fp t990;
fp t968;
fp t989;
fp t988;
fp t928;
fp t960;
fp t958;
fp t987;
fp t912;
fp t930;
fp t986;
fp t985;
fp t984;
fp t876;
fp t983;
fp t869;
fp t982;
fp t981;
fp t980;
fp t979;
fp t966;
fp t970;
fp t963;
fp t921;
fp t959;
fp t957;
fp t953;
fp t929;
fp t914;
fp t952;
fp t904;
fp t951;
fp t916;
fp t950;
fp t949;
fp t948;
fp t947;
fp t946;
fp t945;
fp t944;
fp t943;
fp t942;
fp t941;
fp t918;
fp t917;
fp t915;
      t935 = y*y;
      t933 = t935*t935;
      t894 = RATIONAL(-1.0,288.0);
      t897 = RATIONAL(1.0,288.0);
      t1007 = t897*t935+t894*t933;
      t920 = RATIONAL(-4.0,9.0);
      t934 = x*x;
      t936 = t934*x;
      t1006 = t920*t936+x*RATIONAL(8.0,9.0);
      t937 = y*t935;
      t895 = RATIONAL(1.0,144.0);
      t896 = RATIONAL(-1.0,144.0);
      t1005 = t896*t937+t895*y;
      t1004 = t895*t937+t896*y;
      t932 = RATIONAL(-2.0,3.0);
      t1003 = t932+RATIONAL(5.0,6.0)*t935;
      t906 = RATIONAL(-1.0,48.0);
      t907 = RATIONAL(-1.0,12.0);
      t1002 = t906*t933+t907;
      t899 = RATIONAL(1.0,48.0);
      t956 = t934*t899;
      t913 = RATIONAL(-1.0,36.0);
      t1001 = t913+t956;
      t911 = RATIONAL(-1.0,18.0);
      t873 = t911*y;
      t955 = t934*t937;
      t903 = RATIONAL(-1.0,24.0);
      t1000 = t873+t903*t955;
      t908 = RATIONAL(1.0,18.0);
      t879 = t908*y;
      t898 = RATIONAL(1.0,24.0);
      t999 = t898*t955+t879;
      t919 = RATIONAL(1.0,9.0);
      t964 = t919*t936;
      t909 = RATIONAL(1.0,72.0);
      t870 = t909*t933;
      t965 = t933*t936;
      t871 = t911*t935;
      t971 = x*t871+t913*t965;
      t998 = t935*t964+x*t870+t971;
      t910 = RATIONAL(1.0,36.0);
      t874 = t910*t935;
      t972 = x*t933;
      t997 = t936*t874+t971+t908*t972;
      t927 = RATIONAL(-1.0,9.0);
      t974 = y*t927;
      t954 = t936*t937;
      t889 = t919*t937;
      t996 = t911*t954+t936*t879+(t974+t889)*x;
      t961 = t934*t898;
      t995 = t908*t937+y*t961+t1000;
      t922 = RATIONAL(-2.0,9.0);
      t994 = t922*t972+t933*t964+t1006*t935;
      t923 = RATIONAL(4.0,9.0);
      t973 = y*t934;
      t924 = RATIONAL(-1.0,3.0);
      t892 = t927*t937;
      t905 = RATIONAL(1.0,12.0);
      t993 = t905*t955+t924*t973+t923*y+t892;
      t975 = x*t937;
      t992 = t922*t975+t919*t954+t1006*y;
      t976 = y*t936;
      t978 = x*y;
      t991 = t927*t954+t923*t976+RATIONAL(2.0,9.0)*t975+RATIONAL(-8.0,9.0)*t978
;
      t925 = RATIONAL(1.0,3.0);
      t990 = t907*t955+t925*t973+t920*y+t889;
      t968 = t935*t936;
      t989 = t895*t965+t896*t968+t1007*x;
      t988 = x*t892+t908*t954+t919*t978+t911*t976;
      t928 = RATIONAL(1.0,6.0);
      t960 = t934*t903;
      t958 = t934*t935;
      t987 = t928*t958+t870+t933*t960+t871;
      t912 = RATIONAL(-1.0,72.0);
      t930 = RATIONAL(-1.0,6.0);
      t986 = t908*t935+t930*t958+(t961+t912)*t933;
      t985 = t912*t954+t909*t976+t1004*x;
      t984 = t899*t955+t906*t973+t1005;
      t876 = t909*t937;
      t983 = y*t964+t913*t954+(t873+t876)*x;
      t869 = t912*t937;
      t982 = t910*t954+t936*t974+(t879+t869)*x;
      t981 = y*t960+t911*t937+t999;
      t980 = y*t956+t906*t955+t1004;
      t979 = t909*t954+t912*t976+t1005*x;
      t966 = RATIONAL(5.0,48.0)*t935;
      t970 = t898*t965+x*t966;
      t963 = t933*t934;
      t921 = RATIONAL(-5.0,3.0);
      t959 = t930*t965+t921*x*t935;
      t957 = t934*RATIONAL(1.0,96.0);
      t953 = t934*RATIONAL(-1.0,96.0);
      t929 = RATIONAL(2.0,3.0);
      t914 = RATIONAL(5.0,12.0);
      t952 = t929*t968+t914*t972+t959;
      t904 = RATIONAL(-5.0,48.0);
      t951 = t904*t972+t903*t968+t970;
      t916 = RATIONAL(-5.0,24.0);
      t950 = t970+t916*t968+t928*t936+t1002*x;
      t949 = t925*t972+RATIONAL(4.0,3.0)*x+t959+t1003*t936;
      t948 = t997+t874+t906*t958+t1001*t933;
      t947 = t930*t973+t869+t998+t999;
      t946 = t989+(t957+t894)*t935+(t953+t897)*t933;
      t945 = t997+t910*t933+t906*t963+t1001*t935;
      t944 = t905*t963+t923*t935+t927*t933+t924*t958+t994;
      t943 = t933*t957+t935*t953+t989+t1007;
      t942 = t925*t958+t919*t933+t920*t935+t907*t963+t994;
      t941 = t876+t928*t973+t998+t1000;
      t918 = RATIONAL(1.0,4.0);
      t917 = RATIONAL(-5.0,8.0);
      t915 = RATIONAL(5.0,24.0);
      coeffs_dx->coeff_m2_m2 = t946+t984+t985;
      coeffs_dx->coeff_m1_m2 = t948+t988+t995;
      coeffs_dx->coeff_0_m2 = (t907*t937+t905*y)*t936+(t916*y+t915*t937)*x+t951
;
      coeffs_dx->coeff_p1_m2 = t945+t981+t988;
      coeffs_dx->coeff_p2_m2 = t943+t980+t985;
      coeffs_dx->coeff_m2_m1 = t941+t982+t986;
      coeffs_dx->coeff_m1_m1 = t942+t991+t993;
      coeffs_dx->coeff_0_m1 = (t928*t937+t932*y)*t936+(RATIONAL(-5.0,12.0)*t937
+RATIONAL(5.0,3.0)*y)*x+t952;
      coeffs_dx->coeff_p1_m1 = t944+t990+t991;
      coeffs_dx->coeff_p2_m1 = t947+t982+t987;
      coeffs_dx->coeff_m2_0 = t899*t933+t905+t904*t935+(RATIONAL(-1.0,4.0)+
RATIONAL(-1.0,16.0)*t933+RATIONAL(5.0,16.0)*t935)*t934+t950;
      coeffs_dx->coeff_m1_0 = t930*t933+(t917*t935+RATIONAL(1.0,2.0)+RATIONAL(
1.0,8.0)*t933)*t934+t949+t1003;
      coeffs_dx->coeff_0_0 = (t918*t933+1.0+RATIONAL(-5.0,4.0)*t935)*t936+(
RATIONAL(-5.0,2.0)+t917*t933+RATIONAL(25.0,8.0)*t935)*x;
      coeffs_dx->coeff_p1_0 = t929+t928*t933+RATIONAL(-5.0,6.0)*t935+(RATIONAL(
-1.0,2.0)+RATIONAL(5.0,8.0)*t935+RATIONAL(-1.0,8.0)*t933)*t934+t949;
      coeffs_dx->coeff_p2_0 = t966+(RATIONAL(1.0,16.0)*t933+t918+RATIONAL(-5.0,
16.0)*t935)*t934+t950+t1002;
      coeffs_dx->coeff_m2_p1 = t947+t983+t986;
      coeffs_dx->coeff_m1_p1 = t942+t990+t992;
      coeffs_dx->coeff_0_p1 = (t929*y+t930*t937)*t936+(t921*y+t914*t937)*x+t952
;
      coeffs_dx->coeff_p1_p1 = t944+t992+t993;
      coeffs_dx->coeff_p2_p1 = t941+t983+t987;
      coeffs_dx->coeff_m2_p2 = t946+t979+t980;
      coeffs_dx->coeff_m1_p2 = t948+t981+t996;
      coeffs_dx->coeff_0_p2 = (t907*y+t905*t937)*t936+(t916*t937+t915*y)*x+t951
;
      coeffs_dx->coeff_p1_p2 = t945+t995+t996;
      coeffs_dx->coeff_p2_p2 = t943+t979+t984;