aboutsummaryrefslogtreecommitdiff
path: root/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeff-dy.compute.c
blob: 229c17ad5b047b04319af125f94e1f95f0b7a0fe (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
fp t174,
   t173,
   t189,
   t172,
   t188,
   t187,
   t186,
   t160,
   t185,
   t157,
   t182,
   t140,
   t159,
   t184,
   t183,
   t171,
   t181,
   t180,
   t161,
   t179,
   t145,
   t154,
   t178,
   t177,
   t156,
   t162,
   t176,
   t170,
   t169,
   t168,
   t167,
   t166,
   t165,
   t164,
   t163,
   t158,
   t155,
   t153,
   t152,
   t151,
   t150,
   t148,
   t147,
   t146,
   t144,
   t143,
   t139,
   t138,
   t137,
   t135,
   t133;
      t174 = x*x;
      t173 = t174*x;
      t189 = t173*y;
      t172 = y*y;
      t188 = x*t172;
      t187 = t174*y;
      t186 = x+t173;
      t160 = RATIONAL(9.0,4.0);
      t185 = t160*t172;
      t157 = RATIONAL(-9.0,4.0);
      t182 = t172*t173;
      t140 = t157*t182;
      t159 = RATIONAL(-5.0,1.0);
      t184 = t140+t159*t187;
      t183 = RATIONAL(-27.0,4.0)*t182+RATIONAL(-10.0,1.0)*t187;
      t171 = RATIONAL(-2.0,1.0);
      t181 = t171*t187+t140;
      t180 = t172*t174;
      t161 = RATIONAL(3.0,4.0);
      t179 = t161*t172;
      t145 = t160*t182;
      t154 = RATIONAL(5.0,2.0);
      t178 = t154*t187+t145;
      t177 = y*t186;
      t156 = RATIONAL(4.0,1.0);
      t162 = RATIONAL(-1.0,4.0);
      t176 = t156*t187+t145+t162*x;
      t170 = RATIONAL(-1.0,2.0);
      t169 = RATIONAL(1.0,4.0);
      t168 = RATIONAL(-5.0,2.0);
      t167 = RATIONAL(-3.0,4.0);
      t166 = RATIONAL(3.0,2.0);
      t165 = RATIONAL(2.0,1.0);
      t164 = RATIONAL(-3.0,2.0);
      t163 = RATIONAL(1.0,2.0);
      t158 = RATIONAL(3.0,1.0);
      t155 = RATIONAL(-3.0,1.0);
      t153 = RATIONAL(9.0,2.0);
      t152 = RATIONAL(-9.0,2.0);
      t151 = t167*t173;
      t150 = t161*t173;
      t148 = t162*t173;
      t147 = t169*x;
      t146 = t169*t173;
      t144 = x*t179;
      t143 = t167*t188;
      t139 = t167*t182;
      t138 = t157*t188;
      t137 = x*t185;
      t135 = t173*t179;
      t133 = RATIONAL(27.0,4.0)*t182;
      coeff_dy_m1_m1 = t144+t135+t147+t146+(t164*t172+t170)*t174+(t165*t174-
t186)*y;
      coeff_dy_0_m1 = RATIONAL(5.0,4.0)*t174+t151+t170+(t164+RATIONAL(15.0,4.0)
*t174)*t172+(t158*t173+t165)*y+t184;
      coeff_dy_p1_m1 = t155*t180+t143+t150-t174+(t155*t173+x)*y+t176;
      coeff_dy_p2_m1 = t148+t189+t139+(t179+t169-y)*t174;
      coeff_dy_m1_0 = t153*t180+t138+t154*t177+t184;
      coeff_dy_0_0 = t133+(RATIONAL(-45.0,4.0)*t174+t153)*t172+(t159+RATIONAL(
-15.0,2.0)*t173+RATIONAL(25.0,2.0)*t174)*y;
      coeff_dy_p1_0 = t137+RATIONAL(9.0,1.0)*t180+(RATIONAL(15.0,2.0)*t173+t168
*x)*y+t183;
      coeff_dy_p2_0 = t168*t189+t157*t180+t178;
      coeff_dy_m1_p1 = t148+t137+(t152*t172+t163)*t174+t171*t177+t176;
      coeff_dy_0_p1 = t163+t150+RATIONAL(-5.0,4.0)*t174+(t152+RATIONAL(45.0,4.0
)*t174)*t172+(RATIONAL(6.0,1.0)*t173+t156)*y+t183;
      coeff_dy_p1_p1 = t151+RATIONAL(-9.0,1.0)*t180+t133+t147+t174+t138+(
RATIONAL(-6.0,1.0)*t173+RATIONAL(8.0,1.0)*t174+t165*x)*y;
      coeff_dy_p2_p1 = t146+t165*t189+(t185+t162)*t174+t181;
      coeff_dy_m1_p2 = t166*t180+t139+t143+(-t174+t186*t163)*y;
      coeff_dy_0_p2 = t164*t189-y+(t166+RATIONAL(-15.0,4.0)*t174)*t172+t178;
      coeff_dy_p1_p2 = t158*t180+t144+(t170*x+t166*t173)*y+t181;
      coeff_dy_p2_p2 = t135+t167*t180+(t170*t173+t163*t174)*y;