1 //
2 // BAGEL - Brilliantly Advanced General Electronic Structure Library
3 // Filename: _carsph_44.cc
4 // Copyright (C) 2009 Toru Shiozaki
5 //
6 // Author: Toru Shiozaki <shiozaki@northwestern.edu>
7 // Maintainer: Shiozaki group
8 //
9 // This file is part of the BAGEL package.
10 //
11 // This program is free software: you can redistribute it and/or modify
12 // it under the terms of the GNU General Public License as published by
13 // the Free Software Foundation, either version 3 of the License, or
14 // (at your option) any later version.
15 //
16 // This program is distributed in the hope that it will be useful,
17 // but WITHOUT ANY WARRANTY; without even the implied warranty of
18 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
19 // GNU General Public License for more details.
20 //
21 // You should have received a copy of the GNU General Public License
22 // along with this program. If not, see <http://www.gnu.org/licenses/>.
23 //
24
25 #include <src/integral/carsphlist.h>
26 #include <algorithm>
27
28 using namespace std;
29 using namespace bagel;
30
31
carsph_44(const int nloop,const double * source,double * target)32 void CarSphList::carsph_44(const int nloop, const double* source, double* target) {
33 const double c71 = 45;
34 const double c44 = 42.093645601206838;
35 const double c40 = 39.375;
36 const double c14 = 29.764702249476645;
37 const double c8 = 27.84232950922031;
38 const double c56 = 22.5;
39 const double c74 = 21.213203435596427;
40 const double c42 = 21.046822800603419;
41 const double c79 = 20.124611797498108;
42 const double c32 = 19.843134832984429;
43 const double c2 = 19.6875;
44 const double c50 = 18.824850597016699;
45 const double c29 = 18.561553006146873;
46 const double c75 = 15.90990257669732;
47 const double c11 = 14.882351124738323;
48 const double c17 = 14.031215200402281;
49 const double c24 = 13.311179511974137;
50 const double c4 = 13.125;
51 const double c53 = 11.25;
52 const double c59 = 10.606601717798213;
53 const double c18 = 10.523411400301709;
54 const double c64 = 10.062305898749054;
55 const double c83 = 10;
56 const double c30 = 9.9215674164922145;
57 const double c86 = 9.4868329805051381;
58 const double c34 = 9.354143466934854;
59 const double c7 = 9.2807765030734366;
60 const double c95 = 9;
61 const double c37 = 8.8741196746494246;
62 const double c27 = 8.75;
63 const double c60 = 7.9549512883486599;
64 const double c72 = 7.5;
65 const double c89 = 7.1151247353788536;
66 const double c35 = 7.0156076002011405;
67 const double c78 = 6.7082039324993694;
68 const double c45 = 6.6143782776614763;
69 const double c47 = 6.2749501990055663;
70 const double c28 = 6.1871843353822911;
71 const double c84 = 5.625;
72 const double c80 = 5.0311529493745271;
73 const double c12 = 4.9607837082461073;
74 const double c52 = 4.7062126492541747;
75 const double c6 = 4.6403882515367183;
76 const double c23 = 4.4370598373247123;
77 const double c39 = 4.375;
78 const double c57 = 3.75;
79 const double c76 = 3.5355339059327378;
80 const double c43 = 3.5078038001005702;
81 const double c63 = 3.3541019662496847;
82 const double c26 = 3.3277948779935342;
83 const double c33 = 3.3071891388307382;
84 const double c1 = 3.28125;
85 const double c85 = 3.1622776601683795;
86 const double c92 = 3;
87 const double c36 = 2.9580398915498081;
88 const double c77 = 2.6516504294495533;
89 const double c66 = 2.5155764746872635;
90 const double c9 = 2.4803918541230536;
91 const double c88 = 2.3717082451262845;
92 const double c51 = 2.3531063246270874;
93 const double c15 = 2.3385358667337135;
94 const double c97 = 2.25;
95 const double c20 = 2.2185299186623562;
96 const double c3 = 2.1875;
97 const double c46 = 2.0916500663351889;
98 const double c54 = 1.875;
99 const double c91 = 1.7787811838447134;
100 const double c61 = 1.7677669529663689;
101 const double c16 = 1.7539019000502851;
102 const double c68 = 1.6770509831248424;
103 const double c25 = 1.6638974389967671;
104 const double c31 = 1.6535945694153691;
105 const double c49 = 1.5687375497513916;
106 const double c5 = 1.5467960838455728;
107 const double c62 = 1.3258252147247767;
108 const double c65 = 1.2577882373436318;
109 const double c73 = 1.25;
110 const double c87 = 1.1858541225631423;
111 const double c41 = 1.1692679333668567;
112 const double c96 = 1.125;
113 const double c81 = 1.1180339887498949;
114 const double c38 = 1.1092649593311781;
115 const double c90 = 0.8893905919223567;
116 const double c82 = 0.83852549156242118;
117 const double c13 = 0.82679728470768454;
118 const double c48 = 0.78436877487569578;
119 const double c94 = 0.75;
120 const double c19 = 0.73950997288745202;
121 const double c58 = 0.625;
122 const double c100 = 0.5625;
123 const double c67 = 0.55901699437494745;
124 const double c22 = 0.55463247966558904;
125 const double c0 = 0.546875;
126 const double c70 = 0.41926274578121059;
127 const double c10 = 0.41339864235384227;
128 const double c93 = 0.375;
129 const double c55 = 0.3125;
130 const double c99 = 0.28125;
131 const double c21 = 0.27731623983279452;
132 const double c69 = 0.20963137289060529;
133 const double c98 = 0.140625;
134 for (int iloop = 0; iloop != nloop; ++iloop, target += 81, source += 225) {
135 target[0] = c0 * source[0] - c1 * source[2] + c0 * source[4]
136 - c1 * source[30] + c2 * source[32] - c1 * source[34]
137 + c0 * source[60] - c1 * source[62] + c0 * source[64];
138 target[1] = c3 * source[1] - c3 * source[3] - c4 * source[31]
139 + c4 * source[33] + c3 * source[61] - c3 * source[63];
140 target[2] = c5 * source[5] - c6 * source[7] - c7 * source[35]
141 + c8 * source[37] + c5 * source[65] - c6 * source[67];
142 target[3] = c6 * source[6] - c5 * source[8] - c8 * source[36]
143 + c7 * source[38] + c6 * source[66] - c5 * source[68];
144 target[4] = c9 * source[9] - c9 * source[11] - c10 * source[0]
145 + c10 * source[2] - c10 * source[2] + c10 * source[4]
146 - c11 * source[39] + c11 * source[41] + c9 * source[30]
147 - c9 * source[32] + c9 * source[32] - c9 * source[34]
148 + c9 * source[69] - c9 * source[71] - c10 * source[60]
149 + c10 * source[62] - c10 * source[62] + c10 * source[64];
150 target[5] = c12 * source[10] - c13 * source[1] - c13 * source[3]
151 - c14 * source[40] + c12 * source[31] + c12 * source[33]
152 + c12 * source[70] - c13 * source[61] - c13 * source[63];
153 target[6] = c15 * source[12] - c16 * source[5] - c16 * source[7]
154 - c17 * source[42] + c18 * source[35] + c18 * source[37]
155 + c15 * source[72] - c16 * source[65] - c16 * source[67];
156 target[7] = c15 * source[13] - c16 * source[6] - c16 * source[8]
157 - c17 * source[43] + c18 * source[36] + c18 * source[38]
158 + c15 * source[73] - c16 * source[66] - c16 * source[68];
159 target[8] = c19 * source[14] - c20 * source[9] - c20 * source[11]
160 + c21 * source[0] + c22 * source[2] + c21 * source[4]
161 - c23 * source[44] + c24 * source[39] + c24 * source[41]
162 - c25 * source[30] - c26 * source[32] - c25 * source[34]
163 + c19 * source[74] - c20 * source[69] - c20 * source[71]
164 + c21 * source[60] + c22 * source[62] + c21 * source[64];
165 target[9] = c3 * source[15] - c4 * source[17] + c3 * source[19]
166 - c3 * source[45] + c4 * source[47] - c3 * source[49];
167 target[10] = c27 * source[16] - c27 * source[18] - c27 * source[46]
168 + c27 * source[48];
169 target[11] = c28 * source[20] - c29 * source[22] - c28 * source[50]
170 + c29 * source[52];
171 target[12] = c29 * source[21] - c28 * source[23] - c29 * source[51]
172 + c28 * source[53];
173 target[13] = c30 * source[24] - c30 * source[26] - c31 * source[15]
174 + c31 * source[17] - c31 * source[17] + c31 * source[19]
175 - c30 * source[54] + c30 * source[56] + c31 * source[45]
176 - c31 * source[47] + c31 * source[47] - c31 * source[49];
177 target[14] = c32 * source[25] - c33 * source[16] - c33 * source[18]
178 - c32 * source[55] + c33 * source[46] + c33 * source[48];
179 target[15] = c34 * source[27] - c35 * source[20] - c35 * source[22]
180 - c34 * source[57] + c35 * source[50] + c35 * source[52];
181 target[16] = c34 * source[28] - c35 * source[21] - c35 * source[23]
182 - c34 * source[58] + c35 * source[51] + c35 * source[53];
183 target[17] = c36 * source[29] - c37 * source[24] - c37 * source[26]
184 + c38 * source[15] + c20 * source[17] + c38 * source[19]
185 - c36 * source[59] + c37 * source[54] + c37 * source[56]
186 - c38 * source[45] - c20 * source[47] - c38 * source[49];
187 target[18] = c5 * source[75] - c7 * source[77] + c5 * source[79]
188 - c6 * source[105] + c8 * source[107] - c6 * source[109];
189 target[19] = c28 * source[76] - c28 * source[78] - c29 * source[106]
190 + c29 * source[108];
191 target[20] = c39 * source[80] - c4 * source[82] - c4 * source[110]
192 + c40 * source[112];
193 target[21] = c4 * source[81] - c39 * source[83] - c40 * source[111]
194 + c4 * source[113];
195 target[22] = c35 * source[84] - c35 * source[86] - c41 * source[75]
196 + c41 * source[77] - c41 * source[77] + c41 * source[79]
197 - c42 * source[114] + c42 * source[116] + c43 * source[105]
198 - c43 * source[107] + c43 * source[107] - c43 * source[109];
199 target[23] = c17 * source[85] - c15 * source[76] - c15 * source[78]
200 - c44 * source[115] + c35 * source[106] + c35 * source[108];
201 target[24] = c45 * source[87] - c12 * source[80] - c12 * source[82]
202 - c32 * source[117] + c11 * source[110] + c11 * source[112];
203 target[25] = c45 * source[88] - c12 * source[81] - c12 * source[83]
204 - c32 * source[118] + c11 * source[111] + c11 * source[113];
205 target[26] = c46 * source[89] - c47 * source[84] - c47 * source[86]
206 + c48 * source[75] + c49 * source[77] + c48 * source[79]
207 - c47 * source[119] + c50 * source[114] + c50 * source[116]
208 - c51 * source[105] - c52 * source[107] - c51 * source[109];
209 target[27] = c6 * source[90] - c8 * source[92] + c6 * source[94]
210 - c5 * source[120] + c7 * source[122] - c5 * source[124];
211 target[28] = c29 * source[91] - c29 * source[93] - c28 * source[121]
212 + c28 * source[123];
213 target[29] = c4 * source[95] - c40 * source[97] - c39 * source[125]
214 + c4 * source[127];
215 target[30] = c40 * source[96] - c4 * source[98] - c4 * source[126]
216 + c39 * source[128];
217 target[31] = c42 * source[99] - c42 * source[101] - c43 * source[90]
218 + c43 * source[92] - c43 * source[92] + c43 * source[94]
219 - c35 * source[129] + c35 * source[131] + c41 * source[120]
220 - c41 * source[122] + c41 * source[122] - c41 * source[124];
221 target[32] = c44 * source[100] - c35 * source[91] - c35 * source[93]
222 - c17 * source[130] + c15 * source[121] + c15 * source[123];
223 target[33] = c32 * source[102] - c11 * source[95] - c11 * source[97]
224 - c45 * source[132] + c12 * source[125] + c12 * source[127];
225 target[34] = c32 * source[103] - c11 * source[96] - c11 * source[98]
226 - c45 * source[133] + c12 * source[126] + c12 * source[128];
227 target[35] = c47 * source[104] - c50 * source[99] - c50 * source[101]
228 + c51 * source[90] + c52 * source[92] + c51 * source[94]
229 - c46 * source[134] + c47 * source[129] + c47 * source[131]
230 - c48 * source[120] - c49 * source[122] - c48 * source[124];
231 target[36] = c9 * source[135] - c11 * source[137] + c9 * source[139]
232 - c9 * source[165] + c11 * source[167] - c9 * source[169]
233 - c10 * source[0] + c9 * source[2] - c10 * source[4]
234 + c10 * source[30] - c9 * source[32] + c10 * source[34]
235 - c10 * source[30] + c9 * source[32] - c10 * source[34]
236 + c10 * source[60] - c9 * source[62] + c10 * source[64];
237 target[37] = c30 * source[136] - c30 * source[138] - c30 * source[166]
238 + c30 * source[168] - c31 * source[1] + c31 * source[3]
239 + c31 * source[31] - c31 * source[33] - c31 * source[31]
240 + c31 * source[33] + c31 * source[61] - c31 * source[63];
241 target[38] = c35 * source[140] - c42 * source[142] - c35 * source[170]
242 + c42 * source[172] - c41 * source[5] + c43 * source[7]
243 + c41 * source[35] - c43 * source[37] - c41 * source[35]
244 + c43 * source[37] + c41 * source[65] - c43 * source[67];
245 target[39] = c42 * source[141] - c35 * source[143] - c42 * source[171]
246 + c35 * source[173] - c43 * source[6] + c41 * source[8]
247 + c43 * source[36] - c41 * source[38] - c43 * source[36]
248 + c41 * source[38] + c43 * source[66] - c41 * source[68];
249 target[40] = c53 * source[144] - c53 * source[146] - c54 * source[135]
250 + c54 * source[137] - c54 * source[137] + c54 * source[139]
251 - c53 * source[174] + c53 * source[176] + c54 * source[165]
252 - c54 * source[167] + c54 * source[167] - c54 * source[169]
253 - c54 * source[9] + c54 * source[11] + c55 * source[0]
254 - c55 * source[2] + c55 * source[2] - c55 * source[4]
255 + c54 * source[39] - c54 * source[41] - c55 * source[30]
256 + c55 * source[32] - c55 * source[32] + c55 * source[34]
257 - c54 * source[39] + c54 * source[41] + c55 * source[30]
258 - c55 * source[32] + c55 * source[32] - c55 * source[34]
259 + c54 * source[69] - c54 * source[71] - c55 * source[60]
260 + c55 * source[62] - c55 * source[62] + c55 * source[64];
261 target[41] = c56 * source[145] - c57 * source[136] - c57 * source[138]
262 - c56 * source[175] + c57 * source[166] + c57 * source[168]
263 - c57 * source[10] + c58 * source[1] + c58 * source[3]
264 + c57 * source[40] - c58 * source[31] - c58 * source[33]
265 - c57 * source[40] + c58 * source[31] + c58 * source[33]
266 + c57 * source[70] - c58 * source[61] - c58 * source[63];
267 target[42] = c59 * source[147] - c60 * source[140] - c60 * source[142]
268 - c59 * source[177] + c60 * source[170] + c60 * source[172]
269 - c61 * source[12] + c62 * source[5] + c62 * source[7]
270 + c61 * source[42] - c62 * source[35] - c62 * source[37]
271 - c61 * source[42] + c62 * source[35] + c62 * source[37]
272 + c61 * source[72] - c62 * source[65] - c62 * source[67];
273 target[43] = c59 * source[148] - c60 * source[141] - c60 * source[143]
274 - c59 * source[178] + c60 * source[171] + c60 * source[173]
275 - c61 * source[13] + c62 * source[6] + c62 * source[8]
276 + c61 * source[43] - c62 * source[36] - c62 * source[38]
277 - c61 * source[43] + c62 * source[36] + c62 * source[38]
278 + c61 * source[73] - c62 * source[66] - c62 * source[68];
279 target[44] = c63 * source[149] - c64 * source[144] - c64 * source[146]
280 + c65 * source[135] + c66 * source[137] + c65 * source[139]
281 - c63 * source[179] + c64 * source[174] + c64 * source[176]
282 - c65 * source[165] - c66 * source[167] - c65 * source[169]
283 - c67 * source[14] + c68 * source[9] + c68 * source[11]
284 - c69 * source[0] - c70 * source[2] - c69 * source[4]
285 + c67 * source[44] - c68 * source[39] - c68 * source[41]
286 + c69 * source[30] + c70 * source[32] + c69 * source[34]
287 - c67 * source[44] + c68 * source[39] + c68 * source[41]
288 - c69 * source[30] - c70 * source[32] - c69 * source[34]
289 + c67 * source[74] - c68 * source[69] - c68 * source[71]
290 + c69 * source[60] + c70 * source[62] + c69 * source[64];
291 target[45] = c12 * source[150] - c14 * source[152] + c12 * source[154]
292 - c13 * source[15] + c12 * source[17] - c13 * source[19]
293 - c13 * source[45] + c12 * source[47] - c13 * source[49];
294 target[46] = c32 * source[151] - c32 * source[153] - c33 * source[16]
295 + c33 * source[18] - c33 * source[46] + c33 * source[48];
296 target[47] = c17 * source[155] - c44 * source[157] - c15 * source[20]
297 + c35 * source[22] - c15 * source[50] + c35 * source[52];
298 target[48] = c44 * source[156] - c17 * source[158] - c35 * source[21]
299 + c15 * source[23] - c35 * source[51] + c15 * source[53];
300 target[49] = c56 * source[159] - c56 * source[161] - c57 * source[150]
301 + c57 * source[152] - c57 * source[152] + c57 * source[154]
302 - c57 * source[24] + c57 * source[26] + c58 * source[15]
303 - c58 * source[17] + c58 * source[17] - c58 * source[19]
304 - c57 * source[54] + c57 * source[56] + c58 * source[45]
305 - c58 * source[47] + c58 * source[47] - c58 * source[49];
306 target[50] = c71 * source[160] - c72 * source[151] - c72 * source[153]
307 - c72 * source[25] + c73 * source[16] + c73 * source[18]
308 - c72 * source[55] + c73 * source[46] + c73 * source[48];
309 target[51] = c74 * source[162] - c75 * source[155] - c75 * source[157]
310 - c76 * source[27] + c77 * source[20] + c77 * source[22]
311 - c76 * source[57] + c77 * source[50] + c77 * source[52];
312 target[52] = c74 * source[163] - c75 * source[156] - c75 * source[158]
313 - c76 * source[28] + c77 * source[21] + c77 * source[23]
314 - c76 * source[58] + c77 * source[51] + c77 * source[53];
315 target[53] = c78 * source[164] - c79 * source[159] - c79 * source[161]
316 + c66 * source[150] + c80 * source[152] + c66 * source[154]
317 - c81 * source[29] + c63 * source[24] + c63 * source[26]
318 - c70 * source[15] - c82 * source[17] - c70 * source[19]
319 - c81 * source[59] + c63 * source[54] + c63 * source[56]
320 - c70 * source[45] - c82 * source[47] - c70 * source[49];
321 target[54] = c15 * source[180] - c17 * source[182] + c15 * source[184]
322 - c16 * source[75] + c18 * source[77] - c16 * source[79]
323 - c16 * source[105] + c18 * source[107] - c16 * source[109];
324 target[55] = c34 * source[181] - c34 * source[183] - c35 * source[76]
325 + c35 * source[78] - c35 * source[106] + c35 * source[108];
326 target[56] = c45 * source[185] - c32 * source[187] - c12 * source[80]
327 + c11 * source[82] - c12 * source[110] + c11 * source[112];
328 target[57] = c32 * source[186] - c45 * source[188] - c11 * source[81]
329 + c12 * source[83] - c11 * source[111] + c12 * source[113];
330 target[58] = c59 * source[189] - c59 * source[191] - c61 * source[180]
331 + c61 * source[182] - c61 * source[182] + c61 * source[184]
332 - c60 * source[84] + c60 * source[86] + c62 * source[75]
333 - c62 * source[77] + c62 * source[77] - c62 * source[79]
334 - c60 * source[114] + c60 * source[116] + c62 * source[105]
335 - c62 * source[107] + c62 * source[107] - c62 * source[109];
336 target[59] = c74 * source[190] - c76 * source[181] - c76 * source[183]
337 - c75 * source[85] + c77 * source[76] + c77 * source[78]
338 - c75 * source[115] + c77 * source[106] + c77 * source[108];
339 target[60] = c83 * source[192] - c72 * source[185] - c72 * source[187]
340 - c72 * source[87] + c84 * source[80] + c84 * source[82]
341 - c72 * source[117] + c84 * source[110] + c84 * source[112];
342 target[61] = c83 * source[193] - c72 * source[186] - c72 * source[188]
343 - c72 * source[88] + c84 * source[81] + c84 * source[83]
344 - c72 * source[118] + c84 * source[111] + c84 * source[113];
345 target[62] = c85 * source[194] - c86 * source[189] - c86 * source[191]
346 + c87 * source[180] + c88 * source[182] + c87 * source[184]
347 - c88 * source[89] + c89 * source[84] + c89 * source[86]
348 - c90 * source[75] - c91 * source[77] - c90 * source[79]
349 - c88 * source[119] + c89 * source[114] + c89 * source[116]
350 - c90 * source[105] - c91 * source[107] - c90 * source[109];
351 target[63] = c15 * source[195] - c17 * source[197] + c15 * source[199]
352 - c16 * source[90] + c18 * source[92] - c16 * source[94]
353 - c16 * source[120] + c18 * source[122] - c16 * source[124];
354 target[64] = c34 * source[196] - c34 * source[198] - c35 * source[91]
355 + c35 * source[93] - c35 * source[121] + c35 * source[123];
356 target[65] = c45 * source[200] - c32 * source[202] - c12 * source[95]
357 + c11 * source[97] - c12 * source[125] + c11 * source[127];
358 target[66] = c32 * source[201] - c45 * source[203] - c11 * source[96]
359 + c12 * source[98] - c11 * source[126] + c12 * source[128];
360 target[67] = c59 * source[204] - c59 * source[206] - c61 * source[195]
361 + c61 * source[197] - c61 * source[197] + c61 * source[199]
362 - c60 * source[99] + c60 * source[101] + c62 * source[90]
363 - c62 * source[92] + c62 * source[92] - c62 * source[94]
364 - c60 * source[129] + c60 * source[131] + c62 * source[120]
365 - c62 * source[122] + c62 * source[122] - c62 * source[124];
366 target[68] = c74 * source[205] - c76 * source[196] - c76 * source[198]
367 - c75 * source[100] + c77 * source[91] + c77 * source[93]
368 - c75 * source[130] + c77 * source[121] + c77 * source[123];
369 target[69] = c83 * source[207] - c72 * source[200] - c72 * source[202]
370 - c72 * source[102] + c84 * source[95] + c84 * source[97]
371 - c72 * source[132] + c84 * source[125] + c84 * source[127];
372 target[70] = c83 * source[208] - c72 * source[201] - c72 * source[203]
373 - c72 * source[103] + c84 * source[96] + c84 * source[98]
374 - c72 * source[133] + c84 * source[126] + c84 * source[128];
375 target[71] = c85 * source[209] - c86 * source[204] - c86 * source[206]
376 + c87 * source[195] + c88 * source[197] + c87 * source[199]
377 - c88 * source[104] + c89 * source[99] + c89 * source[101]
378 - c90 * source[90] - c91 * source[92] - c90 * source[94]
379 - c88 * source[134] + c89 * source[129] + c89 * source[131]
380 - c90 * source[120] - c91 * source[122] - c90 * source[124];
381 target[72] = c19 * source[210] - c23 * source[212] + c19 * source[214]
382 - c20 * source[135] + c24 * source[137] - c20 * source[139]
383 - c20 * source[165] + c24 * source[167] - c20 * source[169]
384 + c21 * source[0] - c25 * source[2] + c21 * source[4]
385 + c22 * source[30] - c26 * source[32] + c22 * source[34]
386 + c21 * source[60] - c25 * source[62] + c21 * source[64];
387 target[73] = c36 * source[211] - c36 * source[213] - c37 * source[136]
388 + c37 * source[138] - c37 * source[166] + c37 * source[168]
389 + c38 * source[1] - c38 * source[3] + c20 * source[31]
390 - c20 * source[33] + c38 * source[61] - c38 * source[63];
391 target[74] = c46 * source[215] - c47 * source[217] - c47 * source[140]
392 + c50 * source[142] - c47 * source[170] + c50 * source[172]
393 + c48 * source[5] - c51 * source[7] + c49 * source[35]
394 - c52 * source[37] + c48 * source[65] - c51 * source[67];
395 target[75] = c47 * source[216] - c46 * source[218] - c50 * source[141]
396 + c47 * source[143] - c50 * source[171] + c47 * source[173]
397 + c51 * source[6] - c48 * source[8] + c52 * source[36]
398 - c49 * source[38] + c51 * source[66] - c48 * source[68];
399 target[76] = c63 * source[219] - c63 * source[221] - c67 * source[210]
400 + c67 * source[212] - c67 * source[212] + c67 * source[214]
401 - c64 * source[144] + c64 * source[146] + c68 * source[135]
402 - c68 * source[137] + c68 * source[137] - c68 * source[139]
403 - c64 * source[174] + c64 * source[176] + c68 * source[165]
404 - c68 * source[167] + c68 * source[167] - c68 * source[169]
405 + c65 * source[9] - c65 * source[11] - c69 * source[0]
406 + c69 * source[2] - c69 * source[2] + c69 * source[4]
407 + c66 * source[39] - c66 * source[41] - c70 * source[30]
408 + c70 * source[32] - c70 * source[32] + c70 * source[34]
409 + c65 * source[69] - c65 * source[71] - c69 * source[60]
410 + c69 * source[62] - c69 * source[62] + c69 * source[64];
411 target[77] = c78 * source[220] - c81 * source[211] - c81 * source[213]
412 - c79 * source[145] + c63 * source[136] + c63 * source[138]
413 - c79 * source[175] + c63 * source[166] + c63 * source[168]
414 + c66 * source[10] - c70 * source[1] - c70 * source[3]
415 + c80 * source[40] - c82 * source[31] - c82 * source[33]
416 + c66 * source[70] - c70 * source[61] - c70 * source[63];
417 target[78] = c85 * source[222] - c88 * source[215] - c88 * source[217]
418 - c86 * source[147] + c89 * source[140] + c89 * source[142]
419 - c86 * source[177] + c89 * source[170] + c89 * source[172]
420 + c87 * source[12] - c90 * source[5] - c90 * source[7]
421 + c88 * source[42] - c91 * source[35] - c91 * source[37]
422 + c87 * source[72] - c90 * source[65] - c90 * source[67];
423 target[79] = c85 * source[223] - c88 * source[216] - c88 * source[218]
424 - c86 * source[148] + c89 * source[141] + c89 * source[143]
425 - c86 * source[178] + c89 * source[171] + c89 * source[173]
426 + c87 * source[13] - c90 * source[6] - c90 * source[8]
427 + c88 * source[43] - c91 * source[36] - c91 * source[38]
428 + c87 * source[73] - c90 * source[66] - c90 * source[68];
429 target[80] = source[224] - c92 * source[219] - c92 * source[221]
430 + c93 * source[210] + c94 * source[212] + c93 * source[214]
431 - c92 * source[149] + c95 * source[144] + c95 * source[146]
432 - c96 * source[135] - c97 * source[137] - c96 * source[139]
433 - c92 * source[179] + c95 * source[174] + c95 * source[176]
434 - c96 * source[165] - c97 * source[167] - c96 * source[169]
435 + c93 * source[14] - c96 * source[9] - c96 * source[11]
436 + c98 * source[0] + c99 * source[2] + c98 * source[4]
437 + c94 * source[44] - c97 * source[39] - c97 * source[41]
438 + c99 * source[30] + c100 * source[32] + c99 * source[34]
439 + c93 * source[74] - c96 * source[69] - c96 * source[71]
440 + c98 * source[60] + c99 * source[62] + c98 * source[64];
441 }
442 }
443
carsph_44(const int nloop,const complex<double> * source,complex<double> * target)444 void CCarSphList::carsph_44(const int nloop, const complex<double>* source, complex<double>* target) {
445 const double c71 = 45;
446 const double c44 = 42.093645601206838;
447 const double c40 = 39.375;
448 const double c14 = 29.764702249476645;
449 const double c8 = 27.84232950922031;
450 const double c56 = 22.5;
451 const double c74 = 21.213203435596427;
452 const double c42 = 21.046822800603419;
453 const double c79 = 20.124611797498108;
454 const double c32 = 19.843134832984429;
455 const double c2 = 19.6875;
456 const double c50 = 18.824850597016699;
457 const double c29 = 18.561553006146873;
458 const double c75 = 15.90990257669732;
459 const double c11 = 14.882351124738323;
460 const double c17 = 14.031215200402281;
461 const double c24 = 13.311179511974137;
462 const double c4 = 13.125;
463 const double c53 = 11.25;
464 const double c59 = 10.606601717798213;
465 const double c18 = 10.523411400301709;
466 const double c64 = 10.062305898749054;
467 const double c83 = 10;
468 const double c30 = 9.9215674164922145;
469 const double c86 = 9.4868329805051381;
470 const double c34 = 9.354143466934854;
471 const double c7 = 9.2807765030734366;
472 const double c95 = 9;
473 const double c37 = 8.8741196746494246;
474 const double c27 = 8.75;
475 const double c60 = 7.9549512883486599;
476 const double c72 = 7.5;
477 const double c89 = 7.1151247353788536;
478 const double c35 = 7.0156076002011405;
479 const double c78 = 6.7082039324993694;
480 const double c45 = 6.6143782776614763;
481 const double c47 = 6.2749501990055663;
482 const double c28 = 6.1871843353822911;
483 const double c84 = 5.625;
484 const double c80 = 5.0311529493745271;
485 const double c12 = 4.9607837082461073;
486 const double c52 = 4.7062126492541747;
487 const double c6 = 4.6403882515367183;
488 const double c23 = 4.4370598373247123;
489 const double c39 = 4.375;
490 const double c57 = 3.75;
491 const double c76 = 3.5355339059327378;
492 const double c43 = 3.5078038001005702;
493 const double c63 = 3.3541019662496847;
494 const double c26 = 3.3277948779935342;
495 const double c33 = 3.3071891388307382;
496 const double c1 = 3.28125;
497 const double c85 = 3.1622776601683795;
498 const double c92 = 3;
499 const double c36 = 2.9580398915498081;
500 const double c77 = 2.6516504294495533;
501 const double c66 = 2.5155764746872635;
502 const double c9 = 2.4803918541230536;
503 const double c88 = 2.3717082451262845;
504 const double c51 = 2.3531063246270874;
505 const double c15 = 2.3385358667337135;
506 const double c97 = 2.25;
507 const double c20 = 2.2185299186623562;
508 const double c3 = 2.1875;
509 const double c46 = 2.0916500663351889;
510 const double c54 = 1.875;
511 const double c91 = 1.7787811838447134;
512 const double c61 = 1.7677669529663689;
513 const double c16 = 1.7539019000502851;
514 const double c68 = 1.6770509831248424;
515 const double c25 = 1.6638974389967671;
516 const double c31 = 1.6535945694153691;
517 const double c49 = 1.5687375497513916;
518 const double c5 = 1.5467960838455728;
519 const double c62 = 1.3258252147247767;
520 const double c65 = 1.2577882373436318;
521 const double c73 = 1.25;
522 const double c87 = 1.1858541225631423;
523 const double c41 = 1.1692679333668567;
524 const double c96 = 1.125;
525 const double c81 = 1.1180339887498949;
526 const double c38 = 1.1092649593311781;
527 const double c90 = 0.8893905919223567;
528 const double c82 = 0.83852549156242118;
529 const double c13 = 0.82679728470768454;
530 const double c48 = 0.78436877487569578;
531 const double c94 = 0.75;
532 const double c19 = 0.73950997288745202;
533 const double c58 = 0.625;
534 const double c100 = 0.5625;
535 const double c67 = 0.55901699437494745;
536 const double c22 = 0.55463247966558904;
537 const double c0 = 0.546875;
538 const double c70 = 0.41926274578121059;
539 const double c10 = 0.41339864235384227;
540 const double c93 = 0.375;
541 const double c55 = 0.3125;
542 const double c99 = 0.28125;
543 const double c21 = 0.27731623983279452;
544 const double c69 = 0.20963137289060529;
545 const double c98 = 0.140625;
546 for (int iloop = 0; iloop != nloop; ++iloop, target += 81, source += 225) {
547 target[0] = c0 * source[0] - c1 * source[2] + c0 * source[4]
548 - c1 * source[30] + c2 * source[32] - c1 * source[34]
549 + c0 * source[60] - c1 * source[62] + c0 * source[64];
550 target[1] = c3 * source[1] - c3 * source[3] - c4 * source[31]
551 + c4 * source[33] + c3 * source[61] - c3 * source[63];
552 target[2] = c5 * source[5] - c6 * source[7] - c7 * source[35]
553 + c8 * source[37] + c5 * source[65] - c6 * source[67];
554 target[3] = c6 * source[6] - c5 * source[8] - c8 * source[36]
555 + c7 * source[38] + c6 * source[66] - c5 * source[68];
556 target[4] = c9 * source[9] - c9 * source[11] - c10 * source[0]
557 + c10 * source[2] - c10 * source[2] + c10 * source[4]
558 - c11 * source[39] + c11 * source[41] + c9 * source[30]
559 - c9 * source[32] + c9 * source[32] - c9 * source[34]
560 + c9 * source[69] - c9 * source[71] - c10 * source[60]
561 + c10 * source[62] - c10 * source[62] + c10 * source[64];
562 target[5] = c12 * source[10] - c13 * source[1] - c13 * source[3]
563 - c14 * source[40] + c12 * source[31] + c12 * source[33]
564 + c12 * source[70] - c13 * source[61] - c13 * source[63];
565 target[6] = c15 * source[12] - c16 * source[5] - c16 * source[7]
566 - c17 * source[42] + c18 * source[35] + c18 * source[37]
567 + c15 * source[72] - c16 * source[65] - c16 * source[67];
568 target[7] = c15 * source[13] - c16 * source[6] - c16 * source[8]
569 - c17 * source[43] + c18 * source[36] + c18 * source[38]
570 + c15 * source[73] - c16 * source[66] - c16 * source[68];
571 target[8] = c19 * source[14] - c20 * source[9] - c20 * source[11]
572 + c21 * source[0] + c22 * source[2] + c21 * source[4]
573 - c23 * source[44] + c24 * source[39] + c24 * source[41]
574 - c25 * source[30] - c26 * source[32] - c25 * source[34]
575 + c19 * source[74] - c20 * source[69] - c20 * source[71]
576 + c21 * source[60] + c22 * source[62] + c21 * source[64];
577 target[9] = c3 * source[15] - c4 * source[17] + c3 * source[19]
578 - c3 * source[45] + c4 * source[47] - c3 * source[49];
579 target[10] = c27 * source[16] - c27 * source[18] - c27 * source[46]
580 + c27 * source[48];
581 target[11] = c28 * source[20] - c29 * source[22] - c28 * source[50]
582 + c29 * source[52];
583 target[12] = c29 * source[21] - c28 * source[23] - c29 * source[51]
584 + c28 * source[53];
585 target[13] = c30 * source[24] - c30 * source[26] - c31 * source[15]
586 + c31 * source[17] - c31 * source[17] + c31 * source[19]
587 - c30 * source[54] + c30 * source[56] + c31 * source[45]
588 - c31 * source[47] + c31 * source[47] - c31 * source[49];
589 target[14] = c32 * source[25] - c33 * source[16] - c33 * source[18]
590 - c32 * source[55] + c33 * source[46] + c33 * source[48];
591 target[15] = c34 * source[27] - c35 * source[20] - c35 * source[22]
592 - c34 * source[57] + c35 * source[50] + c35 * source[52];
593 target[16] = c34 * source[28] - c35 * source[21] - c35 * source[23]
594 - c34 * source[58] + c35 * source[51] + c35 * source[53];
595 target[17] = c36 * source[29] - c37 * source[24] - c37 * source[26]
596 + c38 * source[15] + c20 * source[17] + c38 * source[19]
597 - c36 * source[59] + c37 * source[54] + c37 * source[56]
598 - c38 * source[45] - c20 * source[47] - c38 * source[49];
599 target[18] = c5 * source[75] - c7 * source[77] + c5 * source[79]
600 - c6 * source[105] + c8 * source[107] - c6 * source[109];
601 target[19] = c28 * source[76] - c28 * source[78] - c29 * source[106]
602 + c29 * source[108];
603 target[20] = c39 * source[80] - c4 * source[82] - c4 * source[110]
604 + c40 * source[112];
605 target[21] = c4 * source[81] - c39 * source[83] - c40 * source[111]
606 + c4 * source[113];
607 target[22] = c35 * source[84] - c35 * source[86] - c41 * source[75]
608 + c41 * source[77] - c41 * source[77] + c41 * source[79]
609 - c42 * source[114] + c42 * source[116] + c43 * source[105]
610 - c43 * source[107] + c43 * source[107] - c43 * source[109];
611 target[23] = c17 * source[85] - c15 * source[76] - c15 * source[78]
612 - c44 * source[115] + c35 * source[106] + c35 * source[108];
613 target[24] = c45 * source[87] - c12 * source[80] - c12 * source[82]
614 - c32 * source[117] + c11 * source[110] + c11 * source[112];
615 target[25] = c45 * source[88] - c12 * source[81] - c12 * source[83]
616 - c32 * source[118] + c11 * source[111] + c11 * source[113];
617 target[26] = c46 * source[89] - c47 * source[84] - c47 * source[86]
618 + c48 * source[75] + c49 * source[77] + c48 * source[79]
619 - c47 * source[119] + c50 * source[114] + c50 * source[116]
620 - c51 * source[105] - c52 * source[107] - c51 * source[109];
621 target[27] = c6 * source[90] - c8 * source[92] + c6 * source[94]
622 - c5 * source[120] + c7 * source[122] - c5 * source[124];
623 target[28] = c29 * source[91] - c29 * source[93] - c28 * source[121]
624 + c28 * source[123];
625 target[29] = c4 * source[95] - c40 * source[97] - c39 * source[125]
626 + c4 * source[127];
627 target[30] = c40 * source[96] - c4 * source[98] - c4 * source[126]
628 + c39 * source[128];
629 target[31] = c42 * source[99] - c42 * source[101] - c43 * source[90]
630 + c43 * source[92] - c43 * source[92] + c43 * source[94]
631 - c35 * source[129] + c35 * source[131] + c41 * source[120]
632 - c41 * source[122] + c41 * source[122] - c41 * source[124];
633 target[32] = c44 * source[100] - c35 * source[91] - c35 * source[93]
634 - c17 * source[130] + c15 * source[121] + c15 * source[123];
635 target[33] = c32 * source[102] - c11 * source[95] - c11 * source[97]
636 - c45 * source[132] + c12 * source[125] + c12 * source[127];
637 target[34] = c32 * source[103] - c11 * source[96] - c11 * source[98]
638 - c45 * source[133] + c12 * source[126] + c12 * source[128];
639 target[35] = c47 * source[104] - c50 * source[99] - c50 * source[101]
640 + c51 * source[90] + c52 * source[92] + c51 * source[94]
641 - c46 * source[134] + c47 * source[129] + c47 * source[131]
642 - c48 * source[120] - c49 * source[122] - c48 * source[124];
643 target[36] = c9 * source[135] - c11 * source[137] + c9 * source[139]
644 - c9 * source[165] + c11 * source[167] - c9 * source[169]
645 - c10 * source[0] + c9 * source[2] - c10 * source[4]
646 + c10 * source[30] - c9 * source[32] + c10 * source[34]
647 - c10 * source[30] + c9 * source[32] - c10 * source[34]
648 + c10 * source[60] - c9 * source[62] + c10 * source[64];
649 target[37] = c30 * source[136] - c30 * source[138] - c30 * source[166]
650 + c30 * source[168] - c31 * source[1] + c31 * source[3]
651 + c31 * source[31] - c31 * source[33] - c31 * source[31]
652 + c31 * source[33] + c31 * source[61] - c31 * source[63];
653 target[38] = c35 * source[140] - c42 * source[142] - c35 * source[170]
654 + c42 * source[172] - c41 * source[5] + c43 * source[7]
655 + c41 * source[35] - c43 * source[37] - c41 * source[35]
656 + c43 * source[37] + c41 * source[65] - c43 * source[67];
657 target[39] = c42 * source[141] - c35 * source[143] - c42 * source[171]
658 + c35 * source[173] - c43 * source[6] + c41 * source[8]
659 + c43 * source[36] - c41 * source[38] - c43 * source[36]
660 + c41 * source[38] + c43 * source[66] - c41 * source[68];
661 target[40] = c53 * source[144] - c53 * source[146] - c54 * source[135]
662 + c54 * source[137] - c54 * source[137] + c54 * source[139]
663 - c53 * source[174] + c53 * source[176] + c54 * source[165]
664 - c54 * source[167] + c54 * source[167] - c54 * source[169]
665 - c54 * source[9] + c54 * source[11] + c55 * source[0]
666 - c55 * source[2] + c55 * source[2] - c55 * source[4]
667 + c54 * source[39] - c54 * source[41] - c55 * source[30]
668 + c55 * source[32] - c55 * source[32] + c55 * source[34]
669 - c54 * source[39] + c54 * source[41] + c55 * source[30]
670 - c55 * source[32] + c55 * source[32] - c55 * source[34]
671 + c54 * source[69] - c54 * source[71] - c55 * source[60]
672 + c55 * source[62] - c55 * source[62] + c55 * source[64];
673 target[41] = c56 * source[145] - c57 * source[136] - c57 * source[138]
674 - c56 * source[175] + c57 * source[166] + c57 * source[168]
675 - c57 * source[10] + c58 * source[1] + c58 * source[3]
676 + c57 * source[40] - c58 * source[31] - c58 * source[33]
677 - c57 * source[40] + c58 * source[31] + c58 * source[33]
678 + c57 * source[70] - c58 * source[61] - c58 * source[63];
679 target[42] = c59 * source[147] - c60 * source[140] - c60 * source[142]
680 - c59 * source[177] + c60 * source[170] + c60 * source[172]
681 - c61 * source[12] + c62 * source[5] + c62 * source[7]
682 + c61 * source[42] - c62 * source[35] - c62 * source[37]
683 - c61 * source[42] + c62 * source[35] + c62 * source[37]
684 + c61 * source[72] - c62 * source[65] - c62 * source[67];
685 target[43] = c59 * source[148] - c60 * source[141] - c60 * source[143]
686 - c59 * source[178] + c60 * source[171] + c60 * source[173]
687 - c61 * source[13] + c62 * source[6] + c62 * source[8]
688 + c61 * source[43] - c62 * source[36] - c62 * source[38]
689 - c61 * source[43] + c62 * source[36] + c62 * source[38]
690 + c61 * source[73] - c62 * source[66] - c62 * source[68];
691 target[44] = c63 * source[149] - c64 * source[144] - c64 * source[146]
692 + c65 * source[135] + c66 * source[137] + c65 * source[139]
693 - c63 * source[179] + c64 * source[174] + c64 * source[176]
694 - c65 * source[165] - c66 * source[167] - c65 * source[169]
695 - c67 * source[14] + c68 * source[9] + c68 * source[11]
696 - c69 * source[0] - c70 * source[2] - c69 * source[4]
697 + c67 * source[44] - c68 * source[39] - c68 * source[41]
698 + c69 * source[30] + c70 * source[32] + c69 * source[34]
699 - c67 * source[44] + c68 * source[39] + c68 * source[41]
700 - c69 * source[30] - c70 * source[32] - c69 * source[34]
701 + c67 * source[74] - c68 * source[69] - c68 * source[71]
702 + c69 * source[60] + c70 * source[62] + c69 * source[64];
703 target[45] = c12 * source[150] - c14 * source[152] + c12 * source[154]
704 - c13 * source[15] + c12 * source[17] - c13 * source[19]
705 - c13 * source[45] + c12 * source[47] - c13 * source[49];
706 target[46] = c32 * source[151] - c32 * source[153] - c33 * source[16]
707 + c33 * source[18] - c33 * source[46] + c33 * source[48];
708 target[47] = c17 * source[155] - c44 * source[157] - c15 * source[20]
709 + c35 * source[22] - c15 * source[50] + c35 * source[52];
710 target[48] = c44 * source[156] - c17 * source[158] - c35 * source[21]
711 + c15 * source[23] - c35 * source[51] + c15 * source[53];
712 target[49] = c56 * source[159] - c56 * source[161] - c57 * source[150]
713 + c57 * source[152] - c57 * source[152] + c57 * source[154]
714 - c57 * source[24] + c57 * source[26] + c58 * source[15]
715 - c58 * source[17] + c58 * source[17] - c58 * source[19]
716 - c57 * source[54] + c57 * source[56] + c58 * source[45]
717 - c58 * source[47] + c58 * source[47] - c58 * source[49];
718 target[50] = c71 * source[160] - c72 * source[151] - c72 * source[153]
719 - c72 * source[25] + c73 * source[16] + c73 * source[18]
720 - c72 * source[55] + c73 * source[46] + c73 * source[48];
721 target[51] = c74 * source[162] - c75 * source[155] - c75 * source[157]
722 - c76 * source[27] + c77 * source[20] + c77 * source[22]
723 - c76 * source[57] + c77 * source[50] + c77 * source[52];
724 target[52] = c74 * source[163] - c75 * source[156] - c75 * source[158]
725 - c76 * source[28] + c77 * source[21] + c77 * source[23]
726 - c76 * source[58] + c77 * source[51] + c77 * source[53];
727 target[53] = c78 * source[164] - c79 * source[159] - c79 * source[161]
728 + c66 * source[150] + c80 * source[152] + c66 * source[154]
729 - c81 * source[29] + c63 * source[24] + c63 * source[26]
730 - c70 * source[15] - c82 * source[17] - c70 * source[19]
731 - c81 * source[59] + c63 * source[54] + c63 * source[56]
732 - c70 * source[45] - c82 * source[47] - c70 * source[49];
733 target[54] = c15 * source[180] - c17 * source[182] + c15 * source[184]
734 - c16 * source[75] + c18 * source[77] - c16 * source[79]
735 - c16 * source[105] + c18 * source[107] - c16 * source[109];
736 target[55] = c34 * source[181] - c34 * source[183] - c35 * source[76]
737 + c35 * source[78] - c35 * source[106] + c35 * source[108];
738 target[56] = c45 * source[185] - c32 * source[187] - c12 * source[80]
739 + c11 * source[82] - c12 * source[110] + c11 * source[112];
740 target[57] = c32 * source[186] - c45 * source[188] - c11 * source[81]
741 + c12 * source[83] - c11 * source[111] + c12 * source[113];
742 target[58] = c59 * source[189] - c59 * source[191] - c61 * source[180]
743 + c61 * source[182] - c61 * source[182] + c61 * source[184]
744 - c60 * source[84] + c60 * source[86] + c62 * source[75]
745 - c62 * source[77] + c62 * source[77] - c62 * source[79]
746 - c60 * source[114] + c60 * source[116] + c62 * source[105]
747 - c62 * source[107] + c62 * source[107] - c62 * source[109];
748 target[59] = c74 * source[190] - c76 * source[181] - c76 * source[183]
749 - c75 * source[85] + c77 * source[76] + c77 * source[78]
750 - c75 * source[115] + c77 * source[106] + c77 * source[108];
751 target[60] = c83 * source[192] - c72 * source[185] - c72 * source[187]
752 - c72 * source[87] + c84 * source[80] + c84 * source[82]
753 - c72 * source[117] + c84 * source[110] + c84 * source[112];
754 target[61] = c83 * source[193] - c72 * source[186] - c72 * source[188]
755 - c72 * source[88] + c84 * source[81] + c84 * source[83]
756 - c72 * source[118] + c84 * source[111] + c84 * source[113];
757 target[62] = c85 * source[194] - c86 * source[189] - c86 * source[191]
758 + c87 * source[180] + c88 * source[182] + c87 * source[184]
759 - c88 * source[89] + c89 * source[84] + c89 * source[86]
760 - c90 * source[75] - c91 * source[77] - c90 * source[79]
761 - c88 * source[119] + c89 * source[114] + c89 * source[116]
762 - c90 * source[105] - c91 * source[107] - c90 * source[109];
763 target[63] = c15 * source[195] - c17 * source[197] + c15 * source[199]
764 - c16 * source[90] + c18 * source[92] - c16 * source[94]
765 - c16 * source[120] + c18 * source[122] - c16 * source[124];
766 target[64] = c34 * source[196] - c34 * source[198] - c35 * source[91]
767 + c35 * source[93] - c35 * source[121] + c35 * source[123];
768 target[65] = c45 * source[200] - c32 * source[202] - c12 * source[95]
769 + c11 * source[97] - c12 * source[125] + c11 * source[127];
770 target[66] = c32 * source[201] - c45 * source[203] - c11 * source[96]
771 + c12 * source[98] - c11 * source[126] + c12 * source[128];
772 target[67] = c59 * source[204] - c59 * source[206] - c61 * source[195]
773 + c61 * source[197] - c61 * source[197] + c61 * source[199]
774 - c60 * source[99] + c60 * source[101] + c62 * source[90]
775 - c62 * source[92] + c62 * source[92] - c62 * source[94]
776 - c60 * source[129] + c60 * source[131] + c62 * source[120]
777 - c62 * source[122] + c62 * source[122] - c62 * source[124];
778 target[68] = c74 * source[205] - c76 * source[196] - c76 * source[198]
779 - c75 * source[100] + c77 * source[91] + c77 * source[93]
780 - c75 * source[130] + c77 * source[121] + c77 * source[123];
781 target[69] = c83 * source[207] - c72 * source[200] - c72 * source[202]
782 - c72 * source[102] + c84 * source[95] + c84 * source[97]
783 - c72 * source[132] + c84 * source[125] + c84 * source[127];
784 target[70] = c83 * source[208] - c72 * source[201] - c72 * source[203]
785 - c72 * source[103] + c84 * source[96] + c84 * source[98]
786 - c72 * source[133] + c84 * source[126] + c84 * source[128];
787 target[71] = c85 * source[209] - c86 * source[204] - c86 * source[206]
788 + c87 * source[195] + c88 * source[197] + c87 * source[199]
789 - c88 * source[104] + c89 * source[99] + c89 * source[101]
790 - c90 * source[90] - c91 * source[92] - c90 * source[94]
791 - c88 * source[134] + c89 * source[129] + c89 * source[131]
792 - c90 * source[120] - c91 * source[122] - c90 * source[124];
793 target[72] = c19 * source[210] - c23 * source[212] + c19 * source[214]
794 - c20 * source[135] + c24 * source[137] - c20 * source[139]
795 - c20 * source[165] + c24 * source[167] - c20 * source[169]
796 + c21 * source[0] - c25 * source[2] + c21 * source[4]
797 + c22 * source[30] - c26 * source[32] + c22 * source[34]
798 + c21 * source[60] - c25 * source[62] + c21 * source[64];
799 target[73] = c36 * source[211] - c36 * source[213] - c37 * source[136]
800 + c37 * source[138] - c37 * source[166] + c37 * source[168]
801 + c38 * source[1] - c38 * source[3] + c20 * source[31]
802 - c20 * source[33] + c38 * source[61] - c38 * source[63];
803 target[74] = c46 * source[215] - c47 * source[217] - c47 * source[140]
804 + c50 * source[142] - c47 * source[170] + c50 * source[172]
805 + c48 * source[5] - c51 * source[7] + c49 * source[35]
806 - c52 * source[37] + c48 * source[65] - c51 * source[67];
807 target[75] = c47 * source[216] - c46 * source[218] - c50 * source[141]
808 + c47 * source[143] - c50 * source[171] + c47 * source[173]
809 + c51 * source[6] - c48 * source[8] + c52 * source[36]
810 - c49 * source[38] + c51 * source[66] - c48 * source[68];
811 target[76] = c63 * source[219] - c63 * source[221] - c67 * source[210]
812 + c67 * source[212] - c67 * source[212] + c67 * source[214]
813 - c64 * source[144] + c64 * source[146] + c68 * source[135]
814 - c68 * source[137] + c68 * source[137] - c68 * source[139]
815 - c64 * source[174] + c64 * source[176] + c68 * source[165]
816 - c68 * source[167] + c68 * source[167] - c68 * source[169]
817 + c65 * source[9] - c65 * source[11] - c69 * source[0]
818 + c69 * source[2] - c69 * source[2] + c69 * source[4]
819 + c66 * source[39] - c66 * source[41] - c70 * source[30]
820 + c70 * source[32] - c70 * source[32] + c70 * source[34]
821 + c65 * source[69] - c65 * source[71] - c69 * source[60]
822 + c69 * source[62] - c69 * source[62] + c69 * source[64];
823 target[77] = c78 * source[220] - c81 * source[211] - c81 * source[213]
824 - c79 * source[145] + c63 * source[136] + c63 * source[138]
825 - c79 * source[175] + c63 * source[166] + c63 * source[168]
826 + c66 * source[10] - c70 * source[1] - c70 * source[3]
827 + c80 * source[40] - c82 * source[31] - c82 * source[33]
828 + c66 * source[70] - c70 * source[61] - c70 * source[63];
829 target[78] = c85 * source[222] - c88 * source[215] - c88 * source[217]
830 - c86 * source[147] + c89 * source[140] + c89 * source[142]
831 - c86 * source[177] + c89 * source[170] + c89 * source[172]
832 + c87 * source[12] - c90 * source[5] - c90 * source[7]
833 + c88 * source[42] - c91 * source[35] - c91 * source[37]
834 + c87 * source[72] - c90 * source[65] - c90 * source[67];
835 target[79] = c85 * source[223] - c88 * source[216] - c88 * source[218]
836 - c86 * source[148] + c89 * source[141] + c89 * source[143]
837 - c86 * source[178] + c89 * source[171] + c89 * source[173]
838 + c87 * source[13] - c90 * source[6] - c90 * source[8]
839 + c88 * source[43] - c91 * source[36] - c91 * source[38]
840 + c87 * source[73] - c90 * source[66] - c90 * source[68];
841 target[80] = source[224] - c92 * source[219] - c92 * source[221]
842 + c93 * source[210] + c94 * source[212] + c93 * source[214]
843 - c92 * source[149] + c95 * source[144] + c95 * source[146]
844 - c96 * source[135] - c97 * source[137] - c96 * source[139]
845 - c92 * source[179] + c95 * source[174] + c95 * source[176]
846 - c96 * source[165] - c97 * source[167] - c96 * source[169]
847 + c93 * source[14] - c96 * source[9] - c96 * source[11]
848 + c98 * source[0] + c99 * source[2] + c98 * source[4]
849 + c94 * source[44] - c97 * source[39] - c97 * source[41]
850 + c99 * source[30] + c100 * source[32] + c99 * source[34]
851 + c93 * source[74] - c96 * source[69] - c96 * source[71]
852 + c98 * source[60] + c99 * source[62] + c98 * source[64];
853 }
854 }
855
856