v0.14.0
Tensor3_antisymmetric_times_Tensor3.hpp
Go to the documentation of this file.
1 /* Fully contracts a Tensor3_antisymmetric with a Tensor3, yielding a
2  typename promote<T,U>::V. I removed the identically zero components of
3  Tensor3_antisymmetric.*/
4 
5 #pragma once
6 
7 namespace FTensor
8 {
9  /* A(i,j,k)*B(i,j,k) */
10 
11  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
12  char j, char k, int Current_Dim0, int Current_Dim1,
13  int Current_Dim2>
14  typename promote<T, U>::V T3as_times_T3_012(
15  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
16  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, i, j, k> &b,
17  const Number<Current_Dim0> &, const Number<Current_Dim1> &,
18  const Number<Current_Dim2> &)
19  {
20  return a(Current_Dim0 - 1, Current_Dim1 - 1, Current_Dim2 - 1)
21  * b(Current_Dim0 - 1, Current_Dim1 - 1, Current_Dim2 - 1)
22  + T3as_times_T3_012(a, b, Number<Current_Dim0>(),
23  Number<Current_Dim1 - 1>(),
24  Number<Current_Dim2>());
25  }
26 
27  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
28  char j, char k, int Current_Dim0, int Current_Dim2>
29  typename promote<T, U>::V T3as_times_T3_012(
30  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
31  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, i, j, k> &b,
32  const Number<Current_Dim0> &, const Number<1> &,
33  const Number<Current_Dim2> &)
34  {
35  return a(Current_Dim0 - 1, 0, Current_Dim2 - 1)
36  * b(Current_Dim0 - 1, 0, Current_Dim2 - 1)
37  + T3as_times_T3_012(a, b, Number<Current_Dim0>(), Number<Dim12>(),
38  Number<Current_Dim2 - 1>());
39  }
40 
41  /* A special case for when the last two indices are equal. */
42 
43  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
44  char j, char k, int Current_Dim0, int Current_Dim2>
45  typename promote<T, U>::V T3as_times_T3_012(
46  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
47  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, i, j, k> &b,
48  const Number<Current_Dim0> &, const Number<Current_Dim2> &,
49  const Number<Current_Dim2> &)
50  {
51  return T3as_times_T3_012(a, b, Number<Current_Dim0>(),
52  Number<Current_Dim2 - 1>(),
53  Number<Current_Dim2>());
54  }
55 
56  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
57  char j, char k, int Current_Dim0>
58  typename promote<T, U>::V T3as_times_T3_012(
59  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
60  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, i, j, k> &b,
61  const Number<Current_Dim0> &, const Number<2> &, const Number<1> &)
62  {
63  return a(Current_Dim0 - 1, 1, 0) * b(Current_Dim0 - 1, 1, 0)
64  + T3as_times_T3_012(a, b, Number<Current_Dim0 - 1>(),
65  Number<Dim12 - 1>(), Number<Dim12>());
66  }
67 
68  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
69  char j, char k>
70  typename promote<T, U>::V T3as_times_T3_012(
71  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
72  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, i, j, k> &b,
73  const Number<1> &, const Number<2> &, const Number<1> &)
74  {
75  return a(0, 1, 0) * b(0, 1, 0);
76  }
77 
78  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
79  char j, char k>
80  typename promote<T, U>::V
81  operator*(const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
82  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, i, j, k> &b)
83  {
84  return T3as_times_T3_012(a, b, Number<Dim0>(), Number<Dim12 - 1>(),
85  Number<Dim12>());
86  }
87 
88  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
89  char j, char k>
90  typename promote<T, U>::V
91  operator*(const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, i, j, k> &b,
92  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a)
93  {
94  return a * b;
95  }
96 
97  /* A(i,j,k)*B(k,i,j) */
98 
99  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
100  char j, char k, int Current_Dim0, int Current_Dim1,
101  int Current_Dim2>
102  typename promote<T, U>::V T3as_times_T3_201(
103  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
104  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, k, i, j> &b,
105  const Number<Current_Dim0> &, const Number<Current_Dim1> &,
106  const Number<Current_Dim2> &)
107  {
108  return a(Current_Dim0 - 1, Current_Dim1 - 1, Current_Dim2 - 1)
109  * b(Current_Dim2 - 1, Current_Dim0 - 1, Current_Dim1 - 1)
110  + T3as_times_T3_201(a, b, Number<Current_Dim0>(),
111  Number<Current_Dim1 - 1>(),
112  Number<Current_Dim2>());
113  }
114 
115  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
116  char j, char k, int Current_Dim0, int Current_Dim2>
117  typename promote<T, U>::V T3as_times_T3_201(
118  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
119  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, k, i, j> &b,
120  const Number<Current_Dim0> &, const Number<1> &,
121  const Number<Current_Dim2> &)
122  {
123  return a(Current_Dim0 - 1, 0, Current_Dim2 - 1)
124  * b(Current_Dim2 - 1, Current_Dim0 - 1, 0)
125  + T3as_times_T3_201(a, b, Number<Current_Dim0>(), Number<Dim12>(),
126  Number<Current_Dim2 - 1>());
127  }
128 
129  /* A special case for when the last two indices are equal. */
130 
131  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
132  char j, char k, int Current_Dim0, int Current_Dim2>
133  typename promote<T, U>::V T3as_times_T3_201(
134  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
135  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, k, i, j> &b,
136  const Number<Current_Dim0> &, const Number<Current_Dim2> &,
137  const Number<Current_Dim2> &)
138  {
139  return T3as_times_T3_201(a, b, Number<Current_Dim0>(),
140  Number<Current_Dim2 - 1>(),
141  Number<Current_Dim2>());
142  }
143 
144  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
145  char j, char k, int Current_Dim0>
146  typename promote<T, U>::V T3as_times_T3_201(
147  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
148  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, k, i, j> &b,
149  const Number<Current_Dim0> &, const Number<2> &, const Number<1> &)
150  {
151  return a(Current_Dim0 - 1, 1, 0) * b(0, Current_Dim0 - 1, 1)
152  + T3as_times_T3_201(a, b, Number<Current_Dim0 - 1>(),
153  Number<Dim12 - 1>(), Number<Dim12>());
154  }
155 
156  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
157  char j, char k>
158  typename promote<T, U>::V T3as_times_T3_201(
159  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
160  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, k, i, j> &b,
161  const Number<1> &, const Number<2> &, const Number<1> &)
162  {
163  return a(0, 1, 0) * b(0, 0, 1);
164  }
165 
166  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
167  char j, char k>
168  typename promote<T, U>::V
169  operator*(const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
170  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, k, i, j> &b)
171  {
172  return T3as_times_T3_201(a, b, Number<Dim0>(), Number<Dim12 - 1>(),
173  Number<Dim12>());
174  }
175 
176  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
177  char j, char k>
178  typename promote<T, U>::V
179  operator*(const Tensor3_Expr<B, U, Dim12, Dim0, Dim12, k, i, j> &b,
180  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a)
181  {
182  return a * b;
183  }
184 
185  /* A(i,j,k)*B(j,k,i) */
186 
187  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
188  char j, char k, int Current_Dim0, int Current_Dim1,
189  int Current_Dim2>
190  typename promote<T, U>::V T3as_times_T3_120(
191  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
192  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, j, k, i> &b,
193  const Number<Current_Dim0> &, const Number<Current_Dim1> &,
194  const Number<Current_Dim2> &)
195  {
196  return a(Current_Dim0 - 1, Current_Dim1 - 1, Current_Dim2 - 1)
197  * b(Current_Dim1 - 1, Current_Dim2 - 1, Current_Dim0 - 1)
198  + T3as_times_T3_120(a, b, Number<Current_Dim0>(),
199  Number<Current_Dim1 - 1>(),
200  Number<Current_Dim2>());
201  }
202 
203  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
204  char j, char k, int Current_Dim0, int Current_Dim2>
205  typename promote<T, U>::V T3as_times_T3_120(
206  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
207  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, j, k, i> &b,
208  const Number<Current_Dim0> &, const Number<1> &,
209  const Number<Current_Dim2> &)
210  {
211  return a(Current_Dim0 - 1, 0, Current_Dim2 - 1)
212  * b(0, Current_Dim2 - 1, Current_Dim0 - 1)
213  + T3as_times_T3_120(a, b, Number<Current_Dim0>(), Number<Dim12>(),
214  Number<Current_Dim2 - 1>());
215  }
216 
217  /* A special case for when the last two indices are equal. */
218 
219  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
220  char j, char k, int Current_Dim0, int Current_Dim2>
221  typename promote<T, U>::V T3as_times_T3_120(
222  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
223  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, j, k, i> &b,
224  const Number<Current_Dim0> &, const Number<Current_Dim2> &,
225  const Number<Current_Dim2> &)
226  {
227  return T3as_times_T3_120(a, b, Number<Current_Dim0>(),
228  Number<Current_Dim2 - 1>(),
229  Number<Current_Dim2>());
230  }
231 
232  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
233  char j, char k, int Current_Dim0>
234  typename promote<T, U>::V T3as_times_T3_120(
235  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
236  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, j, k, i> &b,
237  const Number<Current_Dim0> &, const Number<2> &, const Number<1> &)
238  {
239  return a(Current_Dim0 - 1, 1, 0) * b(1, 0, Current_Dim0 - 1)
240  + T3as_times_T3_120(a, b, Number<Current_Dim0 - 1>(),
241  Number<Dim12 - 1>(), Number<Dim12>());
242  }
243 
244  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
245  char j, char k>
246  typename promote<T, U>::V T3as_times_T3_120(
247  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
248  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, j, k, i> &b,
249  const Number<1> &, const Number<2> &, const Number<1> &)
250  {
251  return a(0, 1, 0) * b(1, 0, 0);
252  }
253 
254  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
255  char j, char k>
256  typename promote<T, U>::V
257  operator*(const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
258  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, j, k, i> &b)
259  {
260  return T3as_times_T3_120(a, b, Number<Dim0>(), Number<Dim12 - 1>(),
261  Number<Dim12>());
262  }
263 
264  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
265  char j, char k>
266  typename promote<T, U>::V
267  operator*(const Tensor3_Expr<B, U, Dim12, Dim12, Dim0, j, k, i> &b,
268  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a)
269  {
270  return a * b;
271  }
272 
273  /* A(i,j,k)*B(j,i,k) */
274 
275  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
276  char j, char k, int Current_Dim0, int Current_Dim1,
277  int Current_Dim2>
278  typename promote<T, U>::V T3as_times_T3_102(
279  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
280  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, j, i, k> &b,
281  const Number<Current_Dim0> &, const Number<Current_Dim1> &,
282  const Number<Current_Dim2> &)
283  {
284  return a(Current_Dim0 - 1, Current_Dim1 - 1, Current_Dim2 - 1)
285  * b(Current_Dim1 - 1, Current_Dim0 - 1, Current_Dim2 - 1)
286  + T3as_times_T3_102(a, b, Number<Current_Dim0>(),
287  Number<Current_Dim1 - 1>(),
288  Number<Current_Dim2>());
289  }
290 
291  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
292  char j, char k, int Current_Dim0, int Current_Dim2>
293  typename promote<T, U>::V T3as_times_T3_102(
294  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
295  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, j, i, k> &b,
296  const Number<Current_Dim0> &, const Number<1> &,
297  const Number<Current_Dim2> &)
298  {
299  return a(Current_Dim0 - 1, 0, Current_Dim2 - 1)
300  * b(0, Current_Dim0 - 1, Current_Dim2 - 1)
301  + T3as_times_T3_102(a, b, Number<Current_Dim0>(), Number<Dim12>(),
302  Number<Current_Dim2 - 1>());
303  }
304 
305  /* A special case for when the last two indices are equal. */
306 
307  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
308  char j, char k, int Current_Dim0, int Current_Dim2>
309  typename promote<T, U>::V T3as_times_T3_102(
310  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
311  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, j, i, k> &b,
312  const Number<Current_Dim0> &, const Number<Current_Dim2> &,
313  const Number<Current_Dim2> &)
314  {
315  return T3as_times_T3_102(a, b, Number<Current_Dim0>(),
316  Number<Current_Dim2 - 1>(),
317  Number<Current_Dim2>());
318  }
319 
320  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
321  char j, char k, int Current_Dim0>
322  typename promote<T, U>::V T3as_times_T3_102(
323  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
324  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, j, i, k> &b,
325  const Number<Current_Dim0> &, const Number<2> &, const Number<1> &)
326  {
327  return a(Current_Dim0 - 1, 1, 0) * b(1, Current_Dim0 - 1, 0)
328  + T3as_times_T3_102(a, b, Number<Current_Dim0 - 1>(),
329  Number<Dim12 - 1>(), Number<Dim12>());
330  }
331 
332  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
333  char j, char k>
334  typename promote<T, U>::V T3as_times_T3_102(
335  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
336  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, j, i, k> &b,
337  const Number<1> &, const Number<2> &, const Number<1> &)
338  {
339  return a(0, 1, 0) * b(1, 0, 0);
340  }
341 
342  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
343  char j, char k>
344  typename promote<T, U>::V
345  operator*(const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
346  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, j, i, k> &b)
347  {
348  return T3as_times_T3_102(a, b, Number<Dim0>(), Number<Dim12 - 1>(),
349  Number<Dim12>());
350  }
351 
352  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
353  char j, char k>
354  typename promote<T, U>::V
355  operator*(const Tensor3_Expr<B, U, Dim12, Dim0, Dim12, j, i, k> &b,
356  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a)
357  {
358  return a * b;
359  }
360 
361  /* A(i,j,k)*B(k,j,i) */
362 
363  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
364  char j, char k, int Current_Dim0, int Current_Dim1,
365  int Current_Dim2>
366  typename promote<T, U>::V T3as_times_T3_210(
367  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
368  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, k, j, i> &b,
369  const Number<Current_Dim0> &, const Number<Current_Dim1> &,
370  const Number<Current_Dim2> &)
371  {
372  return a(Current_Dim0 - 1, Current_Dim1 - 1, Current_Dim2 - 1)
373  * b(Current_Dim2 - 1, Current_Dim1 - 1, Current_Dim0 - 1)
374  + T3as_times_T3_210(a, b, Number<Current_Dim0>(),
375  Number<Current_Dim1 - 1>(),
376  Number<Current_Dim2>());
377  }
378 
379  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
380  char j, char k, int Current_Dim0, int Current_Dim2>
381  typename promote<T, U>::V T3as_times_T3_210(
382  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
383  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, k, j, i> &b,
384  const Number<Current_Dim0> &, const Number<1> &,
385  const Number<Current_Dim2> &)
386  {
387  return a(Current_Dim0 - 1, 0, Current_Dim2 - 1)
388  * b(Current_Dim2 - 1, 0, Current_Dim0 - 1)
389  + T3as_times_T3_210(a, b, Number<Current_Dim0>(), Number<Dim12>(),
390  Number<Current_Dim2 - 1>());
391  }
392 
393  /* A special case for when the last two indices are equal. */
394 
395  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
396  char j, char k, int Current_Dim0, int Current_Dim2>
397  typename promote<T, U>::V T3as_times_T3_210(
398  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
399  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, k, j, i> &b,
400  const Number<Current_Dim0> &, const Number<Current_Dim2> &,
401  const Number<Current_Dim2> &)
402  {
403  return T3as_times_T3_210(a, b, Number<Current_Dim0>(),
404  Number<Current_Dim2 - 1>(),
405  Number<Current_Dim2>());
406  }
407 
408  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
409  char j, char k, int Current_Dim0>
410  typename promote<T, U>::V T3as_times_T3_210(
411  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
412  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, k, j, i> &b,
413  const Number<Current_Dim0> &, const Number<2> &, const Number<1> &)
414  {
415  return a(Current_Dim0 - 1, 1, 0) * b(0, 1, Current_Dim0 - 1)
416  + T3as_times_T3_210(a, b, Number<Current_Dim0 - 1>(),
417  Number<Dim12 - 1>(), Number<Dim12>());
418  }
419 
420  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
421  char j, char k>
422  typename promote<T, U>::V T3as_times_T3_210(
423  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
424  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, k, j, i> &b,
425  const Number<1> &, const Number<2> &, const Number<1> &)
426  {
427  return a(0, 1, 0) * b(0, 1, 0);
428  }
429 
430  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
431  char j, char k>
432  typename promote<T, U>::V
433  operator*(const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
434  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, k, j, i> &b)
435  {
436  return T3as_times_T3_210(a, b, Number<Dim0>(), Number<Dim12 - 1>(),
437  Number<Dim12>());
438  }
439 
440  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
441  char j, char k>
442  typename promote<T, U>::V
443  operator*(const Tensor3_Expr<B, U, Dim12, Dim12, Dim0, k, j, i> &b,
444  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a)
445  {
446  return a * b;
447  }
448 
449  /* A(i,j,k)*B(i,k,j) */
450 
451  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
452  char j, char k, int Current_Dim0, int Current_Dim1,
453  int Current_Dim2>
454  typename promote<T, U>::V T3as_times_T3_012(
455  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
456  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, i, k, j> &b,
457  const Number<Current_Dim0> &, const Number<Current_Dim1> &,
458  const Number<Current_Dim2> &)
459  {
460  return a(Current_Dim0 - 1, Current_Dim1 - 1, Current_Dim2 - 1)
461  * b(Current_Dim0 - 1, Current_Dim2 - 1, Current_Dim1 - 1)
462  + T3as_times_T3_012(a, b, Number<Current_Dim0>(),
463  Number<Current_Dim1 - 1>(),
464  Number<Current_Dim2>());
465  }
466 
467  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
468  char j, char k, int Current_Dim0, int Current_Dim2>
469  typename promote<T, U>::V T3as_times_T3_012(
470  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
471  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, i, k, j> &b,
472  const Number<Current_Dim0> &, const Number<1> &,
473  const Number<Current_Dim2> &)
474  {
475  return a(Current_Dim0 - 1, 0, Current_Dim2 - 1)
476  * b(Current_Dim0 - 1, Current_Dim2 - 1, 0)
477  + T3as_times_T3_012(a, b, Number<Current_Dim0>(), Number<Dim12>(),
478  Number<Current_Dim2 - 1>());
479  }
480 
481  /* A special case for when the last two indices are equal. */
482 
483  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
484  char j, char k, int Current_Dim0, int Current_Dim2>
485  typename promote<T, U>::V T3as_times_T3_012(
486  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
487  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, i, k, j> &b,
488  const Number<Current_Dim0> &, const Number<Current_Dim2> &,
489  const Number<Current_Dim2> &)
490  {
491  return T3as_times_T3_012(a, b, Number<Current_Dim0>(),
492  Number<Current_Dim2 - 1>(),
493  Number<Current_Dim2>());
494  }
495 
496  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
497  char j, char k, int Current_Dim0>
498  typename promote<T, U>::V T3as_times_T3_012(
499  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
500  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, i, k, j> &b,
501  const Number<Current_Dim0> &, const Number<2> &, const Number<1> &)
502  {
503  return a(Current_Dim0 - 1, 1, 0) * b(Current_Dim0 - 1, 0, 1)
504  + T3as_times_T3_012(a, b, Number<Current_Dim0 - 1>(),
505  Number<Dim12 - 1>(), Number<Dim12>());
506  }
507 
508  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
509  char j, char k>
510  typename promote<T, U>::V T3as_times_T3_012(
511  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
512  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, i, k, j> &b,
513  const Number<1> &, const Number<2> &, const Number<1> &)
514  {
515  return a(0, 1, 0) * b(0, 0, 1);
516  }
517 
518  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
519  char j, char k>
520  typename promote<T, U>::V
521  operator*(const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a,
522  const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, i, k, j> &b)
523  {
524  return T3as_times_T3_012(a, b, Number<Dim0>(), Number<Dim12 - 1>(),
525  Number<Dim12>());
526  }
527 
528  template <class A, class B, class T, class U, int Dim0, int Dim12, char i,
529  char j, char k>
530  typename promote<T, U>::V
531  operator*(const Tensor3_Expr<B, U, Dim0, Dim12, Dim12, i, k, j> &b,
532  const Tensor3_antisymmetric_Expr<A, T, Dim0, Dim12, i, j, k> &a)
533  {
534  return a * b;
535  }
536 }
FTensor
JSON compatible output.
Definition: Christof_constructor.hpp:6
FTensor::T3as_times_T3_012
promote< T, U >::V T3as_times_T3_012(const Tensor3_antisymmetric_Expr< A, T, Dim0, Dim12, i, j, k > &a, const Tensor3_Expr< B, U, Dim0, Dim12, Dim12, i, j, k > &b, const Number< Current_Dim0 > &, const Number< Current_Dim1 > &, const Number< Current_Dim2 > &)
Definition: Tensor3_antisymmetric_times_Tensor3.hpp:14
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
FTensor::T3as_times_T3_120
promote< T, U >::V T3as_times_T3_120(const Tensor3_antisymmetric_Expr< A, T, Dim0, Dim12, i, j, k > &a, const Tensor3_Expr< B, U, Dim0, Dim12, Dim12, j, k, i > &b, const Number< Current_Dim0 > &, const Number< Current_Dim1 > &, const Number< Current_Dim2 > &)
Definition: Tensor3_antisymmetric_times_Tensor3.hpp:190
A
constexpr AssemblyType A
Definition: operators_tests.cpp:30
FTensor::Tensor3_antisymmetric_Expr
Definition: Tensor3_antisymmetric_Expr.hpp:18
FTensor::Number
Definition: Number.hpp:11
a
constexpr double a
Definition: approx_sphere.cpp:30
FTensor::T3as_times_T3_102
promote< T, U >::V T3as_times_T3_102(const Tensor3_antisymmetric_Expr< A, T, Dim0, Dim12, i, j, k > &a, const Tensor3_Expr< B, U, Dim0, Dim12, Dim12, j, i, k > &b, const Number< Current_Dim0 > &, const Number< Current_Dim1 > &, const Number< Current_Dim2 > &)
Definition: Tensor3_antisymmetric_times_Tensor3.hpp:278
FTensor::promote::V
T1 V
Definition: promote.hpp:17
FTensor::Tensor3_Expr
Definition: Tensor3_Expr.hpp:24
i
FTensor::Index< 'i', SPACE_DIM > i
Definition: hcurl_divergence_operator_2d.cpp:27
FTensor::T3as_times_T3_201
promote< T, U >::V T3as_times_T3_201(const Tensor3_antisymmetric_Expr< A, T, Dim0, Dim12, i, j, k > &a, const Tensor3_Expr< B, U, Dim0, Dim12, Dim12, k, i, j > &b, const Number< Current_Dim0 > &, const Number< Current_Dim1 > &, const Number< Current_Dim2 > &)
Definition: Tensor3_antisymmetric_times_Tensor3.hpp:102
j
FTensor::Index< 'j', 3 > j
Definition: matrix_function.cpp:19
k
FTensor::Index< 'k', 3 > k
Definition: matrix_function.cpp:20
FTensor::T3as_times_T3_210
promote< T, U >::V T3as_times_T3_210(const Tensor3_antisymmetric_Expr< A, T, Dim0, Dim12, i, j, k > &a, const Tensor3_Expr< B, U, Dim0, Dim12, Dim12, k, j, i > &b, const Number< Current_Dim0 > &, const Number< Current_Dim1 > &, const Number< Current_Dim2 > &)
Definition: Tensor3_antisymmetric_times_Tensor3.hpp:366
EshelbianPlasticity::U
@ U
Definition: EshelbianContact.cpp:193
Generated by Doxygen 1.8.17 and hosted at University of Glasgow