aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dx.compute.c
blob: df1341fd9c8cd2dbbcf09a4a19f03181fdf5f07d (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 t187;
fp t194;
fp t195;
fp t208;
fp t161;
fp t143;
fp t218;
fp t210;
fp t238;
fp t190;
fp t237;
fp t179;
fp t162;
fp t144;
fp t177;
fp t163;
fp t145;
fp t236;
fp t181;
fp t174;
fp t151;
fp t235;
fp t183;
fp t173;
fp t150;
fp t234;
fp t188;
fp t209;
fp t166;
fp t147;
fp t191;
fp t216;
fp t175;
fp t233;
fp t217;
fp t164;
fp t148;
fp t176;
fp t232;
fp t231;
fp t193;
fp t153;
fp t230;
fp t229;
fp t154;
fp t165;
fp t228;
fp t155;
fp t227;
fp t186;
fp t157;
fp t160;
fp t226;
fp t215;
fp t158;
fp t189;
fp t214;
fp t170;
fp t225;
fp t159;
fp t169;
fp t224;
fp t223;
fp t222;
fp t168;
fp t171;
fp t221;
fp t172;
fp t220;
fp t219;
fp t213;
fp t212;
fp t211;
fp t207;
fp t206;
fp t205;
fp t204;
fp t203;
fp t202;
fp t201;
fp t200;
fp t199;
fp t198;
fp t197;
fp t196;
      t187 = RATIONAL(1.0,4.0);
      t194 = z*z;
      t195 = y*y;
      t208 = t195*t194;
      t161 = t187*t208;
      t143 = x*t161;
      t218 = y*z;
      t210 = x*t218;
      t238 = t143+t187*t210;
      t190 = RATIONAL(-1.0,4.0);
      t237 = t143+t190*t210;
      t179 = t187*t195;
      t162 = z*t179;
      t144 = x*t162;
      t177 = t187*t194;
      t163 = y*t177;
      t145 = x*t163;
      t236 = t144+t145;
      t181 = t190*t194;
      t174 = y*t181;
      t151 = x*t174;
      t235 = t144+t151;
      t183 = t190*t195;
      t173 = z*t183;
      t150 = x*t173;
      t234 = t145+t150;
      t188 = RATIONAL(1.0,2.0);
      t209 = t188*t195;
      t166 = x*t209;
      t147 = z*t166;
      t191 = RATIONAL(-1.0,2.0);
      t216 = t191*x;
      t175 = z*t216;
      t233 = t147+t175;
      t217 = t188*x;
      t164 = t194*t217;
      t148 = y*t164;
      t176 = y*t216;
      t232 = t148+t176;
      t231 = t150+t151;
      t193 = x*t208;
      t153 = t191*t193;
      t230 = t153+t164;
      t229 = t153+t166;
      t154 = t195*t175;
      t165 = z*t217;
      t228 = t154+t165;
      t155 = t194*t176;
      t227 = t155+y*t217;
      t186 = RATIONAL(1.0,8.0);
      t157 = t186*t208;
      t160 = t186*t218;
      t226 = t157+t160;
      t215 = t195*z;
      t158 = t186*t215;
      t189 = RATIONAL(-1.0,8.0);
      t214 = y*t194;
      t170 = t189*t214;
      t225 = t158+t170;
      t159 = t186*t214;
      t169 = t189*t215;
      t224 = t159+t169;
      t223 = t161+t181;
      t222 = t161+t183;
      t168 = t189*t208;
      t171 = t189*t218;
      t221 = t168+t171;
      t172 = t190*t208;
      t220 = t172+t177;
      t219 = t172+t179;
      t213 = t191-x;
      t212 = t193+x;
      t211 = -x+t188;
      t207 = y*t165+t153;
      t206 = t153+y*t175;
      t205 = z-t215;
      t204 = t214-y;
      t203 = t158+t159+t238;
      t202 = t169+t170+t238;
      t201 = t157+t171+t237;
      t200 = t160+t168+t237;
      t199 = t162+t190*z+t230;
      t198 = t163+t190*y+t229;
      t197 = t173+t187*z+t230;
      t196 = t174+t187*y+t229;
      coeffs_dx->coeff_m1_m1_m1 = t203+t221+t231;
      coeffs_dx->coeff_0_m1_m1 = t148+t147+t206;
      coeffs_dx->coeff_p1_m1_m1 = t202+t226+t231;
      coeffs_dx->coeff_m1_0_m1 = t197+t223+t233;
      coeffs_dx->coeff_0_0_m1 = t193+(-t194+t205)*x;
      coeffs_dx->coeff_p1_0_m1 = t199+t220+t233;
      coeffs_dx->coeff_m1_p1_m1 = t200+t225+t234;
      coeffs_dx->coeff_0_p1_m1 = t147+t155+t207;
      coeffs_dx->coeff_p1_p1_m1 = t201+t224+t234;
      coeffs_dx->coeff_m1_m1_0 = t196+t222+t232;
      coeffs_dx->coeff_0_m1_0 = t193+(-t195-t204)*x;
      coeffs_dx->coeff_p1_m1_0 = t198+t219+t232;
      coeffs_dx->coeff_m1_0_0 = t191+t211*t195+(t191*t195+t211)*t194+t212;
      coeffs_dx->coeff_0_0_0 = ((t208+1.0)*RATIONAL(-2.0,1.0)+(t195+t194)*
RATIONAL(2.0,1.0))*x;
      coeffs_dx->coeff_p1_0_0 = t188+t213*t195+(t209+t213)*t194+t212;
      coeffs_dx->coeff_m1_p1_0 = t198+t222+t227;
      coeffs_dx->coeff_0_p1_0 = t193+(-t195+t204)*x;
      coeffs_dx->coeff_p1_p1_0 = t196+t219+t227;
      coeffs_dx->coeff_m1_m1_p1 = t200+t224+t235;
      coeffs_dx->coeff_0_m1_p1 = t154+t148+t207;
      coeffs_dx->coeff_p1_m1_p1 = t201+t225+t235;
      coeffs_dx->coeff_m1_0_p1 = t199+t223+t228;
      coeffs_dx->coeff_0_0_p1 = t193+(-t194-t205)*x;
      coeffs_dx->coeff_p1_0_p1 = t197+t220+t228;
      coeffs_dx->coeff_m1_p1_p1 = t202+t221+t236;
      coeffs_dx->coeff_0_p1_p1 = t154+t155+t206;
      coeffs_dx->coeff_p1_p1_p1 = t203+t226+t236;