aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dyy.compute.c
blob: b9f99e14d09b04a8b853356ee85cd289a3e9e8d9 (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
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 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;