v0.14.0
Loading...
Searching...
No Matches
Tensor4_times_Tensor2.hpp
Go to the documentation of this file.
1/* This file has all of the declarations for expressions like
2 Tensor4*Tensor2 and Tensor2*Tensor4, yielding a
3 Tensor2. */
4
5#pragma once
6
7namespace FTensor
8{
9 /* A(i,j,k,l)*B(k,l) */
10
11 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
12 int Dim3, char i, char j, char k, char l>
13 class Tensor4_times_Tensor2_23
14 {
15 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> iterA;
16 const Tensor2_Expr<B, U, Dim2, Dim3, k, l> iterB;
17
18 template <int Current_Dim0, int Current_Dim1>
19 typename promote<T, U>::V
20 eval(const int N1, const int N2, const Number<Current_Dim0> &,
21 const Number<Current_Dim1> &) const
22 {
23 return iterA(N1, N2, Current_Dim0 - 1, Current_Dim1 - 1)
24 * iterB(Current_Dim0 - 1, Current_Dim1 - 1)
25 + eval(N1, N2, Number<Current_Dim0 - 1>(),
26 Number<Current_Dim1>());
27 }
28 template <int Current_Dim1>
29 typename promote<T, U>::V
30 eval(const int N1, const int N2, const Number<1> &,
31 const Number<Current_Dim1> &) const
32 {
33 return iterA(N1, N2, 0, Current_Dim1 - 1) * iterB(0, Current_Dim1 - 1)
34 + eval(N1, N2, Number<Dim2>(), Number<Current_Dim1 - 1>());
35 }
36 typename promote<T, U>::V eval(const int N1, const int N2,
37 const Number<1> &, const Number<1> &) const
38 {
39 return iterA(N1, N2, 0, 0) * iterB(0, 0);
40 }
41
42 public:
43 typename promote<T, U>::V operator()(const int N1, const int N2) const
44 {
45 return eval(N1, N2, Number<Dim2>(), Number<Dim3>());
46 }
47
48 Tensor4_times_Tensor2_23(
49 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a,
50 const Tensor2_Expr<B, U, Dim2, Dim3, k, l> &b)
51 : iterA(a), iterB(b)
52 {}
53 };
54
55 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
56 int Dim3, char i, char j, char k, char l>
57 inline Tensor2_Expr<const Tensor4_times_Tensor2_23<A, B, T, U, Dim0, Dim1,
58 Dim2, Dim3, i, j, k, l>,
59 typename promote<T, U>::V, Dim0, Dim1, i, j>
60 operator*(const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a,
61 const Tensor2_Expr<B, U, Dim2, Dim3, k, l> &b)
62 {
63 typedef const Tensor4_times_Tensor2_23<A, B, T, U, Dim0, Dim1, Dim2, Dim3,
64 i, j, k, l>
65 TensorExpr;
66 return Tensor2_Expr<TensorExpr, typename promote<T, U>::V, Dim0, Dim1, i,
67 j>(TensorExpr(a, b));
68 }
69
70 /* B(k,l)*A(i,j,k,l) */
71
72 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
73 int Dim3, char i, char j, char k, char l>
74 inline Tensor2_Expr<const Tensor4_times_Tensor2_23<A, B, T, U, Dim0, Dim1,
75 Dim2, Dim3, i, j, k, l>,
76 typename promote<T, U>::V, Dim0, Dim1, i, j>
77 operator*(const Tensor2_Expr<B, U, Dim2, Dim3, k, l> &b,
78 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a)
79 {
80 typedef const Tensor4_times_Tensor2_23<A, B, T, U, Dim0, Dim1, Dim2, Dim3,
81 i, j, k, l>
82 TensorExpr;
83 return Tensor2_Expr<TensorExpr, typename promote<T, U>::V, Dim0, Dim1, i,
84 j>(TensorExpr(a, b));
85 }
86
87 /* A(i,j,k,l)*B(l,k) */
88
89 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
90 int Dim3, char i, char j, char k, char l>
91 class Tensor4_times_Tensor2_32
92 {
93 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> iterA;
94 const Tensor2_Expr<B, U, Dim3, Dim2, l, k> iterB;
95
96 template <int Current_Dim0, int Current_Dim1>
97 typename promote<T, U>::V
98 eval(const int N1, const int N2, const Number<Current_Dim0> &,
99 const Number<Current_Dim1> &) const
100 {
101 return iterA(N1, N2, Current_Dim0 - 1, Current_Dim1 - 1)
102 * iterB(Current_Dim1 - 1, Current_Dim0 - 1)
103 + eval(N1, N2, Number<Current_Dim0 - 1>(),
104 Number<Current_Dim1>());
105 }
106 template <int Current_Dim1>
107 typename promote<T, U>::V
108 eval(const int N1, const int N2, const Number<1> &,
109 const Number<Current_Dim1> &) const
110 {
111 return iterA(N1, N2, 0, Current_Dim1 - 1) * iterB(Current_Dim1 - 1, 0)
112 + eval(N1, N2, Number<Dim2>(), Number<Current_Dim1 - 1>());
113 }
114 typename promote<T, U>::V eval(const int N1, const int N2,
115 const Number<1> &, const Number<1> &) const
116 {
117 return iterA(N1, N2, 0, 0) * iterB(0, 0);
118 }
119
120 public:
121 typename promote<T, U>::V operator()(const int N1, const int N2) const
122 {
123 return eval(N1, N2, Number<Dim2>(), Number<Dim3>());
124 }
125
126 Tensor4_times_Tensor2_32(
127 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a,
128 const Tensor2_Expr<B, U, Dim3, Dim2, l, k> &b)
129 : iterA(a), iterB(b)
130 {}
131 };
132
133 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
134 int Dim3, char i, char j, char k, char l>
135 inline Tensor2_Expr<const Tensor4_times_Tensor2_32<A, B, T, U, Dim0, Dim1,
136 Dim2, Dim3, i, j, k, l>,
137 typename promote<T, U>::V, Dim0, Dim1, i, j>
138 operator*(const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a,
139 const Tensor2_Expr<B, U, Dim3, Dim2, l, k> &b)
140 {
141 typedef const Tensor4_times_Tensor2_32<A, B, T, U, Dim0, Dim1, Dim2, Dim3,
142 i, j, k, l>
143 TensorExpr;
144 return Tensor2_Expr<TensorExpr, typename promote<T, U>::V, Dim0, Dim1, i,
145 j>(TensorExpr(a, b));
146 }
147
148 /* B(l,k)*A(i,j,k,l) */
149
150 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
151 int Dim3, char i, char j, char k, char l>
152 inline Tensor2_Expr<const Tensor4_times_Tensor2_32<A, B, T, U, Dim0, Dim1,
153 Dim2, Dim3, i, j, k, l>,
154 typename promote<T, U>::V, Dim0, Dim1, i, j>
155 operator*(const Tensor2_Expr<B, U, Dim3, Dim2, l, k> &b,
156 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a)
157 {
158 typedef const Tensor4_times_Tensor2_32<A, B, T, U, Dim0, Dim1, Dim2, Dim3,
159 i, j, k, l>
160 TensorExpr;
161 return Tensor2_Expr<TensorExpr, typename promote<T, U>::V, Dim0, Dim1, i,
162 j>(TensorExpr(a, b));
163 }
164
165 /* A(i,j,k,l)*B(i,l) */
166
167 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
168 int Dim3, char i, char j, char k, char l>
169 class Tensor4_times_Tensor2_03
170 {
171 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> iterA;
172 const Tensor2_Expr<B, U, Dim0, Dim3, i, l> iterB;
173
174 template <int Current_Dim0, int Current_Dim1>
175 typename promote<T, U>::V
176 eval(const int N1, const int N2, const Number<Current_Dim0> &,
177 const Number<Current_Dim1> &) const
178 {
179 return iterA(Current_Dim0 - 1, N1, N2, Current_Dim1 - 1)
180 * iterB(Current_Dim0 - 1, Current_Dim1 - 1)
181 + eval(N1, N2, Number<Current_Dim0 - 1>(),
182 Number<Current_Dim1>());
183 }
184 template <int Current_Dim1>
185 typename promote<T, U>::V
186 eval(const int N1, const int N2, const Number<1> &,
187 const Number<Current_Dim1> &) const
188 {
189 return iterA(0, N1, N2, Current_Dim1 - 1) * iterB(0, Current_Dim1 - 1)
190 + eval(N1, N2, Number<Dim0>(), Number<Current_Dim1 - 1>());
191 }
192 typename promote<T, U>::V eval(const int N1, const int N2,
193 const Number<1> &, const Number<1> &) const
194 {
195 return iterA(0, N1, N2, 0) * iterB(0, 0);
196 }
197
198 public:
199 typename promote<T, U>::V operator()(const int N1, const int N2) const
200 {
201 return eval(N1, N2, Number<Dim0>(), Number<Dim3>());
202 }
203
204 Tensor4_times_Tensor2_03(
205 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a,
206 const Tensor2_Expr<B, U, Dim0, Dim3, i, l> &b)
207 : iterA(a), iterB(b)
208 {}
209 };
210
211 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
212 int Dim3, char i, char j, char k, char l>
213 inline Tensor2_Expr<const Tensor4_times_Tensor2_03<A, B, T, U, Dim0, Dim1,
214 Dim2, Dim3, i, j, k, l>,
215 typename promote<T, U>::V, Dim1, Dim2, j, k>
216 operator*(const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a,
217 const Tensor2_Expr<B, U, Dim0, Dim3, i, l> &b)
218 {
219 typedef const Tensor4_times_Tensor2_03<A, B, T, U, Dim0, Dim1, Dim2, Dim3,
220 i, j, k, l>
221 TensorExpr;
222 return Tensor2_Expr<TensorExpr, typename promote<T, U>::V, Dim1, Dim2, j,
223 k>(TensorExpr(a, b));
224 }
225
226 /* B(i,l)*A(i,j,k,l) */
227
228 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
229 int Dim3, char i, char j, char k, char l>
230 inline Tensor2_Expr<const Tensor4_times_Tensor2_03<A, B, T, U, Dim0, Dim1,
231 Dim2, Dim3, i, j, k, l>,
232 typename promote<T, U>::V, Dim1, Dim2, j, k>
233 operator*(const Tensor2_Expr<B, U, Dim0, Dim3, i, l> &b,
234 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a)
235 {
236 typedef const Tensor4_times_Tensor2_03<A, B, T, U, Dim0, Dim1, Dim2, Dim3,
237 i, j, k, l>
238 TensorExpr;
239 return Tensor2_Expr<TensorExpr, typename promote<T, U>::V, Dim1, Dim2, j,
240 k>(TensorExpr(a, b));
241 }
242
243 /* A(i,j,k,l)*B(l,i) */
244
245 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
246 int Dim3, char i, char j, char k, char l>
247 class Tensor4_times_Tensor2_30
248 {
249 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> iterA;
250 const Tensor2_Expr<B, U, Dim3, Dim0, l, i> iterB;
251
252 template <int Current_Dim0, int Current_Dim1>
253 typename promote<T, U>::V
254 eval(const int N1, const int N2, const Number<Current_Dim0> &,
255 const Number<Current_Dim1> &) const
256 {
257 return iterA(Current_Dim0 - 1, N1, N2, Current_Dim1 - 1)
258 * iterB(Current_Dim1 - 1, Current_Dim0 - 1)
259 + eval(N1, N2, Number<Current_Dim0 - 1>(),
260 Number<Current_Dim1>());
261 }
262 template <int Current_Dim1>
263 typename promote<T, U>::V
264 eval(const int N1, const int N2, const Number<1> &,
265 const Number<Current_Dim1> &) const
266 {
267 return iterA(0, N1, N2, Current_Dim1 - 1) * iterB(Current_Dim1 - 1, 0)
268 + eval(N1, N2, Number<Dim0>(), Number<Current_Dim1 - 1>());
269 }
270 typename promote<T, U>::V eval(const int N1, const int N2,
271 const Number<1> &, const Number<1> &) const
272 {
273 return iterA(0, N1, N2, 0) * iterB(0, 0);
274 }
275
276 public:
277 typename promote<T, U>::V operator()(const int N1, const int N2) const
278 {
279 return eval(N1, N2, Number<Dim0>(), Number<Dim3>());
280 }
281
282 Tensor4_times_Tensor2_30(
283 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a,
284 const Tensor2_Expr<B, U, Dim3, Dim0, l, i> &b)
285 : iterA(a), iterB(b)
286 {}
287 };
288
289 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
290 int Dim3, char i, char j, char k, char l>
291 inline Tensor2_Expr<const Tensor4_times_Tensor2_30<A, B, T, U, Dim0, Dim1,
292 Dim2, Dim3, i, j, k, l>,
293 typename promote<T, U>::V, Dim1, Dim2, j, k>
294 operator*(const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a,
295 const Tensor2_Expr<B, U, Dim3, Dim0, l, i> &b)
296 {
297 typedef const Tensor4_times_Tensor2_30<A, B, T, U, Dim0, Dim1, Dim2, Dim3,
298 i, j, k, l>
299 TensorExpr;
300 return Tensor2_Expr<TensorExpr, typename promote<T, U>::V, Dim1, Dim2, j,
301 k>(TensorExpr(a, b));
302 }
303
304 /* B(l,i)*A(i,j,k,l) */
305
306 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
307 int Dim3, char i, char j, char k, char l>
308 inline Tensor2_Expr<const Tensor4_times_Tensor2_30<A, B, T, U, Dim0, Dim1,
309 Dim2, Dim3, i, j, k, l>,
310 typename promote<T, U>::V, Dim1, Dim2, j, k>
311 operator*(const Tensor2_Expr<B, U, Dim3, Dim0, l, i> &b,
312 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a)
313 {
314 typedef const Tensor4_times_Tensor2_30<A, B, T, U, Dim0, Dim1, Dim2, Dim3,
315 i, j, k, l>
316 TensorExpr;
317 return Tensor2_Expr<TensorExpr, typename promote<T, U>::V, Dim1, Dim2, j,
318 k>(TensorExpr(a, b));
319 }
320
321 /* A(i,j,k,l)*B(j,l) */
322
323 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
324 int Dim3, char i, char j, char k, char l>
325 class Tensor4_times_Tensor2_13
326 {
327 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> iterA;
328 const Tensor2_Expr<B, U, Dim2, Dim3, j, l> iterB;
329
330 template <int Current_Dim0, int Current_Dim1>
331 typename promote<T, U>::V
332 eval(const int N1, const int N2, const Number<Current_Dim0> &,
333 const Number<Current_Dim1> &) const
334 {
335 return iterA(N1, Current_Dim0 - 1, N2, Current_Dim1 - 1)
336 * iterB(Current_Dim0 - 1, Current_Dim1 - 1)
337 + eval(N1, N2, Number<Current_Dim0 - 1>(),
338 Number<Current_Dim1>());
339 }
340 template <int Current_Dim1>
341 typename promote<T, U>::V
342 eval(const int N1, const int N2, const Number<1> &,
343 const Number<Current_Dim1> &) const
344 {
345 return iterA(N1, 0, N2, Current_Dim1 - 1) * iterB(0, Current_Dim1 - 1)
346 + eval(N1, N2, Number<Dim0>(), Number<Current_Dim1 - 1>());
347 }
348 typename promote<T, U>::V eval(const int N1, const int N2,
349 const Number<1> &, const Number<1> &) const
350 {
351 return iterA(N1, 0, N2, 0) * iterB(0, 0);
352 }
353
354 public:
355 typename promote<T, U>::V operator()(const int N1, const int N2) const
356 {
357 return eval(N1, N2, Number<Dim0>(), Number<Dim3>());
358 }
359
360 Tensor4_times_Tensor2_13(
361 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a,
362 const Tensor2_Expr<B, U, Dim2, Dim3, j, l> &b)
363 : iterA(a), iterB(b)
364 {}
365 };
366
367 /* B(j,l)*A(i,j,k,l) */
368
369 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
370 int Dim3, char i, char j, char k, char l>
371 inline Tensor2_Expr<const Tensor4_times_Tensor2_13<A, B, T, U, Dim0, Dim1,
372 Dim2, Dim3, i, j, k, l>,
373 typename promote<T, U>::V, Dim0, Dim2, i, k>
374 operator*(const Tensor2_Expr<B, U, Dim1, Dim3, j, l> &b,
375 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a)
376 {
377 typedef const Tensor4_times_Tensor2_13<A, B, T, U, Dim0, Dim1, Dim2, Dim3,
378 i, j, k, l>
379 TensorExpr;
380 return Tensor2_Expr<TensorExpr, typename promote<T, U>::V, Dim1, Dim2, i,
381 k>(TensorExpr(a, b));
382 }
383
384 /* This file has all of the declarations for expressions like
385 Tensor4*Tensor2 and Tensor2*Tensor4, yielding a
386 Tensor4. */
387
388 // TODO: Check dimensions could be errors
389
390 /* A(i,j,k,l)*B(l,m) */
391
392 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
393 int Dim3, int Dim4, char i, char j, char k, char l, char m>
394 class Tensor4_times_Tensor2_3_1
395 {
396 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> iterA;
397 const Tensor2_Expr<B, U, Dim3, Dim4, l, m> iterB;
398
399 template <int Current_Dim0>
400 typename promote<T, U>::V
401 eval(const int N1, const int N2, const int N3, const int N4,
402 const Number<Current_Dim0> &) const
403 {
404 return iterA(N1, N2, N3, Current_Dim0 - 1) * iterB(Current_Dim0 - 1, N4)
405 + eval(N1, N2, N3, N4, Number<Current_Dim0 - 1>());
406 }
407 typename promote<T, U>::V eval(const int N1, const int N2, const int N3,
408 const int N4, const Number<1> &) const
409 {
410 return iterA(N1, N2, N3, 0) * iterB(0, N4);
411 }
412
413 public:
414 Tensor4_times_Tensor2_3_1(
415 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a,
416 const Tensor2_Expr<B, U, Dim3, Dim4, l, m> &b)
417 : iterA(a), iterB(b)
418 {}
419 typename promote<T, U>::V
420 operator()(const int N1, const int N2, const int N3, const int N4) const
421 {
422 return eval(N1, N2, N3, N4, Number<Dim3>());
423 }
424 };
425
426 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
427 int Dim3, int Dim4, char i, char j, char k, char l, char m>
428 inline const Tensor4_Expr<
429 const Tensor4_times_Tensor2_3_1<A, B, T, U, Dim0, Dim1, Dim2, Dim3, Dim4,
430 i, j, k, l, m>,
431 typename promote<T, U>::V, Dim0, Dim1, Dim2, Dim4, i, j, k, m>
432 operator*(const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a,
433 const Tensor2_Expr<B, U, Dim3, Dim4, l, m> &b)
434 {
435 typedef const Tensor4_times_Tensor2_3_1<A, B, T, U, Dim0, Dim1, Dim2, Dim3,
436 Dim4, i, j, k, l, m>
437 TensorExpr;
438 return Tensor4_Expr<TensorExpr, typename promote<T, U>::V, Dim0, Dim1,
439 Dim2, Dim4, i, j, k, m>(TensorExpr(a, b));
440 }
441
442 /* B(l,m)*A(i,j,k,l) */
443
444 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
445 int Dim3, int Dim4, char i, char j, char k, char l, char m>
446 inline const Tensor4_Expr<
447 const Tensor4_times_Tensor2_3_1<A, B, T, U, Dim0, Dim1, Dim2, Dim3, Dim4,
448 i, j, k, l, m>,
449 typename promote<T, U>::V, Dim0, Dim1, Dim2, Dim4, i, j, k, m>
450 operator*(const Tensor2_Expr<B, U, Dim3, Dim4, l, m> &b,
451 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a)
452 {
453 typedef const Tensor4_times_Tensor2_3_1<A, B, T, U, Dim0, Dim1, Dim2, Dim3,
454 Dim4, i, j, k, l, m>
455 TensorExpr;
456 return Tensor4_Expr<TensorExpr, typename promote<T, U>::V, Dim0, Dim1,
457 Dim2, Dim4, i, j, k, m>(TensorExpr(a, b));
458 }
459
460 /* A(i,j,k,l)*B(m,l) */
461
462 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
463 int Dim3, int Dim4, char i, char j, char k, char l, char m>
464 class Tensor4_times_Tensor2_3_0
465 {
466 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> iterA;
467 const Tensor2_Expr<B, U, Dim4, Dim3, m, l> iterB;
468
469 template <int Current_Dim0>
470 typename promote<T, U>::V
471 eval(const int N1, const int N2, const int N3, const int N4,
472 const Number<Current_Dim0> &) const
473 {
474 return iterA(N1, N2, N3, Current_Dim0 - 1) * iterB(N4, Current_Dim0 - 1)
475 + eval(N1, N2, N3, N4, Number<Current_Dim0 - 1>());
476 }
477 typename promote<T, U>::V eval(const int N1, const int N2, const int N3,
478 const int N4, const Number<1> &) const
479 {
480 return iterA(N1, N2, N3, 0) * iterB(N4, 0);
481 }
482
483 public:
484 Tensor4_times_Tensor2_3_0(
485 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a,
486 const Tensor2_Expr<B, U, Dim4, Dim3, m, l> &b)
487 : iterA(a), iterB(b)
488 {}
489 typename promote<T, U>::V
490 operator()(const int N1, const int N2, const int N3, const int N4) const
491 {
492 return eval(N1, N2, N3, N4, Number<Dim3>());
493 }
494 };
495
496 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
497 int Dim3, int Dim4, char i, char j, char k, char l, char m>
498 inline const Tensor4_Expr<
499 const Tensor4_times_Tensor2_3_0<A, B, T, U, Dim0, Dim1, Dim2, Dim3, Dim4,
500 i, j, k, l, m>,
501 typename promote<T, U>::V, Dim0, Dim1, Dim2, Dim4, i, j, k, m>
502 operator*(const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a,
503 const Tensor2_Expr<B, U, Dim4, Dim3, m, l> &b)
504 {
505 typedef const Tensor4_times_Tensor2_3_0<A, B, T, U, Dim0, Dim1, Dim2, Dim3,
506 Dim4, i, j, k, l, m>
507 TensorExpr;
508 return Tensor4_Expr<TensorExpr, typename promote<T, U>::V, Dim0, Dim1,
509 Dim2, Dim4, i, j, k, m>(TensorExpr(a, b));
510 }
511
512 /* B(m,l)*A(i,j,k,l) */
513
514 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
515 int Dim3, int Dim4, char i, char j, char k, char l, char m>
516 inline const Tensor4_Expr<
517 const Tensor4_times_Tensor2_3_0<A, B, T, U, Dim0, Dim1, Dim2, Dim3, Dim4,
518 i, j, k, l, m>,
519 typename promote<T, U>::V, Dim0, Dim1, Dim2, Dim4, i, j, k, m>
520 operator*(const Tensor2_Expr<B, U, Dim4, Dim3, m, l> &b,
521 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a)
522 {
523 typedef const Tensor4_times_Tensor2_3_0<A, B, T, U, Dim0, Dim1, Dim2, Dim3,
524 Dim4, i, j, k, l, m>
525 TensorExpr;
526 return Tensor4_Expr<TensorExpr, typename promote<T, U>::V, Dim0, Dim1,
527 Dim2, Dim4, i, j, k, m>(TensorExpr(a, b));
528 }
529
530 /* A(i,j,k,l)*B(j,m) */
531
532 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
533 int Dim3, int Dim4, char i, char j, char k, char l, char m>
534 class Tensor4_times_Tensor2_1_0
535 {
536 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> iterA;
537 const Tensor2_Expr<B, U, Dim1, Dim4, j, m> iterB;
538
539 template <int Current_Dim0>
540 typename promote<T, U>::V
541 eval(const int N1, const int N2, const int N3, const int N4,
542 const Number<Current_Dim0> &) const
543 {
544 return iterA(N1, Current_Dim0 - 1, N3, N4) * iterB(Current_Dim0 - 1, N2)
545 + eval(N1, N2, N3, N4, Number<Current_Dim0 - 1>());
546 }
547 typename promote<T, U>::V eval(const int N1, const int N2, const int N3,
548 const int N4, const Number<1> &) const
549 {
550 return iterA(N1, 0, N3, N4) * iterB(0, N2);
551 }
552
553 public:
554 Tensor4_times_Tensor2_1_0(
555 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a,
556 const Tensor2_Expr<B, U, Dim1, Dim4, j, m> &b)
557 : iterA(a), iterB(b)
558 {}
559 typename promote<T, U>::V
560 operator()(const int N1, const int N2, const int N3, const int N4) const
561 {
562 return eval(N1, N2, N3, N4, Number<Dim1>());
563 }
564 };
565
566 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
567 int Dim3, int Dim4, char i, char j, char k, char l, char m>
568 inline const Tensor4_Expr<
569 const Tensor4_times_Tensor2_1_0<A, B, T, U, Dim0, Dim1, Dim2, Dim3, Dim4,
570 i, j, k, l, m>,
571 typename promote<T, U>::V, Dim0, Dim4, Dim2, Dim3, i, m, k, l>
572 operator*(const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a,
573 const Tensor2_Expr<B, U, Dim1, Dim4, j, m> &b)
574 {
575 typedef const Tensor4_times_Tensor2_1_0<A, B, T, U, Dim0, Dim1, Dim2, Dim3,
576 Dim4, i, j, k, l, m>
577 TensorExpr;
578 return Tensor4_Expr<TensorExpr, typename promote<T, U>::V, Dim0, Dim1,
579 Dim4, Dim3, i, m, k, l>(TensorExpr(a, b));
580 }
581
582 /* B(j,m)*A(i,j,k,l) */
583
584 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
585 int Dim3, int Dim4, char i, char j, char k, char l, char m>
586 inline const Tensor4_Expr<
587 const Tensor4_times_Tensor2_1_0<A, B, T, U, Dim0, Dim1, Dim2, Dim3, Dim4,
588 i, j, k, l, m>,
589 typename promote<T, U>::V, Dim0, Dim4, Dim2, Dim3, i, m, k, l>
590 operator*(const Tensor2_Expr<B, U, Dim1, Dim4, j, m> &b,
591 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a)
592 {
593 typedef const Tensor4_times_Tensor2_1_0<A, B, T, U, Dim0, Dim1, Dim2, Dim3,
594 Dim4, i, j, k, l, m>
595 TensorExpr;
596 return Tensor4_Expr<TensorExpr, typename promote<T, U>::V, Dim0, Dim1,
597 Dim4, Dim3, i, m, k, l>(TensorExpr(a, b));
598 }
599
600 /* A(i,j,k,l)*B(m,j) */
601
602 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
603 int Dim3, int Dim4, char i, char j, char k, char l, char m>
604 class Tensor4_times_Tensor2_1_1
605 {
606 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> iterA;
607 const Tensor2_Expr<B, U, Dim4, Dim1, m, j> iterB;
608
609 template <int Current_Dim0>
610 typename promote<T, U>::V
611 eval(const int N1, const int N2, const int N3, const int N4,
612 const Number<Current_Dim0> &) const
613 {
614 return iterA(N1, Current_Dim0 - 1, N3, N4) * iterB(N2, Current_Dim0 - 1)
615 + eval(N1, N2, N3, N4, Number<Current_Dim0 - 1>());
616 }
617 typename promote<T, U>::V eval(const int N1, const int N2, const int N3,
618 const int N4, const Number<1> &) const
619 {
620 return iterA(N1, 0, N3, N4) * iterB(N2, 0);
621 }
622
623 public:
624 Tensor4_times_Tensor2_1_1(
625 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a,
626 const Tensor2_Expr<B, U, Dim4, Dim1, m, j> &b)
627 : iterA(a), iterB(b)
628 {}
629 typename promote<T, U>::V
630 operator()(const int N1, const int N2, const int N3, const int N4) const
631 {
632 return eval(N1, N2, N3, N4, Number<Dim1>());
633 }
634 };
635
636 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
637 int Dim3, int Dim4, char i, char j, char k, char l, char m>
638 inline const Tensor4_Expr<
639 const Tensor4_times_Tensor2_1_1<A, B, T, U, Dim0, Dim1, Dim2, Dim3, Dim4,
640 i, j, k, l, m>,
641 typename promote<T, U>::V, Dim0, Dim4, Dim2, Dim3, i, m, k, l>
642 operator*(const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a,
643 const Tensor2_Expr<B, U, Dim4, Dim1, m, j> &b)
644 {
645 typedef const Tensor4_times_Tensor2_1_1<A, B, T, U, Dim0, Dim1, Dim2, Dim3,
646 Dim4, i, j, k, l, m>
647 TensorExpr;
648 return Tensor4_Expr<TensorExpr, typename promote<T, U>::V, Dim0, Dim4,
649 Dim2, Dim3, i, m, k, l>(TensorExpr(a, b));
650 }
651
652 /* B(m,j)*A(i,j,k,l) */
653
654 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
655 int Dim3, int Dim4, char i, char j, char k, char l, char m>
656 inline const Tensor4_Expr<
657 const Tensor4_times_Tensor2_1_1<A, B, T, U, Dim0, Dim1, Dim2, Dim3, Dim4,
658 i, j, k, l, m>,
659 typename promote<T, U>::V, Dim0, Dim4, Dim2, Dim3, i, m, k, l>
660 operator*(const Tensor2_Expr<B, U, Dim4, Dim1, m, j> &b,
661 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a)
662 {
663 typedef const Tensor4_times_Tensor2_1_1<A, B, T, U, Dim0, Dim1, Dim2, Dim3,
664 Dim4, i, j, k, l, m>
665 TensorExpr;
666 return Tensor4_Expr<TensorExpr, typename promote<T, U>::V, Dim0, Dim4,
667 Dim2, Dim3, i, m, k, l>(TensorExpr(a, b));
668 }
669
670 /* A(i,j,k,l)*B(i,m) */
671
672 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
673 int Dim3, int Dim4, char i, char j, char k, char l, char m>
674 class Tensor4_times_Tensor2_0_0
675 {
676 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> iterA;
677 const Tensor2_Expr<B, U, Dim0, Dim4, i, m> iterB;
678
679 template <int Current_Dim0>
680 typename promote<T, U>::V
681 eval(const int N1, const int N2, const int N3, const int N4,
682 const Number<Current_Dim0> &) const
683 {
684 return iterA(Current_Dim0 - 1, N2, N3, N4) * iterB(Current_Dim0 - 1, N1)
685 + eval(N1, N2, N3, N4, Number<Current_Dim0 - 1>());
686 }
687 typename promote<T, U>::V eval(const int N1, const int N2, const int N3,
688 const int N4, const Number<1> &) const
689 {
690 return iterA(0, N2, N3, N4) * iterB(0, N1);
691 }
692
693 public:
694 Tensor4_times_Tensor2_0_0(
695 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a,
696 const Tensor2_Expr<B, U, Dim0, Dim4, i, m> &b)
697 : iterA(a), iterB(b)
698 {}
699 typename promote<T, U>::V
700 operator()(const int N1, const int N2, const int N3, const int N4) const
701 {
702 return eval(N1, N2, N3, N4, Number<Dim0>());
703 }
704 };
705
706 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
707 int Dim3, int Dim4, char i, char j, char k, char l, char m>
708 inline const Tensor4_Expr<
709 const Tensor4_times_Tensor2_0_0<A, B, T, U, Dim0, Dim1, Dim2, Dim3, Dim4,
710 i, j, k, l, m>,
711 typename promote<T, U>::V, Dim1, Dim1, Dim2, Dim3, m, j, k, l>
712 operator*(const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a,
713 const Tensor2_Expr<B, U, Dim0, Dim4, i, m> &b)
714 {
715 typedef const Tensor4_times_Tensor2_0_0<A, B, T, U, Dim0, Dim1, Dim2, Dim3,
716 Dim4, i, j, k, l, m>
717 TensorExpr;
718 return Tensor4_Expr<TensorExpr, typename promote<T, U>::V, Dim4, Dim1,
719 Dim2, Dim3, m, j, k, l>(TensorExpr(a, b));
720 }
721
722 /* B(i,m)*A(i,j,k,l) */
723
724 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
725 int Dim3, int Dim4, char i, char j, char k, char l, char m>
726 inline const Tensor4_Expr<
727 const Tensor4_times_Tensor2_0_0<A, B, T, U, Dim0, Dim1, Dim2, Dim3, Dim4,
728 i, j, k, l, m>,
729 typename promote<T, U>::V, Dim1, Dim1, Dim2, Dim3, m, j, k, l>
730 operator*(const Tensor2_Expr<B, U, Dim0, Dim4, i, m> &b,
731 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a)
732 {
733 typedef const Tensor4_times_Tensor2_0_0<A, B, T, U, Dim0, Dim1, Dim2, Dim3,
734 Dim4, i, j, k, l, m>
735 TensorExpr;
736 return Tensor4_Expr<TensorExpr, typename promote<T, U>::V, Dim4, Dim1,
737 Dim2, Dim3, m, j, k, l>(TensorExpr(a, b));
738 }
739
740 /* A(i,j,k,l)*B(m,i) */
741
742 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
743 int Dim3, int Dim4, char i, char j, char k, char l, char m>
744 class Tensor4_times_Tensor2_0_1
745 {
746 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> iterA;
747 const Tensor2_Expr<B, U, Dim4, Dim0, m, i> iterB;
748
749 template <int Current_Dim0>
750 typename promote<T, U>::V
751 eval(const int N1, const int N2, const int N3, const int N4,
752 const Number<Current_Dim0> &) const
753 {
754 return iterA(Current_Dim0 - 1, N2, N3, N4) * iterB(N1, Current_Dim0 - 1)
755 + eval(N1, N2, N3, N4, Number<Current_Dim0 - 1>());
756 }
757 typename promote<T, U>::V eval(const int N1, const int N2, const int N3,
758 const int N4, const Number<1> &) const
759 {
760 return iterA(0, N2, N3, N4) * iterB(N1, 0);
761 }
762
763 public:
764 Tensor4_times_Tensor2_0_1(
765 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a,
766 const Tensor2_Expr<B, U, Dim4, Dim1, m, i> &b)
767 : iterA(a), iterB(b)
768 {}
769 typename promote<T, U>::V
770 operator()(const int N1, const int N2, const int N3, const int N4) const
771 {
772 return eval(N1, N2, N3, N4, Number<Dim0>());
773 }
774 };
775
776 // FIXME: Template is ambiguous
777 // template<
778 // class A, class B, class T, class U, int Dim0,int Dim1, int Dim2,int Dim3,
779 // int Dim4,char i, char j, char k, char l,char m> inline const Tensor4_Expr
780 // <const
781 // Tensor4_times_Tensor2_0_1<A,B,T,U,Dim0,Dim1,Dim2,Dim3,Dim4,i,j,k,l,m>,typename
782 // promote<T,U>::V,Dim4,Dim1,Dim2,Dim3,m,j,k,l> operator*(
783 // const Tensor4_Expr<A,T,Dim0,Dim1,Dim2,Dim3,i,j,k,l> &a,
784 // const Tensor2_Expr<B,U,Dim4,Dim1,m,i> &b
785 // ) {
786 // typedef const
787 // Tensor4_times_Tensor2_0_1<A,B,T,U,Dim0,Dim1,Dim2,Dim3,Dim4,i,j,k,l,m>
788 // TensorExpr; return Tensor4_Expr<TensorExpr,typename
789 // promote<T,U>::V,Dim4,Dim1,Dim2,Dim3,m,j,k,l> (TensorExpr(a,b));
790 // }
791
792 /* B(i,m)*A(i,j,k,l) */
793
794 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
795 int Dim3, int Dim4, char i, char j, char k, char l, char m>
796 inline const Tensor4_Expr<
797 const Tensor4_times_Tensor2_0_1<A, B, T, U, Dim0, Dim1, Dim2, Dim3, Dim4,
798 i, j, k, l, m>,
799 typename promote<T, U>::V, Dim4, Dim1, Dim2, Dim3, m, j, k, l>
800 operator*(const Tensor2_Expr<B, U, Dim4, Dim1, m, i> &b,
801 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a)
802 {
803 typedef const Tensor4_times_Tensor2_0_1<A, B, T, U, Dim0, Dim1, Dim2, Dim3,
804 Dim4, i, j, k, l, m>
805 TensorExpr;
806 return Tensor4_Expr<TensorExpr, typename promote<T, U>::V, Dim4, Dim1,
807 Dim2, Dim3, m, j, k, l>(TensorExpr(a, b));
808 }
809
810 /* A(i,j,k,l)*B(k,m) */
811
812 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
813 int Dim3, int Dim4, char i, char j, char k, char l, char m>
814 class Tensor4_times_Tensor2_2_0
815 {
816 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> iterA;
817 const Tensor2_Expr<B, U, Dim2, Dim4, k, m> iterB;
818
819 template <int Current_Dim0>
820 typename promote<T, U>::V
821 eval(const int N1, const int N2, const int N3, const int N4,
822 const Number<Current_Dim0> &) const
823 {
824 return iterA(N1, N2, Current_Dim0 - 1, N4) * iterB(Current_Dim0 - 1, N3)
825 + eval(N1, N2, N3, N4, Number<Current_Dim0 - 1>());
826 }
827 typename promote<T, U>::V eval(const int N1, const int N2, const int N3,
828 const int N4, const Number<1> &) const
829 {
830 return iterA(N1, N2, 0, N4) * iterB(0, N3);
831 }
832
833 public:
834 Tensor4_times_Tensor2_2_0(
835 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a,
836 const Tensor2_Expr<B, U, Dim2, Dim4, k, m> &b)
837 : iterA(a), iterB(b)
838 {}
839 typename promote<T, U>::V
840 operator()(const int N1, const int N2, const int N3, const int N4) const
841 {
842 return eval(N1, N2, N3, N4, Number<Dim3>());
843 }
844 };
845
846 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
847 int Dim3, int Dim4, char i, char j, char k, char l, char m>
848 inline const Tensor4_Expr<
849 const Tensor4_times_Tensor2_2_0<A, B, T, U, Dim0, Dim1, Dim2, Dim3, Dim4,
850 i, j, k, l, m>,
851 typename promote<T, U>::V, Dim0, Dim1, Dim4, Dim3, i, j, m, l>
852 operator*(const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a,
853 const Tensor2_Expr<B, U, Dim2, Dim4, k, m> &b)
854 {
855 typedef const Tensor4_times_Tensor2_2_0<A, B, T, U, Dim0, Dim1, Dim2, Dim3,
856 Dim4, i, j, k, l, m>
857 TensorExpr;
858 return Tensor4_Expr<TensorExpr, typename promote<T, U>::V, Dim0, Dim1,
859 Dim4, Dim3, i, j, m, l>(TensorExpr(a, b));
860 }
861
862 /* B(k,m)*A(i,j,k,l) */
863
864 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
865 int Dim3, int Dim4, char i, char j, char k, char l, char m>
866 inline const Tensor4_Expr<
867 const Tensor4_times_Tensor2_2_0<A, B, T, U, Dim0, Dim1, Dim2, Dim3, Dim4,
868 i, j, k, l, m>,
869 typename promote<T, U>::V, Dim0, Dim1, Dim4, Dim3, i, j, m, l>
870 operator*(const Tensor2_Expr<B, U, Dim2, Dim4, k, m> &b,
871 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a)
872 {
873 typedef const Tensor4_times_Tensor2_2_0<A, B, T, U, Dim0, Dim1, Dim2, Dim3,
874 Dim4, i, j, k, l, m>
875 TensorExpr;
876 return Tensor4_Expr<TensorExpr, typename promote<T, U>::V, Dim0, Dim1,
877 Dim4, Dim3, i, j, m, l>(TensorExpr(a, b));
878 }
879
880 /* A(i,j,k,l)*B(m,k) */
881
882 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
883 int Dim3, int Dim4, char i, char j, char k, char l, char m>
884 class Tensor4_times_Tensor2_2_1
885 {
886 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> iterA;
887 const Tensor2_Expr<B, U, Dim4, Dim2, m, k> iterB;
888
889 template <int Current_Dim0>
890 typename promote<T, U>::V
891 eval(const int N1, const int N2, const int N3, const int N4,
892 const Number<Current_Dim0> &) const
893 {
894 return iterA(N1, N2, Current_Dim0 - 1, N4) * iterB(N3, Current_Dim0 - 1)
895 + eval(N1, N2, N3, N4, Number<Current_Dim0 - 1>());
896 }
897 typename promote<T, U>::V eval(const int N1, const int N2, const int N3,
898 const int N4, const Number<1> &) const
899 {
900 return iterA(N1, N2, 0, N4) * iterB(N3, 0);
901 }
902
903 public:
904 Tensor4_times_Tensor2_2_1(
905 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a,
906 const Tensor2_Expr<B, U, Dim4, Dim2, m, k> &b)
907 : iterA(a), iterB(b)
908 {}
909 typename promote<T, U>::V
910 operator()(const int N1, const int N2, const int N3, const int N4) const
911 {
912 return eval(N1, N2, N3, N4, Number<Dim3>());
913 }
914 };
915
916 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
917 int Dim3, int Dim4, char i, char j, char k, char l, char m>
918 inline const Tensor4_Expr<
919 const Tensor4_times_Tensor2_2_1<A, B, T, U, Dim0, Dim1, Dim2, Dim3, Dim4,
920 i, j, k, l, m>,
921 typename promote<T, U>::V, Dim0, Dim1, Dim4, Dim3, i, j, m, l>
922 operator*(const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a,
923 const Tensor2_Expr<B, U, Dim4, Dim2, m, k> &b)
924 {
925 typedef const Tensor4_times_Tensor2_2_1<A, B, T, U, Dim0, Dim1, Dim2, Dim3,
926 Dim4, i, j, k, l, m>
927 TensorExpr;
928 return Tensor4_Expr<TensorExpr, typename promote<T, U>::V, Dim0, Dim1,
929 Dim4, Dim3, i, j, m, l>(TensorExpr(a, b));
930 }
931
932 /* B(m,k)*A(i,j,k,l) */
933
934 template <class A, class B, class T, class U, int Dim0, int Dim1, int Dim2,
935 int Dim3, int Dim4, char i, char j, char k, char l, char m>
936 inline const Tensor4_Expr<
937 const Tensor4_times_Tensor2_2_1<A, B, T, U, Dim0, Dim1, Dim2, Dim3, Dim4,
938 i, j, k, l, m>,
939 typename promote<T, U>::V, Dim0, Dim1, Dim4, Dim3, i, j, m, l>
940 operator*(const Tensor2_Expr<B, U, Dim4, Dim2, m, k> &b,
941 const Tensor4_Expr<A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l> &a)
942 {
943 typedef const Tensor4_times_Tensor2_2_1<A, B, T, U, Dim0, Dim1, Dim2, Dim3,
944 Dim4, i, j, k, l, m>
945 TensorExpr;
946 return Tensor4_Expr<TensorExpr, typename promote<T, U>::V, Dim0, Dim1,
947 Dim4, Dim3, i, j, m, l>(TensorExpr(a, b));
948 }
949
950} // namespace FTensor
N2
static Number< 2 > N2
Definition: BasicFeTools.hpp:101
N1
static Number< 1 > N1
Definition: BasicFeTools.hpp:100
a
constexpr double a
Definition: approx_sphere.cpp:30
FTensor::Tensor2_Expr
Definition: Tensor2_Expr.hpp:27
FTensor::Tensor4_Expr
Definition: Tensor4_Expr.hpp:26
FTensor::Tensor4_times_Tensor2_03
Definition: Tensor4_times_Tensor2.hpp:170
FTensor::Tensor4_times_Tensor2_03::operator()
promote< T, U >::V operator()(const int N1, const int N2) const
Definition: Tensor4_times_Tensor2.hpp:199
FTensor::Tensor4_times_Tensor2_03::eval
promote< T, U >::V eval(const int N1, const int N2, const Number< 1 > &, const Number< 1 > &) const
Definition: Tensor4_times_Tensor2.hpp:192
FTensor::Tensor4_times_Tensor2_03::eval
promote< T, U >::V eval(const int N1, const int N2, const Number< 1 > &, const Number< Current_Dim1 > &) const
Definition: Tensor4_times_Tensor2.hpp:186
FTensor::Tensor4_times_Tensor2_03::Tensor4_times_Tensor2_03
Tensor4_times_Tensor2_03(const Tensor4_Expr< A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l > &a, const Tensor2_Expr< B, U, Dim0, Dim3, i, l > &b)
Definition: Tensor4_times_Tensor2.hpp:204
FTensor::Tensor4_times_Tensor2_03::iterB
const Tensor2_Expr< B, U, Dim0, Dim3, i, l > iterB
Definition: Tensor4_times_Tensor2.hpp:172
FTensor::Tensor4_times_Tensor2_03::iterA
const Tensor4_Expr< A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l > iterA
Definition: Tensor4_times_Tensor2.hpp:171
FTensor::Tensor4_times_Tensor2_03::eval
promote< T, U >::V eval(const int N1, const int N2, const Number< Current_Dim0 > &, const Number< Current_Dim1 > &) const
Definition: Tensor4_times_Tensor2.hpp:176
FTensor::Tensor4_times_Tensor2_0_0
Definition: Tensor4_times_Tensor2.hpp:675
FTensor::Tensor4_times_Tensor2_0_0::operator()
promote< T, U >::V operator()(const int N1, const int N2, const int N3, const int N4) const
Definition: Tensor4_times_Tensor2.hpp:700
FTensor::Tensor4_times_Tensor2_0_0::eval
promote< T, U >::V eval(const int N1, const int N2, const int N3, const int N4, const Number< 1 > &) const
Definition: Tensor4_times_Tensor2.hpp:687
FTensor::Tensor4_times_Tensor2_0_0::eval
promote< T, U >::V eval(const int N1, const int N2, const int N3, const int N4, const Number< Current_Dim0 > &) const
Definition: Tensor4_times_Tensor2.hpp:681
FTensor::Tensor4_times_Tensor2_0_0::Tensor4_times_Tensor2_0_0
Tensor4_times_Tensor2_0_0(const Tensor4_Expr< A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l > &a, const Tensor2_Expr< B, U, Dim0, Dim4, i, m > &b)
Definition: Tensor4_times_Tensor2.hpp:694
FTensor::Tensor4_times_Tensor2_0_0::iterA
const Tensor4_Expr< A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l > iterA
Definition: Tensor4_times_Tensor2.hpp:676
FTensor::Tensor4_times_Tensor2_0_0::iterB
const Tensor2_Expr< B, U, Dim0, Dim4, i, m > iterB
Definition: Tensor4_times_Tensor2.hpp:677
FTensor::Tensor4_times_Tensor2_0_1
Definition: Tensor4_times_Tensor2.hpp:745
FTensor::Tensor4_times_Tensor2_0_1::Tensor4_times_Tensor2_0_1
Tensor4_times_Tensor2_0_1(const Tensor4_Expr< A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l > &a, const Tensor2_Expr< B, U, Dim4, Dim1, m, i > &b)
Definition: Tensor4_times_Tensor2.hpp:764
FTensor::Tensor4_times_Tensor2_0_1::operator()
promote< T, U >::V operator()(const int N1, const int N2, const int N3, const int N4) const
Definition: Tensor4_times_Tensor2.hpp:770
FTensor::Tensor4_times_Tensor2_0_1::eval
promote< T, U >::V eval(const int N1, const int N2, const int N3, const int N4, const Number< 1 > &) const
Definition: Tensor4_times_Tensor2.hpp:757
FTensor::Tensor4_times_Tensor2_0_1::eval
promote< T, U >::V eval(const int N1, const int N2, const int N3, const int N4, const Number< Current_Dim0 > &) const
Definition: Tensor4_times_Tensor2.hpp:751
FTensor::Tensor4_times_Tensor2_0_1::iterA
const Tensor4_Expr< A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l > iterA
Definition: Tensor4_times_Tensor2.hpp:746
FTensor::Tensor4_times_Tensor2_0_1::iterB
const Tensor2_Expr< B, U, Dim4, Dim0, m, i > iterB
Definition: Tensor4_times_Tensor2.hpp:747
FTensor::Tensor4_times_Tensor2_13
Definition: Tensor4_times_Tensor2.hpp:326
FTensor::Tensor4_times_Tensor2_13::operator()
promote< T, U >::V operator()(const int N1, const int N2) const
Definition: Tensor4_times_Tensor2.hpp:355
FTensor::Tensor4_times_Tensor2_13::iterA
const Tensor4_Expr< A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l > iterA
Definition: Tensor4_times_Tensor2.hpp:327
FTensor::Tensor4_times_Tensor2_13::iterB
const Tensor2_Expr< B, U, Dim2, Dim3, j, l > iterB
Definition: Tensor4_times_Tensor2.hpp:328
FTensor::Tensor4_times_Tensor2_13::eval
promote< T, U >::V eval(const int N1, const int N2, const Number< Current_Dim0 > &, const Number< Current_Dim1 > &) const
Definition: Tensor4_times_Tensor2.hpp:332
FTensor::Tensor4_times_Tensor2_13::Tensor4_times_Tensor2_13
Tensor4_times_Tensor2_13(const Tensor4_Expr< A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l > &a, const Tensor2_Expr< B, U, Dim2, Dim3, j, l > &b)
Definition: Tensor4_times_Tensor2.hpp:360
FTensor::Tensor4_times_Tensor2_13::eval
promote< T, U >::V eval(const int N1, const int N2, const Number< 1 > &, const Number< 1 > &) const
Definition: Tensor4_times_Tensor2.hpp:348
FTensor::Tensor4_times_Tensor2_13::eval
promote< T, U >::V eval(const int N1, const int N2, const Number< 1 > &, const Number< Current_Dim1 > &) const
Definition: Tensor4_times_Tensor2.hpp:342
FTensor::Tensor4_times_Tensor2_1_0
Definition: Tensor4_times_Tensor2.hpp:535
FTensor::Tensor4_times_Tensor2_1_0::iterA
const Tensor4_Expr< A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l > iterA
Definition: Tensor4_times_Tensor2.hpp:536
FTensor::Tensor4_times_Tensor2_1_0::eval
promote< T, U >::V eval(const int N1, const int N2, const int N3, const int N4, const Number< 1 > &) const
Definition: Tensor4_times_Tensor2.hpp:547
FTensor::Tensor4_times_Tensor2_1_0::operator()
promote< T, U >::V operator()(const int N1, const int N2, const int N3, const int N4) const
Definition: Tensor4_times_Tensor2.hpp:560
FTensor::Tensor4_times_Tensor2_1_0::iterB
const Tensor2_Expr< B, U, Dim1, Dim4, j, m > iterB
Definition: Tensor4_times_Tensor2.hpp:537
FTensor::Tensor4_times_Tensor2_1_0::Tensor4_times_Tensor2_1_0
Tensor4_times_Tensor2_1_0(const Tensor4_Expr< A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l > &a, const Tensor2_Expr< B, U, Dim1, Dim4, j, m > &b)
Definition: Tensor4_times_Tensor2.hpp:554
FTensor::Tensor4_times_Tensor2_1_0::eval
promote< T, U >::V eval(const int N1, const int N2, const int N3, const int N4, const Number< Current_Dim0 > &) const
Definition: Tensor4_times_Tensor2.hpp:541
FTensor::Tensor4_times_Tensor2_1_1
Definition: Tensor4_times_Tensor2.hpp:605
FTensor::Tensor4_times_Tensor2_1_1::eval
promote< T, U >::V eval(const int N1, const int N2, const int N3, const int N4, const Number< 1 > &) const
Definition: Tensor4_times_Tensor2.hpp:617
FTensor::Tensor4_times_Tensor2_1_1::iterA
const Tensor4_Expr< A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l > iterA
Definition: Tensor4_times_Tensor2.hpp:606
FTensor::Tensor4_times_Tensor2_1_1::iterB
const Tensor2_Expr< B, U, Dim4, Dim1, m, j > iterB
Definition: Tensor4_times_Tensor2.hpp:607
FTensor::Tensor4_times_Tensor2_1_1::eval
promote< T, U >::V eval(const int N1, const int N2, const int N3, const int N4, const Number< Current_Dim0 > &) const
Definition: Tensor4_times_Tensor2.hpp:611
FTensor::Tensor4_times_Tensor2_1_1::Tensor4_times_Tensor2_1_1
Tensor4_times_Tensor2_1_1(const Tensor4_Expr< A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l > &a, const Tensor2_Expr< B, U, Dim4, Dim1, m, j > &b)
Definition: Tensor4_times_Tensor2.hpp:624
FTensor::Tensor4_times_Tensor2_1_1::operator()
promote< T, U >::V operator()(const int N1, const int N2, const int N3, const int N4) const
Definition: Tensor4_times_Tensor2.hpp:630
FTensor::Tensor4_times_Tensor2_23
Definition: Tensor4_times_Tensor2.hpp:14
FTensor::Tensor4_times_Tensor2_23::eval
promote< T, U >::V eval(const int N1, const int N2, const Number< 1 > &, const Number< 1 > &) const
Definition: Tensor4_times_Tensor2.hpp:36
FTensor::Tensor4_times_Tensor2_23::eval
promote< T, U >::V eval(const int N1, const int N2, const Number< 1 > &, const Number< Current_Dim1 > &) const
Definition: Tensor4_times_Tensor2.hpp:30
FTensor::Tensor4_times_Tensor2_23::eval
promote< T, U >::V eval(const int N1, const int N2, const Number< Current_Dim0 > &, const Number< Current_Dim1 > &) const
Definition: Tensor4_times_Tensor2.hpp:20
FTensor::Tensor4_times_Tensor2_23::Tensor4_times_Tensor2_23
Tensor4_times_Tensor2_23(const Tensor4_Expr< A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l > &a, const Tensor2_Expr< B, U, Dim2, Dim3, k, l > &b)
Definition: Tensor4_times_Tensor2.hpp:48
FTensor::Tensor4_times_Tensor2_23::operator()
promote< T, U >::V operator()(const int N1, const int N2) const
Definition: Tensor4_times_Tensor2.hpp:43
FTensor::Tensor4_times_Tensor2_23::iterB
const Tensor2_Expr< B, U, Dim2, Dim3, k, l > iterB
Definition: Tensor4_times_Tensor2.hpp:16
FTensor::Tensor4_times_Tensor2_23::iterA
const Tensor4_Expr< A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l > iterA
Definition: Tensor4_times_Tensor2.hpp:15
FTensor::Tensor4_times_Tensor2_2_0
Definition: Tensor4_times_Tensor2.hpp:815
FTensor::Tensor4_times_Tensor2_2_0::operator()
promote< T, U >::V operator()(const int N1, const int N2, const int N3, const int N4) const
Definition: Tensor4_times_Tensor2.hpp:840
FTensor::Tensor4_times_Tensor2_2_0::iterA
const Tensor4_Expr< A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l > iterA
Definition: Tensor4_times_Tensor2.hpp:816
FTensor::Tensor4_times_Tensor2_2_0::iterB
const Tensor2_Expr< B, U, Dim2, Dim4, k, m > iterB
Definition: Tensor4_times_Tensor2.hpp:817
FTensor::Tensor4_times_Tensor2_2_0::eval
promote< T, U >::V eval(const int N1, const int N2, const int N3, const int N4, const Number< 1 > &) const
Definition: Tensor4_times_Tensor2.hpp:827
FTensor::Tensor4_times_Tensor2_2_0::Tensor4_times_Tensor2_2_0
Tensor4_times_Tensor2_2_0(const Tensor4_Expr< A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l > &a, const Tensor2_Expr< B, U, Dim2, Dim4, k, m > &b)
Definition: Tensor4_times_Tensor2.hpp:834
FTensor::Tensor4_times_Tensor2_2_0::eval
promote< T, U >::V eval(const int N1, const int N2, const int N3, const int N4, const Number< Current_Dim0 > &) const
Definition: Tensor4_times_Tensor2.hpp:821
FTensor::Tensor4_times_Tensor2_2_1
Definition: Tensor4_times_Tensor2.hpp:885
FTensor::Tensor4_times_Tensor2_2_1::Tensor4_times_Tensor2_2_1
Tensor4_times_Tensor2_2_1(const Tensor4_Expr< A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l > &a, const Tensor2_Expr< B, U, Dim4, Dim2, m, k > &b)
Definition: Tensor4_times_Tensor2.hpp:904
FTensor::Tensor4_times_Tensor2_2_1::operator()
promote< T, U >::V operator()(const int N1, const int N2, const int N3, const int N4) const
Definition: Tensor4_times_Tensor2.hpp:910
FTensor::Tensor4_times_Tensor2_2_1::eval
promote< T, U >::V eval(const int N1, const int N2, const int N3, const int N4, const Number< Current_Dim0 > &) const
Definition: Tensor4_times_Tensor2.hpp:891
FTensor::Tensor4_times_Tensor2_2_1::eval
promote< T, U >::V eval(const int N1, const int N2, const int N3, const int N4, const Number< 1 > &) const
Definition: Tensor4_times_Tensor2.hpp:897
FTensor::Tensor4_times_Tensor2_2_1::iterB
const Tensor2_Expr< B, U, Dim4, Dim2, m, k > iterB
Definition: Tensor4_times_Tensor2.hpp:887
FTensor::Tensor4_times_Tensor2_2_1::iterA
const Tensor4_Expr< A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l > iterA
Definition: Tensor4_times_Tensor2.hpp:886
FTensor::Tensor4_times_Tensor2_30
Definition: Tensor4_times_Tensor2.hpp:248
FTensor::Tensor4_times_Tensor2_30::Tensor4_times_Tensor2_30
Tensor4_times_Tensor2_30(const Tensor4_Expr< A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l > &a, const Tensor2_Expr< B, U, Dim3, Dim0, l, i > &b)
Definition: Tensor4_times_Tensor2.hpp:282
FTensor::Tensor4_times_Tensor2_30::operator()
promote< T, U >::V operator()(const int N1, const int N2) const
Definition: Tensor4_times_Tensor2.hpp:277
FTensor::Tensor4_times_Tensor2_30::eval
promote< T, U >::V eval(const int N1, const int N2, const Number< 1 > &, const Number< Current_Dim1 > &) const
Definition: Tensor4_times_Tensor2.hpp:264
FTensor::Tensor4_times_Tensor2_30::iterA
const Tensor4_Expr< A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l > iterA
Definition: Tensor4_times_Tensor2.hpp:249
FTensor::Tensor4_times_Tensor2_30::eval
promote< T, U >::V eval(const int N1, const int N2, const Number< 1 > &, const Number< 1 > &) const
Definition: Tensor4_times_Tensor2.hpp:270
FTensor::Tensor4_times_Tensor2_30::eval
promote< T, U >::V eval(const int N1, const int N2, const Number< Current_Dim0 > &, const Number< Current_Dim1 > &) const
Definition: Tensor4_times_Tensor2.hpp:254
FTensor::Tensor4_times_Tensor2_30::iterB
const Tensor2_Expr< B, U, Dim3, Dim0, l, i > iterB
Definition: Tensor4_times_Tensor2.hpp:250
FTensor::Tensor4_times_Tensor2_32
Definition: Tensor4_times_Tensor2.hpp:92
FTensor::Tensor4_times_Tensor2_32::iterB
const Tensor2_Expr< B, U, Dim3, Dim2, l, k > iterB
Definition: Tensor4_times_Tensor2.hpp:94
FTensor::Tensor4_times_Tensor2_32::Tensor4_times_Tensor2_32
Tensor4_times_Tensor2_32(const Tensor4_Expr< A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l > &a, const Tensor2_Expr< B, U, Dim3, Dim2, l, k > &b)
Definition: Tensor4_times_Tensor2.hpp:126
FTensor::Tensor4_times_Tensor2_32::eval
promote< T, U >::V eval(const int N1, const int N2, const Number< 1 > &, const Number< Current_Dim1 > &) const
Definition: Tensor4_times_Tensor2.hpp:108
FTensor::Tensor4_times_Tensor2_32::iterA
const Tensor4_Expr< A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l > iterA
Definition: Tensor4_times_Tensor2.hpp:93
FTensor::Tensor4_times_Tensor2_32::operator()
promote< T, U >::V operator()(const int N1, const int N2) const
Definition: Tensor4_times_Tensor2.hpp:121
FTensor::Tensor4_times_Tensor2_32::eval
promote< T, U >::V eval(const int N1, const int N2, const Number< Current_Dim0 > &, const Number< Current_Dim1 > &) const
Definition: Tensor4_times_Tensor2.hpp:98
FTensor::Tensor4_times_Tensor2_32::eval
promote< T, U >::V eval(const int N1, const int N2, const Number< 1 > &, const Number< 1 > &) const
Definition: Tensor4_times_Tensor2.hpp:114
FTensor::Tensor4_times_Tensor2_3_0
Definition: Tensor4_times_Tensor2.hpp:465
FTensor::Tensor4_times_Tensor2_3_0::Tensor4_times_Tensor2_3_0
Tensor4_times_Tensor2_3_0(const Tensor4_Expr< A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l > &a, const Tensor2_Expr< B, U, Dim4, Dim3, m, l > &b)
Definition: Tensor4_times_Tensor2.hpp:484
FTensor::Tensor4_times_Tensor2_3_0::operator()
promote< T, U >::V operator()(const int N1, const int N2, const int N3, const int N4) const
Definition: Tensor4_times_Tensor2.hpp:490
FTensor::Tensor4_times_Tensor2_3_0::eval
promote< T, U >::V eval(const int N1, const int N2, const int N3, const int N4, const Number< 1 > &) const
Definition: Tensor4_times_Tensor2.hpp:477
FTensor::Tensor4_times_Tensor2_3_0::iterB
const Tensor2_Expr< B, U, Dim4, Dim3, m, l > iterB
Definition: Tensor4_times_Tensor2.hpp:467
FTensor::Tensor4_times_Tensor2_3_0::eval
promote< T, U >::V eval(const int N1, const int N2, const int N3, const int N4, const Number< Current_Dim0 > &) const
Definition: Tensor4_times_Tensor2.hpp:471
FTensor::Tensor4_times_Tensor2_3_0::iterA
const Tensor4_Expr< A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l > iterA
Definition: Tensor4_times_Tensor2.hpp:466
FTensor::Tensor4_times_Tensor2_3_1
Definition: Tensor4_times_Tensor2.hpp:395
FTensor::Tensor4_times_Tensor2_3_1::operator()
promote< T, U >::V operator()(const int N1, const int N2, const int N3, const int N4) const
Definition: Tensor4_times_Tensor2.hpp:420
FTensor::Tensor4_times_Tensor2_3_1::iterA
const Tensor4_Expr< A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l > iterA
Definition: Tensor4_times_Tensor2.hpp:396
FTensor::Tensor4_times_Tensor2_3_1::iterB
const Tensor2_Expr< B, U, Dim3, Dim4, l, m > iterB
Definition: Tensor4_times_Tensor2.hpp:397
FTensor::Tensor4_times_Tensor2_3_1::eval
promote< T, U >::V eval(const int N1, const int N2, const int N3, const int N4, const Number< Current_Dim0 > &) const
Definition: Tensor4_times_Tensor2.hpp:401
FTensor::Tensor4_times_Tensor2_3_1::Tensor4_times_Tensor2_3_1
Tensor4_times_Tensor2_3_1(const Tensor4_Expr< A, T, Dim0, Dim1, Dim2, Dim3, i, j, k, l > &a, const Tensor2_Expr< B, U, Dim3, Dim4, l, m > &b)
Definition: Tensor4_times_Tensor2.hpp:414
FTensor::Tensor4_times_Tensor2_3_1::eval
promote< T, U >::V eval(const int N1, const int N2, const int N3, const int N4, const Number< 1 > &) const
Definition: Tensor4_times_Tensor2.hpp:407
m
FTensor::Index< 'm', SPACE_DIM > m
Definition: fluid_structure_eigenproblem.cpp:65
i
FTensor::Index< 'i', SPACE_DIM > i
Definition: hcurl_divergence_operator_2d.cpp:27
l
FTensor::Index< 'l', 3 > l
Definition: matrix_function.cpp:21
j
FTensor::Index< 'j', 3 > j
Definition: matrix_function.cpp:19
k
FTensor::Index< 'k', 3 > k
Definition: matrix_function.cpp:20
T
const double T
Definition: mixed_reac_diff.cpp:158
FTensor
Tensors class implemented by Walter Landry.
Definition: FTensor.hpp:51
FTensor::operator*
promote< T, U >::V operator*(const Ddg_Expr< A, T, Dim, Dim, i, j, k, l > &a, const Ddg_Expr< B, U, Dim, Dim, i, k, j, l > &b)
Definition: Ddg_times_Ddg.hpp:79
A
constexpr AssemblyType A
Definition: operators_tests.cpp:30
Generated by Doxygen 1.9.5 and hosted at University of Glasgow