1 files changed, 265 insertions, 256 deletions
diff --git a/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxy.compute.c b/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxy.compute.c
index c0d0630..f27238b 100644
--- a/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxy.compute.c
+++ b/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxy.compute.c
@@ -1,256 +1,265 @@
-fp t1387;
-fp t1428;
-fp t1430;
-fp t1390;
-fp t1485;
-fp t1396;
-fp t1484;
-fp t1397;
-fp t1388;
-fp t1483;
-fp t1426;
-fp t1403;
-fp t1482;
-fp t1458;
-fp t1427;
-fp t1429;
-fp t1389;
-fp t1361;
-fp t1425;
-fp t1481;
-fp t1343;
-fp t1460;
-fp t1480;
-fp t1450;
-fp t1418;
-fp t1398;
-fp t1395;
-fp t1344;
-fp t1479;
-fp t1448;
-fp t1414;
-fp t1417;
-fp t1423;
-fp t1369;
-fp t1401;
-fp t1478;
-fp t1459;
-fp t1416;
-fp t1422;
-fp t1421;
-fp t1402;
-fp t1477;
-fp t1366;
-fp t1367;
-fp t1413;
-fp t1452;
-fp t1476;
-fp t1391;
-fp t1475;
-fp t1350;
-fp t1474;
-fp t1356;
-fp t1400;
-fp t1473;
-fp t1373;
-fp t1472;
-fp t1362;
-fp t1471;
-fp t1351;
-fp t1380;
-fp t1470;
-fp t1455;
-fp t1469;
-fp t1399;
-fp t1355;
-fp t1468;
-fp t1386;
-fp t1415;
-fp t1467;
-fp t1466;
-fp t1449;
-fp t1347;
-fp t1465;
-fp t1461;
-fp t1464;
-fp t1412;
-fp t1463;
-fp t1411;
-fp t1447;
-fp t1462;
-fp t1457;
-fp t1456;
-fp t1454;
-fp t1453;
-fp t1393;
-fp t1451;
-fp t1446;
-fp t1445;
-fp t1409;
-fp t1420;
-fp t1444;
-fp t1394;
-fp t1443;
-fp t1442;
-fp t1441;
-fp t1440;
-fp t1439;
-fp t1438;
-fp t1437;
-fp t1436;
-fp t1435;
-fp t1434;
-fp t1433;
-fp t1424;
-fp t1419;
-fp t1410;
-fp t1408;
-fp t1407;
-fp t1406;
-fp t1405;
-fp t1404;
-fp t1392;
-      t1387 = RATIONAL(1.0,144.0);
-      t1428 = y*y;
-      t1430 = y*t1428;
-      t1390 = RATIONAL(-1.0,72.0);
-      t1485 = t1387*y+t1390*t1430;
-      t1396 = RATIONAL(-1.0,48.0);
-      t1484 = t1396*t1428+t1387;
-      t1397 = RATIONAL(1.0,48.0);
-      t1388 = RATIONAL(-1.0,144.0);
-      t1483 = t1397*t1428+t1388;
-      t1426 = RATIONAL(-8.0,9.0);
-      t1403 = RATIONAL(4.0,9.0);
-      t1482 = t1426*y+t1403*t1430;
-      t1458 = x*t1430;
-      t1427 = x*x;
-      t1429 = t1427*x;
-      t1389 = RATIONAL(1.0,18.0);
-      t1361 = t1389*t1429;
-      t1425 = RATIONAL(2.0,9.0);
-      t1481 = y*t1361+t1425*t1458;
-      t1343 = t1389*t1430;
-      t1460 = y*t1429;
-      t1480 = x*t1343+t1425*t1460;
-      t1450 = t1428*t1429;
-      t1418 = RATIONAL(-1.0,9.0);
-      t1398 = RATIONAL(1.0,12.0);
-      t1395 = RATIONAL(-1.0,24.0);
-      t1344 = t1395*t1428;
-      t1479 = t1398*t1450+t1418*t1429+(t1389+t1344)*x;
-      t1448 = t1427*t1430;
-      t1414 = RATIONAL(1.0,6.0);
-      t1417 = RATIONAL(1.0,9.0);
-      t1423 = RATIONAL(-1.0,3.0);
-      t1369 = t1423*t1427;
-      t1401 = RATIONAL(-1.0,18.0);
-      t1478 = t1401*t1430+t1414*t1448+(t1369+t1417)*y;
-      t1459 = x*t1428;
-      t1416 = RATIONAL(-2.0,3.0);
-      t1422 = RATIONAL(8.0,9.0);
-      t1421 = RATIONAL(1.0,3.0);
-      t1402 = RATIONAL(-4.0,9.0);
-      t1477 = t1416*t1459+t1421*t1450+t1422*x+t1402*t1429;
-      t1366 = t1421*t1427;
-      t1367 = t1421*t1428;
-      t1413 = RATIONAL(-1.0,4.0);
-      t1452 = t1413*t1427;
-      t1476 = t1366+t1367+t1428*t1452+t1402;
-      t1391 = RATIONAL(1.0,72.0);
-      t1475 = t1391*t1429+t1395*t1450+t1483*x;
-      t1350 = t1397*t1427;
-      t1474 = t1391*t1430+t1395*t1448+(t1350+t1388)*y;
-      t1356 = t1396*t1427;
-      t1400 = RATIONAL(1.0,24.0);
-      t1473 = t1400*t1448+y*t1356+t1485;
-      t1373 = t1423*t1428;
-      t1472 = t1401*t1429+t1414*t1450+(t1373+t1417)*x;
-      t1362 = t1390*t1429;
-      t1471 = t1362+t1400*t1450+t1484*x;
-      t1351 = t1395*t1427;
-      t1380 = t1418*t1430;
-      t1470 = t1398*t1448+t1380+(t1351+t1389)*y;
-      t1455 = t1429*t1430;
-      t1469 = y*t1362+RATIONAL(1.0,36.0)*t1455+t1485*x;
-      t1399 = RATIONAL(-1.0,12.0);
-      t1355 = t1400*t1427;
-      t1468 = t1417*t1430+t1399*t1448+(t1355+t1401)*y;
-      t1386 = t1418*x;
-      t1415 = RATIONAL(-1.0,6.0);
-      t1467 = t1386+t1361+t1415*t1450+x*t1367;
-      t1466 = t1343+t1415*t1448+(t1366+t1418)*y;
-      t1449 = t1399*t1429;
-      t1347 = t1400*t1428;
-      t1465 = t1417*t1429+t1428*t1449+(t1401+t1347)*x;
-      t1461 = x*y;
-      t1464 = RATIONAL(16.0,9.0)*t1461+t1426*t1458+t1482*t1429;
-      t1412 = RATIONAL(2.0,3.0);
-      t1463 = t1412*t1459+t1423*t1450+t1426*x+t1403*t1429;
-      t1411 = RATIONAL(1.0,4.0);
-      t1447 = t1411*t1428;
-      t1462 = t1427*t1447+t1369+t1373+t1403;
-      t1457 = y*t1427;
-      t1456 = t1429*t1380+y*t1386;
-      t1454 = t1428*t1427;
-      t1453 = t1416*t1455+RATIONAL(-10.0,3.0)*t1461;
-      t1393 = RATIONAL(5.0,24.0);
-      t1451 = t1393*t1461+t1414*t1455;
-      t1446 = RATIONAL(1.0,8.0)*t1454+t1389+t1456;
-      t1445 = t1401+RATIONAL(-1.0,8.0)*t1454+t1456;
-      t1409 = RATIONAL(4.0,3.0);
-      t1420 = RATIONAL(5.0,3.0);
-      t1444 = t1420*t1458+t1409*t1460+t1453;
-      t1394 = RATIONAL(-5.0,12.0);
-      t1443 = t1394*t1458+y*t1449+t1451;
-      t1442 = t1420*t1460+t1409*t1458+t1453;
-      t1441 = t1399*t1458+t1394*t1460+t1451;
-      t1440 = t1350+RATIONAL(-1.0,16.0)*t1454+t1469+t1483;
-      t1439 = RATIONAL(1.0,16.0)*t1454+t1356+t1469+t1484;
-      t1438 = t1351+t1415*t1428+t1446+t1481;
-      t1437 = t1344+t1415*t1427+t1446+t1480;
-      t1436 = t1347+t1414*t1427+t1445+t1480;
-      t1435 = t1402*t1430+t1422*y+t1416*t1457+t1421*t1448+t1464;
-      t1434 = t1414*t1428+t1355+t1445+t1481;
-      t1433 = t1423*t1448+t1412*t1457+t1464+t1482;
-      t1424 = RATIONAL(5.0,4.0);
-      t1419 = RATIONAL(-5.0,3.0);
-      t1410 = RATIONAL(1.0,2.0);
-      t1408 = RATIONAL(-5.0,8.0);
-      t1407 = RATIONAL(-5.0,4.0);
-      t1406 = RATIONAL(-5.0,2.0);
-      t1405 = RATIONAL(-1.0,2.0);
-      t1404 = RATIONAL(5.0,8.0);
-      t1392 = RATIONAL(-5.0,24.0);
-      coeffs_dxy->coeff_m2_m2 = t1439+t1474+t1475;
-      coeffs_dxy->coeff_m1_m2 = t1434+t1470+t1472;
-      coeffs_dxy->coeff_0_m2 = (t1413*t1428+t1398)*t1429+(t1404*t1428+t1392)*x+
-t1443;
-      coeffs_dxy->coeff_p1_m2 = t1438+t1468+t1472;
-      coeffs_dxy->coeff_p2_m2 = t1440+t1473+t1475;
-      coeffs_dxy->coeff_m2_m1 = t1436+t1478+t1479;
-      coeffs_dxy->coeff_m1_m1 = t1433+t1462+t1463;
-      coeffs_dxy->coeff_0_m1 = (t1410*t1428+t1416)*t1429+(t1420+t1407*t1428)*x+
-t1444;
-      coeffs_dxy->coeff_p1_m1 = t1435+t1463+t1476;
-      coeffs_dxy->coeff_p2_m1 = t1437+t1466+t1479;
-      coeffs_dxy->coeff_m2_0 = (t1452+t1398)*t1430+(t1392+t1404*t1427)*y+t1441;
-      coeffs_dxy->coeff_m1_0 = (t1410*t1427+t1416)*t1430+(t1407*t1427+t1420)*y+
-t1442;
-      coeffs_dxy->coeff_0_0 = (t1406*y+t1430)*t1429+(RATIONAL(25.0,4.0)*y+t1406
-*t1430)*x;
-      coeffs_dxy->coeff_p1_0 = (t1405*t1427+t1412)*t1430+(t1424*t1427+t1419)*y+
-t1442;
-      coeffs_dxy->coeff_p2_0 = (t1399+t1411*t1427)*t1430+(t1393+t1408*t1427)*y+
-t1441;
-      coeffs_dxy->coeff_m2_p1 = t1437+t1465+t1478;
-      coeffs_dxy->coeff_m1_p1 = t1433+t1476+t1477;
-      coeffs_dxy->coeff_0_p1 = (t1405*t1428+t1412)*t1429+(t1419+t1424*t1428)*x+
-t1444;
-      coeffs_dxy->coeff_p1_p1 = t1435+t1462+t1477;
-      coeffs_dxy->coeff_p2_p1 = t1436+t1465+t1466;
-      coeffs_dxy->coeff_m2_p2 = t1440+t1471+t1474;
-      coeffs_dxy->coeff_m1_p2 = t1438+t1467+t1470;
-      coeffs_dxy->coeff_0_p2 = (t1447+t1399)*t1429+(t1393+t1408*t1428)*x+t1443;
-      coeffs_dxy->coeff_p1_p2 = t1434+t1467+t1468;
-      coeffs_dxy->coeff_p2_p2 = t1439+t1471+t1473;
+fp t1564;
+fp t1553;
+fp t1593;
+fp t1595;
+fp t1653;
+fp t1569;
+fp t1586;
+fp t1594;
+fp t1596;
+fp t1652;
+fp t1559;
+fp t1651;
+fp t1568;
+fp t1583;
+fp t1531;
+fp t1650;
+fp t1582;
+fp t1538;
+fp t1649;
+fp t1554;
+fp t1556;
+fp t1648;
+fp t1624;
+fp t1585;
+fp t1561;
+fp t1528;
+fp t1647;
+fp t1525;
+fp t1623;
+fp t1646;
+fp t1584;
+fp t1625;
+fp t1591;
+fp t1645;
+fp t1589;
+fp t1644;
+fp t1626;
+fp t1617;
+fp t1643;
+fp t1642;
+fp t1558;
+fp t1557;
+fp t1515;
+fp t1587;
+fp t1545;
+fp t1641;
+fp t1563;
+fp t1579;
+fp t1588;
+fp t1544;
+fp t1640;
+fp t1560;
+fp t1613;
+fp t1555;
+fp t1517;
+fp t1639;
+fp t1546;
+fp t1590;
+fp t1542;
+fp t1638;
+fp t1526;
+fp t1637;
+fp t1615;
+fp t1536;
+fp t1636;
+fp t1532;
+fp t1635;
+fp t1619;
+fp t1514;
+fp t1634;
+fp t1513;
+fp t1562;
+fp t1633;
+fp t1632;
+fp t1511;
+fp t1631;
+fp t1621;
+fp t1630;
+fp t1627;
+fp t1629;
+fp t1520;
+fp t1628;
+fp t1566;
+fp t1622;
+fp t1620;
+fp t1618;
+fp t1616;
+fp t1576;
+fp t1614;
+fp t1565;
+fp t1612;
+fp t1611;
+fp t1610;
+fp t1573;
+fp t1580;
+fp t1609;
+fp t1608;
+fp t1607;
+fp t1606;
+fp t1605;
+fp t1604;
+fp t1603;
+fp t1602;
+fp t1601;
+fp t1577;
+fp t1600;
+fp t1599;
+fp t1592;
+fp t1581;
+fp t1578;
+fp t1575;
+fp t1574;
+fp t1572;
+fp t1571;
+fp t1570;
+fp t1567;
+      t1564 = RATIONAL(-1.0,72.0);
+      t1553 = RATIONAL(1.0,144.0);
+      t1593 = x*x;
+      t1595 = t1593*x;
+      t1653 = t1553*x+t1564*t1595;
+      t1569 = RATIONAL(4.0,9.0);
+      t1586 = RATIONAL(-8.0,9.0);
+      t1594 = y*y;
+      t1596 = y*t1594;
+      t1652 = t1569*t1596+t1586*y;
+      t1559 = RATIONAL(-1.0,48.0);
+      t1651 = t1553+t1559*t1593;
+      t1568 = RATIONAL(-4.0,9.0);
+      t1583 = RATIONAL(1.0,3.0);
+      t1531 = t1583*t1594;
+      t1650 = t1531+t1568;
+      t1582 = RATIONAL(-1.0,3.0);
+      t1538 = t1582*t1594;
+      t1649 = t1538+t1569;
+      t1554 = RATIONAL(-1.0,144.0);
+      t1556 = RATIONAL(1.0,48.0);
+      t1648 = t1554+t1556*t1593;
+      t1624 = y*t1595;
+      t1585 = RATIONAL(2.0,9.0);
+      t1561 = RATIONAL(1.0,18.0);
+      t1528 = t1561*t1596;
+      t1647 = x*t1528+t1585*t1624;
+      t1525 = t1561*t1595;
+      t1623 = x*t1596;
+      t1646 = t1585*t1623+y*t1525;
+      t1584 = RATIONAL(8.0,9.0);
+      t1625 = x*t1594;
+      t1591 = RATIONAL(-2.0,3.0);
+      t1645 = t1584*x+t1591*t1625+t1650*t1595;
+      t1589 = RATIONAL(2.0,3.0);
+      t1644 = t1589*t1625+t1586*x+t1649*t1595;
+      t1626 = y*t1593;
+      t1617 = t1593*t1596;
+      t1643 = t1582*t1617+t1589*t1626+t1652;
+      t1642 = t1583*t1617+t1568*t1596+t1584*y+t1591*t1626;
+      t1558 = RATIONAL(1.0,12.0);
+      t1557 = RATIONAL(-1.0,24.0);
+      t1515 = t1557*t1593;
+      t1587 = RATIONAL(-1.0,9.0);
+      t1545 = t1587*t1596;
+      t1641 = t1558*t1617+t1545+(t1515+t1561)*y;
+      t1563 = RATIONAL(-1.0,18.0);
+      t1579 = RATIONAL(1.0,9.0);
+      t1588 = RATIONAL(1.0,6.0);
+      t1544 = t1588*t1594;
+      t1640 = (t1563+t1544)*t1595+(t1579+t1538)*x;
+      t1560 = RATIONAL(-1.0,12.0);
+      t1613 = t1560*t1596;
+      t1555 = RATIONAL(1.0,24.0);
+      t1517 = t1555*t1593;
+      t1639 = t1593*t1613+t1579*t1596+(t1517+t1563)*y;
+      t1546 = t1587*x;
+      t1590 = RATIONAL(-1.0,6.0);
+      t1542 = t1590*t1594;
+      t1638 = t1546+t1525+t1595*t1542+x*t1531;
+      t1526 = t1557*t1594;
+      t1637 = (t1558*t1594+t1587)*t1595+(t1561+t1526)*x;
+      t1615 = t1588*t1596;
+      t1536 = t1582*t1593;
+      t1636 = t1593*t1615+t1563*t1596+(t1536+t1579)*y;
+      t1532 = t1583*t1593;
+      t1635 = t1528+t1590*t1617+(t1587+t1532)*y;
+      t1619 = t1560*t1595;
+      t1514 = t1555*t1594;
+      t1634 = t1594*t1619+t1579*t1595+(t1563+t1514)*x;
+      t1513 = t1556*t1594;
+      t1562 = RATIONAL(1.0,72.0);
+      t1633 = (t1562+t1526)*t1595+(t1554+t1513)*x;
+      t1632 = t1562*t1596+t1557*t1617+t1648*y;
+      t1511 = t1564*t1596;
+      t1631 = t1511+t1555*t1617+t1651*y;
+      t1621 = t1595*t1596;
+      t1630 = x*t1511+RATIONAL(1.0,36.0)*t1621+t1653*y;
+      t1627 = x*y;
+      t1629 = RATIONAL(16.0,9.0)*t1627+t1586*t1623+t1652*t1595;
+      t1520 = t1559*t1594;
+      t1628 = t1595*t1514+x*t1520+t1653;
+      t1566 = RATIONAL(5.0,24.0);
+      t1622 = t1566*t1627+t1595*t1615;
+      t1620 = t1595*t1545+y*t1546;
+      t1618 = t1593*t1594;
+      t1616 = RATIONAL(-10.0,3.0)*t1627+t1591*t1621;
+      t1576 = RATIONAL(-1.0,4.0);
+      t1614 = t1576*t1594;
+      t1565 = RATIONAL(-5.0,12.0);
+      t1612 = x*t1613+t1565*t1624+t1622;
+      t1611 = RATIONAL(1.0,8.0)*t1618+t1561+t1620;
+      t1610 = RATIONAL(-1.0,8.0)*t1618+t1563+t1620;
+      t1573 = RATIONAL(4.0,3.0);
+      t1580 = RATIONAL(5.0,3.0);
+      t1609 = t1573*t1623+t1580*t1624+t1616;
+      t1608 = t1580*t1623+t1573*t1624+t1616;
+      t1607 = y*t1619+t1565*t1623+t1622;
+      t1606 = t1513+RATIONAL(-1.0,16.0)*t1618+t1630+t1648;
+      t1605 = t1520+RATIONAL(1.0,16.0)*t1618+t1630+t1651;
+      t1604 = t1515+t1542+t1611+t1646;
+      t1603 = t1590*t1593+t1526+t1611+t1647;
+      t1602 = t1514+t1588*t1593+t1610+t1647;
+      t1601 = t1544+t1517+t1610+t1646;
+      t1577 = RATIONAL(1.0,4.0);
+      t1600 = t1577*t1618+t1536+t1629+t1649;
+      t1599 = t1532+t1593*t1614+t1629+t1650;
+      t1592 = RATIONAL(5.0,8.0);
+      t1581 = RATIONAL(-5.0,3.0);
+      t1578 = RATIONAL(-5.0,4.0);
+      t1575 = RATIONAL(-5.0,8.0);
+      t1574 = RATIONAL(1.0,2.0);
+      t1572 = RATIONAL(5.0,4.0);
+      t1571 = RATIONAL(-5.0,2.0);
+      t1570 = RATIONAL(-1.0,2.0);
+      t1567 = RATIONAL(-5.0,24.0);
+      coeffs_dxy->coeff_m2_m2 = t1605+t1632+t1633;
+      coeffs_dxy->coeff_m1_m2 = t1601+t1640+t1641;
+      coeffs_dxy->coeff_0_m2 = (t1558+t1614)*t1595+(t1592*t1594+t1567)*x+t1607;
+      coeffs_dxy->coeff_p1_m2 = t1604+t1639+t1640;
+      coeffs_dxy->coeff_p2_m2 = t1606+t1631+t1633;
+      coeffs_dxy->coeff_m2_m1 = t1602+t1636+t1637;
+      coeffs_dxy->coeff_m1_m1 = t1600+t1643+t1644;
+      coeffs_dxy->coeff_0_m1 = (t1591+t1574*t1594)*t1595+(t1578*t1594+t1580)*x+
+t1608;
+      coeffs_dxy->coeff_p1_m1 = t1599+t1642+t1644;
+      coeffs_dxy->coeff_p2_m1 = t1603+t1635+t1637;
+      coeffs_dxy->coeff_m2_0 = (t1558+t1576*t1593)*t1596+(t1592*t1593+t1567)*y+
+t1612;
+      coeffs_dxy->coeff_m1_0 = (t1591+t1574*t1593)*t1596+(t1578*t1593+t1580)*y+
+t1609;
+      coeffs_dxy->coeff_0_0 = (t1571*y+t1596)*t1595+(t1571*t1596+RATIONAL(25.0,
+4.0)*y)*x;
+      coeffs_dxy->coeff_p1_0 = (t1570*t1593+t1589)*t1596+(t1581+t1572*t1593)*y+
+t1609;
+      coeffs_dxy->coeff_p2_0 = (t1560+t1577*t1593)*t1596+(t1575*t1593+t1566)*y+
+t1612;
+      coeffs_dxy->coeff_m2_p1 = t1603+t1634+t1636;
+      coeffs_dxy->coeff_m1_p1 = t1599+t1643+t1645;
+      coeffs_dxy->coeff_0_p1 = (t1589+t1570*t1594)*t1595+(t1572*t1594+t1581)*x+
+t1608;
+      coeffs_dxy->coeff_p1_p1 = t1600+t1642+t1645;
+      coeffs_dxy->coeff_p2_p1 = t1602+t1634+t1635;
+      coeffs_dxy->coeff_m2_p2 = t1606+t1628+t1632;
+      coeffs_dxy->coeff_m1_p2 = t1604+t1638+t1641;
+      coeffs_dxy->coeff_0_p2 = (t1577*t1594+t1560)*t1595+(t1575*t1594+t1566)*x+
+t1607;
+      coeffs_dxy->coeff_p1_p2 = t1601+t1638+t1639;
+      coeffs_dxy->coeff_p2_p2 = t1605+t1628+t1631;