diff options
Diffstat (limited to 'src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dyy.compute.c')
-rw-r--r-- | src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dyy.compute.c | 441 |
1 files changed, 219 insertions, 222 deletions
diff --git a/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dyy.compute.c b/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dyy.compute.c index b1a740e..b9f99e1 100644 --- a/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dyy.compute.c +++ b/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dyy.compute.c @@ -1,223 +1,220 @@ -fp t1613; -fp t1675; -fp t1611; -fp t1645; -fp t1612; -fp t1614; -fp t1586; -fp t1644; -fp t1504; -fp t1674; -fp t1585; -fp t1615; -fp t1649; -fp t1584; -fp t1635; -fp t1673; -fp t1583; -fp t1646; -fp t1648; -fp t1588; -fp t1672; -fp t1650; -fp t1671; -fp t1555; -fp t1654; -fp t1670; -fp t1651; -fp t1604; -fp t1647; -fp t1669; -fp t1603; -fp t1605; -fp t1668; -fp t1587; -fp t1643; -fp t1494; -fp t1667; -fp t1666; -fp t1633; -fp t1665; -fp t1664; +fp t1747; +fp t1781; +fp t1780; +fp t1750; +fp t1850; +fp t1782; +fp t1849; +fp t1772; +fp t1811; +fp t1848; +fp t1775; +fp t1847; +fp t1746; +fp t1846; +fp t1748; +fp t1845; +fp t1771; +fp t1783; +fp t1749; +fp t1815; fp t1663; -fp t1609; -fp t1652; -fp t1662; -fp t1661; -fp t1660; -fp t1659; -fp t1589; -fp t1632; -fp t1636; -fp t1642; -fp t1579; -fp t1658; -fp t1580; -fp t1657; -fp t1601; -fp t1656; -fp t1655; -fp t1653; -fp t1641; -fp t1598; -fp t1640; -fp t1608; -fp t1591; -fp t1639; -fp t1568; -fp t1597; -fp t1573; -fp t1595; -fp t1637; -fp t1634; -fp t1599; -fp t1631; -fp t1602; -fp t1630; -fp t1592; -fp t1629; -fp t1578; -fp t1628; -fp t1596; -fp t1627; -fp t1626; -fp t1576; -fp t1577; -fp t1625; -fp t1624; -fp t1623; -fp t1622; -fp t1621; -fp t1620; -fp t1606; -fp t1607; -fp t1619; -fp t1618; - t1613 = y*y; - t1675 = RATIONAL(8.0,9.0)+t1613*RATIONAL(-4.0,3.0); - t1611 = RATIONAL(-1.0,3.0); - t1645 = t1611*x; - t1612 = x*x; - t1614 = t1612*x; - t1586 = RATIONAL(1.0,12.0); - t1644 = y*t1586; - t1504 = t1614*t1644; - t1674 = y*t1645+t1504; - t1585 = RATIONAL(1.0,48.0); - t1615 = t1612*t1612; - t1649 = y*t1615; - t1584 = RATIONAL(-1.0,48.0); - t1635 = t1612*t1584; - t1673 = t1585*t1649+y*t1635; - t1583 = RATIONAL(-1.0,24.0); - t1646 = x*t1583; - t1648 = y*t1614; - t1588 = RATIONAL(1.0,24.0); - t1672 = y*t1646+t1588*t1648; - t1650 = x*t1588; - t1671 = y*t1650+t1583*t1648; - t1555 = t1584*t1615; - t1654 = y*t1612; - t1670 = y*t1555+t1585*t1654; - t1651 = RATIONAL(2.0,3.0)*y; - t1604 = RATIONAL(-1.0,6.0); - t1647 = y*t1604; - t1669 = t1614*t1647+x*t1651; - t1603 = RATIONAL(1.0,6.0); - t1605 = RATIONAL(-2.0,3.0); - t1668 = t1605*t1654+t1603*t1649; - t1587 = RATIONAL(-1.0,12.0); - t1643 = y*t1587; - t1494 = t1614*t1643; - t1667 = x*t1644+t1494; - t1666 = t1583*t1654+t1588*t1649; - t1633 = t1603*t1614; - t1665 = (t1605*x+t1633)*y; - t1664 = t1615*t1647+t1612*t1651; - t1663 = t1615*t1644+t1611*t1654; - t1609 = RATIONAL(1.0,3.0); - t1652 = y*t1609; - t1662 = t1494+x*t1652; - t1661 = t1583*t1649+t1588*t1654; - t1660 = x*t1643+t1504; - t1659 = t1612*t1652+t1615*t1643; - t1589 = RATIONAL(-1.0,18.0); - t1632 = t1613*t1615; - t1636 = t1589*t1612+t1587*t1632; - t1642 = t1613*t1612; - t1579 = RATIONAL(1.0,18.0); - t1658 = t1636+t1579*t1615+t1586*t1642; - t1580 = RATIONAL(1.0,72.0); - t1657 = t1636+t1580*t1615+t1609*t1642; - t1601 = RATIONAL(-2.0,9.0); - t1656 = t1601*t1615+t1609*t1632+t1675*t1612; - t1655 = RATIONAL(-1.0,288.0)*t1615+t1613*t1635+RATIONAL(1.0,288.0)*t1612+ -t1585*t1632; - t1653 = x*t1613; - t1641 = RATIONAL(2.0,1.0)*t1613; - t1598 = RATIONAL(1.0,8.0); - t1640 = RATIONAL(5.0,48.0)*t1612+t1598*t1632; - t1608 = RATIONAL(-5.0,3.0); - t1591 = RATIONAL(-1.0,2.0); - t1639 = t1608*t1612+t1591*t1632; - t1568 = RATIONAL(-2.0,1.0)*t1613; - t1597 = RATIONAL(1.0,2.0); - t1573 = t1597*t1613; - t1595 = RATIONAL(-5.0,8.0); - t1637 = t1595*t1612; - t1634 = t1613*t1614; - t1599 = RATIONAL(1.0,4.0); - t1631 = t1599*t1613+RATIONAL(-5.0,24.0); - t1602 = RATIONAL(-1.0,4.0); - t1630 = t1602*t1613+RATIONAL(5.0,24.0); - t1592 = RATIONAL(-1.0,8.0); - t1629 = t1592*t1642+RATIONAL(-5.0,48.0)*t1615+t1640; - t1578 = RATIONAL(5.0,12.0); - t1628 = t1612*t1641+t1578*t1615+t1639; - t1596 = RATIONAL(4.0,3.0); - t1627 = t1609*t1615+RATIONAL(5.0,2.0)*t1642+t1568+t1596+t1639; - t1626 = t1555+t1587+t1573+t1613*t1637+t1640; - t1576 = RATIONAL(1.0,144.0); - t1577 = RATIONAL(-1.0,144.0); - t1625 = t1613*t1646+t1577*t1614+t1576*x+t1588*t1634+t1655; - t1624 = t1583*t1634+t1577*x+t1613*t1650+t1576*t1614+t1655; - t1623 = t1609*t1653+t1589*x+t1580*t1614+t1587*t1634+t1657; - t1622 = t1579*x+RATIONAL(-1.0,72.0)*t1614+t1586*t1634+t1613*t1645+t1657; - t1621 = t1601*t1614+t1609*t1634+t1656+t1675*x; - t1620 = RATIONAL(-8.0,9.0)*x+RATIONAL(2.0,9.0)*t1614+t1596*t1653+t1611* -t1634+t1656; - t1606 = RATIONAL(1.0,9.0); - t1607 = RATIONAL(-1.0,9.0); - t1619 = t1607*x+t1604*t1634+t1606*t1614+t1603*t1653+t1658; - t1618 = t1606*x+t1607*t1614+t1613*t1633+t1604*t1653+t1658; - coeffs_dyy->coeff_m2_m2 = t1624+t1670+t1672; - coeffs_dyy->coeff_m1_m2 = t1622+t1662+t1663; - coeffs_dyy->coeff_0_m2 = (RATIONAL(5.0,8.0)*t1612+t1592*t1615+t1591)*y+ -t1626; - coeffs_dyy->coeff_p1_m2 = t1623+t1663+t1674; - coeffs_dyy->coeff_p2_m2 = t1625+t1670+t1671; - coeffs_dyy->coeff_m2_m1 = t1618+t1666+t1667; - coeffs_dyy->coeff_m1_m1 = t1620+t1664+t1665; - coeffs_dyy->coeff_0_m1 = (1.0+t1599*t1615+RATIONAL(-5.0,4.0)*t1612)*y+ -t1627; - coeffs_dyy->coeff_p1_m1 = t1621+t1664+t1669; - coeffs_dyy->coeff_p2_m1 = t1619+t1660+t1666; - coeffs_dyy->coeff_m2_0 = t1630*t1614+t1631*x+t1629; - coeffs_dyy->coeff_m1_0 = (t1573+RATIONAL(-5.0,12.0))*t1614+(RATIONAL(5.0, -3.0)+t1568)*x+t1628; - coeffs_dyy->coeff_0_0 = RATIONAL(25.0,8.0)*t1612+RATIONAL(-5.0,2.0)+t1595 -*t1615+(RATIONAL(3.0,1.0)+RATIONAL(3.0,4.0)*t1615+RATIONAL(-15.0,4.0)*t1612)* -t1613; - coeffs_dyy->coeff_p1_0 = (t1591*t1613+t1578)*t1614+(t1641+t1608)*x+t1628; - coeffs_dyy->coeff_p2_0 = t1631*t1614+t1630*x+t1629; - coeffs_dyy->coeff_m2_p1 = t1618+t1660+t1661; - coeffs_dyy->coeff_m1_p1 = t1620+t1668+t1669; - coeffs_dyy->coeff_0_p1 = (t1602*t1615-1.0+RATIONAL(5.0,4.0)*t1612)*y+ -t1627; - coeffs_dyy->coeff_p1_p1 = t1621+t1665+t1668; - coeffs_dyy->coeff_p2_p1 = t1619+t1661+t1667; - coeffs_dyy->coeff_m2_p2 = t1624+t1671+t1673; - coeffs_dyy->coeff_m1_p2 = t1622+t1659+t1674; - coeffs_dyy->coeff_0_p2 = (t1637+t1597+t1598*t1615)*y+t1626; - coeffs_dyy->coeff_p1_p2 = t1623+t1659+t1662; - coeffs_dyy->coeff_p2_p2 = t1625+t1672+t1673; +fp t1824; +fp t1844; +fp t1819; +fp t1778; +fp t1843; +fp t1751; +fp t1822; +fp t1818; +fp t1842; +fp t1841; +fp t1670; +fp t1840; +fp t1821; +fp t1839; +fp t1776; +fp t1777; +fp t1814; +fp t1838; +fp t1837; +fp t1836; +fp t1835; +fp t1820; +fp t1823; +fp t1834; +fp t1833; +fp t1832; +fp t1724; +fp t1831; +fp t1830; +fp t1829; +fp t1805; +fp t1809; +fp t1754; +fp t1804; +fp t1753; +fp t1828; +fp t1752; +fp t1827; +fp t1826; +fp t1825; +fp t1817; +fp t1770; +fp t1761; +fp t1807; +fp t1816; +fp t1765; +fp t1740; +fp t1732; +fp t1813; +fp t1808; +fp t1810; +fp t1767; +fp t1799; +fp t1768; +fp t1798; +fp t1779; +fp t1757; +fp t1797; +fp t1762; +fp t1796; +fp t1763; +fp t1795; +fp t1766; +fp t1794; +fp t1793; +fp t1792; +fp t1791; +fp t1769; +fp t1790; +fp t1789; +fp t1788; +fp t1787; +fp t1786; + t1747 = RATIONAL(1.0,48.0); + t1781 = x*x; + t1780 = t1781*t1781; + t1750 = RATIONAL(-1.0,48.0); + t1850 = t1747*t1780+t1781*t1750; + t1782 = y*y; + t1849 = RATIONAL(8.0,9.0)+t1782*RATIONAL(-4.0,3.0); + t1772 = RATIONAL(1.0,3.0); + t1811 = t1772*t1782; + t1848 = t1811+RATIONAL(-2.0,9.0); + t1775 = RATIONAL(1.0,6.0); + t1847 = RATIONAL(-1.0,9.0)+t1775*t1782; + t1746 = RATIONAL(1.0,24.0); + t1846 = RATIONAL(-1.0,144.0)+t1746*t1782; + t1748 = RATIONAL(-1.0,24.0); + t1845 = RATIONAL(1.0,144.0)+t1782*t1748; + t1771 = RATIONAL(-1.0,3.0); + t1783 = x*t1781; + t1749 = RATIONAL(1.0,12.0); + t1815 = t1783*t1749; + t1663 = y*t1815; + t1824 = x*y; + t1844 = t1771*t1824+t1663; + t1819 = y*t1775; + t1778 = RATIONAL(-2.0,3.0); + t1843 = t1778*t1824+t1783*t1819; + t1751 = RATIONAL(-1.0,12.0); + t1822 = y*t1751; + t1818 = y*t1781; + t1842 = t1780*t1822+t1772*t1818; + t1841 = t1778*t1818+t1780*t1819; + t1670 = t1783*t1822; + t1840 = t1670+t1772*t1824; + t1821 = y*t1749; + t1839 = t1780*t1821+t1771*t1818; + t1776 = RATIONAL(2.0,3.0); + t1777 = RATIONAL(-1.0,6.0); + t1814 = t1777*t1783; + t1838 = t1776*t1824+y*t1814; + t1837 = t1776*t1818+t1777*y*t1780; + t1836 = t1850*y; + t1835 = x*t1822+t1663; + t1820 = y*t1748; + t1823 = y*t1746; + t1834 = t1783*t1820+x*t1823; + t1833 = t1746*t1818+t1780*t1820; + t1832 = x*t1820+t1783*t1823; + t1724 = t1750*t1780; + t1831 = t1747*t1818+y*t1724; + t1830 = t1748*t1818+t1780*t1823; + t1829 = x*t1821+t1670; + t1805 = t1781*t1782; + t1809 = t1782*t1751; + t1754 = RATIONAL(-1.0,18.0); + t1804 = t1780*t1809+t1754*t1781; + t1753 = RATIONAL(1.0,72.0); + t1828 = t1772*t1805+t1804+t1753*t1780; + t1752 = RATIONAL(1.0,18.0); + t1827 = t1749*t1805+t1752*t1780+t1804; + t1826 = t1848*t1780+t1849*t1781; + t1825 = RATIONAL(-1.0,288.0)*t1780+RATIONAL(1.0,288.0)*t1781+t1850*t1782; + t1817 = x*t1782; + t1770 = RATIONAL(-5.0,3.0); + t1761 = RATIONAL(-1.0,2.0); + t1807 = t1761*t1782; + t1816 = t1780*t1807+t1770*t1781; + t1765 = RATIONAL(1.0,2.0); + t1740 = t1765*t1782; + t1732 = RATIONAL(-2.0,1.0)*t1782; + t1813 = t1771*t1782; + t1808 = RATIONAL(1.0,8.0)*t1780; + t1810 = t1782*t1808+RATIONAL(5.0,48.0)*t1781; + t1767 = RATIONAL(-1.0,4.0); + t1799 = t1767*t1782+RATIONAL(5.0,24.0); + t1768 = RATIONAL(1.0,4.0); + t1798 = RATIONAL(-5.0,24.0)+t1768*t1782; + t1779 = RATIONAL(2.0,1.0); + t1757 = RATIONAL(5.0,12.0); + t1797 = t1779*t1805+t1757*t1780+t1816; + t1762 = RATIONAL(-1.0,8.0); + t1796 = t1762*t1805+RATIONAL(-5.0,48.0)*t1780+t1810; + t1763 = RATIONAL(4.0,3.0); + t1795 = t1732+t1763+RATIONAL(5.0,2.0)*t1805+t1772*t1780+t1816; + t1766 = RATIONAL(-5.0,8.0); + t1794 = t1766*t1805+t1751+t1724+t1740+t1810; + t1793 = t1782*t1815+RATIONAL(-1.0,72.0)*t1783+t1828+(t1813+t1752)*x; + t1792 = RATIONAL(-8.0,9.0)*x+t1763*t1817+t1826+(t1813+RATIONAL(2.0,9.0))* +t1783; + t1791 = t1825+t1845*t1783+t1846*x; + t1769 = RATIONAL(1.0,9.0); + t1790 = t1769*x+t1777*t1817+t1827+t1847*t1783; + t1789 = t1782*t1814+t1769*t1783+t1827+t1847*x; + t1788 = t1826+t1848*t1783+t1849*x; + t1787 = t1825+t1846*t1783+t1845*x; + t1786 = t1828+(t1753+t1809)*t1783+(t1754+t1811)*x; + coeffs_dyy->coeff_m2_m2 = t1791+t1831+t1832; + coeffs_dyy->coeff_m1_m2 = t1793+t1839+t1840; + coeffs_dyy->coeff_0_m2 = (RATIONAL(5.0,8.0)*t1781+t1761+t1762*t1780)*y+ +t1794; + coeffs_dyy->coeff_p1_m2 = t1786+t1839+t1844; + coeffs_dyy->coeff_p2_m2 = t1787+t1831+t1834; + coeffs_dyy->coeff_m2_m1 = t1790+t1829+t1830; + coeffs_dyy->coeff_m1_m1 = t1792+t1837+t1843; + coeffs_dyy->coeff_0_m1 = (RATIONAL(-5.0,4.0)*t1781+1.0+t1768*t1780)*y+ +t1795; + coeffs_dyy->coeff_p1_m1 = t1788+t1837+t1838; + coeffs_dyy->coeff_p2_m1 = t1789+t1830+t1835; + coeffs_dyy->coeff_m2_0 = t1799*t1783+t1798*x+t1796; + coeffs_dyy->coeff_m1_0 = (t1740+RATIONAL(-5.0,12.0))*t1783+(t1732+ +RATIONAL(5.0,3.0))*x+t1797; + coeffs_dyy->coeff_0_0 = RATIONAL(25.0,8.0)*t1781+t1766*t1780+RATIONAL( +-5.0,2.0)+(RATIONAL(3.0,1.0)+RATIONAL(3.0,4.0)*t1780+RATIONAL(-15.0,4.0)*t1781) +*t1782; + coeffs_dyy->coeff_p1_0 = (t1807+t1757)*t1783+(t1770+t1779*t1782)*x+t1797; + coeffs_dyy->coeff_p2_0 = t1798*t1783+t1799*x+t1796; + coeffs_dyy->coeff_m2_p1 = t1790+t1833+t1835; + coeffs_dyy->coeff_m1_p1 = t1792+t1838+t1841; + coeffs_dyy->coeff_0_p1 = (RATIONAL(5.0,4.0)*t1781-1.0+t1767*t1780)*y+ +t1795; + coeffs_dyy->coeff_p1_p1 = t1788+t1841+t1843; + coeffs_dyy->coeff_p2_p1 = t1789+t1829+t1833; + coeffs_dyy->coeff_m2_p2 = t1791+t1834+t1836; + coeffs_dyy->coeff_m1_p2 = t1793+t1842+t1844; + coeffs_dyy->coeff_0_p2 = (t1765+t1808+t1766*t1781)*y+t1794; + coeffs_dyy->coeff_p1_p2 = t1786+t1840+t1842; + coeffs_dyy->coeff_p2_p2 = t1787+t1832+t1836; |