aboutsummaryrefslogtreecommitdiff
path: root/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dyy.compute.c
blob: b1a740e18df55b8f84f24c4e56ca806e7cc362d3 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
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;