aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dy.compute.c
blob: 763a422220cc70134c34ea29a89d0fb1f759d1e0 (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
fp t1155;
fp t1170;
fp t1172;
fp t1120;
fp t1156;
fp t1245;
fp t1132;
fp t1131;
fp t1244;
fp t1243;
fp t1163;
fp t1122;
fp t1159;
fp t1242;
fp t1169;
fp t1241;
fp t1130;
fp t1133;
fp t1240;
fp t1143;
fp t1142;
fp t1196;
fp t1239;
fp t1171;
fp t1173;
fp t1213;
fp t1144;
fp t1146;
fp t1105;
fp t1238;
fp t1200;
fp t1145;
fp t1117;
fp t1237;
fp t1166;
fp t1205;
fp t1134;
fp t1148;
fp t1201;
fp t1236;
fp t1215;
fp t1135;
fp t1198;
fp t1235;
fp t1112;
fp t1217;
fp t1234;
fp t1214;
fp t1212;
fp t1233;
fp t1161;
fp t1232;
fp t1231;
fp t1207;
fp t1199;
fp t1230;
fp t1195;
fp t1210;
fp t1149;
fp t1106;
fp t1229;
fp t1203;
fp t1228;
fp t1218;
fp t1147;
fp t1111;
fp t1104;
fp t1227;
fp t1226;
fp t1141;
fp t1160;
fp t1225;
fp t1114;
fp t1108;
fp t1224;
fp t1204;
fp t1223;
fp t1118;
fp t1139;
fp t1164;
fp t1208;
fp t1222;
fp t1221;
fp t1107;
fp t1220;
fp t1219;
fp t1216;
fp t1168;
fp t1128;
fp t1197;
fp t1211;
fp t1119;
fp t1193;
fp t1209;
fp t1206;
fp t1165;
fp t1202;
fp t1194;
fp t1192;
fp t1140;
fp t1191;
fp t1190;
fp t1150;
fp t1189;
fp t1152;
fp t1188;
fp t1187;
fp t1186;
fp t1185;
fp t1136;
fp t1137;
fp t1184;
fp t1183;
fp t1182;
fp t1181;
fp t1180;
fp t1179;
fp t1178;
fp t1177;
fp t1154;
fp t1153;
      t1155 = RATIONAL(1.0,9.0);
      t1170 = x*x;
      t1172 = x*t1170;
      t1120 = t1155*t1172;
      t1156 = RATIONAL(-4.0,9.0);
      t1245 = t1156*x+t1120;
      t1132 = RATIONAL(-1.0,144.0);
      t1131 = RATIONAL(1.0,144.0);
      t1244 = t1132*t1172+t1131*x;
      t1243 = t1132*x+t1131*t1172;
      t1163 = RATIONAL(-1.0,9.0);
      t1122 = t1163*t1172;
      t1159 = RATIONAL(4.0,9.0);
      t1242 = t1122+t1159*x;
      t1169 = t1170*t1170;
      t1241 = t1156*t1170+t1155*t1169;
      t1130 = RATIONAL(-1.0,288.0);
      t1133 = RATIONAL(1.0,288.0);
      t1240 = t1130*t1169+t1133*t1170;
      t1143 = RATIONAL(-1.0,12.0);
      t1142 = RATIONAL(-1.0,48.0);
      t1196 = t1142*t1169;
      t1239 = t1196+t1143;
      t1171 = y*y;
      t1173 = y*t1171;
      t1213 = y*t1169;
      t1144 = RATIONAL(1.0,18.0);
      t1146 = RATIONAL(1.0,36.0);
      t1105 = t1146*t1170;
      t1238 = t1173*t1105+t1144*t1213;
      t1200 = t1173*t1155;
      t1145 = RATIONAL(1.0,72.0);
      t1117 = t1145*t1169;
      t1237 = t1170*t1200+y*t1117;
      t1166 = RATIONAL(-1.0,6.0);
      t1205 = t1166*t1171;
      t1134 = RATIONAL(1.0,24.0);
      t1148 = RATIONAL(-1.0,72.0);
      t1201 = t1169*t1171;
      t1236 = t1148*t1169+t1134*t1201+(t1144+t1205)*t1170;
      t1215 = x*t1171;
      t1135 = RATIONAL(1.0,48.0);
      t1198 = t1171*t1172;
      t1235 = t1142*t1215+t1135*t1198+t1244;
      t1112 = t1145*t1172;
      t1217 = x*t1173;
      t1234 = t1173*t1112+t1148*t1217+t1244*y;
      t1214 = y*RATIONAL(8.0,9.0);
      t1212 = y*RATIONAL(-2.0,9.0);
      t1233 = t1170*t1214+t1169*t1212+t1241*t1173;
      t1161 = RATIONAL(1.0,3.0);
      t1232 = t1143*t1198+t1161*t1215+t1245;
      t1231 = x*t1214+t1172*t1212+t1245*t1173;
      t1207 = t1173*t1170;
      t1199 = t1173*t1169;
      t1230 = t1131*t1199+t1132*t1207+t1240*y;
      t1195 = t1135*t1169;
      t1210 = t1170*t1171;
      t1149 = RATIONAL(-1.0,36.0);
      t1106 = t1149*t1169;
      t1229 = t1142*t1210+t1171*t1195+t1105+t1106;
      t1203 = t1135*t1171;
      t1228 = t1171*t1196+t1146*t1169+(t1203+t1149)*t1170;
      t1218 = x*y;
      t1147 = RATIONAL(-1.0,18.0);
      t1111 = t1147*x;
      t1104 = t1144*t1172;
      t1227 = y*t1122+t1155*t1218+(t1111+t1104)*t1173;
      t1226 = RATIONAL(2.0,9.0)*y*t1172+RATIONAL(-8.0,9.0)*t1218+t1242*t1173;
      t1141 = RATIONAL(1.0,12.0);
      t1160 = RATIONAL(-1.0,3.0);
      t1225 = t1160*t1215+t1141*t1198+t1242;
      t1114 = t1147*t1172;
      t1108 = t1144*x;
      t1224 = y*t1120+t1163*t1218+(t1108+t1114)*t1173;
      t1204 = t1173*t1172;
      t1223 = t1149*t1204+x*t1200+(t1112+t1111)*y;
      t1118 = t1147*t1170;
      t1139 = RATIONAL(-1.0,24.0);
      t1164 = RATIONAL(1.0,6.0);
      t1208 = t1164*t1171;
      t1222 = t1118+t1139*t1201+t1170*t1208+t1117;
      t1221 = t1142*t1198+x*t1203+t1243;
      t1107 = t1148*t1172;
      t1220 = t1146*t1204+t1163*t1217+(t1108+t1107)*y;
      t1219 = t1173*t1107+t1145*t1217+t1243*y;
      t1216 = RATIONAL(-5.0,3.0)*y;
      t1168 = RATIONAL(-2.0,3.0);
      t1128 = t1168*t1173;
      t1197 = t1166*t1169;
      t1211 = t1170*t1216+t1173*t1197;
      t1119 = t1164*t1173;
      t1193 = RATIONAL(5.0,48.0)*t1170;
      t1209 = t1134*t1199+y*t1193;
      t1206 = y*t1118+t1173*t1106;
      t1165 = RATIONAL(2.0,3.0);
      t1202 = t1165*t1173;
      t1194 = RATIONAL(5.0,6.0)*t1170;
      t1192 = t1111+t1139*t1198+t1206;
      t1140 = RATIONAL(-5.0,48.0);
      t1191 = t1139*t1207+t1140*t1213+t1209;
      t1190 = t1134*t1198+t1108+t1206;
      t1150 = RATIONAL(5.0,12.0);
      t1189 = t1170*t1202+t1150*t1213+t1211;
      t1152 = RATIONAL(-5.0,24.0);
      t1188 = t1141*t1173+t1152*y;
      t1187 = RATIONAL(5.0,24.0)*y+t1143*t1173;
      t1186 = t1128+t1173*t1194+RATIONAL(4.0,3.0)*y+t1161*t1213+t1211;
      t1185 = t1119+t1152*t1207+t1209+t1239*y;
      t1136 = RATIONAL(1.0,96.0);
      t1137 = RATIONAL(-1.0,96.0);
      t1184 = t1137*t1201+t1130*t1170+t1133*t1169+t1136*t1210+t1230;
      t1183 = t1134*t1215+t1104+t1192+t1238;
      t1182 = t1136*t1201+t1137*t1210+t1230+t1240;
      t1181 = t1139*t1215+t1114+t1190+t1238;
      t1180 = t1112+x*t1208+t1192+t1237;
      t1179 = t1159*t1170+t1163*t1169+t1160*t1210+t1141*t1201+t1233;
      t1178 = x*t1205+t1107+t1190+t1237;
      t1177 = t1161*t1210+t1143*t1201+t1233+t1241;
      t1154 = RATIONAL(1.0,4.0);
      t1153 = RATIONAL(-5.0,8.0);
      coeffs_dy->coeff_m2_m2 = t1184+t1219+t1235;
      coeffs_dy->coeff_m1_m2 = t1180+t1220+t1236;
      coeffs_dy->coeff_0_m2 = t1140*t1170+t1195+t1141+(RATIONAL(-1.0,16.0)*
t1169+RATIONAL(5.0,16.0)*t1170+RATIONAL(-1.0,4.0))*t1171+t1185;
      coeffs_dy->coeff_p1_m2 = t1178+t1223+t1236;
      coeffs_dy->coeff_p2_m2 = t1184+t1221+t1234;
      coeffs_dy->coeff_m2_m1 = t1183+t1227+t1229;
      coeffs_dy->coeff_m1_m1 = t1177+t1225+t1226;
      coeffs_dy->coeff_0_m1 = t1194+t1197+t1168+(RATIONAL(1.0,2.0)+t1153*t1170+
RATIONAL(1.0,8.0)*t1169)*t1171+t1186;
      coeffs_dy->coeff_p1_m1 = t1177+t1231+t1232;
      coeffs_dy->coeff_p2_m1 = t1181+t1224+t1229;
      coeffs_dy->coeff_m2_0 = t1187*t1172+t1188*x+t1191;
      coeffs_dy->coeff_m1_0 = (RATIONAL(-5.0,12.0)*y+t1119)*t1172+(t1128+
RATIONAL(5.0,3.0)*y)*x+t1189;
      coeffs_dy->coeff_0_0 = (RATIONAL(-5.0,4.0)*t1170+t1154*t1169+1.0)*t1173+(
RATIONAL(-5.0,2.0)+RATIONAL(25.0,8.0)*t1170+t1153*t1169)*y;
      coeffs_dy->coeff_p1_0 = (t1166*t1173+t1150*y)*t1172+(t1202+t1216)*x+t1189
;
      coeffs_dy->coeff_p2_0 = t1188*t1172+t1187*x+t1191;
      coeffs_dy->coeff_m2_p1 = t1181+t1227+t1228;
      coeffs_dy->coeff_m1_p1 = t1179+t1226+t1232;
      coeffs_dy->coeff_0_p1 = t1164*t1169+t1165+RATIONAL(-5.0,6.0)*t1170+(
RATIONAL(5.0,8.0)*t1170+RATIONAL(-1.0,8.0)*t1169+RATIONAL(-1.0,2.0))*t1171+
t1186;
      coeffs_dy->coeff_p1_p1 = t1179+t1225+t1231;
      coeffs_dy->coeff_p2_p1 = t1183+t1224+t1228;
      coeffs_dy->coeff_m2_p2 = t1182+t1219+t1221;
      coeffs_dy->coeff_m1_p2 = t1178+t1220+t1222;
      coeffs_dy->coeff_0_p2 = t1193+(RATIONAL(1.0,16.0)*t1169+t1154+RATIONAL(
-5.0,16.0)*t1170)*t1171+t1185+t1239;
      coeffs_dy->coeff_p1_p2 = t1180+t1222+t1223;
      coeffs_dy->coeff_p2_p2 = t1182+t1234+t1235;