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 | 223 |
1 files changed, 223 insertions, 0 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 new file mode 100644 index 0000000..b1a740e --- /dev/null +++ b/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dyy.compute.c @@ -0,0 +1,223 @@ +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 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; |