v0.15.0
Loading...
Searching...
No Matches
test_T4R.cpp
Go to the documentation of this file.
1#include "../../src/FTensor.hpp"
2#include "test_for_zero.hpp"
3#include <iostream>
4using namespace FTensor;
5using namespace std;
6
7void 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::Tensor2_symmetric
Definition Tensor2_symmetric_value.hpp:14
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
FTensor
Tensors class implemented by Walter Landry.
Definition FTensor.hpp:51
test_T4R
void test_T4R()
Definition test_T4R.cpp:7
test_for_zero.hpp
test_for_zero
void test_for_zero(const T &t, const std::string &s)
Definition test_for_zero.hpp:7
Generated by Doxygen 1.9.8 and hosted at University of Glasgow