aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxx.compute.c
blob: 15d553b69ff43c97223204bf7a3195d8c4cc96c2 (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
fp t1210;
fp t1181;
fp t1212;
fp t1182;
fp t1274;
fp t1185;
fp t1249;
fp t1180;
fp t1245;
fp t1273;
fp t1184;
fp t1244;
fp t1211;
fp t1248;
fp t1183;
fp t1106;
fp t1272;
fp t1271;
fp t1270;
fp t1250;
fp t1197;
fp t1252;
fp t1269;
fp t1268;
fp t1155;
fp t1267;
fp t1098;
fp t1253;
fp t1203;
fp t1266;
fp t1265;
fp t1264;
fp t1205;
fp t1246;
fp t1263;
fp t1262;
fp t1236;
fp t1240;
fp t1261;
fp t1247;
fp t1198;
fp t1260;
fp t1259;
fp t1258;
fp t1176;
fp t1186;
fp t1209;
fp t1242;
fp t1257;
fp t1175;
fp t1229;
fp t1256;
fp t1255;
fp t1237;
fp t1238;
fp t1193;
fp t1204;
fp t1206;
fp t1254;
fp t1251;
fp t1233;
fp t1243;
fp t1187;
fp t1232;
fp t1202;
fp t1239;
fp t1189;
fp t1235;
fp t1231;
fp t1230;
fp t1191;
fp t1171;
fp t1169;
fp t1194;
fp t1228;
fp t1196;
fp t1227;
fp t1188;
fp t1226;
fp t1179;
fp t1225;
fp t1190;
fp t1224;
fp t1223;
fp t1222;
fp t1200;
fp t1201;
fp t1221;
fp t1220;
fp t1219;
fp t1218;
fp t1217;
fp t1174;
fp t1173;
fp t1216;
fp t1215;
      t1210 = y*y;
      t1181 = RATIONAL(-1.0,48.0);
      t1212 = t1210*t1210;
      t1182 = RATIONAL(1.0,48.0);
      t1274 = t1212*t1182+t1210*t1181;
      t1185 = RATIONAL(1.0,24.0);
      t1249 = x*t1185;
      t1180 = RATIONAL(-1.0,24.0);
      t1245 = x*t1180;
      t1273 = t1212*t1249+t1210*t1245;
      t1184 = RATIONAL(-1.0,12.0);
      t1244 = x*t1184;
      t1211 = y*t1210;
      t1248 = x*t1211;
      t1183 = RATIONAL(1.0,12.0);
      t1106 = t1183*t1248;
      t1272 = t1106+y*t1244;
      t1271 = t1274*x;
      t1270 = y*t1249+t1211*t1245;
      t1250 = x*RATIONAL(-2.0,3.0);
      t1197 = RATIONAL(1.0,6.0);
      t1252 = x*t1212;
      t1269 = t1210*t1250+t1197*t1252;
      t1268 = y*t1250+t1197*t1248;
      t1155 = t1181*t1212;
      t1267 = (t1182*t1210+t1155)*x;
      t1098 = t1211*t1244;
      t1253 = x*y;
      t1203 = RATIONAL(1.0,3.0);
      t1266 = t1203*t1253+t1098;
      t1265 = t1185*t1248+y*t1245;
      t1264 = t1183*t1253+t1098;
      t1205 = RATIONAL(-1.0,3.0);
      t1246 = x*t1205;
      t1263 = t1210*t1246+t1183*t1252;
      t1262 = t1212*t1245+t1210*t1249;
      t1236 = t1210*t1203;
      t1240 = t1212*t1184;
      t1261 = (t1236+t1240)*x;
      t1247 = x*RATIONAL(2.0,3.0);
      t1198 = RATIONAL(-1.0,6.0);
      t1260 = t1210*t1247+t1198*t1252;
      t1259 = y*t1247+t1198*t1248;
      t1258 = y*t1246+t1106;
      t1176 = RATIONAL(1.0,72.0);
      t1186 = RATIONAL(-1.0,18.0);
      t1209 = x*x;
      t1242 = t1209*t1240+t1186*t1210;
      t1257 = t1242+t1209*t1236+t1176*t1212;
      t1175 = RATIONAL(1.0,18.0);
      t1229 = t1183*t1209;
      t1256 = t1242+t1175*t1212+t1210*t1229;
      t1255 = RATIONAL(-1.0,288.0)*t1212+RATIONAL(1.0,288.0)*t1210+t1274*t1209;
      t1237 = t1203*t1209;
      t1238 = t1210*t1209;
      t1193 = RATIONAL(-4.0,3.0);
      t1204 = RATIONAL(8.0,9.0);
      t1206 = RATIONAL(-2.0,9.0);
      t1254 = t1193*t1238+t1204*t1210+(t1237+t1206)*t1212;
      t1251 = y*t1209;
      t1233 = RATIONAL(1.0,8.0)*t1212;
      t1243 = RATIONAL(5.0,48.0)*t1210+t1209*t1233;
      t1187 = RATIONAL(-1.0,2.0);
      t1232 = t1187*t1209;
      t1202 = RATIONAL(-5.0,3.0);
      t1239 = t1202*t1210+t1212*t1232;
      t1189 = RATIONAL(-5.0,8.0);
      t1235 = t1189*t1210;
      t1231 = t1211*t1209;
      t1230 = RATIONAL(2.0,1.0)*t1209;
      t1191 = RATIONAL(1.0,2.0);
      t1171 = t1191*t1209;
      t1169 = RATIONAL(-2.0,1.0)*t1209;
      t1194 = RATIONAL(1.0,4.0);
      t1228 = RATIONAL(-5.0,24.0)+t1194*t1209;
      t1196 = RATIONAL(-1.0,4.0);
      t1227 = t1196*t1209+RATIONAL(5.0,24.0);
      t1188 = RATIONAL(-1.0,8.0);
      t1226 = RATIONAL(-5.0,48.0)*t1212+t1188*t1238+t1243;
      t1179 = RATIONAL(5.0,12.0);
      t1225 = t1179*t1212+t1210*t1230+t1239;
      t1190 = RATIONAL(4.0,3.0);
      t1224 = t1169+t1203*t1212+RATIONAL(5.0,2.0)*t1238+t1190+t1239;
      t1223 = t1171+t1209*t1235+t1184+t1155+t1243;
      t1222 = t1204*y+t1206*t1211+t1193*t1251+t1203*t1231+t1254;
      t1200 = RATIONAL(1.0,9.0);
      t1201 = RATIONAL(-1.0,9.0);
      t1221 = t1197*t1251+t1201*y+t1200*t1211+t1198*t1231+t1256;
      t1220 = RATIONAL(2.0,9.0)*t1211+t1205*t1231+t1190*t1251+RATIONAL(-8.0,9.0
)*y+t1254;
      t1219 = t1201*t1211+t1198*t1251+t1197*t1231+t1200*y+t1256;
      t1218 = t1205*t1251+t1175*y+t1257+(RATIONAL(-1.0,72.0)+t1229)*t1211;
      t1217 = t1257+t1184*t1231+t1176*t1211+(t1186+t1237)*y;
      t1174 = RATIONAL(-1.0,144.0);
      t1173 = RATIONAL(1.0,144.0);
      t1216 = t1174*y+t1185*t1251+t1180*t1231+t1173*t1211+t1255;
      t1215 = t1185*t1231+t1180*t1251+t1173*y+t1174*t1211+t1255;
      coeffs_dxx->coeff_m2_m2 = t1216+t1265+t1267;
      coeffs_dxx->coeff_m1_m2 = t1219+t1264+t1273;
      coeffs_dxx->coeff_0_m2 = t1227*t1211+t1228*y+t1226;
      coeffs_dxx->coeff_p1_m2 = t1219+t1262+t1272;
      coeffs_dxx->coeff_p2_m2 = t1216+t1270+t1271;
      coeffs_dxx->coeff_m2_m1 = t1218+t1263+t1266;
      coeffs_dxx->coeff_m1_m1 = t1220+t1260+t1268;
      coeffs_dxx->coeff_0_m1 = (t1171+RATIONAL(-5.0,12.0))*t1211+(t1169+
RATIONAL(5.0,3.0))*y+t1225;
      coeffs_dxx->coeff_p1_m1 = t1220+t1259+t1269;
      coeffs_dxx->coeff_p2_m1 = t1218+t1258+t1261;
      coeffs_dxx->coeff_m2_0 = (RATIONAL(5.0,8.0)*t1210+t1187+t1188*t1212)*x+
t1223;
      coeffs_dxx->coeff_m1_0 = (t1194*t1212+1.0+RATIONAL(-5.0,4.0)*t1210)*x+
t1224;
      coeffs_dxx->coeff_0_0 = t1189*t1212+RATIONAL(-5.0,2.0)+RATIONAL(25.0,8.0)
*t1210+(RATIONAL(3.0,1.0)+RATIONAL(3.0,4.0)*t1212+RATIONAL(-15.0,4.0)*t1210)*
t1209;
      coeffs_dxx->coeff_p1_0 = (t1196*t1212+RATIONAL(5.0,4.0)*t1210-1.0)*x+
t1224;
      coeffs_dxx->coeff_p2_0 = (t1191+t1233+t1235)*x+t1223;
      coeffs_dxx->coeff_m2_p1 = t1217+t1258+t1263;
      coeffs_dxx->coeff_m1_p1 = t1222+t1259+t1260;
      coeffs_dxx->coeff_0_p1 = (t1179+t1232)*t1211+(t1202+t1230)*y+t1225;
      coeffs_dxx->coeff_p1_p1 = t1222+t1268+t1269;
      coeffs_dxx->coeff_p2_p1 = t1217+t1261+t1266;
      coeffs_dxx->coeff_m2_p2 = t1215+t1267+t1270;
      coeffs_dxx->coeff_m1_p2 = t1221+t1272+t1273;
      coeffs_dxx->coeff_0_p2 = t1228*t1211+t1227*y+t1226;
      coeffs_dxx->coeff_p1_p2 = t1221+t1262+t1264;
      coeffs_dxx->coeff_p2_p2 = t1215+t1265+t1271;