aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dy.compute.c
blob: 68bb28c3452e3fbc47701dc1b2cac8a62e5ce482 (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
fp t283;
fp t290;
fp t291;
fp t304;
fp t257;
fp t239;
fp t314;
fp t306;
fp t334;
fp t286;
fp t333;
fp t275;
fp t258;
fp t240;
fp t273;
fp t259;
fp t241;
fp t332;
fp t277;
fp t270;
fp t247;
fp t331;
fp t279;
fp t269;
fp t246;
fp t330;
fp t284;
fp t305;
fp t262;
fp t243;
fp t287;
fp t312;
fp t271;
fp t329;
fp t313;
fp t260;
fp t244;
fp t272;
fp t328;
fp t327;
fp t289;
fp t249;
fp t326;
fp t325;
fp t250;
fp t261;
fp t324;
fp t251;
fp t323;
fp t282;
fp t253;
fp t256;
fp t322;
fp t311;
fp t254;
fp t285;
fp t310;
fp t266;
fp t321;
fp t255;
fp t265;
fp t320;
fp t319;
fp t318;
fp t264;
fp t267;
fp t317;
fp t268;
fp t316;
fp t315;
fp t309;
fp t308;
fp t307;
fp t303;
fp t302;
fp t301;
fp t300;
fp t299;
fp t298;
fp t297;
fp t296;
fp t295;
fp t294;
fp t293;
fp t292;
      t283 = RATIONAL(1.0,4.0);
      t290 = z*z;
      t291 = x*x;
      t304 = t291*t290;
      t257 = t283*t304;
      t239 = y*t257;
      t314 = x*z;
      t306 = y*t314;
      t334 = t239+t283*t306;
      t286 = RATIONAL(-1.0,4.0);
      t333 = t239+t286*t306;
      t275 = t283*t291;
      t258 = z*t275;
      t240 = y*t258;
      t273 = t283*t290;
      t259 = x*t273;
      t241 = y*t259;
      t332 = t240+t241;
      t277 = t286*t290;
      t270 = x*t277;
      t247 = y*t270;
      t331 = t240+t247;
      t279 = t286*t291;
      t269 = z*t279;
      t246 = y*t269;
      t330 = t241+t246;
      t284 = RATIONAL(1.0,2.0);
      t305 = t284*t291;
      t262 = y*t305;
      t243 = z*t262;
      t287 = RATIONAL(-1.0,2.0);
      t312 = t287*y;
      t271 = z*t312;
      t329 = t243+t271;
      t313 = t284*y;
      t260 = t290*t313;
      t244 = x*t260;
      t272 = x*t312;
      t328 = t244+t272;
      t327 = t246+t247;
      t289 = y*t304;
      t249 = t287*t289;
      t326 = t249+t260;
      t325 = t249+t262;
      t250 = t291*t271;
      t261 = z*t313;
      t324 = t250+t261;
      t251 = t290*t272;
      t323 = t251+x*t313;
      t282 = RATIONAL(1.0,8.0);
      t253 = t282*t304;
      t256 = t282*t314;
      t322 = t253+t256;
      t311 = t291*z;
      t254 = t282*t311;
      t285 = RATIONAL(-1.0,8.0);
      t310 = x*t290;
      t266 = t285*t310;
      t321 = t254+t266;
      t255 = t282*t310;
      t265 = t285*t311;
      t320 = t255+t265;
      t319 = t257+t277;
      t318 = t257+t279;
      t264 = t285*t304;
      t267 = t285*t314;
      t317 = t264+t267;
      t268 = t286*t304;
      t316 = t268+t273;
      t315 = t268+t275;
      t309 = t289+y;
      t308 = -y+t284;
      t307 = -y+t287;
      t303 = x*t261+t249;
      t302 = t249+x*t271;
      t301 = x-t310;
      t300 = z-t311;
      t299 = t254+t255+t334;
      t298 = t265+t266+t334;
      t297 = t253+t267+t333;
      t296 = t256+t264+t333;
      t295 = t258+t286*z+t326;
      t294 = t259+t286*x+t325;
      t293 = t269+t283*z+t326;
      t292 = t270+t283*x+t325;
      coeffs_dy->coeff_m1_m1_m1 = t299+t317+t327;
      coeffs_dy->coeff_0_m1_m1 = t293+t319+t329;
      coeffs_dy->coeff_p1_m1_m1 = t296+t321+t330;
      coeffs_dy->coeff_m1_0_m1 = t244+t243+t302;
      coeffs_dy->coeff_0_0_m1 = t289+(-t290+t300)*y;
      coeffs_dy->coeff_p1_0_m1 = t243+t251+t303;
      coeffs_dy->coeff_m1_p1_m1 = t298+t322+t327;
      coeffs_dy->coeff_0_p1_m1 = t295+t316+t329;
      coeffs_dy->coeff_p1_p1_m1 = t297+t320+t330;
      coeffs_dy->coeff_m1_m1_0 = t292+t318+t328;
      coeffs_dy->coeff_0_m1_0 = t287+t308*t291+(t287*t291+t308)*t290+t309;
      coeffs_dy->coeff_p1_m1_0 = t294+t318+t323;
      coeffs_dy->coeff_m1_0_0 = t289+(-t291+t301)*y;
      coeffs_dy->coeff_0_0_0 = ((t304+1.0)*RATIONAL(-2.0,1.0)+(t291+t290)*
RATIONAL(2.0,1.0))*y;
      coeffs_dy->coeff_p1_0_0 = t289+(-t291-t301)*y;
      coeffs_dy->coeff_m1_p1_0 = t294+t315+t328;
      coeffs_dy->coeff_0_p1_0 = t284+t307*t291+(t305+t307)*t290+t309;
      coeffs_dy->coeff_p1_p1_0 = t292+t315+t323;
      coeffs_dy->coeff_m1_m1_p1 = t296+t320+t331;
      coeffs_dy->coeff_0_m1_p1 = t295+t319+t324;
      coeffs_dy->coeff_p1_m1_p1 = t298+t317+t332;
      coeffs_dy->coeff_m1_0_p1 = t244+t250+t303;
      coeffs_dy->coeff_0_0_p1 = t289+(-t290-t300)*y;
      coeffs_dy->coeff_p1_0_p1 = t250+t251+t302;
      coeffs_dy->coeff_m1_p1_p1 = t297+t321+t331;
      coeffs_dy->coeff_0_p1_p1 = t293+t316+t324;
      coeffs_dy->coeff_p1_p1_p1 = t299+t322+t332;