aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxy.compute.c
blob: f27238b31e4ff8abb1022e034e1e62446332ddef (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
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
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;