aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dy.compute.c
blob: 59006b6f31b3f41d20d06e5ad4c1d9726f1c0a4a (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
fp t1007;
fp t1006;
fp t1013;
fp t1015;
fp t964;
fp t1082;
fp t974;
fp t976;
fp t1081;
fp t1080;
fp t975;
fp t977;
fp t1017;
fp t1079;
fp t1004;
fp t1005;
fp t1078;
fp t987;
fp t993;
fp t1077;
fp t985;
fp t1014;
fp t1040;
fp t995;
fp t957;
fp t1076;
fp t994;
fp t978;
fp t952;
fp t1075;
fp t1052;
fp t996;
fp t1016;
fp t1039;
fp t983;
fp t1042;
fp t958;
fp t1048;
fp t1074;
fp t980;
fp t959;
fp t1038;
fp t1073;
fp t972;
fp t1054;
fp t1050;
fp t1043;
fp t1072;
fp t1071;
fp t1049;
fp t1010;
fp t1070;
fp t1053;
fp t1009;
fp t1069;
fp t992;
fp t1011;
fp t1068;
fp t1067;
fp t979;
fp t1066;
fp t989;
fp t1065;
fp t955;
fp t1064;
fp t1001;
fp t951;
fp t1063;
fp t1062;
fp t1061;
fp t1060;
fp t1059;
fp t1058;
fp t1000;
fp t1057;
fp t1056;
fp t1051;
fp t1055;
fp t1035;
fp t1008;
fp t1047;
fp t1045;
fp t1046;
fp t1041;
fp t1037;
fp t1036;
fp t1003;
fp t973;
fp t968;
fp t984;
fp t999;
fp t1034;
fp t986;
fp t1033;
fp t1032;
fp t981;
fp t1031;
fp t1030;
fp t1029;
fp t1028;
fp t988;
fp t990;
fp t1027;
fp t1026;
fp t1025;
fp t1024;
fp t1023;
fp t1022;
fp t1021;
fp t998;
fp t997;
      t1007 = RATIONAL(4.0,9.0);
      t1006 = RATIONAL(-1.0,9.0);
      t1013 = x*x;
      t1015 = t1013*x;
      t964 = t1006*t1015;
      t1082 = t964+t1007*x;
      t974 = RATIONAL(1.0,144.0);
      t976 = RATIONAL(-1.0,144.0);
      t1081 = t974*x+t976*t1015;
      t1080 = t974*t1015+t976*x;
      t975 = RATIONAL(-1.0,288.0);
      t977 = RATIONAL(1.0,288.0);
      t1017 = t1013*t1013;
      t1079 = t975*t1017+t977*t1013;
      t1004 = RATIONAL(-4.0,9.0);
      t1005 = RATIONAL(1.0,9.0);
      t1078 = t1004*t1013+t1005*t1017;
      t987 = RATIONAL(-1.0,48.0);
      t993 = RATIONAL(-1.0,12.0);
      t1077 = t987*t1017+t993;
      t985 = RATIONAL(-1.0,24.0);
      t1014 = y*y;
      t1040 = t1014*t1015;
      t995 = RATIONAL(-1.0,18.0);
      t957 = t995*x;
      t1076 = t985*t1040+t957;
      t994 = RATIONAL(1.0,24.0);
      t978 = RATIONAL(1.0,18.0);
      t952 = t978*x;
      t1075 = t994*t1040+t952;
      t1052 = y*t1017;
      t996 = RATIONAL(1.0,36.0);
      t1016 = t1014*y;
      t1039 = t1013*t1016;
      t983 = RATIONAL(-1.0,36.0);
      t1042 = t1017*t1016;
      t958 = t995*t1013;
      t1048 = t983*t1042+y*t958;
      t1074 = t978*t1052+t1048+t996*t1039;
      t980 = RATIONAL(1.0,72.0);
      t959 = t980*t1017;
      t1038 = t1005*t1016;
      t1073 = y*t959+t1048+t1013*t1038;
      t972 = t1005*t1015;
      t1054 = x*y;
      t1050 = x*t1016;
      t1043 = t1016*t1015;
      t1072 = t1006*t1054+y*t972+t978*t1050+t995*t1043;
      t1071 = RATIONAL(2.0,9.0)*y*t1015+RATIONAL(-8.0,9.0)*t1054+t1082*t1016;
      t1049 = y*RATIONAL(-2.0,9.0);
      t1010 = RATIONAL(8.0,9.0);
      t1070 = t1004*t1050+t1016*t972+t1010*t1054+t1015*t1049;
      t1053 = x*t1014;
      t1009 = RATIONAL(1.0,3.0);
      t1069 = t972+t1004*x+t993*t1040+t1009*t1053;
      t992 = RATIONAL(1.0,12.0);
      t1011 = RATIONAL(-1.0,3.0);
      t1068 = t1011*t1053+t992*t1040+t1082;
      t1067 = t995*t1050+y*t964+t978*t1043+t1005*t1054;
      t979 = RATIONAL(-1.0,72.0);
      t1066 = t980*t1043+t979*t1050+t1081*y;
      t989 = RATIONAL(1.0,48.0);
      t1065 = t987*t1053+t989*t1040+t1081;
      t955 = t980*t1015;
      t1064 = t983*t1043+x*t1038+(t957+t955)*y;
      t1001 = RATIONAL(-1.0,6.0);
      t951 = t979*t1015;
      t1063 = t1001*t1053+t951+t1075;
      t1062 = t996*t1043+t1006*t1050+(t952+t951)*y;
      t1061 = t979*t1043+t980*t1050+t1080*y;
      t1060 = t989*t1053+t987*t1040+t1080;
      t1059 = t978*t1015+t994*t1053+t1076;
      t1058 = t995*t1015+t985*t1053+t1075;
      t1000 = RATIONAL(1.0,6.0);
      t1057 = t955+t1000*t1053+t1076;
      t1056 = t974*t1042+t976*t1039+t1079*y;
      t1051 = y*t1013;
      t1055 = t1010*t1051+t1017*t1049+t1078*t1016;
      t1035 = t1001*t1017;
      t1008 = RATIONAL(-5.0,3.0);
      t1047 = t1008*t1051+t1016*t1035;
      t1045 = RATIONAL(5.0,48.0)*t1013;
      t1046 = t994*t1042+y*t1045;
      t1041 = t1017*t1014;
      t1037 = t1013*t1014;
      t1036 = RATIONAL(5.0,6.0)*t1013;
      t1003 = RATIONAL(-2.0,3.0);
      t973 = t1003*t1016;
      t968 = t1000*t1016;
      t984 = RATIONAL(5.0,12.0);
      t999 = RATIONAL(2.0,3.0);
      t1034 = t984*t1052+t999*t1039+t1047;
      t986 = RATIONAL(-5.0,48.0);
      t1033 = t986*t1052+t985*t1039+t1046;
      t1032 = RATIONAL(5.0,24.0)*y+t993*t1016;
      t981 = RATIONAL(-5.0,24.0);
      t1031 = t992*t1016+t981*y;
      t1030 = t973+t1009*t1052+t1016*t1036+RATIONAL(4.0,3.0)*y+t1047;
      t1029 = t968+t981*t1039+t1046+t1077*y;
      t1028 = t994*t1041+t1001*t1037+t978*t1013+t979*t1017+t1073;
      t988 = RATIONAL(1.0,96.0);
      t990 = RATIONAL(-1.0,96.0);
      t1027 = t988*t1037+t977*t1017+t975*t1013+t990*t1041+t1056;
      t1026 = t1007*t1013+t992*t1041+t1006*t1017+t1011*t1037+t1055;
      t1025 = t987*t1041+t983*t1013+t996*t1017+t989*t1037+t1074;
      t1024 = t993*t1041+t1009*t1037+t1055+t1078;
      t1023 = t987*t1037+t983*t1017+t996*t1013+t989*t1041+t1074;
      t1022 = t990*t1037+t988*t1041+t1056+t1079;
      t1021 = t959+t985*t1041+t1000*t1037+t958+t1073;
      t998 = RATIONAL(1.0,4.0);
      t997 = RATIONAL(-5.0,8.0);
      coeffs_dy->coeff_m2_m2 = t1027+t1061+t1065;
      coeffs_dy->coeff_m1_m2 = t1028+t1057+t1062;
      coeffs_dy->coeff_0_m2 = t986*t1013+t989*t1017+t992+(RATIONAL(5.0,16.0)*
t1013+RATIONAL(-1.0,4.0)+RATIONAL(-1.0,16.0)*t1017)*t1014+t1029;
      coeffs_dy->coeff_p1_m2 = t1028+t1063+t1064;
      coeffs_dy->coeff_p2_m2 = t1027+t1060+t1066;
      coeffs_dy->coeff_m2_m1 = t1023+t1059+t1067;
      coeffs_dy->coeff_m1_m1 = t1024+t1068+t1071;
      coeffs_dy->coeff_0_m1 = t1036+t1003+t1035+(t997*t1013+RATIONAL(1.0,2.0)+
RATIONAL(1.0,8.0)*t1017)*t1014+t1030;
      coeffs_dy->coeff_p1_m1 = t1024+t1069+t1070;
      coeffs_dy->coeff_p2_m1 = t1023+t1058+t1072;
      coeffs_dy->coeff_m2_0 = t1032*t1015+t1031*x+t1033;
      coeffs_dy->coeff_m1_0 = (RATIONAL(-5.0,12.0)*y+t968)*t1015+(RATIONAL(5.0,
3.0)*y+t973)*x+t1034;
      coeffs_dy->coeff_0_0 = (RATIONAL(-5.0,4.0)*t1013+t998*t1017+1.0)*t1016+(
t997*t1017+RATIONAL(25.0,8.0)*t1013+RATIONAL(-5.0,2.0))*y;
      coeffs_dy->coeff_p1_0 = (t1001*t1016+t984*y)*t1015+(t1008*y+t999*t1016)*x
+t1034;
      coeffs_dy->coeff_p2_0 = t1031*t1015+t1032*x+t1033;
      coeffs_dy->coeff_m2_p1 = t1025+t1058+t1067;
      coeffs_dy->coeff_m1_p1 = t1026+t1069+t1071;
      coeffs_dy->coeff_0_p1 = t1000*t1017+t999+RATIONAL(-5.0,6.0)*t1013+(
RATIONAL(5.0,8.0)*t1013+RATIONAL(-1.0,8.0)*t1017+RATIONAL(-1.0,2.0))*t1014+
t1030;
      coeffs_dy->coeff_p1_p1 = t1026+t1068+t1070;
      coeffs_dy->coeff_p2_p1 = t1025+t1059+t1072;
      coeffs_dy->coeff_m2_p2 = t1022+t1060+t1061;
      coeffs_dy->coeff_m1_p2 = t1021+t1062+t1063;
      coeffs_dy->coeff_0_p2 = t1045+(RATIONAL(1.0,16.0)*t1017+RATIONAL(-5.0,
16.0)*t1013+t998)*t1014+t1029+t1077;
      coeffs_dy->coeff_p1_p2 = t1021+t1057+t1064;
      coeffs_dy->coeff_p2_p2 = t1022+t1065+t1066;