v0.14.0
test_T4R.cpp
Go to the documentation of this file.
1 #include "../../src/FTensor.hpp"
2 #include "test_for_zero.hpp"
3 #include <iostream>
4 using namespace FTensor;
5 using namespace std;
6 
7 void test_T4R()
8 {
9  Index<'i', 3> i;
10  Index<'j', 3> j;
11  Index<'k', 3> k;
12  Index<'l', 3> l;
13 
14  Number<0> N0;
15  Number<1> N1;
16  Number<2> N2;
17 
18  Riemann<double, 3> t4R_1, t4R_2, t4R_3;
19 
20  Tensor2_symmetric<double, 3> delta_3(3., 0., 0., 3., 0., 3.),
21  delta_5(5., 0., 0., 5., 0., 5.), delta_2(2., 0., 0., 2., 0., 2.),
22  delta_7(7., 0., 0., 7., 0., 7.);
23 
24  t4R_1(i, j, k, l)
25  = (delta_3(i, k) * delta_5(j, l) && delta_3(i, l) * delta_5(k, j));
26 
27  test_for_zero(
28  t4R_1(0, 0, 0, 0)
29  - (delta_3(0, 0) * delta_5(0, 0) - delta_3(0, 0) * delta_5(0, 0)),
30  "T4ddg && T4ddg (0,0,0,0)");
31  test_for_zero(
32  t4R_1(0, 0, 0, 1)
33  - (delta_3(0, 1) * delta_5(0, 0) - delta_3(0, 1) * delta_5(0, 0)),
34  "T4ddg && T4ddg (0,0,0,1)");
35  test_for_zero(
36  t4R_1(0, 0, 0, 2)
37  - (delta_3(0, 2) * delta_5(0, 0) - delta_3(0, 2) * delta_5(0, 0)),
38  "T4ddg && T4ddg (0,0,0,2)");
39  test_for_zero(
40  t4R_1(0, 0, 1, 0)
41  - (delta_3(0, 0) * delta_5(0, 1) - delta_3(0, 0) * delta_5(1, 0)),
42  "T4ddg && T4ddg (0,0,1,0)");
43  test_for_zero(
44  t4R_1(0, 0, 1, 1)
45  - (delta_3(0, 1) * delta_5(0, 1) - delta_3(0, 1) * delta_5(1, 0)),
46  "T4ddg && T4ddg (0,0,1,1)");
47  test_for_zero(
48  t4R_1(0, 0, 1, 2)
49  - (delta_3(0, 2) * delta_5(0, 1) - delta_3(0, 2) * delta_5(1, 0)),
50  "T4ddg && T4ddg (0,0,1,2)");
51  test_for_zero(
52  t4R_1(0, 0, 2, 0)
53  - (delta_3(0, 0) * delta_5(0, 2) - delta_3(0, 0) * delta_5(2, 0)),
54  "T4ddg && T4ddg (0,0,2,0)");
55  test_for_zero(
56  t4R_1(0, 0, 2, 1)
57  - (delta_3(0, 1) * delta_5(0, 2) - delta_3(0, 1) * delta_5(2, 0)),
58  "T4ddg && T4ddg (0,0,2,1)");
59  test_for_zero(
60  t4R_1(0, 0, 2, 2)
61  - (delta_3(0, 2) * delta_5(0, 2) - delta_3(0, 2) * delta_5(2, 0)),
62  "T4ddg && T4ddg (0,0,2,2)");
63  test_for_zero(
64  t4R_1(0, 1, 0, 0)
65  - (delta_3(0, 0) * delta_5(1, 0) - delta_3(0, 0) * delta_5(0, 1)),
66  "T4ddg && T4ddg (0,1,0,0)");
67  test_for_zero(
68  t4R_1(0, 1, 0, 1)
69  - (delta_3(0, 0) * delta_5(1, 1) - delta_3(0, 1) * delta_5(0, 1)),
70  "T4ddg && T4ddg (0,1,0,1)");
71  test_for_zero(
72  t4R_1(0, 1, 0, 2)
73  - (delta_3(0, 2) * delta_5(1, 0) - delta_3(0, 2) * delta_5(0, 1)),
74  "T4ddg && T4ddg (0,1,0,2)");
75  test_for_zero(
76  t4R_1(0, 1, 1, 0)
77  - (delta_3(0, 0) * delta_5(1, 1) - delta_3(0, 0) * delta_5(1, 1)),
78  "T4ddg && T4ddg (0,1,1,0)");
79  test_for_zero(
80  t4R_1(0, 1, 1, 1)
81  - (delta_3(0, 1) * delta_5(1, 1) - delta_3(0, 1) * delta_5(1, 1)),
82  "T4ddg && T4ddg (0,1,1,1)");
83  test_for_zero(
84  t4R_1(0, 1, 1, 2)
85  - (delta_3(0, 2) * delta_5(1, 1) - delta_3(0, 2) * delta_5(1, 1)),
86  "T4ddg && T4ddg (0,1,1,2)");
87  test_for_zero(
88  t4R_1(0, 1, 2, 0)
89  - (delta_3(0, 0) * delta_5(1, 2) - delta_3(0, 0) * delta_5(2, 1)),
90  "T4ddg && T4ddg (0,1,2,0)");
91  test_for_zero(
92  t4R_1(0, 1, 2, 1)
93  - (delta_3(0, 1) * delta_5(1, 2) - delta_3(0, 1) * delta_5(2, 1)),
94  "T4ddg && T4ddg (0,1,2,1)");
95  test_for_zero(
96  t4R_1(0, 1, 2, 2)
97  - (delta_3(0, 2) * delta_5(1, 2) - delta_3(0, 2) * delta_5(2, 1)),
98  "T4ddg && T4ddg (0,1,2,2)");
99  test_for_zero(
100  t4R_1(0, 2, 0, 0)
101  - (delta_3(0, 0) * delta_5(2, 0) - delta_3(0, 0) * delta_5(0, 2)),
102  "T4ddg && T4ddg (0,2,0,0)");
103  test_for_zero(
104  t4R_1(0, 2, 0, 1)
105  - (delta_3(0, 1) * delta_5(2, 0) - delta_3(0, 1) * delta_5(0, 2)),
106  "T4ddg && T4ddg (0,2,0,1)");
107  test_for_zero(
108  t4R_1(0, 2, 0, 2)
109  - (delta_3(0, 0) * delta_5(2, 2) - delta_3(0, 2) * delta_5(0, 2)),
110  "T4ddg && T4ddg (0,2,0,2)");
111  test_for_zero(
112  t4R_1(0, 2, 1, 0)
113  - (delta_3(0, 0) * delta_5(2, 1) - delta_3(0, 0) * delta_5(1, 2)),
114  "T4ddg && T4ddg (0,2,1,0)");
115  test_for_zero(
116  t4R_1(0, 2, 1, 1)
117  - (delta_3(0, 1) * delta_5(2, 1) - delta_3(0, 1) * delta_5(1, 2)),
118  "T4ddg && T4ddg (0,2,1,1)");
119  test_for_zero(
120  t4R_1(0, 2, 1, 2)
121  - (delta_3(0, 2) * delta_5(2, 1) - delta_3(0, 2) * delta_5(1, 2)),
122  "T4ddg && T4ddg (0,2,1,2)");
123  test_for_zero(
124  t4R_1(0, 2, 2, 0)
125  - (delta_3(0, 0) * delta_5(2, 2) - delta_3(0, 0) * delta_5(2, 2)),
126  "T4ddg && T4ddg (0,2,2,0)");
127  test_for_zero(
128  t4R_1(0, 2, 2, 1)
129  - (delta_3(0, 1) * delta_5(2, 2) - delta_3(0, 1) * delta_5(2, 2)),
130  "T4ddg && T4ddg (0,2,2,1)");
131  test_for_zero(
132  t4R_1(0, 2, 2, 2)
133  - (delta_3(0, 2) * delta_5(2, 2) - delta_3(0, 2) * delta_5(2, 2)),
134  "T4ddg && T4ddg (0,2,2,2)");
135  test_for_zero(
136  t4R_1(1, 0, 0, 0)
137  - (delta_3(1, 0) * delta_5(0, 0) - delta_3(1, 0) * delta_5(0, 0)),
138  "T4ddg && T4ddg (1,0,0,0)");
139  test_for_zero(
140  t4R_1(1, 0, 0, 1)
141  - (delta_3(1, 1) * delta_5(0, 0) - delta_3(1, 1) * delta_5(0, 0)),
142  "T4ddg && T4ddg (1,0,0,1)");
143  test_for_zero(
144  t4R_1(1, 0, 0, 2)
145  - (delta_3(1, 2) * delta_5(0, 0) - delta_3(1, 2) * delta_5(0, 0)),
146  "T4ddg && T4ddg (1,0,0,2)");
147  test_for_zero(
148  t4R_1(1, 0, 1, 0)
149  - (delta_3(1, 0) * delta_5(0, 1) - delta_3(1, 0) * delta_5(1, 0)),
150  "T4ddg && T4ddg (1,0,1,0)");
151  test_for_zero(
152  t4R_1(1, 0, 1, 1)
153  - (delta_3(1, 1) * delta_5(0, 1) - delta_3(1, 1) * delta_5(1, 0)),
154  "T4ddg && T4ddg (1,0,1,1)");
155  test_for_zero(
156  t4R_1(1, 0, 1, 2)
157  - (delta_3(1, 2) * delta_5(0, 1) - delta_3(1, 2) * delta_5(1, 0)),
158  "T4ddg && T4ddg (1,0,1,2)");
159  test_for_zero(
160  t4R_1(1, 0, 2, 0)
161  - (delta_3(1, 0) * delta_5(0, 2) - delta_3(1, 0) * delta_5(2, 0)),
162  "T4ddg && T4ddg (1,0,2,0)");
163  test_for_zero(
164  t4R_1(1, 0, 2, 1)
165  - (delta_3(1, 1) * delta_5(0, 2) - delta_3(1, 1) * delta_5(2, 0)),
166  "T4ddg && T4ddg (1,0,2,1)");
167  test_for_zero(
168  t4R_1(1, 0, 2, 2)
169  - (delta_3(1, 2) * delta_5(0, 2) - delta_3(1, 2) * delta_5(2, 0)),
170  "T4ddg && T4ddg (1,0,2,2)");
171  test_for_zero(
172  t4R_1(1, 1, 0, 0)
173  - (delta_3(1, 0) * delta_5(1, 0) - delta_3(1, 0) * delta_5(0, 1)),
174  "T4ddg && T4ddg (1,1,0,0)");
175  test_for_zero(
176  t4R_1(1, 1, 0, 1)
177  - (delta_3(1, 1) * delta_5(1, 0) - delta_3(1, 1) * delta_5(0, 1)),
178  "T4ddg && T4ddg (1,1,0,1)");
179  test_for_zero(
180  t4R_1(1, 1, 0, 2)
181  - (delta_3(1, 2) * delta_5(1, 0) - delta_3(1, 2) * delta_5(0, 1)),
182  "T4ddg && T4ddg (1,1,0,2)");
183  test_for_zero(
184  t4R_1(1, 1, 1, 0)
185  - (delta_3(1, 0) * delta_5(1, 1) - delta_3(1, 0) * delta_5(1, 1)),
186  "T4ddg && T4ddg (1,1,1,0)");
187  test_for_zero(
188  t4R_1(1, 1, 1, 1)
189  - (delta_3(1, 1) * delta_5(1, 1) - delta_3(1, 1) * delta_5(1, 1)),
190  "T4ddg && T4ddg (1,1,1,1)");
191  test_for_zero(
192  t4R_1(1, 1, 1, 2)
193  - (delta_3(1, 2) * delta_5(1, 1) - delta_3(1, 2) * delta_5(1, 1)),
194  "T4ddg && T4ddg (1,1,1,2)");
195  test_for_zero(
196  t4R_1(1, 1, 2, 0)
197  - (delta_3(1, 0) * delta_5(1, 2) - delta_3(1, 0) * delta_5(2, 1)),
198  "T4ddg && T4ddg (1,1,2,0)");
199  test_for_zero(
200  t4R_1(1, 1, 2, 1)
201  - (delta_3(1, 1) * delta_5(1, 2) - delta_3(1, 1) * delta_5(2, 1)),
202  "T4ddg && T4ddg (1,1,2,1)");
203  test_for_zero(
204  t4R_1(1, 1, 2, 2)
205  - (delta_3(1, 2) * delta_5(1, 2) - delta_3(1, 2) * delta_5(2, 1)),
206  "T4ddg && T4ddg (1,1,2,2)");
207  test_for_zero(
208  t4R_1(1, 2, 0, 0)
209  - (delta_3(1, 0) * delta_5(2, 0) - delta_3(1, 0) * delta_5(0, 2)),
210  "T4ddg && T4ddg (1,2,0,0)");
211  test_for_zero(
212  t4R_1(1, 2, 0, 1)
213  - (delta_3(1, 1) * delta_5(2, 0) - delta_3(1, 1) * delta_5(0, 2)),
214  "T4ddg && T4ddg (1,2,0,1)");
215  test_for_zero(
216  t4R_1(1, 2, 0, 2)
217  - (delta_3(1, 2) * delta_5(2, 0) - delta_3(1, 2) * delta_5(0, 2)),
218  "T4ddg && T4ddg (1,2,0,2)");
219  test_for_zero(
220  t4R_1(1, 2, 1, 0)
221  - (delta_3(1, 0) * delta_5(2, 1) - delta_3(1, 0) * delta_5(1, 2)),
222  "T4ddg && T4ddg (1,2,1,0)");
223  test_for_zero(
224  t4R_1(1, 2, 1, 1)
225  - (delta_3(1, 1) * delta_5(2, 1) - delta_3(1, 1) * delta_5(1, 2)),
226  "T4ddg && T4ddg (1,2,1,1)");
227  test_for_zero(
228  t4R_1(1, 2, 1, 2)
229  - (delta_3(1, 1) * delta_5(2, 2) - delta_3(1, 2) * delta_5(1, 2)),
230  "T4ddg && T4ddg (1,2,1,2)");
231  test_for_zero(
232  t4R_1(1, 2, 2, 0)
233  - (delta_3(1, 0) * delta_5(2, 2) - delta_3(1, 0) * delta_5(2, 2)),
234  "T4ddg && T4ddg (1,2,2,0)");
235  test_for_zero(
236  t4R_1(1, 2, 2, 1)
237  - (delta_3(1, 1) * delta_5(2, 2) - delta_3(1, 1) * delta_5(2, 2)),
238  "T4ddg && T4ddg (1,2,2,1)");
239  test_for_zero(
240  t4R_1(1, 2, 2, 2)
241  - (delta_3(1, 2) * delta_5(2, 2) - delta_3(1, 2) * delta_5(2, 2)),
242  "T4ddg && T4ddg (1,2,2,2)");
243  test_for_zero(
244  t4R_1(2, 0, 0, 0)
245  - (delta_3(2, 0) * delta_5(0, 0) - delta_3(2, 0) * delta_5(0, 0)),
246  "T4ddg && T4ddg (2,0,0,0)");
247  test_for_zero(
248  t4R_1(2, 0, 0, 1)
249  - (delta_3(2, 1) * delta_5(0, 0) - delta_3(2, 1) * delta_5(0, 0)),
250  "T4ddg && T4ddg (2,0,0,1)");
251  test_for_zero(
252  t4R_1(2, 0, 0, 2)
253  - (delta_3(2, 2) * delta_5(0, 0) - delta_3(2, 2) * delta_5(0, 0)),
254  "T4ddg && T4ddg (2,0,0,2)");
255  test_for_zero(
256  t4R_1(2, 0, 1, 0)
257  - (delta_3(2, 0) * delta_5(0, 1) - delta_3(2, 0) * delta_5(1, 0)),
258  "T4ddg && T4ddg (2,0,1,0)");
259  test_for_zero(
260  t4R_1(2, 0, 1, 1)
261  - (delta_3(2, 1) * delta_5(0, 1) - delta_3(2, 1) * delta_5(1, 0)),
262  "T4ddg && T4ddg (2,0,1,1)");
263  test_for_zero(
264  t4R_1(2, 0, 1, 2)
265  - (delta_3(2, 2) * delta_5(0, 1) - delta_3(2, 2) * delta_5(1, 0)),
266  "T4ddg && T4ddg (2,0,1,2)");
267  test_for_zero(
268  t4R_1(2, 0, 2, 0)
269  - (delta_3(2, 0) * delta_5(0, 2) - delta_3(2, 0) * delta_5(2, 0)),
270  "T4ddg && T4ddg (2,0,2,0)");
271  test_for_zero(
272  t4R_1(2, 0, 2, 1)
273  - (delta_3(2, 1) * delta_5(0, 2) - delta_3(2, 1) * delta_5(2, 0)),
274  "T4ddg && T4ddg (2,0,2,1)");
275  test_for_zero(
276  t4R_1(2, 0, 2, 2)
277  - (delta_3(2, 2) * delta_5(0, 2) - delta_3(2, 2) * delta_5(2, 0)),
278  "T4ddg && T4ddg (2,0,2,2)");
279  test_for_zero(
280  t4R_1(2, 1, 0, 0)
281  - (delta_3(2, 0) * delta_5(1, 0) - delta_3(2, 0) * delta_5(0, 1)),
282  "T4ddg && T4ddg (2,1,0,0)");
283  test_for_zero(
284  t4R_1(2, 1, 0, 1)
285  - (delta_3(2, 1) * delta_5(1, 0) - delta_3(2, 1) * delta_5(0, 1)),
286  "T4ddg && T4ddg (2,1,0,1)");
287  test_for_zero(
288  t4R_1(2, 1, 0, 2)
289  - (delta_3(2, 2) * delta_5(1, 0) - delta_3(2, 2) * delta_5(0, 1)),
290  "T4ddg && T4ddg (2,1,0,2)");
291  test_for_zero(
292  t4R_1(2, 1, 1, 0)
293  - (delta_3(2, 0) * delta_5(1, 1) - delta_3(2, 0) * delta_5(1, 1)),
294  "T4ddg && T4ddg (2,1,1,0)");
295  test_for_zero(
296  t4R_1(2, 1, 1, 1)
297  - (delta_3(2, 1) * delta_5(1, 1) - delta_3(2, 1) * delta_5(1, 1)),
298  "T4ddg && T4ddg (2,1,1,1)");
299  test_for_zero(
300  t4R_1(2, 1, 1, 2)
301  - (delta_3(2, 2) * delta_5(1, 1) - delta_3(2, 2) * delta_5(1, 1)),
302  "T4ddg && T4ddg (2,1,1,2)");
303  test_for_zero(
304  t4R_1(2, 1, 2, 0)
305  - (delta_3(2, 0) * delta_5(1, 2) - delta_3(2, 0) * delta_5(2, 1)),
306  "T4ddg && T4ddg (2,1,2,0)");
307  test_for_zero(
308  t4R_1(2, 1, 2, 1)
309  - (delta_3(2, 1) * delta_5(1, 2) - delta_3(2, 1) * delta_5(2, 1)),
310  "T4ddg && T4ddg (2,1,2,1)");
311  test_for_zero(
312  t4R_1(2, 1, 2, 2)
313  - (delta_3(2, 2) * delta_5(1, 2) - delta_3(2, 2) * delta_5(2, 1)),
314  "T4ddg && T4ddg (2,1,2,2)");
315  test_for_zero(
316  t4R_1(2, 2, 0, 0)
317  - (delta_3(2, 0) * delta_5(2, 0) - delta_3(2, 0) * delta_5(0, 2)),
318  "T4ddg && T4ddg (2,2,0,0)");
319  test_for_zero(
320  t4R_1(2, 2, 0, 1)
321  - (delta_3(2, 1) * delta_5(2, 0) - delta_3(2, 1) * delta_5(0, 2)),
322  "T4ddg && T4ddg (2,2,0,1)");
323  test_for_zero(
324  t4R_1(2, 2, 0, 2)
325  - (delta_3(2, 2) * delta_5(2, 0) - delta_3(2, 2) * delta_5(0, 2)),
326  "T4ddg && T4ddg (2,2,0,2)");
327  test_for_zero(
328  t4R_1(2, 2, 1, 0)
329  - (delta_3(2, 0) * delta_5(2, 1) - delta_3(2, 0) * delta_5(1, 2)),
330  "T4ddg && T4ddg (2,2,1,0)");
331  test_for_zero(
332  t4R_1(2, 2, 1, 1)
333  - (delta_3(2, 1) * delta_5(2, 1) - delta_3(2, 1) * delta_5(1, 2)),
334  "T4ddg && T4ddg (2,2,1,1)");
335  test_for_zero(
336  t4R_1(2, 2, 1, 2)
337  - (delta_3(2, 2) * delta_5(2, 1) - delta_3(2, 2) * delta_5(1, 2)),
338  "T4ddg && T4ddg (2,2,1,2)");
339  test_for_zero(
340  t4R_1(2, 2, 2, 0)
341  - (delta_3(2, 0) * delta_5(2, 2) - delta_3(2, 0) * delta_5(2, 2)),
342  "T4ddg && T4ddg (2,2,2,0)");
343  test_for_zero(
344  t4R_1(2, 2, 2, 1)
345  - (delta_3(2, 1) * delta_5(2, 2) - delta_3(2, 1) * delta_5(2, 2)),
346  "T4ddg && T4ddg (2,2,2,1)");
347  test_for_zero(
348  t4R_1(2, 2, 2, 2)
349  - (delta_3(2, 2) * delta_5(2, 2) - delta_3(2, 2) * delta_5(2, 2)),
350  "T4ddg && T4ddg (2,2,2,2)");
351 
352  t4R_2(i, j, k, l)
353  = (delta_2(i, k) * delta_7(j, l) && delta_2(i, l) * delta_7(k, j));
354 
355  t4R_3(i, j, k, l) = t4R_1(i, j, k, l) + t4R_2(i, j, k, l);
356  test_for_zero(t4R_3(0, 0, 0, 0) - (t4R_1(0, 0, 0, 0) + t4R_2(0, 0, 0, 0)),
357  "T4R + T4R (0,0,0,0)");
358  test_for_zero(t4R_3(0, 0, 0, 1) - (t4R_1(0, 0, 0, 1) + t4R_2(0, 0, 0, 1)),
359  "T4R + T4R (0,0,0,1)");
360  test_for_zero(t4R_3(0, 0, 0, 2) - (t4R_1(0, 0, 0, 2) + t4R_2(0, 0, 0, 2)),
361  "T4R + T4R (0,0,0,2)");
362  test_for_zero(t4R_3(0, 0, 1, 0) - (t4R_1(0, 0, 1, 0) + t4R_2(0, 0, 1, 0)),
363  "T4R + T4R (0,0,1,0)");
364  test_for_zero(t4R_3(0, 0, 1, 1) - (t4R_1(0, 0, 1, 1) + t4R_2(0, 0, 1, 1)),
365  "T4R + T4R (0,0,1,1)");
366  test_for_zero(t4R_3(0, 0, 1, 2) - (t4R_1(0, 0, 1, 2) + t4R_2(0, 0, 1, 2)),
367  "T4R + T4R (0,0,1,2)");
368  test_for_zero(t4R_3(0, 0, 2, 0) - (t4R_1(0, 0, 2, 0) + t4R_2(0, 0, 2, 0)),
369  "T4R + T4R (0,0,2,0)");
370  test_for_zero(t4R_3(0, 0, 2, 1) - (t4R_1(0, 0, 2, 1) + t4R_2(0, 0, 2, 1)),
371  "T4R + T4R (0,0,2,1)");
372  test_for_zero(t4R_3(0, 0, 2, 2) - (t4R_1(0, 0, 2, 2) + t4R_2(0, 0, 2, 2)),
373  "T4R + T4R (0,0,2,2)");
374  test_for_zero(t4R_3(0, 1, 0, 0) - (t4R_1(0, 1, 0, 0) + t4R_2(0, 1, 0, 0)),
375  "T4R + T4R (0,1,0,0)");
376  test_for_zero(t4R_3(0, 1, 0, 1) - (t4R_1(0, 1, 0, 1) + t4R_2(0, 1, 0, 1)),
377  "T4R + T4R (0,1,0,1)");
378  test_for_zero(t4R_3(0, 1, 0, 2) - (t4R_1(0, 1, 0, 2) + t4R_2(0, 1, 0, 2)),
379  "T4R + T4R (0,1,0,2)");
380  test_for_zero(t4R_3(0, 1, 1, 0) - (t4R_1(0, 1, 1, 0) + t4R_2(0, 1, 1, 0)),
381  "T4R + T4R (0,1,1,0)");
382  test_for_zero(t4R_3(0, 1, 1, 1) - (t4R_1(0, 1, 1, 1) + t4R_2(0, 1, 1, 1)),
383  "T4R + T4R (0,1,1,1)");
384  test_for_zero(t4R_3(0, 1, 1, 2) - (t4R_1(0, 1, 1, 2) + t4R_2(0, 1, 1, 2)),
385  "T4R + T4R (0,1,1,2)");
386  test_for_zero(t4R_3(0, 1, 2, 0) - (t4R_1(0, 1, 2, 0) + t4R_2(0, 1, 2, 0)),
387  "T4R + T4R (0,1,2,0)");
388  test_for_zero(t4R_3(0, 1, 2, 1) - (t4R_1(0, 1, 2, 1) + t4R_2(0, 1, 2, 1)),
389  "T4R + T4R (0,1,2,1)");
390  test_for_zero(t4R_3(0, 1, 2, 2) - (t4R_1(0, 1, 2, 2) + t4R_2(0, 1, 2, 2)),
391  "T4R + T4R (0,1,2,2)");
392  test_for_zero(t4R_3(0, 2, 0, 0) - (t4R_1(0, 2, 0, 0) + t4R_2(0, 2, 0, 0)),
393  "T4R + T4R (0,2,0,0)");
394  test_for_zero(t4R_3(0, 2, 0, 1) - (t4R_1(0, 2, 0, 1) + t4R_2(0, 2, 0, 1)),
395  "T4R + T4R (0,2,0,1)");
396  test_for_zero(t4R_3(0, 2, 0, 2) - (t4R_1(0, 2, 0, 2) + t4R_2(0, 2, 0, 2)),
397  "T4R + T4R (0,2,0,2)");
398  test_for_zero(t4R_3(0, 2, 1, 0) - (t4R_1(0, 2, 1, 0) + t4R_2(0, 2, 1, 0)),
399  "T4R + T4R (0,2,1,0)");
400  test_for_zero(t4R_3(0, 2, 1, 1) - (t4R_1(0, 2, 1, 1) + t4R_2(0, 2, 1, 1)),
401  "T4R + T4R (0,2,1,1)");
402  test_for_zero(t4R_3(0, 2, 1, 2) - (t4R_1(0, 2, 1, 2) + t4R_2(0, 2, 1, 2)),
403  "T4R + T4R (0,2,1,2)");
404  test_for_zero(t4R_3(0, 2, 2, 0) - (t4R_1(0, 2, 2, 0) + t4R_2(0, 2, 2, 0)),
405  "T4R + T4R (0,2,2,0)");
406  test_for_zero(t4R_3(0, 2, 2, 1) - (t4R_1(0, 2, 2, 1) + t4R_2(0, 2, 2, 1)),
407  "T4R + T4R (0,2,2,1)");
408  test_for_zero(t4R_3(0, 2, 2, 2) - (t4R_1(0, 2, 2, 2) + t4R_2(0, 2, 2, 2)),
409  "T4R + T4R (0,2,2,2)");
410  test_for_zero(t4R_3(1, 0, 0, 0) - (t4R_1(1, 0, 0, 0) + t4R_2(1, 0, 0, 0)),
411  "T4R + T4R (1,0,0,0)");
412  test_for_zero(t4R_3(1, 0, 0, 1) - (t4R_1(1, 0, 0, 1) + t4R_2(1, 0, 0, 1)),
413  "T4R + T4R (1,0,0,1)");
414  test_for_zero(t4R_3(1, 0, 0, 2) - (t4R_1(1, 0, 0, 2) + t4R_2(1, 0, 0, 2)),
415  "T4R + T4R (1,0,0,2)");
416  test_for_zero(t4R_3(1, 0, 1, 0) - (t4R_1(1, 0, 1, 0) + t4R_2(1, 0, 1, 0)),
417  "T4R + T4R (1,0,1,0)");
418  test_for_zero(t4R_3(1, 0, 1, 1) - (t4R_1(1, 0, 1, 1) + t4R_2(1, 0, 1, 1)),
419  "T4R + T4R (1,0,1,1)");
420  test_for_zero(t4R_3(1, 0, 1, 2) - (t4R_1(1, 0, 1, 2) + t4R_2(1, 0, 1, 2)),
421  "T4R + T4R (1,0,1,2)");
422  test_for_zero(t4R_3(1, 0, 2, 0) - (t4R_1(1, 0, 2, 0) + t4R_2(1, 0, 2, 0)),
423  "T4R + T4R (1,0,2,0)");
424  test_for_zero(t4R_3(1, 0, 2, 1) - (t4R_1(1, 0, 2, 1) + t4R_2(1, 0, 2, 1)),
425  "T4R + T4R (1,0,2,1)");
426  test_for_zero(t4R_3(1, 0, 2, 2) - (t4R_1(1, 0, 2, 2) + t4R_2(1, 0, 2, 2)),
427  "T4R + T4R (1,0,2,2)");
428  test_for_zero(t4R_3(1, 1, 0, 0) - (t4R_1(1, 1, 0, 0) + t4R_2(1, 1, 0, 0)),
429  "T4R + T4R (1,1,0,0)");
430  test_for_zero(t4R_3(1, 1, 0, 1) - (t4R_1(1, 1, 0, 1) + t4R_2(1, 1, 0, 1)),
431  "T4R + T4R (1,1,0,1)");
432  test_for_zero(t4R_3(1, 1, 0, 2) - (t4R_1(1, 1, 0, 2) + t4R_2(1, 1, 0, 2)),
433  "T4R + T4R (1,1,0,2)");
434  test_for_zero(t4R_3(1, 1, 1, 0) - (t4R_1(1, 1, 1, 0) + t4R_2(1, 1, 1, 0)),
435  "T4R + T4R (1,1,1,0)");
436  test_for_zero(t4R_3(1, 1, 1, 1) - (t4R_1(1, 1, 1, 1) + t4R_2(1, 1, 1, 1)),
437  "T4R + T4R (1,1,1,1)");
438  test_for_zero(t4R_3(1, 1, 1, 2) - (t4R_1(1, 1, 1, 2) + t4R_2(1, 1, 1, 2)),
439  "T4R + T4R (1,1,1,2)");
440  test_for_zero(t4R_3(1, 1, 2, 0) - (t4R_1(1, 1, 2, 0) + t4R_2(1, 1, 2, 0)),
441  "T4R + T4R (1,1,2,0)");
442  test_for_zero(t4R_3(1, 1, 2, 1) - (t4R_1(1, 1, 2, 1) + t4R_2(1, 1, 2, 1)),
443  "T4R + T4R (1,1,2,1)");
444  test_for_zero(t4R_3(1, 1, 2, 2) - (t4R_1(1, 1, 2, 2) + t4R_2(1, 1, 2, 2)),
445  "T4R + T4R (1,1,2,2)");
446  test_for_zero(t4R_3(1, 2, 0, 0) - (t4R_1(1, 2, 0, 0) + t4R_2(1, 2, 0, 0)),
447  "T4R + T4R (1,2,0,0)");
448  test_for_zero(t4R_3(1, 2, 0, 1) - (t4R_1(1, 2, 0, 1) + t4R_2(1, 2, 0, 1)),
449  "T4R + T4R (1,2,0,1)");
450  test_for_zero(t4R_3(1, 2, 0, 2) - (t4R_1(1, 2, 0, 2) + t4R_2(1, 2, 0, 2)),
451  "T4R + T4R (1,2,0,2)");
452  test_for_zero(t4R_3(1, 2, 1, 0) - (t4R_1(1, 2, 1, 0) + t4R_2(1, 2, 1, 0)),
453  "T4R + T4R (1,2,1,0)");
454  test_for_zero(t4R_3(1, 2, 1, 1) - (t4R_1(1, 2, 1, 1) + t4R_2(1, 2, 1, 1)),
455  "T4R + T4R (1,2,1,1)");
456  test_for_zero(t4R_3(1, 2, 1, 2) - (t4R_1(1, 2, 1, 2) + t4R_2(1, 2, 1, 2)),
457  "T4R + T4R (1,2,1,2)");
458  test_for_zero(t4R_3(1, 2, 2, 0) - (t4R_1(1, 2, 2, 0) + t4R_2(1, 2, 2, 0)),
459  "T4R + T4R (1,2,2,0)");
460  test_for_zero(t4R_3(1, 2, 2, 1) - (t4R_1(1, 2, 2, 1) + t4R_2(1, 2, 2, 1)),
461  "T4R + T4R (1,2,2,1)");
462  test_for_zero(t4R_3(1, 2, 2, 2) - (t4R_1(1, 2, 2, 2) + t4R_2(1, 2, 2, 2)),
463  "T4R + T4R (1,2,2,2)");
464  test_for_zero(t4R_3(2, 0, 0, 0) - (t4R_1(2, 0, 0, 0) + t4R_2(2, 0, 0, 0)),
465  "T4R + T4R (2,0,0,0)");
466  test_for_zero(t4R_3(2, 0, 0, 1) - (t4R_1(2, 0, 0, 1) + t4R_2(2, 0, 0, 1)),
467  "T4R + T4R (2,0,0,1)");
468  test_for_zero(t4R_3(2, 0, 0, 2) - (t4R_1(2, 0, 0, 2) + t4R_2(2, 0, 0, 2)),
469  "T4R + T4R (2,0,0,2)");
470  test_for_zero(t4R_3(2, 0, 1, 0) - (t4R_1(2, 0, 1, 0) + t4R_2(2, 0, 1, 0)),
471  "T4R + T4R (2,0,1,0)");
472  test_for_zero(t4R_3(2, 0, 1, 1) - (t4R_1(2, 0, 1, 1) + t4R_2(2, 0, 1, 1)),
473  "T4R + T4R (2,0,1,1)");
474  test_for_zero(t4R_3(2, 0, 1, 2) - (t4R_1(2, 0, 1, 2) + t4R_2(2, 0, 1, 2)),
475  "T4R + T4R (2,0,1,2)");
476  test_for_zero(t4R_3(2, 0, 2, 0) - (t4R_1(2, 0, 2, 0) + t4R_2(2, 0, 2, 0)),
477  "T4R + T4R (2,0,2,0)");
478  test_for_zero(t4R_3(2, 0, 2, 1) - (t4R_1(2, 0, 2, 1) + t4R_2(2, 0, 2, 1)),
479  "T4R + T4R (2,0,2,1)");
480  test_for_zero(t4R_3(2, 0, 2, 2) - (t4R_1(2, 0, 2, 2) + t4R_2(2, 0, 2, 2)),
481  "T4R + T4R (2,0,2,2)");
482  test_for_zero(t4R_3(2, 1, 0, 0) - (t4R_1(2, 1, 0, 0) + t4R_2(2, 1, 0, 0)),
483  "T4R + T4R (2,1,0,0)");
484  test_for_zero(t4R_3(2, 1, 0, 1) - (t4R_1(2, 1, 0, 1) + t4R_2(2, 1, 0, 1)),
485  "T4R + T4R (2,1,0,1)");
486  test_for_zero(t4R_3(2, 1, 0, 2) - (t4R_1(2, 1, 0, 2) + t4R_2(2, 1, 0, 2)),
487  "T4R + T4R (2,1,0,2)");
488  test_for_zero(t4R_3(2, 1, 1, 0) - (t4R_1(2, 1, 1, 0) + t4R_2(2, 1, 1, 0)),
489  "T4R + T4R (2,1,1,0)");
490  test_for_zero(t4R_3(2, 1, 1, 1) - (t4R_1(2, 1, 1, 1) + t4R_2(2, 1, 1, 1)),
491  "T4R + T4R (2,1,1,1)");
492  test_for_zero(t4R_3(2, 1, 1, 2) - (t4R_1(2, 1, 1, 2) + t4R_2(2, 1, 1, 2)),
493  "T4R + T4R (2,1,1,2)");
494  test_for_zero(t4R_3(2, 1, 2, 0) - (t4R_1(2, 1, 2, 0) + t4R_2(2, 1, 2, 0)),
495  "T4R + T4R (2,1,2,0)");
496  test_for_zero(t4R_3(2, 1, 2, 1) - (t4R_1(2, 1, 2, 1) + t4R_2(2, 1, 2, 1)),
497  "T4R + T4R (2,1,2,1)");
498  test_for_zero(t4R_3(2, 1, 2, 2) - (t4R_1(2, 1, 2, 2) + t4R_2(2, 1, 2, 2)),
499  "T4R + T4R (2,1,2,2)");
500  test_for_zero(t4R_3(2, 2, 0, 0) - (t4R_1(2, 2, 0, 0) + t4R_2(2, 2, 0, 0)),
501  "T4R + T4R (2,2,0,0)");
502  test_for_zero(t4R_3(2, 2, 0, 1) - (t4R_1(2, 2, 0, 1) + t4R_2(2, 2, 0, 1)),
503  "T4R + T4R (2,2,0,1)");
504  test_for_zero(t4R_3(2, 2, 0, 2) - (t4R_1(2, 2, 0, 2) + t4R_2(2, 2, 0, 2)),
505  "T4R + T4R (2,2,0,2)");
506  test_for_zero(t4R_3(2, 2, 1, 0) - (t4R_1(2, 2, 1, 0) + t4R_2(2, 2, 1, 0)),
507  "T4R + T4R (2,2,1,0)");
508  test_for_zero(t4R_3(2, 2, 1, 1) - (t4R_1(2, 2, 1, 1) + t4R_2(2, 2, 1, 1)),
509  "T4R + T4R (2,2,1,1)");
510  test_for_zero(t4R_3(2, 2, 1, 2) - (t4R_1(2, 2, 1, 2) + t4R_2(2, 2, 1, 2)),
511  "T4R + T4R (2,2,1,2)");
512  test_for_zero(t4R_3(2, 2, 2, 0) - (t4R_1(2, 2, 2, 0) + t4R_2(2, 2, 2, 0)),
513  "T4R + T4R (2,2,2,0)");
514  test_for_zero(t4R_3(2, 2, 2, 1) - (t4R_1(2, 2, 2, 1) + t4R_2(2, 2, 2, 1)),
515  "T4R + T4R (2,2,2,1)");
516  test_for_zero(t4R_3(2, 2, 2, 2) - (t4R_1(2, 2, 2, 2) + t4R_2(2, 2, 2, 2)),
517  "T4R + T4R (2,2,2,2)");
518 
519  t4R_3(i, j, k, l) = t4R_1(i, j, k, l) - t4R_2(i, j, k, l);
520  test_for_zero(t4R_3(0, 0, 0, 0) - (t4R_1(0, 0, 0, 0) - t4R_2(0, 0, 0, 0)),
521  "T4R - T4R (0,0,0,0)");
522  test_for_zero(t4R_3(0, 0, 0, 1) - (t4R_1(0, 0, 0, 1) - t4R_2(0, 0, 0, 1)),
523  "T4R - T4R (0,0,0,1)");
524  test_for_zero(t4R_3(0, 0, 0, 2) - (t4R_1(0, 0, 0, 2) - t4R_2(0, 0, 0, 2)),
525  "T4R - T4R (0,0,0,2)");
526  test_for_zero(t4R_3(0, 0, 1, 0) - (t4R_1(0, 0, 1, 0) - t4R_2(0, 0, 1, 0)),
527  "T4R - T4R (0,0,1,0)");
528  test_for_zero(t4R_3(0, 0, 1, 1) - (t4R_1(0, 0, 1, 1) - t4R_2(0, 0, 1, 1)),
529  "T4R - T4R (0,0,1,1)");
530  test_for_zero(t4R_3(0, 0, 1, 2) - (t4R_1(0, 0, 1, 2) - t4R_2(0, 0, 1, 2)),
531  "T4R - T4R (0,0,1,2)");
532  test_for_zero(t4R_3(0, 0, 2, 0) - (t4R_1(0, 0, 2, 0) - t4R_2(0, 0, 2, 0)),
533  "T4R - T4R (0,0,2,0)");
534  test_for_zero(t4R_3(0, 0, 2, 1) - (t4R_1(0, 0, 2, 1) - t4R_2(0, 0, 2, 1)),
535  "T4R - T4R (0,0,2,1)");
536  test_for_zero(t4R_3(0, 0, 2, 2) - (t4R_1(0, 0, 2, 2) - t4R_2(0, 0, 2, 2)),
537  "T4R - T4R (0,0,2,2)");
538  test_for_zero(t4R_3(0, 1, 0, 0) - (t4R_1(0, 1, 0, 0) - t4R_2(0, 1, 0, 0)),
539  "T4R - T4R (0,1,0,0)");
540  test_for_zero(t4R_3(0, 1, 0, 1) - (t4R_1(0, 1, 0, 1) - t4R_2(0, 1, 0, 1)),
541  "T4R - T4R (0,1,0,1)");
542  test_for_zero(t4R_3(0, 1, 0, 2) - (t4R_1(0, 1, 0, 2) - t4R_2(0, 1, 0, 2)),
543  "T4R - T4R (0,1,0,2)");
544  test_for_zero(t4R_3(0, 1, 1, 0) - (t4R_1(0, 1, 1, 0) - t4R_2(0, 1, 1, 0)),
545  "T4R - T4R (0,1,1,0)");
546  test_for_zero(t4R_3(0, 1, 1, 1) - (t4R_1(0, 1, 1, 1) - t4R_2(0, 1, 1, 1)),
547  "T4R - T4R (0,1,1,1)");
548  test_for_zero(t4R_3(0, 1, 1, 2) - (t4R_1(0, 1, 1, 2) - t4R_2(0, 1, 1, 2)),
549  "T4R - T4R (0,1,1,2)");
550  test_for_zero(t4R_3(0, 1, 2, 0) - (t4R_1(0, 1, 2, 0) - t4R_2(0, 1, 2, 0)),
551  "T4R - T4R (0,1,2,0)");
552  test_for_zero(t4R_3(0, 1, 2, 1) - (t4R_1(0, 1, 2, 1) - t4R_2(0, 1, 2, 1)),
553  "T4R - T4R (0,1,2,1)");
554  test_for_zero(t4R_3(0, 1, 2, 2) - (t4R_1(0, 1, 2, 2) - t4R_2(0, 1, 2, 2)),
555  "T4R - T4R (0,1,2,2)");
556  test_for_zero(t4R_3(0, 2, 0, 0) - (t4R_1(0, 2, 0, 0) - t4R_2(0, 2, 0, 0)),
557  "T4R - T4R (0,2,0,0)");
558  test_for_zero(t4R_3(0, 2, 0, 1) - (t4R_1(0, 2, 0, 1) - t4R_2(0, 2, 0, 1)),
559  "T4R - T4R (0,2,0,1)");
560  test_for_zero(t4R_3(0, 2, 0, 2) - (t4R_1(0, 2, 0, 2) - t4R_2(0, 2, 0, 2)),
561  "T4R - T4R (0,2,0,2)");
562  test_for_zero(t4R_3(0, 2, 1, 0) - (t4R_1(0, 2, 1, 0) - t4R_2(0, 2, 1, 0)),
563  "T4R - T4R (0,2,1,0)");
564  test_for_zero(t4R_3(0, 2, 1, 1) - (t4R_1(0, 2, 1, 1) - t4R_2(0, 2, 1, 1)),
565  "T4R - T4R (0,2,1,1)");
566  test_for_zero(t4R_3(0, 2, 1, 2) - (t4R_1(0, 2, 1, 2) - t4R_2(0, 2, 1, 2)),
567  "T4R - T4R (0,2,1,2)");
568  test_for_zero(t4R_3(0, 2, 2, 0) - (t4R_1(0, 2, 2, 0) - t4R_2(0, 2, 2, 0)),
569  "T4R - T4R (0,2,2,0)");
570  test_for_zero(t4R_3(0, 2, 2, 1) - (t4R_1(0, 2, 2, 1) - t4R_2(0, 2, 2, 1)),
571  "T4R - T4R (0,2,2,1)");
572  test_for_zero(t4R_3(0, 2, 2, 2) - (t4R_1(0, 2, 2, 2) - t4R_2(0, 2, 2, 2)),
573  "T4R - T4R (0,2,2,2)");
574  test_for_zero(t4R_3(1, 0, 0, 0) - (t4R_1(1, 0, 0, 0) - t4R_2(1, 0, 0, 0)),
575  "T4R - T4R (1,0,0,0)");
576  test_for_zero(t4R_3(1, 0, 0, 1) - (t4R_1(1, 0, 0, 1) - t4R_2(1, 0, 0, 1)),
577  "T4R - T4R (1,0,0,1)");
578  test_for_zero(t4R_3(1, 0, 0, 2) - (t4R_1(1, 0, 0, 2) - t4R_2(1, 0, 0, 2)),
579  "T4R - T4R (1,0,0,2)");
580  test_for_zero(t4R_3(1, 0, 1, 0) - (t4R_1(1, 0, 1, 0) - t4R_2(1, 0, 1, 0)),
581  "T4R - T4R (1,0,1,0)");
582  test_for_zero(t4R_3(1, 0, 1, 1) - (t4R_1(1, 0, 1, 1) - t4R_2(1, 0, 1, 1)),
583  "T4R - T4R (1,0,1,1)");
584  test_for_zero(t4R_3(1, 0, 1, 2) - (t4R_1(1, 0, 1, 2) - t4R_2(1, 0, 1, 2)),
585  "T4R - T4R (1,0,1,2)");
586  test_for_zero(t4R_3(1, 0, 2, 0) - (t4R_1(1, 0, 2, 0) - t4R_2(1, 0, 2, 0)),
587  "T4R - T4R (1,0,2,0)");
588  test_for_zero(t4R_3(1, 0, 2, 1) - (t4R_1(1, 0, 2, 1) - t4R_2(1, 0, 2, 1)),
589  "T4R - T4R (1,0,2,1)");
590  test_for_zero(t4R_3(1, 0, 2, 2) - (t4R_1(1, 0, 2, 2) - t4R_2(1, 0, 2, 2)),
591  "T4R - T4R (1,0,2,2)");
592  test_for_zero(t4R_3(1, 1, 0, 0) - (t4R_1(1, 1, 0, 0) - t4R_2(1, 1, 0, 0)),
593  "T4R - T4R (1,1,0,0)");
594  test_for_zero(t4R_3(1, 1, 0, 1) - (t4R_1(1, 1, 0, 1) - t4R_2(1, 1, 0, 1)),
595  "T4R - T4R (1,1,0,1)");
596  test_for_zero(t4R_3(1, 1, 0, 2) - (t4R_1(1, 1, 0, 2) - t4R_2(1, 1, 0, 2)),
597  "T4R - T4R (1,1,0,2)");
598  test_for_zero(t4R_3(1, 1, 1, 0) - (t4R_1(1, 1, 1, 0) - t4R_2(1, 1, 1, 0)),
599  "T4R - T4R (1,1,1,0)");
600  test_for_zero(t4R_3(1, 1, 1, 1) - (t4R_1(1, 1, 1, 1) - t4R_2(1, 1, 1, 1)),
601  "T4R - T4R (1,1,1,1)");
602  test_for_zero(t4R_3(1, 1, 1, 2) - (t4R_1(1, 1, 1, 2) - t4R_2(1, 1, 1, 2)),
603  "T4R - T4R (1,1,1,2)");
604  test_for_zero(t4R_3(1, 1, 2, 0) - (t4R_1(1, 1, 2, 0) - t4R_2(1, 1, 2, 0)),
605  "T4R - T4R (1,1,2,0)");
606  test_for_zero(t4R_3(1, 1, 2, 1) - (t4R_1(1, 1, 2, 1) - t4R_2(1, 1, 2, 1)),
607  "T4R - T4R (1,1,2,1)");
608  test_for_zero(t4R_3(1, 1, 2, 2) - (t4R_1(1, 1, 2, 2) - t4R_2(1, 1, 2, 2)),
609  "T4R - T4R (1,1,2,2)");
610  test_for_zero(t4R_3(1, 2, 0, 0) - (t4R_1(1, 2, 0, 0) - t4R_2(1, 2, 0, 0)),
611  "T4R - T4R (1,2,0,0)");
612  test_for_zero(t4R_3(1, 2, 0, 1) - (t4R_1(1, 2, 0, 1) - t4R_2(1, 2, 0, 1)),
613  "T4R - T4R (1,2,0,1)");
614  test_for_zero(t4R_3(1, 2, 0, 2) - (t4R_1(1, 2, 0, 2) - t4R_2(1, 2, 0, 2)),
615  "T4R - T4R (1,2,0,2)");
616  test_for_zero(t4R_3(1, 2, 1, 0) - (t4R_1(1, 2, 1, 0) - t4R_2(1, 2, 1, 0)),
617  "T4R - T4R (1,2,1,0)");
618  test_for_zero(t4R_3(1, 2, 1, 1) - (t4R_1(1, 2, 1, 1) - t4R_2(1, 2, 1, 1)),
619  "T4R - T4R (1,2,1,1)");
620  test_for_zero(t4R_3(1, 2, 1, 2) - (t4R_1(1, 2, 1, 2) - t4R_2(1, 2, 1, 2)),
621  "T4R - T4R (1,2,1,2)");
622  test_for_zero(t4R_3(1, 2, 2, 0) - (t4R_1(1, 2, 2, 0) - t4R_2(1, 2, 2, 0)),
623  "T4R - T4R (1,2,2,0)");
624  test_for_zero(t4R_3(1, 2, 2, 1) - (t4R_1(1, 2, 2, 1) - t4R_2(1, 2, 2, 1)),
625  "T4R - T4R (1,2,2,1)");
626  test_for_zero(t4R_3(1, 2, 2, 2) - (t4R_1(1, 2, 2, 2) - t4R_2(1, 2, 2, 2)),
627  "T4R - T4R (1,2,2,2)");
628  test_for_zero(t4R_3(2, 0, 0, 0) - (t4R_1(2, 0, 0, 0) - t4R_2(2, 0, 0, 0)),
629  "T4R - T4R (2,0,0,0)");
630  test_for_zero(t4R_3(2, 0, 0, 1) - (t4R_1(2, 0, 0, 1) - t4R_2(2, 0, 0, 1)),
631  "T4R - T4R (2,0,0,1)");
632  test_for_zero(t4R_3(2, 0, 0, 2) - (t4R_1(2, 0, 0, 2) - t4R_2(2, 0, 0, 2)),
633  "T4R - T4R (2,0,0,2)");
634  test_for_zero(t4R_3(2, 0, 1, 0) - (t4R_1(2, 0, 1, 0) - t4R_2(2, 0, 1, 0)),
635  "T4R - T4R (2,0,1,0)");
636  test_for_zero(t4R_3(2, 0, 1, 1) - (t4R_1(2, 0, 1, 1) - t4R_2(2, 0, 1, 1)),
637  "T4R - T4R (2,0,1,1)");
638  test_for_zero(t4R_3(2, 0, 1, 2) - (t4R_1(2, 0, 1, 2) - t4R_2(2, 0, 1, 2)),
639  "T4R - T4R (2,0,1,2)");
640  test_for_zero(t4R_3(2, 0, 2, 0) - (t4R_1(2, 0, 2, 0) - t4R_2(2, 0, 2, 0)),
641  "T4R - T4R (2,0,2,0)");
642  test_for_zero(t4R_3(2, 0, 2, 1) - (t4R_1(2, 0, 2, 1) - t4R_2(2, 0, 2, 1)),
643  "T4R - T4R (2,0,2,1)");
644  test_for_zero(t4R_3(2, 0, 2, 2) - (t4R_1(2, 0, 2, 2) - t4R_2(2, 0, 2, 2)),
645  "T4R - T4R (2,0,2,2)");
646  test_for_zero(t4R_3(2, 1, 0, 0) - (t4R_1(2, 1, 0, 0) - t4R_2(2, 1, 0, 0)),
647  "T4R - T4R (2,1,0,0)");
648  test_for_zero(t4R_3(2, 1, 0, 1) - (t4R_1(2, 1, 0, 1) - t4R_2(2, 1, 0, 1)),
649  "T4R - T4R (2,1,0,1)");
650  test_for_zero(t4R_3(2, 1, 0, 2) - (t4R_1(2, 1, 0, 2) - t4R_2(2, 1, 0, 2)),
651  "T4R - T4R (2,1,0,2)");
652  test_for_zero(t4R_3(2, 1, 1, 0) - (t4R_1(2, 1, 1, 0) - t4R_2(2, 1, 1, 0)),
653  "T4R - T4R (2,1,1,0)");
654  test_for_zero(t4R_3(2, 1, 1, 1) - (t4R_1(2, 1, 1, 1) - t4R_2(2, 1, 1, 1)),
655  "T4R - T4R (2,1,1,1)");
656  test_for_zero(t4R_3(2, 1, 1, 2) - (t4R_1(2, 1, 1, 2) - t4R_2(2, 1, 1, 2)),
657  "T4R - T4R (2,1,1,2)");
658  test_for_zero(t4R_3(2, 1, 2, 0) - (t4R_1(2, 1, 2, 0) - t4R_2(2, 1, 2, 0)),
659  "T4R - T4R (2,1,2,0)");
660  test_for_zero(t4R_3(2, 1, 2, 1) - (t4R_1(2, 1, 2, 1) - t4R_2(2, 1, 2, 1)),
661  "T4R - T4R (2,1,2,1)");
662  test_for_zero(t4R_3(2, 1, 2, 2) - (t4R_1(2, 1, 2, 2) - t4R_2(2, 1, 2, 2)),
663  "T4R - T4R (2,1,2,2)");
664  test_for_zero(t4R_3(2, 2, 0, 0) - (t4R_1(2, 2, 0, 0) - t4R_2(2, 2, 0, 0)),
665  "T4R - T4R (2,2,0,0)");
666  test_for_zero(t4R_3(2, 2, 0, 1) - (t4R_1(2, 2, 0, 1) - t4R_2(2, 2, 0, 1)),
667  "T4R - T4R (2,2,0,1)");
668  test_for_zero(t4R_3(2, 2, 0, 2) - (t4R_1(2, 2, 0, 2) - t4R_2(2, 2, 0, 2)),
669  "T4R - T4R (2,2,0,2)");
670  test_for_zero(t4R_3(2, 2, 1, 0) - (t4R_1(2, 2, 1, 0) - t4R_2(2, 2, 1, 0)),
671  "T4R - T4R (2,2,1,0)");
672  test_for_zero(t4R_3(2, 2, 1, 1) - (t4R_1(2, 2, 1, 1) - t4R_2(2, 2, 1, 1)),
673  "T4R - T4R (2,2,1,1)");
674  test_for_zero(t4R_3(2, 2, 1, 2) - (t4R_1(2, 2, 1, 2) - t4R_2(2, 2, 1, 2)),
675  "T4R - T4R (2,2,1,2)");
676  test_for_zero(t4R_3(2, 2, 2, 0) - (t4R_1(2, 2, 2, 0) - t4R_2(2, 2, 2, 0)),
677  "T4R - T4R (2,2,2,0)");
678  test_for_zero(t4R_3(2, 2, 2, 1) - (t4R_1(2, 2, 2, 1) - t4R_2(2, 2, 2, 1)),
679  "T4R - T4R (2,2,2,1)");
680  test_for_zero(t4R_3(2, 2, 2, 2) - (t4R_1(2, 2, 2, 2) - t4R_2(2, 2, 2, 2)),
681  "T4R - T4R (2,2,2,2)");
682 
683  // t3as_1(i,j,k)=t4R_1(i,j,k,l)*t1_2(l);
684  // t3as_1(i,j,k)=t1_2(l)*t4R_1(i,j,k,l);
685  // t3as_1(i,j,k)=t4R_1(i,j,l,k)*t1_2(l);
686  // t3as_1(i,j,k)=t1_2(l)*t4R_1(i,j,l,k);
687  // t3as_1(i,j,k)=t4R_1(i,l,j,k)*t1_2(l);
688  // t3as_1(i,j,k)=t1_2(l)*t4R_1(i,l,j,k);
689  // t3as_1(i,j,k)=t4R_1(l,i,j,k)*t1_2(l);
690  // t3as_1(i,j,k)=t1_2(l)*t4R_1(l,i,j,k);
691 
692  // cout << '\n';
693 
694  // test_for_zero(t4R_1(i,j,k,l)*t4ddg_2(i,j,k,l)
695  // ); test_for_zero(t4ddg_2(i,j,k,l)*t4R_1(i,j,k,l)
696  // ); test_for_zero(t4R_1(i,j,k,l)*t4ddg_2(i,k,j,l)
697  // ); test_for_zero(t4ddg_2(i,k,j,l)*t4R_1(i,j,k,l)
698  // ); test_for_zero(t4R_3(i,j,k,l)*t4ddg_2(i,j,k,l)
699  // ); test_for_zero(t4ddg_2(i,j,k,l)*t4R_3(i,j,k,l)
700  // ); test_for_zero(t4R_3(i,j,k,l)*t4ddg_2(i,k,j,l)
701  // ); test_for_zero(t4ddg_2(i,k,j,l)*t4R_3(i,j,k,l));
702 
703  // cout << '\n';
704 
705  // t2s_1(j,l)=t4R_1(i,j,k,l)*t2s_2(i,k);
706  // t2s_1(j,l)=t2s_2(i,k)*t4R_1(i,j,k,l);
707 }
FTensor
JSON compatible output.
Definition: Christof_constructor.hpp:6
FTensor::Tensor2_symmetric
Definition: Tensor2_symmetric_value.hpp:13
FTensor::Number< 0 >
test_for_zero
void test_for_zero(const T &t, const std::string &s)
Definition: test_for_zero.hpp:7
FTensor::Riemann
Definition: Riemann.hpp:10
i
FTensor::Index< 'i', SPACE_DIM > i
Definition: hcurl_divergence_operator_2d.cpp:27
FTensor::Index< 'i', 3 >
test_for_zero.hpp
std
Definition: enable_if.hpp:5
j
FTensor::Index< 'j', 3 > j
Definition: matrix_function.cpp:19
test_T4R
void test_T4R()
Definition: test_T4R.cpp:7
k
FTensor::Index< 'k', 3 > k
Definition: matrix_function.cpp:20
l
FTensor::Index< 'l', 3 > l
Definition: matrix_function.cpp:21
Generated by Doxygen 1.8.17 and hosted at University of Glasgow