aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-I.compute.c
blob: fc1b8a5f192d2beaf7817e8f89da6b858d70967a (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
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
fp t81;
fp t84;
fp t119;
fp t85;
fp t115;
fp t78;
fp t56;
fp t29;
fp t142;
fp t118;
fp t65;
fp t141;
fp t83;
fp t116;
fp t54;
fp t30;
fp t120;
fp t63;
fp t140;
fp t139;
fp t117;
fp t50;
fp t32;
fp t138;
fp t69;
fp t137;
fp t103;
fp t43;
fp t136;
fp t135;
fp t134;
fp t122;
fp t106;
fp t44;
fp t133;
fp t121;
fp t55;
fp t132;
fp t105;
fp t45;
fp t123;
fp t52;
fp t131;
fp t130;
fp t104;
fp t47;
fp t129;
fp t60;
fp t128;
fp t127;
fp t126;
fp t125;
fp t124;
fp t114;
fp t79;
fp t35;
fp t82;
fp t113;
fp t112;
fp t111;
fp t110;
fp t109;
fp t108;
fp t107;
fp t102;
fp t101;
fp t100;
fp t77;
fp t21;
fp t22;
fp t23;
fp t24;
fp t99;
fp t25;
fp t80;
fp t37;
fp t41;
fp t98;
fp t27;
fp t36;
fp t42;
fp t97;
fp t28;
fp t38;
fp t39;
fp t96;
fp t95;
fp t94;
fp t93;
fp t72;
fp t92;
fp t91;
fp t90;
fp t89;
fp t71;
fp t88;
fp t87;
fp t86;
fp t40;
fp t26;
      t81 = RATIONAL(-1.0,4.0);
      t84 = y*y;
      t119 = t81*t84;
      t85 = x*x;
      t115 = t85*t84;
      t78 = RATIONAL(1.0,4.0);
      t56 = t78*t115;
      t29 = z*t56;
      t142 = t29+z*t119;
      t118 = t81*t85;
      t65 = z*t118;
      t141 = t29+t65;
      t83 = z*z;
      t116 = t85*t83;
      t54 = t78*t116;
      t30 = y*t54;
      t120 = t81*t83;
      t63 = y*t120;
      t140 = t30+t63;
      t139 = t30+y*t118;
      t117 = t84*t83;
      t50 = t78*t117;
      t32 = x*t50;
      t138 = t32+x*t120;
      t69 = x*t119;
      t137 = t32+t69;
      t103 = t83*t115;
      t43 = t81*t103;
      t136 = t43+t50;
      t135 = t43+t54;
      t134 = t43+t56;
      t122 = t78*t84;
      t106 = z*t115;
      t44 = t81*t106;
      t133 = t44+z*t122;
      t121 = t78*t85;
      t55 = z*t121;
      t132 = t44+t55;
      t105 = y*t116;
      t45 = t81*t105;
      t123 = t78*t83;
      t52 = y*t123;
      t131 = t45+t52;
      t130 = t45+y*t121;
      t104 = x*t117;
      t47 = t81*t104;
      t129 = t47+x*t123;
      t60 = x*t122;
      t128 = t47+t60;
      t127 = x*y;
      t126 = y*z;
      t125 = t78*z;
      t124 = t81*z;
      t114 = x*t126;
      t79 = RATIONAL(1.0,2.0);
      t35 = t79*t103;
      t82 = RATIONAL(-1.0,2.0);
      t113 = t35+t82*t115;
      t112 = t83+t84;
      t111 = t83+t85;
      t110 = t84+t85;
      t109 = t83*t127;
      t108 = t84*x*z;
      t107 = t85*t126;
      t102 = 1.0+t116;
      t101 = 1.0+t115;
      t100 = t117+1.0;
      t77 = RATIONAL(1.0,8.0);
      t21 = t77*t103;
      t22 = t77*t106;
      t23 = t77*t105;
      t24 = t77*t107;
      t99 = t21+t22+t23+t24;
      t25 = t77*t104;
      t80 = RATIONAL(-1.0,8.0);
      t37 = t80*t105;
      t41 = t80*t109;
      t98 = t21+t25+t37+t41;
      t27 = t77*t109;
      t36 = t80*t106;
      t42 = t80*t114;
      t97 = t21+t27+t36+t42;
      t28 = t77*t114;
      t38 = t80*t107;
      t39 = t80*t104;
      t96 = t21+t28+t38+t39;
      t95 = t136+(t55+t124)*y;
      t94 = z*t60+x*t124+t135;
      t93 = x*t52+t81*t127+t134;
      t72 = t82*t116;
      t92 = t79*t85+t72+t113;
      t91 = t136+(t65+t125)*y;
      t90 = z*t69+x*t125+t135;
      t89 = x*t63+t78*t127+t134;
      t71 = t82*t117;
      t88 = t72+t71+t35+t79*t83;
      t87 = t79*t84+t71+t113;
      t86 = -t115+t110;
      t40 = t80*t108;
      t26 = t77*t108;
      coeffs_I->coeff_m1_m1_m1 = t26+t39+t24+t37+t97;
      coeffs_I->coeff_0_m1_m1 = t91+t140+t142;
      coeffs_I->coeff_p1_m1_m1 = t28+t24+t40+t36+t98;
      coeffs_I->coeff_m1_0_m1 = t90+t138+t141;
      coeffs_I->coeff_0_0_m1 = (t101*t82+t110*t79)*z+t88;
      coeffs_I->coeff_p1_0_m1 = t94+t129+t141;
      coeffs_I->coeff_m1_p1_m1 = t23+t41+t26+t36+t96;
      coeffs_I->coeff_0_p1_m1 = t95+t131+t142;
      coeffs_I->coeff_p1_p1_m1 = t38+t40+t25+t23+t97;
      coeffs_I->coeff_m1_m1_0 = t89+t137+t139;
      coeffs_I->coeff_0_m1_0 = (t102*t82+t111*t79)*y+t87;
      coeffs_I->coeff_p1_m1_0 = t93+t128+t139;
      coeffs_I->coeff_m1_0_0 = (t100*t82+t112*t79)*x+t92;
      coeffs_I->coeff_0_0_0 = RATIONAL(1.0,1.0)+(-1.0+t86)*t83-t86;
      coeffs_I->coeff_p1_0_0 = (t100*t79+t112*t82)*x+t92;
      coeffs_I->coeff_m1_p1_0 = t93+t130+t137;
      coeffs_I->coeff_0_p1_0 = (t102*t79+t111*t82)*y+t87;
      coeffs_I->coeff_p1_p1_0 = t89+t128+t130;
      coeffs_I->coeff_m1_m1_p1 = t37+t27+t22+t40+t96;
      coeffs_I->coeff_0_m1_p1 = t95+t133+t140;
      coeffs_I->coeff_p1_m1_p1 = t42+t22+t26+t38+t98;
      coeffs_I->coeff_m1_0_p1 = t94+t132+t138;
      coeffs_I->coeff_0_0_p1 = (t101*t79+t110*t82)*z+t88;
      coeffs_I->coeff_p1_0_p1 = t90+t129+t132;
      coeffs_I->coeff_m1_p1_p1 = t41+t42+t40+t39+t99;
      coeffs_I->coeff_0_p1_p1 = t91+t131+t133;
      coeffs_I->coeff_p1_p1_p1 = t25+t26+t27+t28+t99;