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