From 99f3fe1eceecccb7d9487e8ecdb171c7d6865428 Mon Sep 17 00:00:00 2001 From: jthorn Date: Wed, 28 Aug 2002 11:42:25 +0000 Subject: add Maple code and C skeletons for Hermite polynomial interpolation git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@101 df1f8a13-aa1d-4dd4-9681-27ded5b42416 --- .../3d.coeffs/3d.cube.order2/coeff-dx.compute.c | 1142 ++++++++++++++++++++ 1 file changed, 1142 insertions(+) create mode 100644 src/GeneralizedPolynomial-Uniform/Hermite/3d.coeffs/3d.cube.order2/coeff-dx.compute.c (limited to 'src/GeneralizedPolynomial-Uniform/Hermite/3d.coeffs/3d.cube.order2/coeff-dx.compute.c') diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/3d.coeffs/3d.cube.order2/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/3d.coeffs/3d.cube.order2/coeff-dx.compute.c new file mode 100644 index 0000000..68acced --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Hermite/3d.coeffs/3d.cube.order2/coeff-dx.compute.c @@ -0,0 +1,1142 @@ +fp t1274, + t1441, + t1260, + t1271, + t1270, + t1422, + t1267, + t1542, + t1442, + t1263, + t1273, + t1419, + t1265, + t1541, + t1449, + t1375, + t1446, + t1372, + t1269, + t1268, + t1225, + t1272, + t1379, + t1216, + t1322, + t886, + t1540, + t1417, + t1368, + t1240, + t1254, + t1408, + t1326, + t998, + t1539, + t1237, + t1454, + t1418, + t1239, + t1380, + t1376, + t929, + t1538, + t1445, + t1253, + t1369, + t979, + t1313, + t1251, + t1537, + t1230, + t1536, + t1242, + t1378, + t1318, + t974, + t1535, + t1257, + t920, + t1448, + t1252, + t1363, + t1534, + t1331, + t1226, + t1533, + t1317, + t950, + t1425, + t1377, + t1532, + t871, + t1235, + t1429, + t1362, + t1531, + t1420, + t1261, + t1248, + t1233, + t1530, + t1247, + t1338, + t1211, + t1371, + t846, + t1529, + t1205, + t1213, + t876, + t1528, + t1243, + t1396, + t1131, + t1385, + t1323, + t1527, + t1342, + t821, + t1526, + t1200, + t1525, + t1524, + t1220, + t867, + t1212, + t1314, + t1523, + t930, + t1522, + t1159, + t1521, + t1256, + t1400, + t1370, + t1039, + t1520, + t1222, + t1327, + t1241, + t1519, + t1518, + t1426, + t1321, + t830, + t1231, + t1517, + t1394, + t1218, + t1330, + t1516, + t1196, + t1207, + t1515, + t1384, + t1514, + t1074, + t1513, + t1246, + t1392, + t1158, + t1512, + t883, + t1511, + t1510, + t1133, + t1214, + t1509, + t1508, + t1507, + t965, + t1506, + t1221, + t1405, + t1505, + t1223, + t1388, + t1142, + t1504, + t1503, + t1502, + t928, + t1501, + t1238, + t1185, + t1245, + t1409, + t1500, + t1407, + t1499, + t1249, + t1498, + t1232, + t1415, + t1118, + t1497, + t1037, + t1496, + t1183, + t1495, + t1102, + t1236, + t1406, + t1174, + t1494, + t1250, + t1414, + t1148, + t1493, + t1492, + t1184, + t1491, + t1398, + t1181, + t1447, + t1122, + t1490, + t1109, + t1134, + t1489, + t1264, + t1451, + t1152, + t1488, + t1191, + t1412, + t1487, + t1413, + t1179, + t1486, + t1177, + t1444, + t1107, + t1485, + t1188, + t1484, + t1483, + t951, + t1482, + t1197, + t1199, + t1481, + t1443, + t1095, + t1480, + t1208, + t1130, + t1479, + t1229, + t1399, + t1478, + t1170, + t1138, + t1477, + t1316, + t1393, + t1476, + t1042, + t1143, + t1475, + t1401, + t1104, + t1474, + t1473, + t1228, + t1204, + t1472, + t1471, + t1404, + t1227, + t1156, + t1470, + t1469, + t943, + t1468, + t1019, + t1450, + t1367, + t1467, + t925, + t1386, + t1105, + t1160, + t1466, + t1387, + t1465, + t955, + t1464, + t1195, + t1215, + t1202, + t1463, + t1151, + t1374, + t1462, + t987, + t1162, + t1461, + t1460, + t1459, + t1259, + t1423, + t1123, + t1458, + t1411, + t1116, + t1244, + t1457, + t1176, + t1194, + t1456, + t1258, + t1453, + t1266, + t1455, + t1452, + t1190, + t1187, + t1186, + t1224, + t1440, + t859, + t1439, + t1361, + t999, + t1438, + t1437, + t1436, + t1383, + t1325, + t898, + t1435, + t1434, + t1332, + t1035, + t1433, + t1432, + t1059, + t1364, + t1430, + t1192, + t1424, + t1048, + t1427, + t1262, + t1421, + t1416, + t1410, + t1403, + t1402, + t1091, + t1397, + t1395, + t1391, + t1390, + t1389, + t1382, + t1255, + t1381, + t1373, + t1168, + t1085, + t1365, + t1178, + t1153, + t1124, + t923, + t1360, + t986, + t1340, + t1359, + t1358, + t1357, + t1356, + t1355, + t1354, + t1353, + t1352, + t1351, + t877, + t1344, + t1350, + t1329, + t1349, + t1348, + t1347, + t1346, + t1345, + t1013, + t1343, + t1341, + t1339, + t1093, + t1129, + t1157, + t1096, + t1002, + t1337, + t1092, + t1336, + t1144, + t1335, + t1334, + t1333, + t1328, + t1324, + t1166, + t1320, + t1319, + t1180, + t1315, + t1217, + t1312, + t1311, + t1310, + t1308, + t1307, + t1306, + t942, + t1305, + t1304, + t961, + t1303, + t1301, + t1300, + t1036, + t1030, + t1299, + t1298, + t1297, + t995, + t1296, + t833, + t1295, + t1206, + t852, + t1294, + t1293, + t1210, + t875, + t1292, + t997, + t1291, + t887, + t1290, + t1289, + t1058, + t1288, + t1046, + t1287, + t1003, + t1286, + t1077, + t1285, + t1060, + t1284, + t1283, + t1282, + t1281, + t1234, + t1280, + t1279, + t1278, + t897, + t1277, + t1219, + t1209, + t1203, + t1201, + t1198, + t1167, + t1164, + t1163, + t1155, + t1149, + t1140, + t1137, + t1136, + t1128, + t1121, + t1119, + t1103, + t1098, + t1097, + t1090, + t1075, + t1072, + t1070, + t1066, + t1055, + t1052, + t1050, + t1047, + t1045, + t1040, + t1038, + t1034, + t1028, + t1025, + t1016, + t1010, + t1009, + t1008, + t1006, + t985, + t984, + t978, + t975, + t971, + t968, + t967, + t966, + t963, + t958, + t957, + t956, + t937, + t933, + t932, + t918, + t917, + t912, + t909, + t907, + t905, + t900, + t896, + t894, + t873, + t870, + t868, + t843, + t841, + t838, + t832, + t831, + t827, + t823; + t1274 = z*z; + t1441 = x*t1274; + t1260 = y*t1441; + t1271 = y*y; + t1270 = t1271*y; + t1422 = t1274*t1270; + t1267 = x*t1422; + t1542 = t1267+t1260; + t1442 = z*t1271; + t1263 = x*t1442; + t1273 = t1274*z; + t1419 = t1273*t1271; + t1265 = x*t1419; + t1541 = t1263+t1265; + t1449 = t1273+z; + t1375 = x*t1449; + t1446 = t1270+y; + t1372 = x*t1446; + t1269 = t1271*t1274; + t1268 = x*t1269; + t1225 = RATIONAL(-4.0,1.0); + t1272 = x*x; + t1379 = t1273*t1272; + t1216 = RATIONAL(27.0,8.0); + t1322 = t1216*t1379; + t886 = t1270*t1322; + t1540 = t1225*t1268+t886; + t1417 = t1271*t1272; + t1368 = z*t1417; + t1240 = RATIONAL(-3.0,2.0); + t1254 = RATIONAL(9.0,8.0); + t1408 = t1254*t1270; + t1326 = t1272*t1408; + t998 = z*t1326; + t1539 = t1240*t1368+t998; + t1237 = RATIONAL(1.0,2.0); + t1454 = x*z; + t1418 = t1237*t1454; + t1239 = RATIONAL(-3.0,8.0); + t1380 = t1239*t1272; + t1376 = z*t1380; + t929 = y*t1376; + t1538 = y*t1418+t929; + t1445 = z*t1272; + t1253 = RATIONAL(3.0,8.0); + t1369 = t1253*t1445; + t979 = t1270*t1369; + t1313 = t1272*t1422; + t1251 = RATIONAL(3.0,2.0); + t1537 = t979+t1251*t1313; + t1230 = RATIONAL(-5.0,1.0); + t1536 = t1271*t1322+t1230*t1267; + t1242 = RATIONAL(-5.0,2.0); + t1378 = t1242*t1441; + t1318 = t1254*t1417; + t974 = t1273*t1318; + t1535 = t1270*t1378+t974; + t1257 = RATIONAL(-9.0,8.0); + t920 = t1257*t1313; + t1448 = x*t1273; + t1252 = RATIONAL(5.0,4.0); + t1363 = t1252*t1448; + t1534 = t920+t1271*t1363; + t1331 = t1274*t1417; + t1226 = RATIONAL(-9.0,4.0); + t1533 = t1270*t1363+t1226*t1331; + t1317 = t1257*t1379; + t950 = t1271*t1317; + t1425 = t1274*t1252; + t1377 = x*t1425; + t1532 = t950+t1270*t1377; + t871 = y*t1322; + t1235 = RATIONAL(-9.0,2.0); + t1429 = t1274*t1272; + t1362 = y*t1429; + t1531 = t871+t1235*t1362; + t1420 = t1273*t1270; + t1261 = x*t1420; + t1248 = RATIONAL(-3.0,4.0); + t1233 = RATIONAL(1.0,4.0); + t1530 = t1248*t1331+t1233*t1261; + t1247 = RATIONAL(9.0,2.0); + t1338 = t1247*t1429; + t1211 = RATIONAL(-27.0,8.0); + t1371 = y*t1379; + t846 = t1211*t1371; + t1529 = y*t1338+t846; + t1205 = RATIONAL(-45.0,4.0); + t1213 = RATIONAL(15.0,4.0); + t876 = t1213*t1261; + t1528 = t876+t1205*t1331; + t1243 = RATIONAL(-1.0,2.0); + t1396 = t1274*t1243; + t1131 = t1270*t1396; + t1385 = t1271*t1253; + t1323 = t1272*t1385; + t1527 = t1273*t1323+x*t1131; + t1342 = t1270*t1379; + t821 = RATIONAL(-81.0,8.0)*t1342; + t1526 = t821+RATIONAL(25.0,1.0)*t1268; + t1200 = RATIONAL(45.0,4.0); + t1525 = t1200*t1261+RATIONAL(-45.0,2.0)*t1331; + t1524 = t1230*t1265+t1216*t1313; + t1220 = RATIONAL(-15.0,4.0); + t867 = t1220*t1261; + t1212 = RATIONAL(45.0,8.0); + t1314 = t1212*t1429; + t1523 = t1271*t1314+t867; + t930 = t1254*t1313; + t1522 = t930+t1242*t1265; + t1159 = t1257*t1420; + t1521 = t1247*t1261+t1159; + t1256 = RATIONAL(9.0,4.0); + t1400 = t1257*t1272; + t1370 = z*t1400; + t1039 = t1270*t1370; + t1520 = t1256*t1313+t1039; + t1222 = RATIONAL(-15.0,8.0); + t1327 = t1271*t1379; + t1241 = RATIONAL(3.0,4.0); + t1519 = t1222*t1327+t1241*t1379; + t1518 = t1235*t1313+t1039; + t1426 = t1270*t1272; + t1321 = t1211*t1426; + t830 = t1273*t1321; + t1231 = RATIONAL(8.0,1.0); + t1517 = t830+t1231*t1268; + t1394 = t1248*t1272; + t1218 = RATIONAL(15.0,8.0); + t1330 = t1218*t1417; + t1516 = (t1394+t1330)*z; + t1196 = RATIONAL(-135.0,8.0); + t1207 = RATIONAL(27.0,4.0); + t1515 = t1196*t1313+t1207*t1426; + t1384 = t1241*t1272; + t1514 = t1270*t1384+t1222*t1313; + t1074 = t1241*t1261; + t1513 = t1220*t1331+t1074; + t1246 = RATIONAL(-1.0,8.0); + t1392 = t1246*t1270; + t1158 = t1273*t1392; + t1512 = t1158+t1237*t1261; + t883 = t1216*z*t1426; + t1511 = t883+t1235*t1368; + t1510 = y*t1394+t1218*t1362; + t1133 = t1230*t1441; + t1214 = RATIONAL(-15.0,2.0); + t1509 = t1133+t1214*t1267; + t1508 = t1212*t1327+t1226*t1379; + t1507 = t1214*t1331+t1256*t1261; + t965 = t1226*t1313; + t1506 = t965+t1541; + t1221 = RATIONAL(-45.0,8.0); + t1405 = t1256*t1272; + t1505 = t1221*t1313+t1270*t1405; + t1223 = RATIONAL(5.0,2.0); + t1388 = t1223*t1271; + t1142 = x*t1388; + t1504 = t1142+t1213*t1265; + t1503 = t1207*t1379+t1196*t1327; + t1502 = -t1267+t974; + t928 = t1240*t1331; + t1501 = t1074+t928; + t1238 = RATIONAL(-1.0,4.0); + t1185 = t1238*t1273; + t1245 = RATIONAL(5.0,8.0); + t1409 = t1245*t1271; + t1500 = t1185+t1273*t1409; + t1407 = t1245*t1274; + t1499 = (t1238+t1407)*t1270; + t1249 = RATIONAL(-5.0,8.0); + t1498 = t1249*t1419+t1233*t1273; + t1232 = RATIONAL(-3.0,1.0); + t1415 = t1239*t1273; + t1118 = t1270*t1415; + t1497 = t1232*t1261+t1118; + t1037 = y*t1317; + t1496 = t1235*t1327+t1037; + t1183 = t1248*t1270; + t1495 = t1183+t1218*t1422; + t1102 = t1237*t1422; + t1236 = RATIONAL(1.0,8.0); + t1406 = t1270*t1236; + t1174 = z*t1406; + t1494 = t1102+t1174; + t1250 = RATIONAL(3.0,1.0); + t1414 = t1273*t1253; + t1148 = t1270*t1414; + t1493 = t1250*t1261+t1148; + t1492 = t1037+t1256*t1327; + t1184 = t1233*t1270; + t1491 = t1249*t1422+t1184; + t1398 = t1237*t1271; + t1181 = z*t1398; + t1447 = z*t1270; + t1122 = t1239*t1447; + t1490 = t1181+t1122; + t1109 = t1273*t1408; + t1134 = t1240*t1419; + t1489 = t1109+t1134; + t1264 = y*t1448; + t1451 = y*t1246; + t1152 = t1273*t1451; + t1488 = -t1264+t1152; + t1191 = t1241*t1270; + t1412 = t1222*t1274; + t1487 = t1191+t1270*t1412; + t1413 = t1273*t1236; + t1179 = y*t1413; + t1486 = t1264+t1179; + t1177 = y*t1415; + t1444 = y*t1274; + t1107 = t1237*t1444; + t1485 = t1177+t1107; + t1188 = t1248*t1273; + t1484 = t1218*t1419+t1188; + t1483 = t1223*t1265+t965+z*t1142; + t951 = t1247*t1313; + t1482 = t1220*t1265+t1252*t1263+t951; + t1197 = RATIONAL(-27.0,2.0); + t1199 = RATIONAL(15.0,1.0); + t1481 = t1197*t1327+y*t1133+t1199*t1267; + t1443 = x*t1271; + t1095 = t1230*t1443; + t1480 = t1197*t1313+t1199*t1265+z*t1095; + t1208 = RATIONAL(15.0,2.0); + t1130 = t1274*t1388; + t1479 = t1130+t1208*t1331+t1235*t1261; + t1229 = RATIONAL(2.0,1.0); + t1399 = t1274*t1236; + t1478 = y*t1399+t920+t1229*t1265; + t1170 = t1251*t1422; + t1138 = y*t1396; + t1477 = t1240*t1327+(t1170+t1138)*x; + t1316 = t1274*t1384; + t1393 = t1233*t1271; + t1476 = t1271*t1316+t1274*t1393+t1243*t1261; + t1042 = t1226*t1327; + t1143 = t1223*t1441; + t1475 = t1042+t1223*t1267+y*t1143; + t1401 = t1270*t1240; + t1104 = t1274*t1401; + t1474 = t1248*t1379+x*t1104+t1218*t1327; + t1473 = x*t1134+t1218*t1313+t1270*t1394; + t1228 = RATIONAL(5.0,1.0); + t1204 = RATIONAL(-27.0,4.0); + t1472 = t1204*t1327+t1542*t1228; + t1471 = t1248*t1327+(t1107+t1102)*x; + t1404 = t1238*t1271; + t1227 = RATIONAL(-2.0,1.0); + t1156 = t1227*t1441; + t1470 = z*t1404+y*t1156+t1227*t1267; + t1469 = t1204*t1313+t1541*t1228; + t943 = t1254*t1371; + t1468 = t1042+t943+t1542; + t1019 = t1247*t1327; + t1450 = x*t1270; + t1367 = t1220*t1450; + t1467 = t1019+y*t1377+t1274*t1367; + t925 = t1270*t1376; + t1386 = t1271*t1251; + t1105 = t1273*t1386; + t1160 = t1243*t1442; + t1466 = t925+t1240*t1313+(t1105+t1160)*x; + t1387 = t1242*t1271; + t1465 = t1208*t1265+t1387*t1454+t1518; + t955 = t1250*t1265; + t1464 = -t1263+t955+t1518; + t1195 = RATIONAL(135.0,8.0); + t1215 = RATIONAL(-15.0,1.0); + t1202 = RATIONAL(-10.0,1.0); + t1463 = t1202*t1441+t1215*t1267+t1195*t1327+t1204*t1379; + t1151 = z*t1392; + t1374 = t1229*t1450; + t1462 = t1274*t1374+t1151+t950+t1236*t1442; + t987 = t1247*t1331; + t1162 = t1271*t1396; + t1461 = t987+t1162+t1497; + t1460 = t1221*t1331+t1271*t1407+t1493; + t1459 = y*t1378+t1208*t1267+t1496; + t1259 = RATIONAL(-5.0,4.0); + t1423 = t1274*t1259; + t1123 = t1271*t1423; + t1458 = t1123+t1200*t1331+t1497; + t1411 = t1229*t1271; + t1116 = y*t1414; + t1244 = RATIONAL(-6.0,1.0); + t1457 = t1116+t951+t1138+t1244*t1265+t1411*t1454; + t1176 = t1227*t1443; + t1194 = t1238*y; + t1456 = z*t1176+t1274*t1194+t1227*t1265+t1520; + t1258 = RATIONAL(-9.0,1.0); + t1453 = x*y; + t1266 = z*t1453; + t1455 = -t1266+t1269+t1258*t1331+t1493; + t1452 = t1243*y; + t1190 = t1233*z; + t1187 = t1233*y; + t1186 = t1238*z; + t1224 = RATIONAL(4.0,1.0); + t1440 = t1224*x; + t859 = RATIONAL(81.0,8.0)*t1342; + t1439 = t859+t1205*t1261; + t1361 = x*t1183; + t999 = t1273*t1326; + t1438 = t1273*t1361+t999; + t1437 = t1257*t1368+t998; + t1436 = t943+t1259*t1264; + t1383 = x*t1187; + t1325 = t1239*t1379; + t898 = y*t1325; + t1435 = t1273*t1383+t898; + t1434 = t1259*t1261+t999; + t1332 = t1253*t1379; + t1035 = t1270*t1332; + t1433 = t1239*t1313+t1035; + t1432 = t1270+t1273; + t1059 = y*t1332; + t1364 = x*t1185; + t1430 = t1059+y*t1364; + t1192 = t1241*t1273; + t1424 = x*t1186; + t1048 = y*t1369; + t1427 = y*t1424+t1048; + t1262 = x*t1447; + t1421 = t1232*t1270; + t1416 = t1246*t1271; + t1410 = t1232*t1272; + t1403 = t1258*t1272; + t1402 = t1274*t1233; + t1091 = t1273*t1406; + t1397 = t1091+t1274*t1392; + t1395 = x*t1452; + t1391 = t1268-t1265; + t1390 = t1229*t1274; + t1389 = t1259*t1454; + t1382 = t1226*t1272; + t1255 = RATIONAL(9.0,1.0); + t1381 = t1255*t1272; + t1373 = t1250*t1450; + t1168 = t1271*t1440; + t1085 = t1202*t1443; + t1365 = t1213*t1450; + t1178 = y*t1405; + t1153 = z*t1384; + t1124 = z*t1405; + t923 = t1270*t1325; + t1360 = t925+t1241*t1313+t923; + t986 = t1254*y*t1445; + t1340 = t1271*t1382; + t1359 = t986+z*t1340+t998; + t1358 = t1221*t1368+t1124+t883; + t1357 = t1222*t1368+t1153+t998; + t1356 = t846+(t1382+t1314)*y; + t1355 = t886+t1505; + t1354 = t1221*t1327+t1256*t1379+t886; + t1353 = t1152+t1237*t1264+t898; + t1352 = t871+t1178+t1221*t1362; + t1351 = t830+t1270*t1382+t1212*t1313; + t877 = z*t1321; + t1344 = t1247*t1417; + t1350 = t877+t986+z*t1344; + t1329 = t1272*t1412; + t1349 = t943+(t1384+t1329)*y; + t1348 = t877+z*t1382+t1212*t1368; + t1347 = t830+t1508; + t1346 = t1240*t1261+t999+t1148; + t1345 = t1249*t1444+t1116+t1187; + t1013 = t1270*t1317; + t1343 = t1118+t1251*t1261+t1013; + t1341 = t1271*t1394; + t1339 = t1272*t1386; + t1093 = t1273*t1398; + t1129 = t1243*t1419; + t1157 = t1273*t1393; + t1096 = t1236*y*z; + t1002 = y*t1370; + t1337 = t1096+t1002+t1266; + t1092 = t1274*t1398; + t1336 = t1251*t1429; + t1144 = t1253*t1447; + t1335 = t1249*t1442+t1190+t1144; + t1334 = t1274*t1394; + t1333 = t1274*t1382; + t1328 = t1273*t1395+t1059+t1179; + t1324 = t1109+t1222*t1419+t1192; + t1166 = z*t1451; + t1320 = t1151+t1166+z*t1393; + t1319 = t1013+t1492; + t1180 = t1241*t1419; + t1315 = t1091+t1238*t1422+t1174; + t1217 = RATIONAL(25.0,4.0); + t1312 = t1217*t1271+t1242; + t1311 = t1200*t1272+t1259; + t1310 = t1217*t1274+t1242; + t1308 = t1274*t1372; + t1307 = t1405+t1238; + t1306 = t1252+t1213*t1272; + t942 = x*t1123; + t1305 = t1274*t1330+t942+t1438; + t1304 = t1240*t1272+t1243; + t961 = x*t1092; + t1303 = t961+t1253*t1313+x*t1129+t923; + t1301 = t1237+t1235*t1272; + t1300 = t1271*t1375; + t1036 = t1235*t1331; + t1030 = t1228*t1268; + t1299 = t876+t830+t1036+t1030; + t1298 = t930+t1013+t1391; + t1297 = t1233+t1384-x; + t995 = t1274*t1176; + t1296 = z*t1383+t929+t1250*t1331+t995+t1438; + t833 = RATIONAL(-20.0,1.0)*t1268; + t1295 = RATIONAL(45.0,2.0)*t1331+t1274*t1387+t859+t833+t1258*t1261+t1159; + t1206 = RATIONAL(-25.0,4.0); + t852 = t1206*t1268; + t1294 = t852+t1143+t1274*t1365+t1354; + t1293 = t1195*t1313+t1204*t1426+t1215*t1265+t1085+t1526; + t1210 = RATIONAL(25.0,2.0); + t875 = t1210*t1268; + t1292 = t1214*t1265+t1095+t875+t1351; + t997 = t1274*t1095; + t1291 = t997+t1271*t1425+t1213*t1331+t1346; + t887 = t1274*t1085; + t1290 = t1255*t1331+t887+t1002+t867+t1252*t1266+t886; + t1289 = t1274*t1404-t1261+t999+t1158+t1256*t1331+t995; + t1058 = t1251*t1331; + t1288 = t1058+t1092+t995+t1346; + t1046 = t1274*t1373; + t1287 = -t1260+t1046+t1496+t1540; + t1003 = t1232*t1265; + t1286 = t1176+t1030+t1003+t1351; + t1077 = t1232*t1267; + t1285 = t1156+t1077+t1030+t1347; + t1060 = x*t1130; + t1284 = t1271*t1389+t1259*t1265+t1060+t1013+t1520; + t1283 = t1060+t1319+t1423*t1372; + t1282 = t1035+t979+t1248*t1313-t1268+(t1093+t1181)*x; + t1281 = t1271*t1329+t1249*t1269+t1060+t1343; + t1234 = RATIONAL(6.0,1.0); + t1280 = y*t1407+t887+t1234*t1265+t1177+t1168+t1194+t1355; + t1279 = t887+t1234*t1267+z*t1409+t1186+t1274*t1440+t1122+t1354; + t1278 = t1244*t1267+t943+t1019+t1390*t1453+t1166+t1160+t1144+t1517; + t897 = t1274*t1168; + t1277 = -t1269+t1096+t1232*t1331+t1048+t897+z*t1395+t1343; + t1219 = RATIONAL(75.0,4.0); + t1209 = RATIONAL(-12.0,1.0); + t1203 = RATIONAL(-75.0,4.0); + t1201 = RATIONAL(27.0,2.0); + t1198 = RATIONAL(12.0,1.0); + t1167 = t1273*t1385; + t1164 = y*t1402; + t1163 = t1246*t1444; + t1155 = t1270*t1402; + t1149 = t1248*t1422; + t1140 = t1271*t1415; + t1137 = t1239*t1422; + t1136 = t1271*t1413; + t1128 = t1270*t1399; + t1121 = t1271*t1188; + t1119 = t1241*t1422; + t1103 = z*t1416; + t1098 = t1273*t1416; + t1097 = t1253*t1422; + t1090 = t1271*t1185; + t1075 = t1241*t1262; + t1072 = t1270*t1424; + t1070 = t1251*t1327; + t1066 = t1421*t1454; + t1055 = z*t1339; + t1052 = t1232*t1264; + t1050 = y*t1334; + t1047 = z*t1341; + t1045 = t1401*t1454; + t1040 = t1271*t1124; + t1038 = t1251*t1262; + t1034 = t1270*t1364; + t1028 = t1274*t1178; + t1025 = y*t1363; + t1016 = t1188*t1453; + t1010 = t1271*t1325; + t1009 = t1241*t1264; + t1008 = z*t1373; + t1006 = z*t1361; + t985 = t1243*t1262; + t984 = t1239*t1368; + t978 = t1252*t1262; + t975 = t1255*t1261; + t971 = z*t1318; + t968 = t1226*t1261; + t967 = t1251*t1264; + t966 = z*t1323; + t963 = t1270*t1389; + t958 = t1233*t1262; + t957 = y*t1389; + t956 = t1240*t1264; + t937 = t1239*t1362; + t933 = y*t1336; + t932 = t1254*t1362; + t918 = y*t1316; + t917 = t1257*t1362; + t912 = t1253*t1362; + t909 = t1240*t1362; + t907 = y*t1333; + t905 = t1272*t1180; + t900 = t1250*t1264; + t896 = t1270*t1418; + t894 = t1271*t1153; + t873 = z*t1365; + t870 = t1201*t1313; + t868 = z*t1367; + t843 = t1220*t1264; + t841 = t1213*t1264; + t838 = t1207*t1313; + t832 = t1211*t1327; + t831 = t1207*t1327; + t827 = t1201*t1327; + t823 = t1211*t1313; + coeff_dx_m1_m1_m1 = t918+t928+t1155+t1229*t1268+t1353+t1320+t896+t1157+ +t905+t1162+t894+t1360+t1164-t1541-t1542+t1512+t1538; + coeff_dx_0_m1_m1 = t963+t957+t997+t987+t907+t1359+t1434+t1436+t1475+t1483 +; + coeff_dx_p1_m1_m1 = t1456+t1261+t1036+t1262+t1028+t1319+t1337+t897+t1090+ +t1315+t1092+t1040+t1470+t1486; + coeff_dx_p2_m1_m1 = t1282+t1050+t1047+t1058+t1072+t1034+t1427+t1430+t1471 +; + coeff_dx_m1_0_m1 = t1304*t1274+t1045+t1046+(t1390-t1449)*x+t1291+t1149+ +t1335+t1357+t1483+t1498+t1519; + coeff_dx_0_0_m1 = t838+t875+t1338+t873+(t1223+t1206*t1271)*t1375+t1347+ +t1348+t1509+t1528; + coeff_dx_p1_0_m1 = t1279+t1066+t1301*t1274+t1227*t1375+t1119+t1358+t1458+ +t1469+t1500; + coeff_dx_p2_0_m1 = t1075+t1336+(-t1274+t1449*t1237)*x+t1284+t1474+t1513+ +t1516; + coeff_dx_m1_p1_m1 = t1055+t1093+t1050+t1038+t1260+t1077+t1070+t1277+t1119 ++t1328+t1456+t1490; + coeff_dx_0_p1_m1 = t868+t1028+t1290+t1025+t1459+t1469+t1511; + coeff_dx_p1_p1_m1 = t1225*t1300+t1278+t907+t1008+t1149+t1129+t1164+t838+ +t1350+t1455+t1488; + coeff_dx_p2_p1_m1 = t1006+t918+t1296+t1435+t1477+t1506+t1539; + coeff_dx_m1_p2_m1 = t1282+t1103+t1010+t984+t1267+t985+t1098+t1315+t1476; + coeff_dx_0_p2_m1 = t971+t978+t1284+t1533+t1535; + coeff_dx_p1_p2_m1 = t1155-t1262+t1136+t1289+t1437+t1462+t1506; + coeff_dx_p2_p2_m1 = t958+t961+t966+t1238*t1300+t1360+t1527+t1530; + coeff_dx_m1_m1_0 = t1121+t955+t956+(t1411-t1446)*x+t1304*t1271+t1291+ +t1345+t1349+t1475+t1491+t1514; + coeff_dx_0_m1_0 = t841+t831+t1344+(t1223+t1206*t1274)*t1372+t1292+t1356+ +t1528; + coeff_dx_p1_m1_0 = t1227*t1372+t1280+t1180+t1052+t1301*t1271+t1352+t1458+ +t1472+t1499; + coeff_dx_p2_m1_0 = t1009+t1339+(-t1271+t1446*t1237)*x+t1283+t1473+t1510+ +t1513; + coeff_dx_m1_0_0 = t1306*t1271+(t1229+t1432*t1250)*x+((RATIONAL(-75.0,8.0) +*t1272+RATIONAL(-25.0,8.0))*t1271+t1306)*t1274+t1292+t1304+t1484+t1495+t1508+ +t1509+t1521; + coeff_dx_0_0_0 = (t1247+RATIONAL(225.0,8.0)*t1269+(t1271+t1274)*t1205)* +t1272+(t1214*t1270+t1230+t1210*t1271+(t1214+t1219*t1271)*t1273+(t1219*t1270+ +RATIONAL(-125.0,4.0)*t1271+t1210)*t1274)*x+t1439+t1503+t1515; + coeff_dx_p1_0_0 = (t1224+t1432*t1234)*x+((RATIONAL(-225.0,8.0)*t1272+ +RATIONAL(25.0,8.0))*t1271+t1311)*t1274+t1311*t1271+t975+t1293+t1301+t1324+t1463 ++t1487; + coeff_dx_p2_0_0 = t968+(-1.0+t1432*t1240)*x+(t1220*t1271+t1251+(t1220+ +RATIONAL(75.0,8.0)*t1271)*t1274)*t1272+t1294+t1504+t1505; + coeff_dx_m1_p1_0 = t1271*t1410+t1280+(y+t1421)*x-t1271+t967+t1459+t1479+ +t1487+t1489+t1510; + coeff_dx_0_p1_0 = t1271*t1381+t827+t843+((t1208+t1203*t1274)*t1270+t1310* +y)*x+t1293+t1352+t1525; + coeff_dx_p1_p1_0 = t900+t1271*t1403+t1105+(t1244*t1270+t1229*y+(t1198* +t1273+t1231)*t1271)*x+t1271+t1295+t1345+t1356+t1481+t1495+t1515; + coeff_dx_p2_p1_0 = t1016+t1250*t1417+(t1452+t1251*t1270)*x+t1286+t1349+ +t1467+t1507; + coeff_dx_m1_p2_0 = t1167+t1450+t1297*t1271+t1281+t1473+t1499+t1535; + coeff_dx_0_p2_0 = t1340+t852+t832+t1310*t1450+t1355+t1504+t1523; + coeff_dx_p1_p2_0 = t1374+t1140+t1307*t1271+t1286+t1460+t1491+t1536; + coeff_dx_p2_p2_0 = t1341+(t1243*t1270+(t1237+t1192)*t1271)*x+t1305+t1514+ +t1532; + coeff_dx_m1_m1_p1 = t1263+t1180+t1047+t933+t1003+t985+t967+t1277+t1470+ +t1485+t1492+t1494+t1537; + coeff_dx_0_m1_p1 = t843+t1290+t978+t1040+t1465+t1472+t1531; + coeff_dx_p1_m1_p1 = t1121+t1131+t900-t1262+t831+t1320+t1225*t1308+t1359+ +t1455+t1457+t1517+t1529; + coeff_dx_p2_m1_p1 = t1296+t958+t909+t1016+t894+t1466+t1468; + coeff_dx_m1_0_p1 = t1279+t1104+(z+t1232*t1273)*x-t1274+t1274*t1410+t1038+ +t1324+t1465+t1479+t1516; + coeff_dx_0_0_p1 = t870+t1274*t1381+t868+((t1203*t1271+t1208)*t1273+t1312* +z)*x+t1358+t1463+t1525+t1526; + coeff_dx_p1_0_p1 = t1274+t1274*t1403+t1170+t1008+(t1229*z+t1244*t1273+( +t1231+t1198*t1270)*t1274)*x+t1295+t1335+t1348+t1480+t1484+t1503; + coeff_dx_p2_0_p1 = t1006+t1250*t1429+(t1251*t1273+t1243*z)*x+t1285+t1357+ +t1482+t1507; + coeff_dx_m1_p1_p1 = t1045+(t1227+t1244*t1272)*t1269+t1105+t1278+t956+ +t1170+t909+t1457+t1521+t1538+t1539; + coeff_dx_0_p1_p1 = t841+t957+t873+t833+RATIONAL(18.0,1.0)*t1331+t1350+ +t1439+t1480+t1481+t1529; + coeff_dx_p1_p1_p1 = t821+t1104+(t1229+RATIONAL(-18.0,1.0)*t1272)*t1269+ +t1066+t1337+t975+t1052+t827+((t1224*z+t1209*t1273)*t1271+(t1224*y+RATIONAL(16.0 +,1.0)*t1271+t1209*t1270)*t1274)*x+t870+t1485+t1489+t1490+t1511+t1531; + coeff_dx_p2_p1_p1 = t1009+t968+t1055+t1234*t1331+t1287+t1075+t933+t1427+ +t1464; + coeff_dx_m1_p2_p1 = t1131+t1288+t896+t966+t1140+t1462+t1466; + coeff_dx_0_p2_p1 = t963+t1299+t1437+t1482+t1536; + coeff_dx_p1_p2_p1 = t1262+t832+t971+t1167+t1103+t1224*t1267+t1461+t1464+ +t1494+t1540; + coeff_dx_p2_p2_p1 = t1268+t984+t1013+t1072+(t1188+t1190)*t1443+t1501+ +t1502+t1537; + coeff_dx_m1_m1_p2 = t937+t1163+t1090+t1328-t1391+t1397+t1433+t1471+t1476; + coeff_dx_0_m1_p2 = t1025+t932+t1283+t1522+t1533; + coeff_dx_p1_m1_p2 = t917+t1157+t1289+t1128+t1468+t1478+t1488; + coeff_dx_p2_m1_p2 = t905+t912+t1238*t1308+t1303+t1435+t1530; + coeff_dx_m1_0_p2 = t1097+t1448+t1297*t1274+t1281+t1474+t1500+t1522; + coeff_dx_0_0_p2 = t823+t1333+t1312*t1448+t1294+t1523; + coeff_dx_p1_0_p2 = t1137+t1229*t1448+t1307*t1274+t1285+t1460+t1498+t1524; + coeff_dx_p2_0_p2 = t1334+(t1243*t1273+(t1237+t1191)*t1274)*x+t1305+t1519+ +t1534; + coeff_dx_m1_p1_p2 = t912+t1288+t1137+t1129+t1353+t1477+t1478; + coeff_dx_0_p1_p2 = t917+t1299+t1436+t1467+t1524; + coeff_dx_p1_p1_p2 = t1224*t1265+t1287+t823+t1163+t1093+t932+t1097+t1461+ +t1486; + coeff_dx_p2_p1_p2 = t1070+t937+(t1187+t1183)*t1441+t1298+t1430+t1501; + coeff_dx_m1_p2_p2 = t1128+t1136+(t1246+t1380)*t1269+t1303+t1512+t1527; + coeff_dx_0_p2_p2 = t942+t1274*t1318+t1434+t1532+t1534; + coeff_dx_p1_p2_p2 = t1098+t1261+(t1236+t1400)*t1269+t1298+t1397+t1502; + coeff_dx_p2_p2_p2 = t1010+t1034+t1274*t1323+(t1157+(t1184+t1404)*t1274)*x ++t1433; -- cgit v1.2.3