1 files changed, 228 insertions, 223 deletions
diff --git a/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxx.compute.c b/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxx.compute.c
index 15d553b..3f214f7 100644
--- a/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxx.compute.c
+++ b/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxx.compute.c
@@ -1,224 +1,229 @@
-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 t1370;
+fp t1373;
+fp t1440;
+fp t1340;
+fp t1439;
+fp t1338;
+fp t1438;
+fp t1374;
+fp t1372;
+fp t1412;
+fp t1410;
+fp t1437;
+fp t1364;
+fp t1415;
+fp t1375;
+fp t1341;
+fp t1397;
 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;
+fp t1436;
+fp t1435;
+fp t1409;
+fp t1414;
+fp t1343;
+fp t1265;
+fp t1434;
+fp t1365;
+fp t1394;
+fp t1433;
+fp t1339;
+fp t1393;
+fp t1342;
+fp t1432;
+fp t1413;
+fp t1417;
+fp t1431;
+fp t1311;
+fp t1430;
+fp t1429;
+fp t1428;
+fp t1371;
+fp t1368;
+fp t1427;
+fp t1400;
+fp t1426;
+fp t1408;
+fp t1369;
+fp t1425;
+fp t1424;
+fp t1411;
+fp t1423;
+fp t1422;
+fp t1392;
+fp t1346;
+fp t1404;
+fp t1345;
+fp t1421;
+fp t1403;
+fp t1344;
+fp t1420;
+fp t1419;
+fp t1405;
+fp t1362;
+fp t1363;
+fp t1366;
+fp t1418;
+fp t1416;
+fp t1350;
+fp t1399;
+fp t1361;
+fp t1407;
+fp t1354;
+fp t1322;
+fp t1324;
+fp t1353;
+fp t1406;
+fp t1402;
+fp t1401;
+fp t1356;
+fp t1391;
+fp t1357;
+fp t1390;
+fp t1351;
+fp t1389;
+fp t1347;
+fp t1388;
+fp t1352;
+fp t1387;
+fp t1355;
+fp t1386;
+fp t1385;
+fp t1384;
+fp t1383;
+fp t1382;
+fp t1367;
+fp t1381;
+fp t1380;
+fp t1379;
+fp t1378;
+      t1370 = RATIONAL(-1.0,6.0);
+      t1373 = x*x;
+      t1440 = t1370*t1373+RATIONAL(1.0,9.0);
+      t1340 = RATIONAL(-1.0,24.0);
+      t1439 = t1373*t1340+RATIONAL(1.0,144.0);
+      t1338 = RATIONAL(1.0,24.0);
+      t1438 = t1373*t1338+RATIONAL(-1.0,144.0);
+      t1374 = y*y;
+      t1372 = t1374*t1374;
+      t1412 = x*t1372;
+      t1410 = x*t1374;
+      t1437 = t1338*t1410+t1340*t1412;
+      t1364 = RATIONAL(-1.0,3.0);
+      t1415 = y*t1364;
+      t1375 = t1374*y;
+      t1341 = RATIONAL(1.0,12.0);
+      t1397 = t1375*t1341;
+      t1261 = x*t1397;
+      t1436 = x*t1415+t1261;
+      t1435 = t1340*t1410+t1338*t1412;
+      t1409 = x*t1341;
+      t1414 = x*t1375;
+      t1343 = RATIONAL(-1.0,12.0);
+      t1265 = t1343*t1414;
+      t1434 = t1265+y*t1409;
+      t1365 = RATIONAL(1.0,3.0);
+      t1394 = t1374*t1365;
+      t1433 = x*t1394+t1343*t1412;
+      t1339 = RATIONAL(1.0,48.0);
+      t1393 = t1372*t1339;
+      t1342 = RATIONAL(-1.0,48.0);
+      t1432 = t1342*t1410+x*t1393;
+      t1413 = x*t1338;
+      t1417 = x*y;
+      t1431 = t1375*t1413+t1340*t1417;
+      t1311 = t1342*t1372;
+      t1430 = t1339*t1410+x*t1311;
+      t1429 = t1340*t1414+y*t1413;
+      t1428 = t1372*t1409+t1364*t1410;
+      t1371 = RATIONAL(-2.0,3.0);
+      t1368 = RATIONAL(1.0,6.0);
+      t1427 = t1368*t1412+t1371*t1410;
+      t1400 = t1375*t1368;
+      t1426 = x*t1400+t1371*t1417;
+      t1408 = t1370*x;
+      t1369 = RATIONAL(2.0,3.0);
+      t1425 = t1369*t1417+t1375*t1408;
+      t1424 = t1372*t1408+t1369*t1410;
+      t1411 = y*t1365;
+      t1423 = t1265+x*t1411;
+      t1422 = t1343*t1417+t1261;
+      t1392 = t1373*t1343;
+      t1346 = RATIONAL(-1.0,18.0);
+      t1404 = t1346*t1374+t1372*t1392;
+      t1345 = RATIONAL(1.0,72.0);
+      t1421 = t1404+t1345*t1372+t1373*t1394;
+      t1403 = t1373*t1374;
+      t1344 = RATIONAL(1.0,18.0);
+      t1420 = t1344*t1372+t1404+t1341*t1403;
+      t1419 = t1342*t1403+t1373*t1393+RATIONAL(-1.0,288.0)*t1372+RATIONAL(1.0,
+288.0)*t1374;
+      t1405 = t1372*t1373;
+      t1362 = RATIONAL(-2.0,9.0);
+      t1363 = RATIONAL(-4.0,3.0);
+      t1366 = RATIONAL(8.0,9.0);
+      t1418 = t1362*t1372+t1366*t1374+t1365*t1405+t1363*t1403;
+      t1416 = y*t1373;
+      t1350 = RATIONAL(-1.0,2.0);
+      t1399 = t1350*t1373;
+      t1361 = RATIONAL(-5.0,3.0);
+      t1407 = t1361*t1374+t1372*t1399;
+      t1354 = RATIONAL(1.0,2.0);
+      t1322 = t1354*t1373;
+      t1324 = RATIONAL(-2.0,1.0)*t1373;
+      t1353 = RATIONAL(1.0,8.0);
+      t1406 = t1353*t1405+RATIONAL(5.0,48.0)*t1374;
+      t1402 = t1373*t1375;
+      t1401 = RATIONAL(2.0,1.0)*t1373;
+      t1356 = RATIONAL(-1.0,4.0);
+      t1391 = t1356*t1373+RATIONAL(5.0,24.0);
+      t1357 = RATIONAL(1.0,4.0);
+      t1390 = t1357*t1373+RATIONAL(-5.0,24.0);
+      t1351 = RATIONAL(-1.0,8.0);
+      t1389 = RATIONAL(-5.0,48.0)*t1372+t1351*t1403+t1406;
+      t1347 = RATIONAL(5.0,12.0);
+      t1388 = t1347*t1372+t1374*t1401+t1407;
+      t1352 = RATIONAL(4.0,3.0);
+      t1387 = t1352+RATIONAL(5.0,2.0)*t1403+t1365*t1372+t1324+t1407;
+      t1355 = RATIONAL(-5.0,8.0);
+      t1386 = t1355*t1403+t1322+t1311+t1343+t1406;
+      t1385 = t1419+t1439*t1375+t1438*y;
+      t1384 = t1419+t1438*t1375+t1439*y;
+      t1383 = t1421+t1373*t1411+t1346*y+(t1392+t1345)*t1375;
+      t1382 = t1344*y+RATIONAL(-1.0,72.0)*t1375+t1421+(t1415+t1397)*t1373;
+      t1367 = RATIONAL(-1.0,9.0);
+      t1381 = t1367*y+t1368*t1416+t1420+t1440*t1375;
+      t1380 = RATIONAL(-8.0,9.0)*y+RATIONAL(2.0,9.0)*t1375+t1364*t1402+t1352*
+t1416+t1418;
+      t1379 = t1420+t1373*t1400+t1367*t1375+t1440*y;
+      t1378 = t1363*t1416+t1362*t1375+t1365*t1402+y*t1366+t1418;
+      coeffs_dxx->coeff_m2_m2 = t1385+t1430+t1431;
+      coeffs_dxx->coeff_m1_m2 = t1379+t1434+t1435;
+      coeffs_dxx->coeff_0_m2 = t1391*t1375+t1390*y+t1389;
+      coeffs_dxx->coeff_p1_m2 = t1379+t1422+t1437;
+      coeffs_dxx->coeff_p2_m2 = t1385+t1429+t1432;
+      coeffs_dxx->coeff_m2_m1 = t1382+t1423+t1428;
+      coeffs_dxx->coeff_m1_m1 = t1380+t1424+t1426;
+      coeffs_dxx->coeff_0_m1 = (RATIONAL(-5.0,12.0)+t1322)*t1375+(RATIONAL(5.0,
+3.0)+t1324)*y+t1388;
+      coeffs_dxx->coeff_p1_m1 = t1380+t1425+t1427;
+      coeffs_dxx->coeff_p2_m1 = t1382+t1433+t1436;
+      coeffs_dxx->coeff_m2_0 = (t1351*t1372+RATIONAL(5.0,8.0)*t1374+t1350)*x+
+t1386;
+      coeffs_dxx->coeff_m1_0 = (t1357*t1372+1.0+RATIONAL(-5.0,4.0)*t1374)*x+
+t1387;
+      coeffs_dxx->coeff_0_0 = RATIONAL(25.0,8.0)*t1374+RATIONAL(-5.0,2.0)+t1355
+*t1372+(RATIONAL(3.0,4.0)*t1372+RATIONAL(-15.0,4.0)*t1374+RATIONAL(3.0,1.0))*
+t1373;
+      coeffs_dxx->coeff_p1_0 = (RATIONAL(5.0,4.0)*t1374+t1356*t1372-1.0)*x+
+t1387;
+      coeffs_dxx->coeff_p2_0 = (t1353*t1372+t1355*t1374+t1354)*x+t1386;
+      coeffs_dxx->coeff_m2_p1 = t1383+t1428+t1436;
+      coeffs_dxx->coeff_m1_p1 = t1378+t1424+t1425;
+      coeffs_dxx->coeff_0_p1 = (t1347+t1399)*t1375+(t1401+t1361)*y+t1388;
+      coeffs_dxx->coeff_p1_p1 = t1378+t1426+t1427;
+      coeffs_dxx->coeff_p2_p1 = t1383+t1423+t1433;
+      coeffs_dxx->coeff_m2_p2 = t1384+t1429+t1430;
+      coeffs_dxx->coeff_m1_p2 = t1381+t1422+t1435;
+      coeffs_dxx->coeff_0_p2 = t1390*t1375+t1391*y+t1389;
+      coeffs_dxx->coeff_p1_p2 = t1381+t1434+t1437;
+      coeffs_dxx->coeff_p2_p2 = t1384+t1431+t1432;