aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dz.compute.c
blob: 57847913e6e2b7b8b6a39542ee4058296128953a (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 t379;
fp t386;
fp t387;
fp t400;
fp t353;
fp t335;
fp t410;
fp t402;
fp t430;
fp t382;
fp t429;
fp t371;
fp t354;
fp t336;
fp t369;
fp t355;
fp t337;
fp t428;
fp t373;
fp t366;
fp t343;
fp t427;
fp t375;
fp t365;
fp t342;
fp t426;
fp t380;
fp t401;
fp t358;
fp t339;
fp t383;
fp t408;
fp t367;
fp t425;
fp t409;
fp t356;
fp t340;
fp t368;
fp t424;
fp t423;
fp t385;
fp t345;
fp t422;
fp t421;
fp t346;
fp t357;
fp t420;
fp t347;
fp t419;
fp t378;
fp t349;
fp t352;
fp t418;
fp t407;
fp t350;
fp t381;
fp t406;
fp t362;
fp t417;
fp t351;
fp t361;
fp t416;
fp t415;
fp t414;
fp t360;
fp t363;
fp t413;
fp t364;
fp t412;
fp t411;
fp t405;
fp t404;
fp t403;
fp t399;
fp t398;
fp t397;
fp t396;
fp t395;
fp t394;
fp t393;
fp t392;
fp t391;
fp t390;
fp t389;
fp t388;
      t379 = RATIONAL(1.0,4.0);
      t386 = y*y;
      t387 = x*x;
      t400 = t387*t386;
      t353 = t379*t400;
      t335 = z*t353;
      t410 = x*y;
      t402 = z*t410;
      t430 = t335+t379*t402;
      t382 = RATIONAL(-1.0,4.0);
      t429 = t335+t382*t402;
      t371 = t379*t387;
      t354 = y*t371;
      t336 = z*t354;
      t369 = t379*t386;
      t355 = x*t369;
      t337 = z*t355;
      t428 = t336+t337;
      t373 = t382*t386;
      t366 = x*t373;
      t343 = z*t366;
      t427 = t336+t343;
      t375 = t382*t387;
      t365 = y*t375;
      t342 = z*t365;
      t426 = t337+t342;
      t380 = RATIONAL(1.0,2.0);
      t401 = t380*t387;
      t358 = z*t401;
      t339 = y*t358;
      t383 = RATIONAL(-1.0,2.0);
      t408 = t383*z;
      t367 = y*t408;
      t425 = t339+t367;
      t409 = t380*z;
      t356 = t386*t409;
      t340 = x*t356;
      t368 = x*t408;
      t424 = t340+t368;
      t423 = t342+t343;
      t385 = z*t400;
      t345 = t383*t385;
      t422 = t345+t356;
      t421 = t345+t358;
      t346 = t387*t367;
      t357 = y*t409;
      t420 = t346+t357;
      t347 = t386*t368;
      t419 = t347+x*t409;
      t378 = RATIONAL(1.0,8.0);
      t349 = t378*t400;
      t352 = t378*t410;
      t418 = t349+t352;
      t407 = t387*y;
      t350 = t378*t407;
      t381 = RATIONAL(-1.0,8.0);
      t406 = x*t386;
      t362 = t381*t406;
      t417 = t350+t362;
      t351 = t378*t406;
      t361 = t381*t407;
      t416 = t351+t361;
      t415 = t353+t373;
      t414 = t353+t375;
      t360 = t381*t400;
      t363 = t381*t410;
      t413 = t360+t363;
      t364 = t382*t400;
      t412 = t364+t369;
      t411 = t364+t371;
      t405 = t380-z;
      t404 = t383-z;
      t403 = t385+z;
      t399 = x*t357+t345;
      t398 = t345+x*t367;
      t397 = t407-y;
      t396 = t406-x;
      t395 = t350+t351+t430;
      t394 = t361+t362+t430;
      t393 = t349+t363+t429;
      t392 = t352+t360+t429;
      t391 = t354+t382*y+t422;
      t390 = t355+t382*x+t421;
      t389 = t365+t379*y+t422;
      t388 = t366+t379*x+t421;
      coeffs_dz->coeff_m1_m1_m1 = t395+t413+t423;
      coeffs_dz->coeff_0_m1_m1 = t389+t415+t425;
      coeffs_dz->coeff_p1_m1_m1 = t392+t417+t426;
      coeffs_dz->coeff_m1_0_m1 = t388+t414+t424;
      coeffs_dz->coeff_0_0_m1 = t383+t405*t387+(t383*t387+t405)*t386+t403;
      coeffs_dz->coeff_p1_0_m1 = t390+t414+t419;
      coeffs_dz->coeff_m1_p1_m1 = t392+t416+t427;
      coeffs_dz->coeff_0_p1_m1 = t391+t415+t420;
      coeffs_dz->coeff_p1_p1_m1 = t394+t413+t428;
      coeffs_dz->coeff_m1_m1_0 = t339+t340+t398;
      coeffs_dz->coeff_0_m1_0 = t385+(-t386-t397)*z;
      coeffs_dz->coeff_p1_m1_0 = t339+t347+t399;
      coeffs_dz->coeff_m1_0_0 = t385+(-t387-t396)*z;
      coeffs_dz->coeff_0_0_0 = ((1.0+t400)*RATIONAL(-2.0,1.0)+(t386+t387)*
RATIONAL(2.0,1.0))*z;
      coeffs_dz->coeff_p1_0_0 = t385+(-t387+t396)*z;
      coeffs_dz->coeff_m1_p1_0 = t340+t346+t399;
      coeffs_dz->coeff_0_p1_0 = t385+(-t386+t397)*z;
      coeffs_dz->coeff_p1_p1_0 = t347+t346+t398;
      coeffs_dz->coeff_m1_m1_p1 = t394+t418+t423;
      coeffs_dz->coeff_0_m1_p1 = t391+t412+t425;
      coeffs_dz->coeff_p1_m1_p1 = t393+t416+t426;
      coeffs_dz->coeff_m1_0_p1 = t390+t411+t424;
      coeffs_dz->coeff_0_0_p1 = t380+t404*t387+(t401+t404)*t386+t403;
      coeffs_dz->coeff_p1_0_p1 = t388+t411+t419;
      coeffs_dz->coeff_m1_p1_p1 = t393+t417+t427;
      coeffs_dz->coeff_0_p1_p1 = t389+t412+t420;
      coeffs_dz->coeff_p1_p1_p1 = t395+t418+t428;