v0.14.0
test_T4ddg_01.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_T4ddg_01(const Tensor2_symmetric<double, 3> &t2s_2,
8  const Tensor2_symmetric<double, 3> &t2s_3,
9  const Dg<double, 3, 3> &t3dg_2,
10  const Dg<double, 3, 3> &t3dg_3)
11 {
12  Index<'i', 3> i;
13  Index<'j', 3> j;
14  Index<'k', 3> k;
15  Index<'l', 3> l;
16  Index<'m', 3> m;
17 
18  Number<0> N0;
19  Number<1> N1;
20  Number<2> N2;
21 
22  Ddg<double, 3, 3> t4ddg_1, t4ddg_2, t4ddg_3;
23 
24  t4ddg_1(i, j, l, m) = t3dg_2(i, j, k) * t3dg_3(l, m, k);
25  test_for_zero(t4ddg_1(0, 0, 0, 0)
26  - (t3dg_2(0, 0, 0) * t3dg_3(0, 0, 0)
27  + t3dg_2(0, 0, 1) * t3dg_3(0, 0, 1)
28  + t3dg_2(0, 0, 2) * t3dg_3(0, 0, 2)),
29  "T3dg(i,j,k)*T3dg(l,m,k)(0,0,0,0)");
30  test_for_zero(t4ddg_1(0, 0, 0, 1)
31  - (t3dg_2(0, 0, 0) * t3dg_3(0, 1, 0)
32  + t3dg_2(0, 0, 1) * t3dg_3(0, 1, 1)
33  + t3dg_2(0, 0, 2) * t3dg_3(0, 1, 2)),
34  "T3dg(i,j,k)*T3dg(l,m,k)(0,0,0,1)");
35  test_for_zero(t4ddg_1(0, 0, 0, 2)
36  - (t3dg_2(0, 0, 0) * t3dg_3(0, 2, 0)
37  + t3dg_2(0, 0, 1) * t3dg_3(0, 2, 1)
38  + t3dg_2(0, 0, 2) * t3dg_3(0, 2, 2)),
39  "T3dg(i,j,k)*T3dg(l,m,k)(0,0,0,2)");
40  test_for_zero(t4ddg_1(0, 0, 1, 0)
41  - (t3dg_2(0, 0, 0) * t3dg_3(1, 0, 0)
42  + t3dg_2(0, 0, 1) * t3dg_3(1, 0, 1)
43  + t3dg_2(0, 0, 2) * t3dg_3(1, 0, 2)),
44  "T3dg(i,j,k)*T3dg(l,m,k)(0,0,1,0)");
45  test_for_zero(t4ddg_1(0, 0, 1, 1)
46  - (t3dg_2(0, 0, 0) * t3dg_3(1, 1, 0)
47  + t3dg_2(0, 0, 1) * t3dg_3(1, 1, 1)
48  + t3dg_2(0, 0, 2) * t3dg_3(1, 1, 2)),
49  "T3dg(i,j,k)*T3dg(l,m,k)(0,0,1,1)");
50  test_for_zero(t4ddg_1(0, 0, 1, 2)
51  - (t3dg_2(0, 0, 0) * t3dg_3(1, 2, 0)
52  + t3dg_2(0, 0, 1) * t3dg_3(1, 2, 1)
53  + t3dg_2(0, 0, 2) * t3dg_3(1, 2, 2)),
54  "T3dg(i,j,k)*T3dg(l,m,k)(0,0,1,2)");
55  test_for_zero(t4ddg_1(0, 0, 2, 0)
56  - (t3dg_2(0, 0, 0) * t3dg_3(2, 0, 0)
57  + t3dg_2(0, 0, 1) * t3dg_3(2, 0, 1)
58  + t3dg_2(0, 0, 2) * t3dg_3(2, 0, 2)),
59  "T3dg(i,j,k)*T3dg(l,m,k)(0,0,2,0)");
60  test_for_zero(t4ddg_1(0, 0, 2, 1)
61  - (t3dg_2(0, 0, 0) * t3dg_3(2, 1, 0)
62  + t3dg_2(0, 0, 1) * t3dg_3(2, 1, 1)
63  + t3dg_2(0, 0, 2) * t3dg_3(2, 1, 2)),
64  "T3dg(i,j,k)*T3dg(l,m,k)(0,0,2,1)");
65  test_for_zero(t4ddg_1(0, 0, 2, 2)
66  - (t3dg_2(0, 0, 0) * t3dg_3(2, 2, 0)
67  + t3dg_2(0, 0, 1) * t3dg_3(2, 2, 1)
68  + t3dg_2(0, 0, 2) * t3dg_3(2, 2, 2)),
69  "T3dg(i,j,k)*T3dg(l,m,k)(0,0,2,2)");
70  test_for_zero(t4ddg_1(0, 1, 0, 0)
71  - (t3dg_2(0, 1, 0) * t3dg_3(0, 0, 0)
72  + t3dg_2(0, 1, 1) * t3dg_3(0, 0, 1)
73  + t3dg_2(0, 1, 2) * t3dg_3(0, 0, 2)),
74  "T3dg(i,j,k)*T3dg(l,m,k)(0,1,0,0)");
75  test_for_zero(t4ddg_1(0, 1, 0, 1)
76  - (t3dg_2(0, 1, 0) * t3dg_3(0, 1, 0)
77  + t3dg_2(0, 1, 1) * t3dg_3(0, 1, 1)
78  + t3dg_2(0, 1, 2) * t3dg_3(0, 1, 2)),
79  "T3dg(i,j,k)*T3dg(l,m,k)(0,1,0,1)");
80  test_for_zero(t4ddg_1(0, 1, 0, 2)
81  - (t3dg_2(0, 1, 0) * t3dg_3(0, 2, 0)
82  + t3dg_2(0, 1, 1) * t3dg_3(0, 2, 1)
83  + t3dg_2(0, 1, 2) * t3dg_3(0, 2, 2)),
84  "T3dg(i,j,k)*T3dg(l,m,k)(0,1,0,2)");
85  test_for_zero(t4ddg_1(0, 1, 1, 0)
86  - (t3dg_2(0, 1, 0) * t3dg_3(1, 0, 0)
87  + t3dg_2(0, 1, 1) * t3dg_3(1, 0, 1)
88  + t3dg_2(0, 1, 2) * t3dg_3(1, 0, 2)),
89  "T3dg(i,j,k)*T3dg(l,m,k)(0,1,1,0)");
90  test_for_zero(t4ddg_1(0, 1, 1, 1)
91  - (t3dg_2(0, 1, 0) * t3dg_3(1, 1, 0)
92  + t3dg_2(0, 1, 1) * t3dg_3(1, 1, 1)
93  + t3dg_2(0, 1, 2) * t3dg_3(1, 1, 2)),
94  "T3dg(i,j,k)*T3dg(l,m,k)(0,1,1,1)");
95  test_for_zero(t4ddg_1(0, 1, 1, 2)
96  - (t3dg_2(0, 1, 0) * t3dg_3(1, 2, 0)
97  + t3dg_2(0, 1, 1) * t3dg_3(1, 2, 1)
98  + t3dg_2(0, 1, 2) * t3dg_3(1, 2, 2)),
99  "T3dg(i,j,k)*T3dg(l,m,k)(0,1,1,2)");
100  test_for_zero(t4ddg_1(0, 1, 2, 0)
101  - (t3dg_2(0, 1, 0) * t3dg_3(2, 0, 0)
102  + t3dg_2(0, 1, 1) * t3dg_3(2, 0, 1)
103  + t3dg_2(0, 1, 2) * t3dg_3(2, 0, 2)),
104  "T3dg(i,j,k)*T3dg(l,m,k)(0,1,2,0)");
105  test_for_zero(t4ddg_1(0, 1, 2, 1)
106  - (t3dg_2(0, 1, 0) * t3dg_3(2, 1, 0)
107  + t3dg_2(0, 1, 1) * t3dg_3(2, 1, 1)
108  + t3dg_2(0, 1, 2) * t3dg_3(2, 1, 2)),
109  "T3dg(i,j,k)*T3dg(l,m,k)(0,1,2,1)");
110  test_for_zero(t4ddg_1(0, 1, 2, 2)
111  - (t3dg_2(0, 1, 0) * t3dg_3(2, 2, 0)
112  + t3dg_2(0, 1, 1) * t3dg_3(2, 2, 1)
113  + t3dg_2(0, 1, 2) * t3dg_3(2, 2, 2)),
114  "T3dg(i,j,k)*T3dg(l,m,k)(0,1,2,2)");
115  test_for_zero(t4ddg_1(0, 2, 0, 0)
116  - (t3dg_2(0, 2, 0) * t3dg_3(0, 0, 0)
117  + t3dg_2(0, 2, 1) * t3dg_3(0, 0, 1)
118  + t3dg_2(0, 2, 2) * t3dg_3(0, 0, 2)),
119  "T3dg(i,j,k)*T3dg(l,m,k)(0,2,0,0)");
120  test_for_zero(t4ddg_1(0, 2, 0, 1)
121  - (t3dg_2(0, 2, 0) * t3dg_3(0, 1, 0)
122  + t3dg_2(0, 2, 1) * t3dg_3(0, 1, 1)
123  + t3dg_2(0, 2, 2) * t3dg_3(0, 1, 2)),
124  "T3dg(i,j,k)*T3dg(l,m,k)(0,2,0,1)");
125  test_for_zero(t4ddg_1(0, 2, 0, 2)
126  - (t3dg_2(0, 2, 0) * t3dg_3(0, 2, 0)
127  + t3dg_2(0, 2, 1) * t3dg_3(0, 2, 1)
128  + t3dg_2(0, 2, 2) * t3dg_3(0, 2, 2)),
129  "T3dg(i,j,k)*T3dg(l,m,k)(0,2,0,2)");
130  test_for_zero(t4ddg_1(0, 2, 1, 0)
131  - (t3dg_2(0, 2, 0) * t3dg_3(1, 0, 0)
132  + t3dg_2(0, 2, 1) * t3dg_3(1, 0, 1)
133  + t3dg_2(0, 2, 2) * t3dg_3(1, 0, 2)),
134  "T3dg(i,j,k)*T3dg(l,m,k)(0,2,1,0)");
135  test_for_zero(t4ddg_1(0, 2, 1, 1)
136  - (t3dg_2(0, 2, 0) * t3dg_3(1, 1, 0)
137  + t3dg_2(0, 2, 1) * t3dg_3(1, 1, 1)
138  + t3dg_2(0, 2, 2) * t3dg_3(1, 1, 2)),
139  "T3dg(i,j,k)*T3dg(l,m,k)(0,2,1,1)");
140  test_for_zero(t4ddg_1(0, 2, 1, 2)
141  - (t3dg_2(0, 2, 0) * t3dg_3(1, 2, 0)
142  + t3dg_2(0, 2, 1) * t3dg_3(1, 2, 1)
143  + t3dg_2(0, 2, 2) * t3dg_3(1, 2, 2)),
144  "T3dg(i,j,k)*T3dg(l,m,k)(0,2,1,2)");
145  test_for_zero(t4ddg_1(0, 2, 2, 0)
146  - (t3dg_2(0, 2, 0) * t3dg_3(2, 0, 0)
147  + t3dg_2(0, 2, 1) * t3dg_3(2, 0, 1)
148  + t3dg_2(0, 2, 2) * t3dg_3(2, 0, 2)),
149  "T3dg(i,j,k)*T3dg(l,m,k)(0,2,2,0)");
150  test_for_zero(t4ddg_1(0, 2, 2, 1)
151  - (t3dg_2(0, 2, 0) * t3dg_3(2, 1, 0)
152  + t3dg_2(0, 2, 1) * t3dg_3(2, 1, 1)
153  + t3dg_2(0, 2, 2) * t3dg_3(2, 1, 2)),
154  "T3dg(i,j,k)*T3dg(l,m,k)(0,2,2,1)");
155  test_for_zero(t4ddg_1(0, 2, 2, 2)
156  - (t3dg_2(0, 2, 0) * t3dg_3(2, 2, 0)
157  + t3dg_2(0, 2, 1) * t3dg_3(2, 2, 1)
158  + t3dg_2(0, 2, 2) * t3dg_3(2, 2, 2)),
159  "T3dg(i,j,k)*T3dg(l,m,k)(0,2,2,2)");
160  test_for_zero(t4ddg_1(1, 0, 0, 0)
161  - (t3dg_2(1, 0, 0) * t3dg_3(0, 0, 0)
162  + t3dg_2(1, 0, 1) * t3dg_3(0, 0, 1)
163  + t3dg_2(1, 0, 2) * t3dg_3(0, 0, 2)),
164  "T3dg(i,j,k)*T3dg(l,m,k)(1,0,0,0)");
165  test_for_zero(t4ddg_1(1, 0, 0, 1)
166  - (t3dg_2(1, 0, 0) * t3dg_3(0, 1, 0)
167  + t3dg_2(1, 0, 1) * t3dg_3(0, 1, 1)
168  + t3dg_2(1, 0, 2) * t3dg_3(0, 1, 2)),
169  "T3dg(i,j,k)*T3dg(l,m,k)(1,0,0,1)");
170  test_for_zero(t4ddg_1(1, 0, 0, 2)
171  - (t3dg_2(1, 0, 0) * t3dg_3(0, 2, 0)
172  + t3dg_2(1, 0, 1) * t3dg_3(0, 2, 1)
173  + t3dg_2(1, 0, 2) * t3dg_3(0, 2, 2)),
174  "T3dg(i,j,k)*T3dg(l,m,k)(1,0,0,2)");
175  test_for_zero(t4ddg_1(1, 0, 1, 0)
176  - (t3dg_2(1, 0, 0) * t3dg_3(1, 0, 0)
177  + t3dg_2(1, 0, 1) * t3dg_3(1, 0, 1)
178  + t3dg_2(1, 0, 2) * t3dg_3(1, 0, 2)),
179  "T3dg(i,j,k)*T3dg(l,m,k)(1,0,1,0)");
180  test_for_zero(t4ddg_1(1, 0, 1, 1)
181  - (t3dg_2(1, 0, 0) * t3dg_3(1, 1, 0)
182  + t3dg_2(1, 0, 1) * t3dg_3(1, 1, 1)
183  + t3dg_2(1, 0, 2) * t3dg_3(1, 1, 2)),
184  "T3dg(i,j,k)*T3dg(l,m,k)(1,0,1,1)");
185  test_for_zero(t4ddg_1(1, 0, 1, 2)
186  - (t3dg_2(1, 0, 0) * t3dg_3(1, 2, 0)
187  + t3dg_2(1, 0, 1) * t3dg_3(1, 2, 1)
188  + t3dg_2(1, 0, 2) * t3dg_3(1, 2, 2)),
189  "T3dg(i,j,k)*T3dg(l,m,k)(1,0,1,2)");
190  test_for_zero(t4ddg_1(1, 0, 2, 0)
191  - (t3dg_2(1, 0, 0) * t3dg_3(2, 0, 0)
192  + t3dg_2(1, 0, 1) * t3dg_3(2, 0, 1)
193  + t3dg_2(1, 0, 2) * t3dg_3(2, 0, 2)),
194  "T3dg(i,j,k)*T3dg(l,m,k)(1,0,2,0)");
195  test_for_zero(t4ddg_1(1, 0, 2, 1)
196  - (t3dg_2(1, 0, 0) * t3dg_3(2, 1, 0)
197  + t3dg_2(1, 0, 1) * t3dg_3(2, 1, 1)
198  + t3dg_2(1, 0, 2) * t3dg_3(2, 1, 2)),
199  "T3dg(i,j,k)*T3dg(l,m,k)(1,0,2,1)");
200  test_for_zero(t4ddg_1(1, 0, 2, 2)
201  - (t3dg_2(1, 0, 0) * t3dg_3(2, 2, 0)
202  + t3dg_2(1, 0, 1) * t3dg_3(2, 2, 1)
203  + t3dg_2(1, 0, 2) * t3dg_3(2, 2, 2)),
204  "T3dg(i,j,k)*T3dg(l,m,k)(1,0,2,2)");
205  test_for_zero(t4ddg_1(1, 1, 0, 0)
206  - (t3dg_2(1, 1, 0) * t3dg_3(0, 0, 0)
207  + t3dg_2(1, 1, 1) * t3dg_3(0, 0, 1)
208  + t3dg_2(1, 1, 2) * t3dg_3(0, 0, 2)),
209  "T3dg(i,j,k)*T3dg(l,m,k)(1,1,0,0)");
210  test_for_zero(t4ddg_1(1, 1, 0, 1)
211  - (t3dg_2(1, 1, 0) * t3dg_3(0, 1, 0)
212  + t3dg_2(1, 1, 1) * t3dg_3(0, 1, 1)
213  + t3dg_2(1, 1, 2) * t3dg_3(0, 1, 2)),
214  "T3dg(i,j,k)*T3dg(l,m,k)(1,1,0,1)");
215  test_for_zero(t4ddg_1(1, 1, 0, 2)
216  - (t3dg_2(1, 1, 0) * t3dg_3(0, 2, 0)
217  + t3dg_2(1, 1, 1) * t3dg_3(0, 2, 1)
218  + t3dg_2(1, 1, 2) * t3dg_3(0, 2, 2)),
219  "T3dg(i,j,k)*T3dg(l,m,k)(1,1,0,2)");
220  test_for_zero(t4ddg_1(1, 1, 1, 0)
221  - (t3dg_2(1, 1, 0) * t3dg_3(1, 0, 0)
222  + t3dg_2(1, 1, 1) * t3dg_3(1, 0, 1)
223  + t3dg_2(1, 1, 2) * t3dg_3(1, 0, 2)),
224  "T3dg(i,j,k)*T3dg(l,m,k)(1,1,1,0)");
225  test_for_zero(t4ddg_1(1, 1, 1, 1)
226  - (t3dg_2(1, 1, 0) * t3dg_3(1, 1, 0)
227  + t3dg_2(1, 1, 1) * t3dg_3(1, 1, 1)
228  + t3dg_2(1, 1, 2) * t3dg_3(1, 1, 2)),
229  "T3dg(i,j,k)*T3dg(l,m,k)(1,1,1,1)");
230  test_for_zero(t4ddg_1(1, 1, 1, 2)
231  - (t3dg_2(1, 1, 0) * t3dg_3(1, 2, 0)
232  + t3dg_2(1, 1, 1) * t3dg_3(1, 2, 1)
233  + t3dg_2(1, 1, 2) * t3dg_3(1, 2, 2)),
234  "T3dg(i,j,k)*T3dg(l,m,k)(1,1,1,2)");
235  test_for_zero(t4ddg_1(1, 1, 2, 0)
236  - (t3dg_2(1, 1, 0) * t3dg_3(2, 0, 0)
237  + t3dg_2(1, 1, 1) * t3dg_3(2, 0, 1)
238  + t3dg_2(1, 1, 2) * t3dg_3(2, 0, 2)),
239  "T3dg(i,j,k)*T3dg(l,m,k)(1,1,2,0)");
240  test_for_zero(t4ddg_1(1, 1, 2, 1)
241  - (t3dg_2(1, 1, 0) * t3dg_3(2, 1, 0)
242  + t3dg_2(1, 1, 1) * t3dg_3(2, 1, 1)
243  + t3dg_2(1, 1, 2) * t3dg_3(2, 1, 2)),
244  "T3dg(i,j,k)*T3dg(l,m,k)(1,1,2,1)");
245  test_for_zero(t4ddg_1(1, 1, 2, 2)
246  - (t3dg_2(1, 1, 0) * t3dg_3(2, 2, 0)
247  + t3dg_2(1, 1, 1) * t3dg_3(2, 2, 1)
248  + t3dg_2(1, 1, 2) * t3dg_3(2, 2, 2)),
249  "T3dg(i,j,k)*T3dg(l,m,k)(1,1,2,2)");
250  test_for_zero(t4ddg_1(1, 2, 0, 0)
251  - (t3dg_2(1, 2, 0) * t3dg_3(0, 0, 0)
252  + t3dg_2(1, 2, 1) * t3dg_3(0, 0, 1)
253  + t3dg_2(1, 2, 2) * t3dg_3(0, 0, 2)),
254  "T3dg(i,j,k)*T3dg(l,m,k)(1,2,0,0)");
255  test_for_zero(t4ddg_1(1, 2, 0, 1)
256  - (t3dg_2(1, 2, 0) * t3dg_3(0, 1, 0)
257  + t3dg_2(1, 2, 1) * t3dg_3(0, 1, 1)
258  + t3dg_2(1, 2, 2) * t3dg_3(0, 1, 2)),
259  "T3dg(i,j,k)*T3dg(l,m,k)(1,2,0,1)");
260  test_for_zero(t4ddg_1(1, 2, 0, 2)
261  - (t3dg_2(1, 2, 0) * t3dg_3(0, 2, 0)
262  + t3dg_2(1, 2, 1) * t3dg_3(0, 2, 1)
263  + t3dg_2(1, 2, 2) * t3dg_3(0, 2, 2)),
264  "T3dg(i,j,k)*T3dg(l,m,k)(1,2,0,2)");
265  test_for_zero(t4ddg_1(1, 2, 1, 0)
266  - (t3dg_2(1, 2, 0) * t3dg_3(1, 0, 0)
267  + t3dg_2(1, 2, 1) * t3dg_3(1, 0, 1)
268  + t3dg_2(1, 2, 2) * t3dg_3(1, 0, 2)),
269  "T3dg(i,j,k)*T3dg(l,m,k)(1,2,1,0)");
270  test_for_zero(t4ddg_1(1, 2, 1, 1)
271  - (t3dg_2(1, 2, 0) * t3dg_3(1, 1, 0)
272  + t3dg_2(1, 2, 1) * t3dg_3(1, 1, 1)
273  + t3dg_2(1, 2, 2) * t3dg_3(1, 1, 2)),
274  "T3dg(i,j,k)*T3dg(l,m,k)(1,2,1,1)");
275  test_for_zero(t4ddg_1(1, 2, 1, 2)
276  - (t3dg_2(1, 2, 0) * t3dg_3(1, 2, 0)
277  + t3dg_2(1, 2, 1) * t3dg_3(1, 2, 1)
278  + t3dg_2(1, 2, 2) * t3dg_3(1, 2, 2)),
279  "T3dg(i,j,k)*T3dg(l,m,k)(1,2,1,2)");
280  test_for_zero(t4ddg_1(1, 2, 2, 0)
281  - (t3dg_2(1, 2, 0) * t3dg_3(2, 0, 0)
282  + t3dg_2(1, 2, 1) * t3dg_3(2, 0, 1)
283  + t3dg_2(1, 2, 2) * t3dg_3(2, 0, 2)),
284  "T3dg(i,j,k)*T3dg(l,m,k)(1,2,2,0)");
285  test_for_zero(t4ddg_1(1, 2, 2, 1)
286  - (t3dg_2(1, 2, 0) * t3dg_3(2, 1, 0)
287  + t3dg_2(1, 2, 1) * t3dg_3(2, 1, 1)
288  + t3dg_2(1, 2, 2) * t3dg_3(2, 1, 2)),
289  "T3dg(i,j,k)*T3dg(l,m,k)(1,2,2,1)");
290  test_for_zero(t4ddg_1(1, 2, 2, 2)
291  - (t3dg_2(1, 2, 0) * t3dg_3(2, 2, 0)
292  + t3dg_2(1, 2, 1) * t3dg_3(2, 2, 1)
293  + t3dg_2(1, 2, 2) * t3dg_3(2, 2, 2)),
294  "T3dg(i,j,k)*T3dg(l,m,k)(1,2,2,2)");
295  test_for_zero(t4ddg_1(2, 0, 0, 0)
296  - (t3dg_2(2, 0, 0) * t3dg_3(0, 0, 0)
297  + t3dg_2(2, 0, 1) * t3dg_3(0, 0, 1)
298  + t3dg_2(2, 0, 2) * t3dg_3(0, 0, 2)),
299  "T3dg(i,j,k)*T3dg(l,m,k)(2,0,0,0)");
300  test_for_zero(t4ddg_1(2, 0, 0, 1)
301  - (t3dg_2(2, 0, 0) * t3dg_3(0, 1, 0)
302  + t3dg_2(2, 0, 1) * t3dg_3(0, 1, 1)
303  + t3dg_2(2, 0, 2) * t3dg_3(0, 1, 2)),
304  "T3dg(i,j,k)*T3dg(l,m,k)(2,0,0,1)");
305  test_for_zero(t4ddg_1(2, 0, 0, 2)
306  - (t3dg_2(2, 0, 0) * t3dg_3(0, 2, 0)
307  + t3dg_2(2, 0, 1) * t3dg_3(0, 2, 1)
308  + t3dg_2(2, 0, 2) * t3dg_3(0, 2, 2)),
309  "T3dg(i,j,k)*T3dg(l,m,k)(2,0,0,2)");
310  test_for_zero(t4ddg_1(2, 0, 1, 0)
311  - (t3dg_2(2, 0, 0) * t3dg_3(1, 0, 0)
312  + t3dg_2(2, 0, 1) * t3dg_3(1, 0, 1)
313  + t3dg_2(2, 0, 2) * t3dg_3(1, 0, 2)),
314  "T3dg(i,j,k)*T3dg(l,m,k)(2,0,1,0)");
315  test_for_zero(t4ddg_1(2, 0, 1, 1)
316  - (t3dg_2(2, 0, 0) * t3dg_3(1, 1, 0)
317  + t3dg_2(2, 0, 1) * t3dg_3(1, 1, 1)
318  + t3dg_2(2, 0, 2) * t3dg_3(1, 1, 2)),
319  "T3dg(i,j,k)*T3dg(l,m,k)(2,0,1,1)");
320  test_for_zero(t4ddg_1(2, 0, 1, 2)
321  - (t3dg_2(2, 0, 0) * t3dg_3(1, 2, 0)
322  + t3dg_2(2, 0, 1) * t3dg_3(1, 2, 1)
323  + t3dg_2(2, 0, 2) * t3dg_3(1, 2, 2)),
324  "T3dg(i,j,k)*T3dg(l,m,k)(2,0,1,2)");
325  test_for_zero(t4ddg_1(2, 0, 2, 0)
326  - (t3dg_2(2, 0, 0) * t3dg_3(2, 0, 0)
327  + t3dg_2(2, 0, 1) * t3dg_3(2, 0, 1)
328  + t3dg_2(2, 0, 2) * t3dg_3(2, 0, 2)),
329  "T3dg(i,j,k)*T3dg(l,m,k)(2,0,2,0)");
330  test_for_zero(t4ddg_1(2, 0, 2, 1)
331  - (t3dg_2(2, 0, 0) * t3dg_3(2, 1, 0)
332  + t3dg_2(2, 0, 1) * t3dg_3(2, 1, 1)
333  + t3dg_2(2, 0, 2) * t3dg_3(2, 1, 2)),
334  "T3dg(i,j,k)*T3dg(l,m,k)(2,0,2,1)");
335  test_for_zero(t4ddg_1(2, 0, 2, 2)
336  - (t3dg_2(2, 0, 0) * t3dg_3(2, 2, 0)
337  + t3dg_2(2, 0, 1) * t3dg_3(2, 2, 1)
338  + t3dg_2(2, 0, 2) * t3dg_3(2, 2, 2)),
339  "T3dg(i,j,k)*T3dg(l,m,k)(2,0,2,2)");
340  test_for_zero(t4ddg_1(2, 1, 0, 0)
341  - (t3dg_2(2, 1, 0) * t3dg_3(0, 0, 0)
342  + t3dg_2(2, 1, 1) * t3dg_3(0, 0, 1)
343  + t3dg_2(2, 1, 2) * t3dg_3(0, 0, 2)),
344  "T3dg(i,j,k)*T3dg(l,m,k)(2,1,0,0)");
345  test_for_zero(t4ddg_1(2, 1, 0, 1)
346  - (t3dg_2(2, 1, 0) * t3dg_3(0, 1, 0)
347  + t3dg_2(2, 1, 1) * t3dg_3(0, 1, 1)
348  + t3dg_2(2, 1, 2) * t3dg_3(0, 1, 2)),
349  "T3dg(i,j,k)*T3dg(l,m,k)(2,1,0,1)");
350  test_for_zero(t4ddg_1(2, 1, 0, 2)
351  - (t3dg_2(2, 1, 0) * t3dg_3(0, 2, 0)
352  + t3dg_2(2, 1, 1) * t3dg_3(0, 2, 1)
353  + t3dg_2(2, 1, 2) * t3dg_3(0, 2, 2)),
354  "T3dg(i,j,k)*T3dg(l,m,k)(2,1,0,2)");
355  test_for_zero(t4ddg_1(2, 1, 1, 0)
356  - (t3dg_2(2, 1, 0) * t3dg_3(1, 0, 0)
357  + t3dg_2(2, 1, 1) * t3dg_3(1, 0, 1)
358  + t3dg_2(2, 1, 2) * t3dg_3(1, 0, 2)),
359  "T3dg(i,j,k)*T3dg(l,m,k)(2,1,1,0)");
360  test_for_zero(t4ddg_1(2, 1, 1, 1)
361  - (t3dg_2(2, 1, 0) * t3dg_3(1, 1, 0)
362  + t3dg_2(2, 1, 1) * t3dg_3(1, 1, 1)
363  + t3dg_2(2, 1, 2) * t3dg_3(1, 1, 2)),
364  "T3dg(i,j,k)*T3dg(l,m,k)(2,1,1,1)");
365  test_for_zero(t4ddg_1(2, 1, 1, 2)
366  - (t3dg_2(2, 1, 0) * t3dg_3(1, 2, 0)
367  + t3dg_2(2, 1, 1) * t3dg_3(1, 2, 1)
368  + t3dg_2(2, 1, 2) * t3dg_3(1, 2, 2)),
369  "T3dg(i,j,k)*T3dg(l,m,k)(2,1,1,2)");
370  test_for_zero(t4ddg_1(2, 1, 2, 0)
371  - (t3dg_2(2, 1, 0) * t3dg_3(2, 0, 0)
372  + t3dg_2(2, 1, 1) * t3dg_3(2, 0, 1)
373  + t3dg_2(2, 1, 2) * t3dg_3(2, 0, 2)),
374  "T3dg(i,j,k)*T3dg(l,m,k)(2,1,2,0)");
375  test_for_zero(t4ddg_1(2, 1, 2, 1)
376  - (t3dg_2(2, 1, 0) * t3dg_3(2, 1, 0)
377  + t3dg_2(2, 1, 1) * t3dg_3(2, 1, 1)
378  + t3dg_2(2, 1, 2) * t3dg_3(2, 1, 2)),
379  "T3dg(i,j,k)*T3dg(l,m,k)(2,1,2,1)");
380  test_for_zero(t4ddg_1(2, 1, 2, 2)
381  - (t3dg_2(2, 1, 0) * t3dg_3(2, 2, 0)
382  + t3dg_2(2, 1, 1) * t3dg_3(2, 2, 1)
383  + t3dg_2(2, 1, 2) * t3dg_3(2, 2, 2)),
384  "T3dg(i,j,k)*T3dg(l,m,k)(2,1,2,2)");
385  test_for_zero(t4ddg_1(2, 2, 0, 0)
386  - (t3dg_2(2, 2, 0) * t3dg_3(0, 0, 0)
387  + t3dg_2(2, 2, 1) * t3dg_3(0, 0, 1)
388  + t3dg_2(2, 2, 2) * t3dg_3(0, 0, 2)),
389  "T3dg(i,j,k)*T3dg(l,m,k)(2,2,0,0)");
390  test_for_zero(t4ddg_1(2, 2, 0, 1)
391  - (t3dg_2(2, 2, 0) * t3dg_3(0, 1, 0)
392  + t3dg_2(2, 2, 1) * t3dg_3(0, 1, 1)
393  + t3dg_2(2, 2, 2) * t3dg_3(0, 1, 2)),
394  "T3dg(i,j,k)*T3dg(l,m,k)(2,2,0,1)");
395  test_for_zero(t4ddg_1(2, 2, 0, 2)
396  - (t3dg_2(2, 2, 0) * t3dg_3(0, 2, 0)
397  + t3dg_2(2, 2, 1) * t3dg_3(0, 2, 1)
398  + t3dg_2(2, 2, 2) * t3dg_3(0, 2, 2)),
399  "T3dg(i,j,k)*T3dg(l,m,k)(2,2,0,2)");
400  test_for_zero(t4ddg_1(2, 2, 1, 0)
401  - (t3dg_2(2, 2, 0) * t3dg_3(1, 0, 0)
402  + t3dg_2(2, 2, 1) * t3dg_3(1, 0, 1)
403  + t3dg_2(2, 2, 2) * t3dg_3(1, 0, 2)),
404  "T3dg(i,j,k)*T3dg(l,m,k)(2,2,1,0)");
405  test_for_zero(t4ddg_1(2, 2, 1, 1)
406  - (t3dg_2(2, 2, 0) * t3dg_3(1, 1, 0)
407  + t3dg_2(2, 2, 1) * t3dg_3(1, 1, 1)
408  + t3dg_2(2, 2, 2) * t3dg_3(1, 1, 2)),
409  "T3dg(i,j,k)*T3dg(l,m,k)(2,2,1,1)");
410  test_for_zero(t4ddg_1(2, 2, 1, 2)
411  - (t3dg_2(2, 2, 0) * t3dg_3(1, 2, 0)
412  + t3dg_2(2, 2, 1) * t3dg_3(1, 2, 1)
413  + t3dg_2(2, 2, 2) * t3dg_3(1, 2, 2)),
414  "T3dg(i,j,k)*T3dg(l,m,k)(2,2,1,2)");
415  test_for_zero(t4ddg_1(2, 2, 2, 0)
416  - (t3dg_2(2, 2, 0) * t3dg_3(2, 0, 0)
417  + t3dg_2(2, 2, 1) * t3dg_3(2, 0, 1)
418  + t3dg_2(2, 2, 2) * t3dg_3(2, 0, 2)),
419  "T3dg(i,j,k)*T3dg(l,m,k)(2,2,2,0)");
420  test_for_zero(t4ddg_1(2, 2, 2, 1)
421  - (t3dg_2(2, 2, 0) * t3dg_3(2, 1, 0)
422  + t3dg_2(2, 2, 1) * t3dg_3(2, 1, 1)
423  + t3dg_2(2, 2, 2) * t3dg_3(2, 1, 2)),
424  "T3dg(i,j,k)*T3dg(l,m,k)(2,2,2,1)");
425  test_for_zero(t4ddg_1(2, 2, 2, 2)
426  - (t3dg_2(2, 2, 0) * t3dg_3(2, 2, 0)
427  + t3dg_2(2, 2, 1) * t3dg_3(2, 2, 1)
428  + t3dg_2(2, 2, 2) * t3dg_3(2, 2, 2)),
429  "T3dg(i,j,k)*T3dg(l,m,k)(2,2,2,2)");
430 
431  t4ddg_2(i, j, l, m) = t2s_2(i, j) * t2s_3(l, m);
432  t4ddg_3(i, j, l, m) = t2s_3(i, j) * t2s_2(l, m);
433  test_for_zero(t4ddg_2(0, 0, 0, 0) - (t2s_2(0, 0) * t2s_3(0, 0)),
434  "T2s(i,j)*T2s(l,m)(0,0,0,0)");
435  test_for_zero(t4ddg_2(0, 0, 0, 1) - (t2s_2(0, 0) * t2s_3(0, 1)),
436  "T2s(i,j)*T2s(l,m)(0,0,0,1)");
437  test_for_zero(t4ddg_2(0, 0, 0, 2) - (t2s_2(0, 0) * t2s_3(0, 2)),
438  "T2s(i,j)*T2s(l,m)(0,0,0,2)");
439  test_for_zero(t4ddg_2(0, 0, 1, 0) - (t2s_2(0, 0) * t2s_3(1, 0)),
440  "T2s(i,j)*T2s(l,m)(0,0,1,0)");
441  test_for_zero(t4ddg_2(0, 0, 1, 1) - (t2s_2(0, 0) * t2s_3(1, 1)),
442  "T2s(i,j)*T2s(l,m)(0,0,1,1)");
443  test_for_zero(t4ddg_2(0, 0, 1, 2) - (t2s_2(0, 0) * t2s_3(1, 2)),
444  "T2s(i,j)*T2s(l,m)(0,0,1,2)");
445  test_for_zero(t4ddg_2(0, 0, 2, 0) - (t2s_2(0, 0) * t2s_3(2, 0)),
446  "T2s(i,j)*T2s(l,m)(0,0,2,0)");
447  test_for_zero(t4ddg_2(0, 0, 2, 1) - (t2s_2(0, 0) * t2s_3(2, 1)),
448  "T2s(i,j)*T2s(l,m)(0,0,2,1)");
449  test_for_zero(t4ddg_2(0, 0, 2, 2) - (t2s_2(0, 0) * t2s_3(2, 2)),
450  "T2s(i,j)*T2s(l,m)(0,0,2,2)");
451  test_for_zero(t4ddg_2(0, 1, 0, 0) - (t2s_2(0, 1) * t2s_3(0, 0)),
452  "T2s(i,j)*T2s(l,m)(0,1,0,0)");
453  test_for_zero(t4ddg_2(0, 1, 0, 1) - (t2s_2(0, 1) * t2s_3(0, 1)),
454  "T2s(i,j)*T2s(l,m)(0,1,0,1)");
455  test_for_zero(t4ddg_2(0, 1, 0, 2) - (t2s_2(0, 1) * t2s_3(0, 2)),
456  "T2s(i,j)*T2s(l,m)(0,1,0,2)");
457  test_for_zero(t4ddg_2(0, 1, 1, 0) - (t2s_2(0, 1) * t2s_3(1, 0)),
458  "T2s(i,j)*T2s(l,m)(0,1,1,0)");
459  test_for_zero(t4ddg_2(0, 1, 1, 1) - (t2s_2(0, 1) * t2s_3(1, 1)),
460  "T2s(i,j)*T2s(l,m)(0,1,1,1)");
461  test_for_zero(t4ddg_2(0, 1, 1, 2) - (t2s_2(0, 1) * t2s_3(1, 2)),
462  "T2s(i,j)*T2s(l,m)(0,1,1,2)");
463  test_for_zero(t4ddg_2(0, 1, 2, 0) - (t2s_2(0, 1) * t2s_3(2, 0)),
464  "T2s(i,j)*T2s(l,m)(0,1,2,0)");
465  test_for_zero(t4ddg_2(0, 1, 2, 1) - (t2s_2(0, 1) * t2s_3(2, 1)),
466  "T2s(i,j)*T2s(l,m)(0,1,2,1)");
467  test_for_zero(t4ddg_2(0, 1, 2, 2) - (t2s_2(0, 1) * t2s_3(2, 2)),
468  "T2s(i,j)*T2s(l,m)(0,1,2,2)");
469  test_for_zero(t4ddg_2(0, 2, 0, 0) - (t2s_2(0, 2) * t2s_3(0, 0)),
470  "T2s(i,j)*T2s(l,m)(0,2,0,0)");
471  test_for_zero(t4ddg_2(0, 2, 0, 1) - (t2s_2(0, 2) * t2s_3(0, 1)),
472  "T2s(i,j)*T2s(l,m)(0,2,0,1)");
473  test_for_zero(t4ddg_2(0, 2, 0, 2) - (t2s_2(0, 2) * t2s_3(0, 2)),
474  "T2s(i,j)*T2s(l,m)(0,2,0,2)");
475  test_for_zero(t4ddg_2(0, 2, 1, 0) - (t2s_2(0, 2) * t2s_3(1, 0)),
476  "T2s(i,j)*T2s(l,m)(0,2,1,0)");
477  test_for_zero(t4ddg_2(0, 2, 1, 1) - (t2s_2(0, 2) * t2s_3(1, 1)),
478  "T2s(i,j)*T2s(l,m)(0,2,1,1)");
479  test_for_zero(t4ddg_2(0, 2, 1, 2) - (t2s_2(0, 2) * t2s_3(1, 2)),
480  "T2s(i,j)*T2s(l,m)(0,2,1,2)");
481  test_for_zero(t4ddg_2(0, 2, 2, 0) - (t2s_2(0, 2) * t2s_3(2, 0)),
482  "T2s(i,j)*T2s(l,m)(0,2,2,0)");
483  test_for_zero(t4ddg_2(0, 2, 2, 1) - (t2s_2(0, 2) * t2s_3(2, 1)),
484  "T2s(i,j)*T2s(l,m)(0,2,2,1)");
485  test_for_zero(t4ddg_2(0, 2, 2, 2) - (t2s_2(0, 2) * t2s_3(2, 2)),
486  "T2s(i,j)*T2s(l,m)(0,2,2,2)");
487  test_for_zero(t4ddg_2(1, 0, 0, 0) - (t2s_2(1, 0) * t2s_3(0, 0)),
488  "T2s(i,j)*T2s(l,m)(1,0,0,0)");
489  test_for_zero(t4ddg_2(1, 0, 0, 1) - (t2s_2(1, 0) * t2s_3(0, 1)),
490  "T2s(i,j)*T2s(l,m)(1,0,0,1)");
491  test_for_zero(t4ddg_2(1, 0, 0, 2) - (t2s_2(1, 0) * t2s_3(0, 2)),
492  "T2s(i,j)*T2s(l,m)(1,0,0,2)");
493  test_for_zero(t4ddg_2(1, 0, 1, 0) - (t2s_2(1, 0) * t2s_3(1, 0)),
494  "T2s(i,j)*T2s(l,m)(1,0,1,0)");
495  test_for_zero(t4ddg_2(1, 0, 1, 1) - (t2s_2(1, 0) * t2s_3(1, 1)),
496  "T2s(i,j)*T2s(l,m)(1,0,1,1)");
497  test_for_zero(t4ddg_2(1, 0, 1, 2) - (t2s_2(1, 0) * t2s_3(1, 2)),
498  "T2s(i,j)*T2s(l,m)(1,0,1,2)");
499  test_for_zero(t4ddg_2(1, 0, 2, 0) - (t2s_2(1, 0) * t2s_3(2, 0)),
500  "T2s(i,j)*T2s(l,m)(1,0,2,0)");
501  test_for_zero(t4ddg_2(1, 0, 2, 1) - (t2s_2(1, 0) * t2s_3(2, 1)),
502  "T2s(i,j)*T2s(l,m)(1,0,2,1)");
503  test_for_zero(t4ddg_2(1, 0, 2, 2) - (t2s_2(1, 0) * t2s_3(2, 2)),
504  "T2s(i,j)*T2s(l,m)(1,0,2,2)");
505  test_for_zero(t4ddg_2(1, 1, 0, 0) - (t2s_2(1, 1) * t2s_3(0, 0)),
506  "T2s(i,j)*T2s(l,m)(1,1,0,0)");
507  test_for_zero(t4ddg_2(1, 1, 0, 1) - (t2s_2(1, 1) * t2s_3(0, 1)),
508  "T2s(i,j)*T2s(l,m)(1,1,0,1)");
509  test_for_zero(t4ddg_2(1, 1, 0, 2) - (t2s_2(1, 1) * t2s_3(0, 2)),
510  "T2s(i,j)*T2s(l,m)(1,1,0,2)");
511  test_for_zero(t4ddg_2(1, 1, 1, 0) - (t2s_2(1, 1) * t2s_3(1, 0)),
512  "T2s(i,j)*T2s(l,m)(1,1,1,0)");
513  test_for_zero(t4ddg_2(1, 1, 1, 1) - (t2s_2(1, 1) * t2s_3(1, 1)),
514  "T2s(i,j)*T2s(l,m)(1,1,1,1)");
515  test_for_zero(t4ddg_2(1, 1, 1, 2) - (t2s_2(1, 1) * t2s_3(1, 2)),
516  "T2s(i,j)*T2s(l,m)(1,1,1,2)");
517  test_for_zero(t4ddg_2(1, 1, 2, 0) - (t2s_2(1, 1) * t2s_3(2, 0)),
518  "T2s(i,j)*T2s(l,m)(1,1,2,0)");
519  test_for_zero(t4ddg_2(1, 1, 2, 1) - (t2s_2(1, 1) * t2s_3(2, 1)),
520  "T2s(i,j)*T2s(l,m)(1,1,2,1)");
521  test_for_zero(t4ddg_2(1, 1, 2, 2) - (t2s_2(1, 1) * t2s_3(2, 2)),
522  "T2s(i,j)*T2s(l,m)(1,1,2,2)");
523  test_for_zero(t4ddg_2(1, 2, 0, 0) - (t2s_2(1, 2) * t2s_3(0, 0)),
524  "T2s(i,j)*T2s(l,m)(1,2,0,0)");
525  test_for_zero(t4ddg_2(1, 2, 0, 1) - (t2s_2(1, 2) * t2s_3(0, 1)),
526  "T2s(i,j)*T2s(l,m)(1,2,0,1)");
527  test_for_zero(t4ddg_2(1, 2, 0, 2) - (t2s_2(1, 2) * t2s_3(0, 2)),
528  "T2s(i,j)*T2s(l,m)(1,2,0,2)");
529  test_for_zero(t4ddg_2(1, 2, 1, 0) - (t2s_2(1, 2) * t2s_3(1, 0)),
530  "T2s(i,j)*T2s(l,m)(1,2,1,0)");
531  test_for_zero(t4ddg_2(1, 2, 1, 1) - (t2s_2(1, 2) * t2s_3(1, 1)),
532  "T2s(i,j)*T2s(l,m)(1,2,1,1)");
533  test_for_zero(t4ddg_2(1, 2, 1, 2) - (t2s_2(1, 2) * t2s_3(1, 2)),
534  "T2s(i,j)*T2s(l,m)(1,2,1,2)");
535  test_for_zero(t4ddg_2(1, 2, 2, 0) - (t2s_2(1, 2) * t2s_3(2, 0)),
536  "T2s(i,j)*T2s(l,m)(1,2,2,0)");
537  test_for_zero(t4ddg_2(1, 2, 2, 1) - (t2s_2(1, 2) * t2s_3(2, 1)),
538  "T2s(i,j)*T2s(l,m)(1,2,2,1)");
539  test_for_zero(t4ddg_2(1, 2, 2, 2) - (t2s_2(1, 2) * t2s_3(2, 2)),
540  "T2s(i,j)*T2s(l,m)(1,2,2,2)");
541  test_for_zero(t4ddg_2(2, 0, 0, 0) - (t2s_2(2, 0) * t2s_3(0, 0)),
542  "T2s(i,j)*T2s(l,m)(2,0,0,0)");
543  test_for_zero(t4ddg_2(2, 0, 0, 1) - (t2s_2(2, 0) * t2s_3(0, 1)),
544  "T2s(i,j)*T2s(l,m)(2,0,0,1)");
545  test_for_zero(t4ddg_2(2, 0, 0, 2) - (t2s_2(2, 0) * t2s_3(0, 2)),
546  "T2s(i,j)*T2s(l,m)(2,0,0,2)");
547  test_for_zero(t4ddg_2(2, 0, 1, 0) - (t2s_2(2, 0) * t2s_3(1, 0)),
548  "T2s(i,j)*T2s(l,m)(2,0,1,0)");
549  test_for_zero(t4ddg_2(2, 0, 1, 1) - (t2s_2(2, 0) * t2s_3(1, 1)),
550  "T2s(i,j)*T2s(l,m)(2,0,1,1)");
551  test_for_zero(t4ddg_2(2, 0, 1, 2) - (t2s_2(2, 0) * t2s_3(1, 2)),
552  "T2s(i,j)*T2s(l,m)(2,0,1,2)");
553  test_for_zero(t4ddg_2(2, 0, 2, 0) - (t2s_2(2, 0) * t2s_3(2, 0)),
554  "T2s(i,j)*T2s(l,m)(2,0,2,0)");
555  test_for_zero(t4ddg_2(2, 0, 2, 1) - (t2s_2(2, 0) * t2s_3(2, 1)),
556  "T2s(i,j)*T2s(l,m)(2,0,2,1)");
557  test_for_zero(t4ddg_2(2, 0, 2, 2) - (t2s_2(2, 0) * t2s_3(2, 2)),
558  "T2s(i,j)*T2s(l,m)(2,0,2,2)");
559  test_for_zero(t4ddg_2(2, 1, 0, 0) - (t2s_2(2, 1) * t2s_3(0, 0)),
560  "T2s(i,j)*T2s(l,m)(2,1,0,0)");
561  test_for_zero(t4ddg_2(2, 1, 0, 1) - (t2s_2(2, 1) * t2s_3(0, 1)),
562  "T2s(i,j)*T2s(l,m)(2,1,0,1)");
563  test_for_zero(t4ddg_2(2, 1, 0, 2) - (t2s_2(2, 1) * t2s_3(0, 2)),
564  "T2s(i,j)*T2s(l,m)(2,1,0,2)");
565  test_for_zero(t4ddg_2(2, 1, 1, 0) - (t2s_2(2, 1) * t2s_3(1, 0)),
566  "T2s(i,j)*T2s(l,m)(2,1,1,0)");
567  test_for_zero(t4ddg_2(2, 1, 1, 1) - (t2s_2(2, 1) * t2s_3(1, 1)),
568  "T2s(i,j)*T2s(l,m)(2,1,1,1)");
569  test_for_zero(t4ddg_2(2, 1, 1, 2) - (t2s_2(2, 1) * t2s_3(1, 2)),
570  "T2s(i,j)*T2s(l,m)(2,1,1,2)");
571  test_for_zero(t4ddg_2(2, 1, 2, 0) - (t2s_2(2, 1) * t2s_3(2, 0)),
572  "T2s(i,j)*T2s(l,m)(2,1,2,0)");
573  test_for_zero(t4ddg_2(2, 1, 2, 1) - (t2s_2(2, 1) * t2s_3(2, 1)),
574  "T2s(i,j)*T2s(l,m)(2,1,2,1)");
575  test_for_zero(t4ddg_2(2, 1, 2, 2) - (t2s_2(2, 1) * t2s_3(2, 2)),
576  "T2s(i,j)*T2s(l,m)(2,1,2,2)");
577  test_for_zero(t4ddg_2(2, 2, 0, 0) - (t2s_2(2, 2) * t2s_3(0, 0)),
578  "T2s(i,j)*T2s(l,m)(2,2,0,0)");
579  test_for_zero(t4ddg_2(2, 2, 0, 1) - (t2s_2(2, 2) * t2s_3(0, 1)),
580  "T2s(i,j)*T2s(l,m)(2,2,0,1)");
581  test_for_zero(t4ddg_2(2, 2, 0, 2) - (t2s_2(2, 2) * t2s_3(0, 2)),
582  "T2s(i,j)*T2s(l,m)(2,2,0,2)");
583  test_for_zero(t4ddg_2(2, 2, 1, 0) - (t2s_2(2, 2) * t2s_3(1, 0)),
584  "T2s(i,j)*T2s(l,m)(2,2,1,0)");
585  test_for_zero(t4ddg_2(2, 2, 1, 1) - (t2s_2(2, 2) * t2s_3(1, 1)),
586  "T2s(i,j)*T2s(l,m)(2,2,1,1)");
587  test_for_zero(t4ddg_2(2, 2, 1, 2) - (t2s_2(2, 2) * t2s_3(1, 2)),
588  "T2s(i,j)*T2s(l,m)(2,2,1,2)");
589  test_for_zero(t4ddg_2(2, 2, 2, 0) - (t2s_2(2, 2) * t2s_3(2, 0)),
590  "T2s(i,j)*T2s(l,m)(2,2,2,0)");
591  test_for_zero(t4ddg_2(2, 2, 2, 1) - (t2s_2(2, 2) * t2s_3(2, 1)),
592  "T2s(i,j)*T2s(l,m)(2,2,2,1)");
593  test_for_zero(t4ddg_2(2, 2, 2, 2) - (t2s_2(2, 2) * t2s_3(2, 2)),
594  "T2s(i,j)*T2s(l,m)(2,2,2,2)");
595 
596  test_for_zero(t4ddg_1(i, j, k, l) * t4ddg_2(i, k, j, l)
597  - t4ddg_1(0, 0, 0, 0) * t4ddg_2(0, 0, 0, 0)
598  - t4ddg_1(0, 0, 0, 1) * t4ddg_2(0, 0, 0, 1)
599  - t4ddg_1(0, 0, 0, 2) * t4ddg_2(0, 0, 0, 2)
600  - t4ddg_1(0, 0, 1, 0) * t4ddg_2(0, 1, 0, 0)
601  - t4ddg_1(0, 0, 1, 1) * t4ddg_2(0, 1, 0, 1)
602  - t4ddg_1(0, 0, 1, 2) * t4ddg_2(0, 1, 0, 2)
603  - t4ddg_1(0, 0, 2, 0) * t4ddg_2(0, 2, 0, 0)
604  - t4ddg_1(0, 0, 2, 1) * t4ddg_2(0, 2, 0, 1)
605  - t4ddg_1(0, 0, 2, 2) * t4ddg_2(0, 2, 0, 2)
606  - t4ddg_1(0, 1, 0, 0) * t4ddg_2(0, 0, 1, 0)
607  - t4ddg_1(0, 1, 0, 1) * t4ddg_2(0, 0, 1, 1)
608  - t4ddg_1(0, 1, 0, 2) * t4ddg_2(0, 0, 1, 2)
609  - t4ddg_1(0, 1, 1, 0) * t4ddg_2(0, 1, 1, 0)
610  - t4ddg_1(0, 1, 1, 1) * t4ddg_2(0, 1, 1, 1)
611  - t4ddg_1(0, 1, 1, 2) * t4ddg_2(0, 1, 1, 2)
612  - t4ddg_1(0, 1, 2, 0) * t4ddg_2(0, 2, 1, 0)
613  - t4ddg_1(0, 1, 2, 1) * t4ddg_2(0, 2, 1, 1)
614  - t4ddg_1(0, 1, 2, 2) * t4ddg_2(0, 2, 1, 2)
615  - t4ddg_1(0, 2, 0, 0) * t4ddg_2(0, 0, 2, 0)
616  - t4ddg_1(0, 2, 0, 1) * t4ddg_2(0, 0, 2, 1)
617  - t4ddg_1(0, 2, 0, 2) * t4ddg_2(0, 0, 2, 2)
618  - t4ddg_1(0, 2, 1, 0) * t4ddg_2(0, 1, 2, 0)
619  - t4ddg_1(0, 2, 1, 1) * t4ddg_2(0, 1, 2, 1)
620  - t4ddg_1(0, 2, 1, 2) * t4ddg_2(0, 1, 2, 2)
621  - t4ddg_1(0, 2, 2, 0) * t4ddg_2(0, 2, 2, 0)
622  - t4ddg_1(0, 2, 2, 1) * t4ddg_2(0, 2, 2, 1)
623  - t4ddg_1(0, 2, 2, 2) * t4ddg_2(0, 2, 2, 2)
624 
625  - t4ddg_1(1, 0, 0, 0) * t4ddg_2(1, 0, 0, 0)
626  - t4ddg_1(1, 0, 0, 1) * t4ddg_2(1, 0, 0, 1)
627  - t4ddg_1(1, 0, 0, 2) * t4ddg_2(1, 0, 0, 2)
628  - t4ddg_1(1, 0, 1, 0) * t4ddg_2(1, 1, 0, 0)
629  - t4ddg_1(1, 0, 1, 1) * t4ddg_2(1, 1, 0, 1)
630  - t4ddg_1(1, 0, 1, 2) * t4ddg_2(1, 1, 0, 2)
631  - t4ddg_1(1, 0, 2, 0) * t4ddg_2(1, 2, 0, 0)
632  - t4ddg_1(1, 0, 2, 1) * t4ddg_2(1, 2, 0, 1)
633  - t4ddg_1(1, 0, 2, 2) * t4ddg_2(1, 2, 0, 2)
634  - t4ddg_1(1, 1, 0, 0) * t4ddg_2(1, 0, 1, 0)
635  - t4ddg_1(1, 1, 0, 1) * t4ddg_2(1, 0, 1, 1)
636  - t4ddg_1(1, 1, 0, 2) * t4ddg_2(1, 0, 1, 2)
637  - t4ddg_1(1, 1, 1, 0) * t4ddg_2(1, 1, 1, 0)
638  - t4ddg_1(1, 1, 1, 1) * t4ddg_2(1, 1, 1, 1)
639  - t4ddg_1(1, 1, 1, 2) * t4ddg_2(1, 1, 1, 2)
640  - t4ddg_1(1, 1, 2, 0) * t4ddg_2(1, 2, 1, 0)
641  - t4ddg_1(1, 1, 2, 1) * t4ddg_2(1, 2, 1, 1)
642  - t4ddg_1(1, 1, 2, 2) * t4ddg_2(1, 2, 1, 2)
643  - t4ddg_1(1, 2, 0, 0) * t4ddg_2(1, 0, 2, 0)
644  - t4ddg_1(1, 2, 0, 1) * t4ddg_2(1, 0, 2, 1)
645  - t4ddg_1(1, 2, 0, 2) * t4ddg_2(1, 0, 2, 2)
646  - t4ddg_1(1, 2, 1, 0) * t4ddg_2(1, 1, 2, 0)
647  - t4ddg_1(1, 2, 1, 1) * t4ddg_2(1, 1, 2, 1)
648  - t4ddg_1(1, 2, 1, 2) * t4ddg_2(1, 1, 2, 2)
649  - t4ddg_1(1, 2, 2, 0) * t4ddg_2(1, 2, 2, 0)
650  - t4ddg_1(1, 2, 2, 1) * t4ddg_2(1, 2, 2, 1)
651  - t4ddg_1(1, 2, 2, 2) * t4ddg_2(1, 2, 2, 2)
652 
653  - t4ddg_1(2, 0, 0, 0) * t4ddg_2(2, 0, 0, 0)
654  - t4ddg_1(2, 0, 0, 1) * t4ddg_2(2, 0, 0, 1)
655  - t4ddg_1(2, 0, 0, 2) * t4ddg_2(2, 0, 0, 2)
656  - t4ddg_1(2, 0, 1, 0) * t4ddg_2(2, 1, 0, 0)
657  - t4ddg_1(2, 0, 1, 1) * t4ddg_2(2, 1, 0, 1)
658  - t4ddg_1(2, 0, 1, 2) * t4ddg_2(2, 1, 0, 2)
659  - t4ddg_1(2, 0, 2, 0) * t4ddg_2(2, 2, 0, 0)
660  - t4ddg_1(2, 0, 2, 1) * t4ddg_2(2, 2, 0, 1)
661  - t4ddg_1(2, 0, 2, 2) * t4ddg_2(2, 2, 0, 2)
662  - t4ddg_1(2, 1, 0, 0) * t4ddg_2(2, 0, 1, 0)
663  - t4ddg_1(2, 1, 0, 1) * t4ddg_2(2, 0, 1, 1)
664  - t4ddg_1(2, 1, 0, 2) * t4ddg_2(2, 0, 1, 2)
665  - t4ddg_1(2, 1, 1, 0) * t4ddg_2(2, 1, 1, 0)
666  - t4ddg_1(2, 1, 1, 1) * t4ddg_2(2, 1, 1, 1)
667  - t4ddg_1(2, 1, 1, 2) * t4ddg_2(2, 1, 1, 2)
668  - t4ddg_1(2, 1, 2, 0) * t4ddg_2(2, 2, 1, 0)
669  - t4ddg_1(2, 1, 2, 1) * t4ddg_2(2, 2, 1, 1)
670  - t4ddg_1(2, 1, 2, 2) * t4ddg_2(2, 2, 1, 2)
671  - t4ddg_1(2, 2, 0, 0) * t4ddg_2(2, 0, 2, 0)
672  - t4ddg_1(2, 2, 0, 1) * t4ddg_2(2, 0, 2, 1)
673  - t4ddg_1(2, 2, 0, 2) * t4ddg_2(2, 0, 2, 2)
674  - t4ddg_1(2, 2, 1, 0) * t4ddg_2(2, 1, 2, 0)
675  - t4ddg_1(2, 2, 1, 1) * t4ddg_2(2, 1, 2, 1)
676  - t4ddg_1(2, 2, 1, 2) * t4ddg_2(2, 1, 2, 2)
677  - t4ddg_1(2, 2, 2, 0) * t4ddg_2(2, 2, 2, 0)
678  - t4ddg_1(2, 2, 2, 1) * t4ddg_2(2, 2, 2, 1)
679  - t4ddg_1(2, 2, 2, 2) * t4ddg_2(2, 2, 2, 2),
680  "T4ddg(i,j,k,l)*T4ddg(i,j,k,l)");
681 
682  t4ddg_1(i, j, k, l) = t4ddg_2(i, j, k, l) + t4ddg_3(i, j, k, l);
683  test_for_zero(t4ddg_1(0, 0, 0, 0)
684  - (t4ddg_2(0, 0, 0, 0) + t4ddg_3(0, 0, 0, 0)),
685  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(0,0,0,0)");
686  test_for_zero(t4ddg_1(0, 0, 0, 1)
687  - (t4ddg_2(0, 0, 0, 1) + t4ddg_3(0, 0, 0, 1)),
688  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(0,0,0,1)");
689  test_for_zero(t4ddg_1(0, 0, 0, 2)
690  - (t4ddg_2(0, 0, 0, 2) + t4ddg_3(0, 0, 0, 2)),
691  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(0,0,0,2)");
692  test_for_zero(t4ddg_1(0, 0, 1, 0)
693  - (t4ddg_2(0, 0, 1, 0) + t4ddg_3(0, 0, 1, 0)),
694  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(0,0,1,0)");
695  test_for_zero(t4ddg_1(0, 0, 1, 1)
696  - (t4ddg_2(0, 0, 1, 1) + t4ddg_3(0, 0, 1, 1)),
697  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(0,0,1,1)");
698  test_for_zero(t4ddg_1(0, 0, 1, 2)
699  - (t4ddg_2(0, 0, 1, 2) + t4ddg_3(0, 0, 1, 2)),
700  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(0,0,1,2)");
701  test_for_zero(t4ddg_1(0, 0, 2, 0)
702  - (t4ddg_2(0, 0, 2, 0) + t4ddg_3(0, 0, 2, 0)),
703  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(0,0,2,0)");
704  test_for_zero(t4ddg_1(0, 0, 2, 1)
705  - (t4ddg_2(0, 0, 2, 1) + t4ddg_3(0, 0, 2, 1)),
706  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(0,0,2,1)");
707  test_for_zero(t4ddg_1(0, 0, 2, 2)
708  - (t4ddg_2(0, 0, 2, 2) + t4ddg_3(0, 0, 2, 2)),
709  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(0,0,2,2)");
710  test_for_zero(t4ddg_1(0, 1, 0, 0)
711  - (t4ddg_2(0, 1, 0, 0) + t4ddg_3(0, 1, 0, 0)),
712  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(0,1,0,0)");
713  test_for_zero(t4ddg_1(0, 1, 0, 1)
714  - (t4ddg_2(0, 1, 0, 1) + t4ddg_3(0, 1, 0, 1)),
715  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(0,1,0,1)");
716  test_for_zero(t4ddg_1(0, 1, 0, 2)
717  - (t4ddg_2(0, 1, 0, 2) + t4ddg_3(0, 1, 0, 2)),
718  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(0,1,0,2)");
719  test_for_zero(t4ddg_1(0, 1, 1, 0)
720  - (t4ddg_2(0, 1, 1, 0) + t4ddg_3(0, 1, 1, 0)),
721  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(0,1,1,0)");
722  test_for_zero(t4ddg_1(0, 1, 1, 1)
723  - (t4ddg_2(0, 1, 1, 1) + t4ddg_3(0, 1, 1, 1)),
724  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(0,1,1,1)");
725  test_for_zero(t4ddg_1(0, 1, 1, 2)
726  - (t4ddg_2(0, 1, 1, 2) + t4ddg_3(0, 1, 1, 2)),
727  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(0,1,1,2)");
728  test_for_zero(t4ddg_1(0, 1, 2, 0)
729  - (t4ddg_2(0, 1, 2, 0) + t4ddg_3(0, 1, 2, 0)),
730  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(0,1,2,0)");
731  test_for_zero(t4ddg_1(0, 1, 2, 1)
732  - (t4ddg_2(0, 1, 2, 1) + t4ddg_3(0, 1, 2, 1)),
733  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(0,1,2,1)");
734  test_for_zero(t4ddg_1(0, 1, 2, 2)
735  - (t4ddg_2(0, 1, 2, 2) + t4ddg_3(0, 1, 2, 2)),
736  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(0,1,2,2)");
737  test_for_zero(t4ddg_1(0, 2, 0, 0)
738  - (t4ddg_2(0, 2, 0, 0) + t4ddg_3(0, 2, 0, 0)),
739  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(0,2,0,0)");
740  test_for_zero(t4ddg_1(0, 2, 0, 1)
741  - (t4ddg_2(0, 2, 0, 1) + t4ddg_3(0, 2, 0, 1)),
742  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(0,2,0,1)");
743  test_for_zero(t4ddg_1(0, 2, 0, 2)
744  - (t4ddg_2(0, 2, 0, 2) + t4ddg_3(0, 2, 0, 2)),
745  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(0,2,0,2)");
746  test_for_zero(t4ddg_1(0, 2, 1, 0)
747  - (t4ddg_2(0, 2, 1, 0) + t4ddg_3(0, 2, 1, 0)),
748  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(0,2,1,0)");
749  test_for_zero(t4ddg_1(0, 2, 1, 1)
750  - (t4ddg_2(0, 2, 1, 1) + t4ddg_3(0, 2, 1, 1)),
751  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(0,2,1,1)");
752  test_for_zero(t4ddg_1(0, 2, 1, 2)
753  - (t4ddg_2(0, 2, 1, 2) + t4ddg_3(0, 2, 1, 2)),
754  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(0,2,1,2)");
755  test_for_zero(t4ddg_1(0, 2, 2, 0)
756  - (t4ddg_2(0, 2, 2, 0) + t4ddg_3(0, 2, 2, 0)),
757  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(0,2,2,0)");
758  test_for_zero(t4ddg_1(0, 2, 2, 1)
759  - (t4ddg_2(0, 2, 2, 1) + t4ddg_3(0, 2, 2, 1)),
760  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(0,2,2,1)");
761  test_for_zero(t4ddg_1(0, 2, 2, 2)
762  - (t4ddg_2(0, 2, 2, 2) + t4ddg_3(0, 2, 2, 2)),
763  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(0,2,2,2)");
764  test_for_zero(t4ddg_1(1, 0, 0, 0)
765  - (t4ddg_2(1, 0, 0, 0) + t4ddg_3(1, 0, 0, 0)),
766  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(1,0,0,0)");
767  test_for_zero(t4ddg_1(1, 0, 0, 1)
768  - (t4ddg_2(1, 0, 0, 1) + t4ddg_3(1, 0, 0, 1)),
769  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(1,0,0,1)");
770  test_for_zero(t4ddg_1(1, 0, 0, 2)
771  - (t4ddg_2(1, 0, 0, 2) + t4ddg_3(1, 0, 0, 2)),
772  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(1,0,0,2)");
773  test_for_zero(t4ddg_1(1, 0, 1, 0)
774  - (t4ddg_2(1, 0, 1, 0) + t4ddg_3(1, 0, 1, 0)),
775  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(1,0,1,0)");
776  test_for_zero(t4ddg_1(1, 0, 1, 1)
777  - (t4ddg_2(1, 0, 1, 1) + t4ddg_3(1, 0, 1, 1)),
778  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(1,0,1,1)");
779  test_for_zero(t4ddg_1(1, 0, 1, 2)
780  - (t4ddg_2(1, 0, 1, 2) + t4ddg_3(1, 0, 1, 2)),
781  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(1,0,1,2)");
782  test_for_zero(t4ddg_1(1, 0, 2, 0)
783  - (t4ddg_2(1, 0, 2, 0) + t4ddg_3(1, 0, 2, 0)),
784  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(1,0,2,0)");
785  test_for_zero(t4ddg_1(1, 0, 2, 1)
786  - (t4ddg_2(1, 0, 2, 1) + t4ddg_3(1, 0, 2, 1)),
787  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(1,0,2,1)");
788  test_for_zero(t4ddg_1(1, 0, 2, 2)
789  - (t4ddg_2(1, 0, 2, 2) + t4ddg_3(1, 0, 2, 2)),
790  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(1,0,2,2)");
791  test_for_zero(t4ddg_1(1, 1, 0, 0)
792  - (t4ddg_2(1, 1, 0, 0) + t4ddg_3(1, 1, 0, 0)),
793  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(1,1,0,0)");
794  test_for_zero(t4ddg_1(1, 1, 0, 1)
795  - (t4ddg_2(1, 1, 0, 1) + t4ddg_3(1, 1, 0, 1)),
796  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(1,1,0,1)");
797  test_for_zero(t4ddg_1(1, 1, 0, 2)
798  - (t4ddg_2(1, 1, 0, 2) + t4ddg_3(1, 1, 0, 2)),
799  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(1,1,0,2)");
800  test_for_zero(t4ddg_1(1, 1, 1, 0)
801  - (t4ddg_2(1, 1, 1, 0) + t4ddg_3(1, 1, 1, 0)),
802  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(1,1,1,0)");
803  test_for_zero(t4ddg_1(1, 1, 1, 1)
804  - (t4ddg_2(1, 1, 1, 1) + t4ddg_3(1, 1, 1, 1)),
805  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(1,1,1,1)");
806  test_for_zero(t4ddg_1(1, 1, 1, 2)
807  - (t4ddg_2(1, 1, 1, 2) + t4ddg_3(1, 1, 1, 2)),
808  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(1,1,1,2)");
809  test_for_zero(t4ddg_1(1, 1, 2, 0)
810  - (t4ddg_2(1, 1, 2, 0) + t4ddg_3(1, 1, 2, 0)),
811  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(1,1,2,0)");
812  test_for_zero(t4ddg_1(1, 1, 2, 1)
813  - (t4ddg_2(1, 1, 2, 1) + t4ddg_3(1, 1, 2, 1)),
814  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(1,1,2,1)");
815  test_for_zero(t4ddg_1(1, 1, 2, 2)
816  - (t4ddg_2(1, 1, 2, 2) + t4ddg_3(1, 1, 2, 2)),
817  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(1,1,2,2)");
818  test_for_zero(t4ddg_1(1, 2, 0, 0)
819  - (t4ddg_2(1, 2, 0, 0) + t4ddg_3(1, 2, 0, 0)),
820  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(1,2,0,0)");
821  test_for_zero(t4ddg_1(1, 2, 0, 1)
822  - (t4ddg_2(1, 2, 0, 1) + t4ddg_3(1, 2, 0, 1)),
823  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(1,2,0,1)");
824  test_for_zero(t4ddg_1(1, 2, 0, 2)
825  - (t4ddg_2(1, 2, 0, 2) + t4ddg_3(1, 2, 0, 2)),
826  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(1,2,0,2)");
827  test_for_zero(t4ddg_1(1, 2, 1, 0)
828  - (t4ddg_2(1, 2, 1, 0) + t4ddg_3(1, 2, 1, 0)),
829  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(1,2,1,0)");
830  test_for_zero(t4ddg_1(1, 2, 1, 1)
831  - (t4ddg_2(1, 2, 1, 1) + t4ddg_3(1, 2, 1, 1)),
832  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(1,2,1,1)");
833  test_for_zero(t4ddg_1(1, 2, 1, 2)
834  - (t4ddg_2(1, 2, 1, 2) + t4ddg_3(1, 2, 1, 2)),
835  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(1,2,1,2)");
836  test_for_zero(t4ddg_1(1, 2, 2, 0)
837  - (t4ddg_2(1, 2, 2, 0) + t4ddg_3(1, 2, 2, 0)),
838  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(1,2,2,0)");
839  test_for_zero(t4ddg_1(1, 2, 2, 1)
840  - (t4ddg_2(1, 2, 2, 1) + t4ddg_3(1, 2, 2, 1)),
841  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(1,2,2,1)");
842  test_for_zero(t4ddg_1(1, 2, 2, 2)
843  - (t4ddg_2(1, 2, 2, 2) + t4ddg_3(1, 2, 2, 2)),
844  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(1,2,2,2)");
845  test_for_zero(t4ddg_1(2, 0, 0, 0)
846  - (t4ddg_2(2, 0, 0, 0) + t4ddg_3(2, 0, 0, 0)),
847  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(2,0,0,0)");
848  test_for_zero(t4ddg_1(2, 0, 0, 1)
849  - (t4ddg_2(2, 0, 0, 1) + t4ddg_3(2, 0, 0, 1)),
850  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(2,0,0,1)");
851  test_for_zero(t4ddg_1(2, 0, 0, 2)
852  - (t4ddg_2(2, 0, 0, 2) + t4ddg_3(2, 0, 0, 2)),
853  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(2,0,0,2)");
854  test_for_zero(t4ddg_1(2, 0, 1, 0)
855  - (t4ddg_2(2, 0, 1, 0) + t4ddg_3(2, 0, 1, 0)),
856  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(2,0,1,0)");
857  test_for_zero(t4ddg_1(2, 0, 1, 1)
858  - (t4ddg_2(2, 0, 1, 1) + t4ddg_3(2, 0, 1, 1)),
859  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(2,0,1,1)");
860  test_for_zero(t4ddg_1(2, 0, 1, 2)
861  - (t4ddg_2(2, 0, 1, 2) + t4ddg_3(2, 0, 1, 2)),
862  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(2,0,1,2)");
863  test_for_zero(t4ddg_1(2, 0, 2, 0)
864  - (t4ddg_2(2, 0, 2, 0) + t4ddg_3(2, 0, 2, 0)),
865  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(2,0,2,0)");
866  test_for_zero(t4ddg_1(2, 0, 2, 1)
867  - (t4ddg_2(2, 0, 2, 1) + t4ddg_3(2, 0, 2, 1)),
868  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(2,0,2,1)");
869  test_for_zero(t4ddg_1(2, 0, 2, 2)
870  - (t4ddg_2(2, 0, 2, 2) + t4ddg_3(2, 0, 2, 2)),
871  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(2,0,2,2)");
872  test_for_zero(t4ddg_1(2, 1, 0, 0)
873  - (t4ddg_2(2, 1, 0, 0) + t4ddg_3(2, 1, 0, 0)),
874  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(2,1,0,0)");
875  test_for_zero(t4ddg_1(2, 1, 0, 1)
876  - (t4ddg_2(2, 1, 0, 1) + t4ddg_3(2, 1, 0, 1)),
877  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(2,1,0,1)");
878  test_for_zero(t4ddg_1(2, 1, 0, 2)
879  - (t4ddg_2(2, 1, 0, 2) + t4ddg_3(2, 1, 0, 2)),
880  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(2,1,0,2)");
881  test_for_zero(t4ddg_1(2, 1, 1, 0)
882  - (t4ddg_2(2, 1, 1, 0) + t4ddg_3(2, 1, 1, 0)),
883  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(2,1,1,0)");
884  test_for_zero(t4ddg_1(2, 1, 1, 1)
885  - (t4ddg_2(2, 1, 1, 1) + t4ddg_3(2, 1, 1, 1)),
886  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(2,1,1,1)");
887  test_for_zero(t4ddg_1(2, 1, 1, 2)
888  - (t4ddg_2(2, 1, 1, 2) + t4ddg_3(2, 1, 1, 2)),
889  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(2,1,1,2)");
890  test_for_zero(t4ddg_1(2, 1, 2, 0)
891  - (t4ddg_2(2, 1, 2, 0) + t4ddg_3(2, 1, 2, 0)),
892  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(2,1,2,0)");
893  test_for_zero(t4ddg_1(2, 1, 2, 1)
894  - (t4ddg_2(2, 1, 2, 1) + t4ddg_3(2, 1, 2, 1)),
895  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(2,1,2,1)");
896  test_for_zero(t4ddg_1(2, 1, 2, 2)
897  - (t4ddg_2(2, 1, 2, 2) + t4ddg_3(2, 1, 2, 2)),
898  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(2,1,2,2)");
899  test_for_zero(t4ddg_1(2, 2, 0, 0)
900  - (t4ddg_2(2, 2, 0, 0) + t4ddg_3(2, 2, 0, 0)),
901  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(2,2,0,0)");
902  test_for_zero(t4ddg_1(2, 2, 0, 1)
903  - (t4ddg_2(2, 2, 0, 1) + t4ddg_3(2, 2, 0, 1)),
904  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(2,2,0,1)");
905  test_for_zero(t4ddg_1(2, 2, 0, 2)
906  - (t4ddg_2(2, 2, 0, 2) + t4ddg_3(2, 2, 0, 2)),
907  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(2,2,0,2)");
908  test_for_zero(t4ddg_1(2, 2, 1, 0)
909  - (t4ddg_2(2, 2, 1, 0) + t4ddg_3(2, 2, 1, 0)),
910  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(2,2,1,0)");
911  test_for_zero(t4ddg_1(2, 2, 1, 1)
912  - (t4ddg_2(2, 2, 1, 1) + t4ddg_3(2, 2, 1, 1)),
913  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(2,2,1,1)");
914  test_for_zero(t4ddg_1(2, 2, 1, 2)
915  - (t4ddg_2(2, 2, 1, 2) + t4ddg_3(2, 2, 1, 2)),
916  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(2,2,1,2)");
917  test_for_zero(t4ddg_1(2, 2, 2, 0)
918  - (t4ddg_2(2, 2, 2, 0) + t4ddg_3(2, 2, 2, 0)),
919  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(2,2,2,0)");
920  test_for_zero(t4ddg_1(2, 2, 2, 1)
921  - (t4ddg_2(2, 2, 2, 1) + t4ddg_3(2, 2, 2, 1)),
922  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(2,2,2,1)");
923  test_for_zero(t4ddg_1(2, 2, 2, 2)
924  - (t4ddg_2(2, 2, 2, 2) + t4ddg_3(2, 2, 2, 2)),
925  "T4ddg(i,j,k,l)+T4ddg(i,j,k,l)(2,2,2,2)");
926 
927  t4ddg_1(i, j, k, l) = t4ddg_2(i, j, k, l) + t4ddg_3(k, l, i, j);
928  test_for_zero(t4ddg_1(0, 0, 0, 0)
929  - (t4ddg_2(0, 0, 0, 0) + t4ddg_3(0, 0, 0, 0)),
930  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(0,0,0,0)");
931  test_for_zero(t4ddg_1(0, 0, 0, 1)
932  - (t4ddg_2(0, 0, 0, 1) + t4ddg_3(0, 1, 0, 0)),
933  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(0,0,0,1)");
934  test_for_zero(t4ddg_1(0, 0, 0, 2)
935  - (t4ddg_2(0, 0, 0, 2) + t4ddg_3(0, 2, 0, 0)),
936  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(0,0,0,2)");
937  test_for_zero(t4ddg_1(0, 0, 1, 0)
938  - (t4ddg_2(0, 0, 1, 0) + t4ddg_3(1, 0, 0, 0)),
939  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(0,0,1,0)");
940  test_for_zero(t4ddg_1(0, 0, 1, 1)
941  - (t4ddg_2(0, 0, 1, 1) + t4ddg_3(1, 1, 0, 0)),
942  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(0,0,1,1)");
943  test_for_zero(t4ddg_1(0, 0, 1, 2)
944  - (t4ddg_2(0, 0, 1, 2) + t4ddg_3(1, 2, 0, 0)),
945  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(0,0,1,2)");
946  test_for_zero(t4ddg_1(0, 0, 2, 0)
947  - (t4ddg_2(0, 0, 2, 0) + t4ddg_3(2, 0, 0, 0)),
948  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(0,0,2,0)");
949  test_for_zero(t4ddg_1(0, 0, 2, 1)
950  - (t4ddg_2(0, 0, 2, 1) + t4ddg_3(2, 1, 0, 0)),
951  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(0,0,2,1)");
952  test_for_zero(t4ddg_1(0, 0, 2, 2)
953  - (t4ddg_2(0, 0, 2, 2) + t4ddg_3(2, 2, 0, 0)),
954  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(0,0,2,2)");
955  test_for_zero(t4ddg_1(0, 1, 0, 0)
956  - (t4ddg_2(0, 1, 0, 0) + t4ddg_3(0, 0, 0, 1)),
957  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(0,1,0,0)");
958  test_for_zero(t4ddg_1(0, 1, 0, 1)
959  - (t4ddg_2(0, 1, 0, 1) + t4ddg_3(0, 1, 0, 1)),
960  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(0,1,0,1)");
961  test_for_zero(t4ddg_1(0, 1, 0, 2)
962  - (t4ddg_2(0, 1, 0, 2) + t4ddg_3(0, 2, 0, 1)),
963  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(0,1,0,2)");
964  test_for_zero(t4ddg_1(0, 1, 1, 0)
965  - (t4ddg_2(0, 1, 1, 0) + t4ddg_3(1, 0, 0, 1)),
966  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(0,1,1,0)");
967  test_for_zero(t4ddg_1(0, 1, 1, 1)
968  - (t4ddg_2(0, 1, 1, 1) + t4ddg_3(1, 1, 0, 1)),
969  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(0,1,1,1)");
970  test_for_zero(t4ddg_1(0, 1, 1, 2)
971  - (t4ddg_2(0, 1, 1, 2) + t4ddg_3(1, 2, 0, 1)),
972  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(0,1,1,2)");
973  test_for_zero(t4ddg_1(0, 1, 2, 0)
974  - (t4ddg_2(0, 1, 2, 0) + t4ddg_3(2, 0, 0, 1)),
975  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(0,1,2,0)");
976  test_for_zero(t4ddg_1(0, 1, 2, 1)
977  - (t4ddg_2(0, 1, 2, 1) + t4ddg_3(2, 1, 0, 1)),
978  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(0,1,2,1)");
979  test_for_zero(t4ddg_1(0, 1, 2, 2)
980  - (t4ddg_2(0, 1, 2, 2) + t4ddg_3(2, 2, 0, 1)),
981  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(0,1,2,2)");
982  test_for_zero(t4ddg_1(0, 2, 0, 0)
983  - (t4ddg_2(0, 2, 0, 0) + t4ddg_3(0, 0, 0, 2)),
984  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(0,2,0,0)");
985  test_for_zero(t4ddg_1(0, 2, 0, 1)
986  - (t4ddg_2(0, 2, 0, 1) + t4ddg_3(0, 1, 0, 2)),
987  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(0,2,0,1)");
988  test_for_zero(t4ddg_1(0, 2, 0, 2)
989  - (t4ddg_2(0, 2, 0, 2) + t4ddg_3(0, 2, 0, 2)),
990  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(0,2,0,2)");
991  test_for_zero(t4ddg_1(0, 2, 1, 0)
992  - (t4ddg_2(0, 2, 1, 0) + t4ddg_3(1, 0, 0, 2)),
993  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(0,2,1,0)");
994  test_for_zero(t4ddg_1(0, 2, 1, 1)
995  - (t4ddg_2(0, 2, 1, 1) + t4ddg_3(1, 1, 0, 2)),
996  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(0,2,1,1)");
997  test_for_zero(t4ddg_1(0, 2, 1, 2)
998  - (t4ddg_2(0, 2, 1, 2) + t4ddg_3(1, 2, 0, 2)),
999  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(0,2,1,2)");
1000  test_for_zero(t4ddg_1(0, 2, 2, 0)
1001  - (t4ddg_2(0, 2, 2, 0) + t4ddg_3(2, 0, 0, 2)),
1002  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(0,2,2,0)");
1003  test_for_zero(t4ddg_1(0, 2, 2, 1)
1004  - (t4ddg_2(0, 2, 2, 1) + t4ddg_3(2, 1, 0, 2)),
1005  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(0,2,2,1)");
1006  test_for_zero(t4ddg_1(0, 2, 2, 2)
1007  - (t4ddg_2(0, 2, 2, 2) + t4ddg_3(2, 2, 0, 2)),
1008  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(0,2,2,2)");
1009  test_for_zero(t4ddg_1(1, 0, 0, 0)
1010  - (t4ddg_2(1, 0, 0, 0) + t4ddg_3(0, 0, 1, 0)),
1011  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(1,0,0,0)");
1012  test_for_zero(t4ddg_1(1, 0, 0, 1)
1013  - (t4ddg_2(1, 0, 0, 1) + t4ddg_3(0, 1, 1, 0)),
1014  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(1,0,0,1)");
1015  test_for_zero(t4ddg_1(1, 0, 0, 2)
1016  - (t4ddg_2(1, 0, 0, 2) + t4ddg_3(0, 2, 1, 0)),
1017  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(1,0,0,2)");
1018  test_for_zero(t4ddg_1(1, 0, 1, 0)
1019  - (t4ddg_2(1, 0, 1, 0) + t4ddg_3(1, 0, 1, 0)),
1020  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(1,0,1,0)");
1021  test_for_zero(t4ddg_1(1, 0, 1, 1)
1022  - (t4ddg_2(1, 0, 1, 1) + t4ddg_3(1, 1, 1, 0)),
1023  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(1,0,1,1)");
1024  test_for_zero(t4ddg_1(1, 0, 1, 2)
1025  - (t4ddg_2(1, 0, 1, 2) + t4ddg_3(1, 2, 1, 0)),
1026  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(1,0,1,2)");
1027  test_for_zero(t4ddg_1(1, 0, 2, 0)
1028  - (t4ddg_2(1, 0, 2, 0) + t4ddg_3(2, 0, 1, 0)),
1029  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(1,0,2,0)");
1030  test_for_zero(t4ddg_1(1, 0, 2, 1)
1031  - (t4ddg_2(1, 0, 2, 1) + t4ddg_3(2, 1, 1, 0)),
1032  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(1,0,2,1)");
1033  test_for_zero(t4ddg_1(1, 0, 2, 2)
1034  - (t4ddg_2(1, 0, 2, 2) + t4ddg_3(2, 2, 1, 0)),
1035  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(1,0,2,2)");
1036  test_for_zero(t4ddg_1(1, 1, 0, 0)
1037  - (t4ddg_2(1, 1, 0, 0) + t4ddg_3(0, 0, 1, 1)),
1038  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(1,1,0,0)");
1039  test_for_zero(t4ddg_1(1, 1, 0, 1)
1040  - (t4ddg_2(1, 1, 0, 1) + t4ddg_3(0, 1, 1, 1)),
1041  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(1,1,0,1)");
1042  test_for_zero(t4ddg_1(1, 1, 0, 2)
1043  - (t4ddg_2(1, 1, 0, 2) + t4ddg_3(0, 2, 1, 1)),
1044  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(1,1,0,2)");
1045  test_for_zero(t4ddg_1(1, 1, 1, 0)
1046  - (t4ddg_2(1, 1, 1, 0) + t4ddg_3(1, 0, 1, 1)),
1047  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(1,1,1,0)");
1048  test_for_zero(t4ddg_1(1, 1, 1, 1)
1049  - (t4ddg_2(1, 1, 1, 1) + t4ddg_3(1, 1, 1, 1)),
1050  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(1,1,1,1)");
1051  test_for_zero(t4ddg_1(1, 1, 1, 2)
1052  - (t4ddg_2(1, 1, 1, 2) + t4ddg_3(1, 2, 1, 1)),
1053  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(1,1,1,2)");
1054  test_for_zero(t4ddg_1(1, 1, 2, 0)
1055  - (t4ddg_2(1, 1, 2, 0) + t4ddg_3(2, 0, 1, 1)),
1056  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(1,1,2,0)");
1057  test_for_zero(t4ddg_1(1, 1, 2, 1)
1058  - (t4ddg_2(1, 1, 2, 1) + t4ddg_3(2, 1, 1, 1)),
1059  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(1,1,2,1)");
1060  test_for_zero(t4ddg_1(1, 1, 2, 2)
1061  - (t4ddg_2(1, 1, 2, 2) + t4ddg_3(2, 2, 1, 1)),
1062  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(1,1,2,2)");
1063  test_for_zero(t4ddg_1(1, 2, 0, 0)
1064  - (t4ddg_2(1, 2, 0, 0) + t4ddg_3(0, 0, 1, 2)),
1065  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(1,2,0,0)");
1066  test_for_zero(t4ddg_1(1, 2, 0, 1)
1067  - (t4ddg_2(1, 2, 0, 1) + t4ddg_3(0, 1, 1, 2)),
1068  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(1,2,0,1)");
1069  test_for_zero(t4ddg_1(1, 2, 0, 2)
1070  - (t4ddg_2(1, 2, 0, 2) + t4ddg_3(0, 2, 1, 2)),
1071  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(1,2,0,2)");
1072  test_for_zero(t4ddg_1(1, 2, 1, 0)
1073  - (t4ddg_2(1, 2, 1, 0) + t4ddg_3(1, 0, 1, 2)),
1074  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(1,2,1,0)");
1075  test_for_zero(t4ddg_1(1, 2, 1, 1)
1076  - (t4ddg_2(1, 2, 1, 1) + t4ddg_3(1, 1, 1, 2)),
1077  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(1,2,1,1)");
1078  test_for_zero(t4ddg_1(1, 2, 1, 2)
1079  - (t4ddg_2(1, 2, 1, 2) + t4ddg_3(1, 2, 1, 2)),
1080  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(1,2,1,2)");
1081  test_for_zero(t4ddg_1(1, 2, 2, 0)
1082  - (t4ddg_2(1, 2, 2, 0) + t4ddg_3(2, 0, 1, 2)),
1083  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(1,2,2,0)");
1084  test_for_zero(t4ddg_1(1, 2, 2, 1)
1085  - (t4ddg_2(1, 2, 2, 1) + t4ddg_3(2, 1, 1, 2)),
1086  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(1,2,2,1)");
1087  test_for_zero(t4ddg_1(1, 2, 2, 2)
1088  - (t4ddg_2(1, 2, 2, 2) + t4ddg_3(2, 2, 1, 2)),
1089  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(1,2,2,2)");
1090  test_for_zero(t4ddg_1(2, 0, 0, 0)
1091  - (t4ddg_2(2, 0, 0, 0) + t4ddg_3(0, 0, 2, 0)),
1092  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(2,0,0,0)");
1093  test_for_zero(t4ddg_1(2, 0, 0, 1)
1094  - (t4ddg_2(2, 0, 0, 1) + t4ddg_3(0, 1, 2, 0)),
1095  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(2,0,0,1)");
1096  test_for_zero(t4ddg_1(2, 0, 0, 2)
1097  - (t4ddg_2(2, 0, 0, 2) + t4ddg_3(0, 2, 2, 0)),
1098  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(2,0,0,2)");
1099  test_for_zero(t4ddg_1(2, 0, 1, 0)
1100  - (t4ddg_2(2, 0, 1, 0) + t4ddg_3(1, 0, 2, 0)),
1101  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(2,0,1,0)");
1102  test_for_zero(t4ddg_1(2, 0, 1, 1)
1103  - (t4ddg_2(2, 0, 1, 1) + t4ddg_3(1, 1, 2, 0)),
1104  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(2,0,1,1)");
1105  test_for_zero(t4ddg_1(2, 0, 1, 2)
1106  - (t4ddg_2(2, 0, 1, 2) + t4ddg_3(1, 2, 2, 0)),
1107  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(2,0,1,2)");
1108  test_for_zero(t4ddg_1(2, 0, 2, 0)
1109  - (t4ddg_2(2, 0, 2, 0) + t4ddg_3(2, 0, 2, 0)),
1110  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(2,0,2,0)");
1111  test_for_zero(t4ddg_1(2, 0, 2, 1)
1112  - (t4ddg_2(2, 0, 2, 1) + t4ddg_3(2, 1, 2, 0)),
1113  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(2,0,2,1)");
1114  test_for_zero(t4ddg_1(2, 0, 2, 2)
1115  - (t4ddg_2(2, 0, 2, 2) + t4ddg_3(2, 2, 2, 0)),
1116  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(2,0,2,2)");
1117  test_for_zero(t4ddg_1(2, 1, 0, 0)
1118  - (t4ddg_2(2, 1, 0, 0) + t4ddg_3(0, 0, 2, 1)),
1119  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(2,1,0,0)");
1120  test_for_zero(t4ddg_1(2, 1, 0, 1)
1121  - (t4ddg_2(2, 1, 0, 1) + t4ddg_3(0, 1, 2, 1)),
1122  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(2,1,0,1)");
1123  test_for_zero(t4ddg_1(2, 1, 0, 2)
1124  - (t4ddg_2(2, 1, 0, 2) + t4ddg_3(0, 2, 2, 1)),
1125  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(2,1,0,2)");
1126  test_for_zero(t4ddg_1(2, 1, 1, 0)
1127  - (t4ddg_2(2, 1, 1, 0) + t4ddg_3(1, 0, 2, 1)),
1128  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(2,1,1,0)");
1129  test_for_zero(t4ddg_1(2, 1, 1, 1)
1130  - (t4ddg_2(2, 1, 1, 1) + t4ddg_3(1, 1, 2, 1)),
1131  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(2,1,1,1)");
1132  test_for_zero(t4ddg_1(2, 1, 1, 2)
1133  - (t4ddg_2(2, 1, 1, 2) + t4ddg_3(1, 2, 2, 1)),
1134  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(2,1,1,2)");
1135  test_for_zero(t4ddg_1(2, 1, 2, 0)
1136  - (t4ddg_2(2, 1, 2, 0) + t4ddg_3(2, 0, 2, 1)),
1137  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(2,1,2,0)");
1138  test_for_zero(t4ddg_1(2, 1, 2, 1)
1139  - (t4ddg_2(2, 1, 2, 1) + t4ddg_3(2, 1, 2, 1)),
1140  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(2,1,2,1)");
1141  test_for_zero(t4ddg_1(2, 1, 2, 2)
1142  - (t4ddg_2(2, 1, 2, 2) + t4ddg_3(2, 2, 2, 1)),
1143  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(2,1,2,2)");
1144  test_for_zero(t4ddg_1(2, 2, 0, 0)
1145  - (t4ddg_2(2, 2, 0, 0) + t4ddg_3(0, 0, 2, 2)),
1146  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(2,2,0,0)");
1147  test_for_zero(t4ddg_1(2, 2, 0, 1)
1148  - (t4ddg_2(2, 2, 0, 1) + t4ddg_3(0, 1, 2, 2)),
1149  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(2,2,0,1)");
1150  test_for_zero(t4ddg_1(2, 2, 0, 2)
1151  - (t4ddg_2(2, 2, 0, 2) + t4ddg_3(0, 2, 2, 2)),
1152  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(2,2,0,2)");
1153  test_for_zero(t4ddg_1(2, 2, 1, 0)
1154  - (t4ddg_2(2, 2, 1, 0) + t4ddg_3(1, 0, 2, 2)),
1155  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(2,2,1,0)");
1156  test_for_zero(t4ddg_1(2, 2, 1, 1)
1157  - (t4ddg_2(2, 2, 1, 1) + t4ddg_3(1, 1, 2, 2)),
1158  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(2,2,1,1)");
1159  test_for_zero(t4ddg_1(2, 2, 1, 2)
1160  - (t4ddg_2(2, 2, 1, 2) + t4ddg_3(1, 2, 2, 2)),
1161  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(2,2,1,2)");
1162  test_for_zero(t4ddg_1(2, 2, 2, 0)
1163  - (t4ddg_2(2, 2, 2, 0) + t4ddg_3(2, 0, 2, 2)),
1164  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(2,2,2,0)");
1165  test_for_zero(t4ddg_1(2, 2, 2, 1)
1166  - (t4ddg_2(2, 2, 2, 1) + t4ddg_3(2, 1, 2, 2)),
1167  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(2,2,2,1)");
1168  test_for_zero(t4ddg_1(2, 2, 2, 2)
1169  - (t4ddg_2(2, 2, 2, 2) + t4ddg_3(2, 2, 2, 2)),
1170  "T4ddg(i,j,k,l)+T4ddg(k,l,i,j)(2,2,2,2)");
1171 
1172  t4ddg_1(i, j, k, l) = t4ddg_2(i, j, k, l) - t4ddg_3(i, j, k, l);
1173  test_for_zero(t4ddg_1(0, 0, 0, 0)
1174  - (t4ddg_2(0, 0, 0, 0) - t4ddg_3(0, 0, 0, 0)),
1175  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(0,0,0,0)");
1176  test_for_zero(t4ddg_1(0, 0, 0, 1)
1177  - (t4ddg_2(0, 0, 0, 1) - t4ddg_3(0, 0, 0, 1)),
1178  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(0,0,0,1)");
1179  test_for_zero(t4ddg_1(0, 0, 0, 2)
1180  - (t4ddg_2(0, 0, 0, 2) - t4ddg_3(0, 0, 0, 2)),
1181  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(0,0,0,2)");
1182  test_for_zero(t4ddg_1(0, 0, 1, 0)
1183  - (t4ddg_2(0, 0, 1, 0) - t4ddg_3(0, 0, 1, 0)),
1184  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(0,0,1,0)");
1185  test_for_zero(t4ddg_1(0, 0, 1, 1)
1186  - (t4ddg_2(0, 0, 1, 1) - t4ddg_3(0, 0, 1, 1)),
1187  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(0,0,1,1)");
1188  test_for_zero(t4ddg_1(0, 0, 1, 2)
1189  - (t4ddg_2(0, 0, 1, 2) - t4ddg_3(0, 0, 1, 2)),
1190  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(0,0,1,2)");
1191  test_for_zero(t4ddg_1(0, 0, 2, 0)
1192  - (t4ddg_2(0, 0, 2, 0) - t4ddg_3(0, 0, 2, 0)),
1193  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(0,0,2,0)");
1194  test_for_zero(t4ddg_1(0, 0, 2, 1)
1195  - (t4ddg_2(0, 0, 2, 1) - t4ddg_3(0, 0, 2, 1)),
1196  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(0,0,2,1)");
1197  test_for_zero(t4ddg_1(0, 0, 2, 2)
1198  - (t4ddg_2(0, 0, 2, 2) - t4ddg_3(0, 0, 2, 2)),
1199  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(0,0,2,2)");
1200  test_for_zero(t4ddg_1(0, 1, 0, 0)
1201  - (t4ddg_2(0, 1, 0, 0) - t4ddg_3(0, 1, 0, 0)),
1202  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(0,1,0,0)");
1203  test_for_zero(t4ddg_1(0, 1, 0, 1)
1204  - (t4ddg_2(0, 1, 0, 1) - t4ddg_3(0, 1, 0, 1)),
1205  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(0,1,0,1)");
1206  test_for_zero(t4ddg_1(0, 1, 0, 2)
1207  - (t4ddg_2(0, 1, 0, 2) - t4ddg_3(0, 1, 0, 2)),
1208  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(0,1,0,2)");
1209  test_for_zero(t4ddg_1(0, 1, 1, 0)
1210  - (t4ddg_2(0, 1, 1, 0) - t4ddg_3(0, 1, 1, 0)),
1211  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(0,1,1,0)");
1212  test_for_zero(t4ddg_1(0, 1, 1, 1)
1213  - (t4ddg_2(0, 1, 1, 1) - t4ddg_3(0, 1, 1, 1)),
1214  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(0,1,1,1)");
1215  test_for_zero(t4ddg_1(0, 1, 1, 2)
1216  - (t4ddg_2(0, 1, 1, 2) - t4ddg_3(0, 1, 1, 2)),
1217  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(0,1,1,2)");
1218  test_for_zero(t4ddg_1(0, 1, 2, 0)
1219  - (t4ddg_2(0, 1, 2, 0) - t4ddg_3(0, 1, 2, 0)),
1220  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(0,1,2,0)");
1221  test_for_zero(t4ddg_1(0, 1, 2, 1)
1222  - (t4ddg_2(0, 1, 2, 1) - t4ddg_3(0, 1, 2, 1)),
1223  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(0,1,2,1)");
1224  test_for_zero(t4ddg_1(0, 1, 2, 2)
1225  - (t4ddg_2(0, 1, 2, 2) - t4ddg_3(0, 1, 2, 2)),
1226  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(0,1,2,2)");
1227  test_for_zero(t4ddg_1(0, 2, 0, 0)
1228  - (t4ddg_2(0, 2, 0, 0) - t4ddg_3(0, 2, 0, 0)),
1229  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(0,2,0,0)");
1230  test_for_zero(t4ddg_1(0, 2, 0, 1)
1231  - (t4ddg_2(0, 2, 0, 1) - t4ddg_3(0, 2, 0, 1)),
1232  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(0,2,0,1)");
1233  test_for_zero(t4ddg_1(0, 2, 0, 2)
1234  - (t4ddg_2(0, 2, 0, 2) - t4ddg_3(0, 2, 0, 2)),
1235  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(0,2,0,2)");
1236  test_for_zero(t4ddg_1(0, 2, 1, 0)
1237  - (t4ddg_2(0, 2, 1, 0) - t4ddg_3(0, 2, 1, 0)),
1238  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(0,2,1,0)");
1239  test_for_zero(t4ddg_1(0, 2, 1, 1)
1240  - (t4ddg_2(0, 2, 1, 1) - t4ddg_3(0, 2, 1, 1)),
1241  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(0,2,1,1)");
1242  test_for_zero(t4ddg_1(0, 2, 1, 2)
1243  - (t4ddg_2(0, 2, 1, 2) - t4ddg_3(0, 2, 1, 2)),
1244  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(0,2,1,2)");
1245  test_for_zero(t4ddg_1(0, 2, 2, 0)
1246  - (t4ddg_2(0, 2, 2, 0) - t4ddg_3(0, 2, 2, 0)),
1247  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(0,2,2,0)");
1248  test_for_zero(t4ddg_1(0, 2, 2, 1)
1249  - (t4ddg_2(0, 2, 2, 1) - t4ddg_3(0, 2, 2, 1)),
1250  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(0,2,2,1)");
1251  test_for_zero(t4ddg_1(0, 2, 2, 2)
1252  - (t4ddg_2(0, 2, 2, 2) - t4ddg_3(0, 2, 2, 2)),
1253  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(0,2,2,2)");
1254  test_for_zero(t4ddg_1(1, 0, 0, 0)
1255  - (t4ddg_2(1, 0, 0, 0) - t4ddg_3(1, 0, 0, 0)),
1256  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(1,0,0,0)");
1257  test_for_zero(t4ddg_1(1, 0, 0, 1)
1258  - (t4ddg_2(1, 0, 0, 1) - t4ddg_3(1, 0, 0, 1)),
1259  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(1,0,0,1)");
1260  test_for_zero(t4ddg_1(1, 0, 0, 2)
1261  - (t4ddg_2(1, 0, 0, 2) - t4ddg_3(1, 0, 0, 2)),
1262  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(1,0,0,2)");
1263  test_for_zero(t4ddg_1(1, 0, 1, 0)
1264  - (t4ddg_2(1, 0, 1, 0) - t4ddg_3(1, 0, 1, 0)),
1265  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(1,0,1,0)");
1266  test_for_zero(t4ddg_1(1, 0, 1, 1)
1267  - (t4ddg_2(1, 0, 1, 1) - t4ddg_3(1, 0, 1, 1)),
1268  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(1,0,1,1)");
1269  test_for_zero(t4ddg_1(1, 0, 1, 2)
1270  - (t4ddg_2(1, 0, 1, 2) - t4ddg_3(1, 0, 1, 2)),
1271  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(1,0,1,2)");
1272  test_for_zero(t4ddg_1(1, 0, 2, 0)
1273  - (t4ddg_2(1, 0, 2, 0) - t4ddg_3(1, 0, 2, 0)),
1274  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(1,0,2,0)");
1275  test_for_zero(t4ddg_1(1, 0, 2, 1)
1276  - (t4ddg_2(1, 0, 2, 1) - t4ddg_3(1, 0, 2, 1)),
1277  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(1,0,2,1)");
1278  test_for_zero(t4ddg_1(1, 0, 2, 2)
1279  - (t4ddg_2(1, 0, 2, 2) - t4ddg_3(1, 0, 2, 2)),
1280  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(1,0,2,2)");
1281  test_for_zero(t4ddg_1(1, 1, 0, 0)
1282  - (t4ddg_2(1, 1, 0, 0) - t4ddg_3(1, 1, 0, 0)),
1283  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(1,1,0,0)");
1284  test_for_zero(t4ddg_1(1, 1, 0, 1)
1285  - (t4ddg_2(1, 1, 0, 1) - t4ddg_3(1, 1, 0, 1)),
1286  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(1,1,0,1)");
1287  test_for_zero(t4ddg_1(1, 1, 0, 2)
1288  - (t4ddg_2(1, 1, 0, 2) - t4ddg_3(1, 1, 0, 2)),
1289  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(1,1,0,2)");
1290  test_for_zero(t4ddg_1(1, 1, 1, 0)
1291  - (t4ddg_2(1, 1, 1, 0) - t4ddg_3(1, 1, 1, 0)),
1292  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(1,1,1,0)");
1293  test_for_zero(t4ddg_1(1, 1, 1, 1)
1294  - (t4ddg_2(1, 1, 1, 1) - t4ddg_3(1, 1, 1, 1)),
1295  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(1,1,1,1)");
1296  test_for_zero(t4ddg_1(1, 1, 1, 2)
1297  - (t4ddg_2(1, 1, 1, 2) - t4ddg_3(1, 1, 1, 2)),
1298  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(1,1,1,2)");
1299  test_for_zero(t4ddg_1(1, 1, 2, 0)
1300  - (t4ddg_2(1, 1, 2, 0) - t4ddg_3(1, 1, 2, 0)),
1301  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(1,1,2,0)");
1302  test_for_zero(t4ddg_1(1, 1, 2, 1)
1303  - (t4ddg_2(1, 1, 2, 1) - t4ddg_3(1, 1, 2, 1)),
1304  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(1,1,2,1)");
1305  test_for_zero(t4ddg_1(1, 1, 2, 2)
1306  - (t4ddg_2(1, 1, 2, 2) - t4ddg_3(1, 1, 2, 2)),
1307  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(1,1,2,2)");
1308  test_for_zero(t4ddg_1(1, 2, 0, 0)
1309  - (t4ddg_2(1, 2, 0, 0) - t4ddg_3(1, 2, 0, 0)),
1310  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(1,2,0,0)");
1311  test_for_zero(t4ddg_1(1, 2, 0, 1)
1312  - (t4ddg_2(1, 2, 0, 1) - t4ddg_3(1, 2, 0, 1)),
1313  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(1,2,0,1)");
1314  test_for_zero(t4ddg_1(1, 2, 0, 2)
1315  - (t4ddg_2(1, 2, 0, 2) - t4ddg_3(1, 2, 0, 2)),
1316  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(1,2,0,2)");
1317  test_for_zero(t4ddg_1(1, 2, 1, 0)
1318  - (t4ddg_2(1, 2, 1, 0) - t4ddg_3(1, 2, 1, 0)),
1319  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(1,2,1,0)");
1320  test_for_zero(t4ddg_1(1, 2, 1, 1)
1321  - (t4ddg_2(1, 2, 1, 1) - t4ddg_3(1, 2, 1, 1)),
1322  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(1,2,1,1)");
1323  test_for_zero(t4ddg_1(1, 2, 1, 2)
1324  - (t4ddg_2(1, 2, 1, 2) - t4ddg_3(1, 2, 1, 2)),
1325  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(1,2,1,2)");
1326  test_for_zero(t4ddg_1(1, 2, 2, 0)
1327  - (t4ddg_2(1, 2, 2, 0) - t4ddg_3(1, 2, 2, 0)),
1328  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(1,2,2,0)");
1329  test_for_zero(t4ddg_1(1, 2, 2, 1)
1330  - (t4ddg_2(1, 2, 2, 1) - t4ddg_3(1, 2, 2, 1)),
1331  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(1,2,2,1)");
1332  test_for_zero(t4ddg_1(1, 2, 2, 2)
1333  - (t4ddg_2(1, 2, 2, 2) - t4ddg_3(1, 2, 2, 2)),
1334  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(1,2,2,2)");
1335  test_for_zero(t4ddg_1(2, 0, 0, 0)
1336  - (t4ddg_2(2, 0, 0, 0) - t4ddg_3(2, 0, 0, 0)),
1337  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(2,0,0,0)");
1338  test_for_zero(t4ddg_1(2, 0, 0, 1)
1339  - (t4ddg_2(2, 0, 0, 1) - t4ddg_3(2, 0, 0, 1)),
1340  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(2,0,0,1)");
1341  test_for_zero(t4ddg_1(2, 0, 0, 2)
1342  - (t4ddg_2(2, 0, 0, 2) - t4ddg_3(2, 0, 0, 2)),
1343  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(2,0,0,2)");
1344  test_for_zero(t4ddg_1(2, 0, 1, 0)
1345  - (t4ddg_2(2, 0, 1, 0) - t4ddg_3(2, 0, 1, 0)),
1346  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(2,0,1,0)");
1347  test_for_zero(t4ddg_1(2, 0, 1, 1)
1348  - (t4ddg_2(2, 0, 1, 1) - t4ddg_3(2, 0, 1, 1)),
1349  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(2,0,1,1)");
1350  test_for_zero(t4ddg_1(2, 0, 1, 2)
1351  - (t4ddg_2(2, 0, 1, 2) - t4ddg_3(2, 0, 1, 2)),
1352  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(2,0,1,2)");
1353  test_for_zero(t4ddg_1(2, 0, 2, 0)
1354  - (t4ddg_2(2, 0, 2, 0) - t4ddg_3(2, 0, 2, 0)),
1355  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(2,0,2,0)");
1356  test_for_zero(t4ddg_1(2, 0, 2, 1)
1357  - (t4ddg_2(2, 0, 2, 1) - t4ddg_3(2, 0, 2, 1)),
1358  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(2,0,2,1)");
1359  test_for_zero(t4ddg_1(2, 0, 2, 2)
1360  - (t4ddg_2(2, 0, 2, 2) - t4ddg_3(2, 0, 2, 2)),
1361  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(2,0,2,2)");
1362  test_for_zero(t4ddg_1(2, 1, 0, 0)
1363  - (t4ddg_2(2, 1, 0, 0) - t4ddg_3(2, 1, 0, 0)),
1364  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(2,1,0,0)");
1365  test_for_zero(t4ddg_1(2, 1, 0, 1)
1366  - (t4ddg_2(2, 1, 0, 1) - t4ddg_3(2, 1, 0, 1)),
1367  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(2,1,0,1)");
1368  test_for_zero(t4ddg_1(2, 1, 0, 2)
1369  - (t4ddg_2(2, 1, 0, 2) - t4ddg_3(2, 1, 0, 2)),
1370  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(2,1,0,2)");
1371  test_for_zero(t4ddg_1(2, 1, 1, 0)
1372  - (t4ddg_2(2, 1, 1, 0) - t4ddg_3(2, 1, 1, 0)),
1373  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(2,1,1,0)");
1374  test_for_zero(t4ddg_1(2, 1, 1, 1)
1375  - (t4ddg_2(2, 1, 1, 1) - t4ddg_3(2, 1, 1, 1)),
1376  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(2,1,1,1)");
1377  test_for_zero(t4ddg_1(2, 1, 1, 2)
1378  - (t4ddg_2(2, 1, 1, 2) - t4ddg_3(2, 1, 1, 2)),
1379  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(2,1,1,2)");
1380  test_for_zero(t4ddg_1(2, 1, 2, 0)
1381  - (t4ddg_2(2, 1, 2, 0) - t4ddg_3(2, 1, 2, 0)),
1382  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(2,1,2,0)");
1383  test_for_zero(t4ddg_1(2, 1, 2, 1)
1384  - (t4ddg_2(2, 1, 2, 1) - t4ddg_3(2, 1, 2, 1)),
1385  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(2,1,2,1)");
1386  test_for_zero(t4ddg_1(2, 1, 2, 2)
1387  - (t4ddg_2(2, 1, 2, 2) - t4ddg_3(2, 1, 2, 2)),
1388  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(2,1,2,2)");
1389  test_for_zero(t4ddg_1(2, 2, 0, 0)
1390  - (t4ddg_2(2, 2, 0, 0) - t4ddg_3(2, 2, 0, 0)),
1391  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(2,2,0,0)");
1392  test_for_zero(t4ddg_1(2, 2, 0, 1)
1393  - (t4ddg_2(2, 2, 0, 1) - t4ddg_3(2, 2, 0, 1)),
1394  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(2,2,0,1)");
1395  test_for_zero(t4ddg_1(2, 2, 0, 2)
1396  - (t4ddg_2(2, 2, 0, 2) - t4ddg_3(2, 2, 0, 2)),
1397  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(2,2,0,2)");
1398  test_for_zero(t4ddg_1(2, 2, 1, 0)
1399  - (t4ddg_2(2, 2, 1, 0) - t4ddg_3(2, 2, 1, 0)),
1400  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(2,2,1,0)");
1401  test_for_zero(t4ddg_1(2, 2, 1, 1)
1402  - (t4ddg_2(2, 2, 1, 1) - t4ddg_3(2, 2, 1, 1)),
1403  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(2,2,1,1)");
1404  test_for_zero(t4ddg_1(2, 2, 1, 2)
1405  - (t4ddg_2(2, 2, 1, 2) - t4ddg_3(2, 2, 1, 2)),
1406  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(2,2,1,2)");
1407  test_for_zero(t4ddg_1(2, 2, 2, 0)
1408  - (t4ddg_2(2, 2, 2, 0) - t4ddg_3(2, 2, 2, 0)),
1409  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(2,2,2,0)");
1410  test_for_zero(t4ddg_1(2, 2, 2, 1)
1411  - (t4ddg_2(2, 2, 2, 1) - t4ddg_3(2, 2, 2, 1)),
1412  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(2,2,2,1)");
1413  test_for_zero(t4ddg_1(2, 2, 2, 2)
1414  - (t4ddg_2(2, 2, 2, 2) - t4ddg_3(2, 2, 2, 2)),
1415  "T4ddg(i,j,k,l)-T4ddg(i,j,k,l)(2,2,2,2)");
1416 
1417  t4ddg_2(i, k, j, l) = t2s_2(i, k) * t2s_2(j, l);
1418  t4ddg_3(i, l, j, k) = t2s_2(i, l) * t2s_2(j, k);
1419  t4ddg_1(i, j, k, l) = (t4ddg_2(i, k, j, l) || t4ddg_3(i, l, j, k));
1420  test_for_zero(t4ddg_1(0, 0, 0, 0)
1421  - (t4ddg_2(0, 0, 0, 0) + t4ddg_3(0, 0, 0, 0)),
1422  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(0,0,0,0)");
1423  test_for_zero(t4ddg_1(0, 0, 0, 1)
1424  - (t4ddg_2(0, 0, 0, 1) + t4ddg_3(0, 1, 0, 0)),
1425  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(0,0,0,1)");
1426  test_for_zero(t4ddg_1(0, 0, 0, 2)
1427  - (t4ddg_2(0, 0, 0, 2) + t4ddg_3(0, 2, 0, 0)),
1428  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(0,0,0,2)");
1429  test_for_zero(t4ddg_1(0, 1, 0, 0)
1430  - (t4ddg_2(0, 0, 1, 0) + t4ddg_3(0, 0, 1, 0)),
1431  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(0,0,1,0)");
1432  test_for_zero(t4ddg_1(0, 1, 0, 1)
1433  - (t4ddg_2(0, 0, 1, 1) + t4ddg_3(0, 1, 1, 0)),
1434  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(0,0,1,1)");
1435  test_for_zero(t4ddg_1(0, 1, 0, 2)
1436  - (t4ddg_2(0, 0, 1, 2) + t4ddg_3(0, 2, 1, 0)),
1437  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(0,0,1,2)");
1438  test_for_zero(t4ddg_1(0, 2, 0, 0)
1439  - (t4ddg_2(0, 0, 2, 0) + t4ddg_3(0, 0, 2, 0)),
1440  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(0,0,2,0)");
1441  test_for_zero(t4ddg_1(0, 2, 0, 1)
1442  - (t4ddg_2(0, 0, 2, 1) + t4ddg_3(0, 1, 2, 0)),
1443  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(0,0,2,1)");
1444  test_for_zero(t4ddg_1(0, 2, 0, 2)
1445  - (t4ddg_2(0, 0, 2, 2) + t4ddg_3(0, 2, 2, 0)),
1446  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(0,0,2,2)");
1447  test_for_zero(t4ddg_1(0, 0, 1, 0)
1448  - (t4ddg_2(0, 1, 0, 0) + t4ddg_3(0, 0, 0, 1)),
1449  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(0,1,0,0)");
1450  test_for_zero(t4ddg_1(0, 0, 1, 1)
1451  - (t4ddg_2(0, 1, 0, 1) + t4ddg_3(0, 1, 0, 1)),
1452  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(0,1,0,1)");
1453  test_for_zero(t4ddg_1(0, 0, 1, 2)
1454  - (t4ddg_2(0, 1, 0, 2) + t4ddg_3(0, 2, 0, 1)),
1455  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(0,1,0,2)");
1456  test_for_zero(t4ddg_1(0, 1, 1, 0)
1457  - (t4ddg_2(0, 1, 1, 0) + t4ddg_3(0, 0, 1, 1)),
1458  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(0,1,1,0)");
1459  test_for_zero(t4ddg_1(0, 1, 1, 1)
1460  - (t4ddg_2(0, 1, 1, 1) + t4ddg_3(0, 1, 1, 1)),
1461  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(0,1,1,1)");
1462  test_for_zero(t4ddg_1(0, 1, 1, 2)
1463  - (t4ddg_2(0, 1, 1, 2) + t4ddg_3(0, 2, 1, 1)),
1464  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(0,1,1,2)");
1465  test_for_zero(t4ddg_1(0, 2, 1, 0)
1466  - (t4ddg_2(0, 1, 2, 0) + t4ddg_3(0, 0, 2, 1)),
1467  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(0,1,2,0)");
1468  test_for_zero(t4ddg_1(0, 2, 1, 1)
1469  - (t4ddg_2(0, 1, 2, 1) + t4ddg_3(0, 1, 2, 1)),
1470  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(0,1,2,1)");
1471  test_for_zero(t4ddg_1(0, 2, 1, 2)
1472  - (t4ddg_2(0, 1, 2, 2) + t4ddg_3(0, 2, 2, 1)),
1473  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(0,1,2,2)");
1474  test_for_zero(t4ddg_1(0, 0, 2, 0)
1475  - (t4ddg_2(0, 2, 0, 0) + t4ddg_3(0, 0, 0, 2)),
1476  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(0,2,0,0)");
1477  test_for_zero(t4ddg_1(0, 0, 2, 1)
1478  - (t4ddg_2(0, 2, 0, 1) + t4ddg_3(0, 1, 0, 2)),
1479  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(0,2,0,1)");
1480  test_for_zero(t4ddg_1(0, 0, 2, 2)
1481  - (t4ddg_2(0, 2, 0, 2) + t4ddg_3(0, 2, 0, 2)),
1482  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(0,2,0,2)");
1483  test_for_zero(t4ddg_1(0, 1, 2, 0)
1484  - (t4ddg_2(0, 2, 1, 0) + t4ddg_3(0, 0, 1, 2)),
1485  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(0,2,1,0)");
1486  test_for_zero(t4ddg_1(0, 1, 2, 1)
1487  - (t4ddg_2(0, 2, 1, 1) + t4ddg_3(0, 1, 1, 2)),
1488  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(0,2,1,1)");
1489  test_for_zero(t4ddg_1(0, 1, 2, 2)
1490  - (t4ddg_2(0, 2, 1, 2) + t4ddg_3(0, 2, 1, 2)),
1491  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(0,2,1,2)");
1492  test_for_zero(t4ddg_1(0, 2, 2, 0)
1493  - (t4ddg_2(0, 2, 2, 0) + t4ddg_3(0, 0, 2, 2)),
1494  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(0,2,2,0)");
1495  test_for_zero(t4ddg_1(0, 2, 2, 1)
1496  - (t4ddg_2(0, 2, 2, 1) + t4ddg_3(0, 1, 2, 2)),
1497  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(0,2,2,1)");
1498  test_for_zero(t4ddg_1(0, 2, 2, 2)
1499  - (t4ddg_2(0, 2, 2, 2) + t4ddg_3(0, 2, 2, 2)),
1500  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(0,2,2,2)");
1501  test_for_zero(t4ddg_1(1, 0, 0, 0)
1502  - (t4ddg_2(1, 0, 0, 0) + t4ddg_3(1, 0, 0, 0)),
1503  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(1,0,0,0)");
1504  test_for_zero(t4ddg_1(1, 0, 0, 1)
1505  - (t4ddg_2(1, 0, 0, 1) + t4ddg_3(1, 1, 0, 0)),
1506  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(1,0,0,1)");
1507  test_for_zero(t4ddg_1(1, 0, 0, 2)
1508  - (t4ddg_2(1, 0, 0, 2) + t4ddg_3(1, 2, 0, 0)),
1509  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(1,0,0,2)");
1510  test_for_zero(t4ddg_1(1, 1, 0, 0)
1511  - (t4ddg_2(1, 0, 1, 0) + t4ddg_3(1, 0, 1, 0)),
1512  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(1,0,1,0)");
1513  test_for_zero(t4ddg_1(1, 1, 0, 1)
1514  - (t4ddg_2(1, 0, 1, 1) + t4ddg_3(1, 1, 1, 0)),
1515  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(1,0,1,1)");
1516  test_for_zero(t4ddg_1(1, 1, 0, 2)
1517  - (t4ddg_2(1, 0, 1, 2) + t4ddg_3(1, 2, 1, 0)),
1518  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(1,0,1,2)");
1519  test_for_zero(t4ddg_1(1, 2, 0, 0)
1520  - (t4ddg_2(1, 0, 2, 0) + t4ddg_3(1, 0, 2, 0)),
1521  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(1,0,2,0)");
1522  test_for_zero(t4ddg_1(1, 2, 0, 1)
1523  - (t4ddg_2(1, 0, 2, 1) + t4ddg_3(1, 1, 2, 0)),
1524  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(1,0,2,1)");
1525  test_for_zero(t4ddg_1(1, 2, 0, 2)
1526  - (t4ddg_2(1, 0, 2, 2) + t4ddg_3(1, 2, 2, 0)),
1527  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(1,0,2,2)");
1528  test_for_zero(t4ddg_1(1, 0, 1, 0)
1529  - (t4ddg_2(1, 1, 0, 0) + t4ddg_3(1, 0, 0, 1)),
1530  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(1,1,0,0)");
1531  test_for_zero(t4ddg_1(1, 0, 1, 1)
1532  - (t4ddg_2(1, 1, 0, 1) + t4ddg_3(1, 1, 0, 1)),
1533  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(1,1,0,1)");
1534  test_for_zero(t4ddg_1(1, 0, 1, 2)
1535  - (t4ddg_2(1, 1, 0, 2) + t4ddg_3(1, 2, 0, 1)),
1536  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(1,1,0,2)");
1537  test_for_zero(t4ddg_1(1, 1, 1, 0)
1538  - (t4ddg_2(1, 1, 1, 0) + t4ddg_3(1, 0, 1, 1)),
1539  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(1,1,1,0)");
1540  test_for_zero(t4ddg_1(1, 1, 1, 1)
1541  - (t4ddg_2(1, 1, 1, 1) + t4ddg_3(1, 1, 1, 1)),
1542  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(1,1,1,1)");
1543  test_for_zero(t4ddg_1(1, 1, 1, 2)
1544  - (t4ddg_2(1, 1, 1, 2) + t4ddg_3(1, 2, 1, 1)),
1545  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(1,1,1,2)");
1546  test_for_zero(t4ddg_1(1, 2, 1, 0)
1547  - (t4ddg_2(1, 1, 2, 0) + t4ddg_3(1, 0, 2, 1)),
1548  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(1,1,2,0)");
1549  test_for_zero(t4ddg_1(1, 2, 1, 1)
1550  - (t4ddg_2(1, 1, 2, 1) + t4ddg_3(1, 1, 2, 1)),
1551  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(1,1,2,1)");
1552  test_for_zero(t4ddg_1(1, 2, 1, 2)
1553  - (t4ddg_2(1, 1, 2, 2) + t4ddg_3(1, 2, 2, 1)),
1554  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(1,1,2,2)");
1555  test_for_zero(t4ddg_1(1, 0, 2, 0)
1556  - (t4ddg_2(1, 2, 0, 0) + t4ddg_3(1, 0, 0, 2)),
1557  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(1,2,0,0)");
1558  test_for_zero(t4ddg_1(1, 0, 2, 1)
1559  - (t4ddg_2(1, 2, 0, 1) + t4ddg_3(1, 1, 0, 2)),
1560  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(1,2,0,1)");
1561  test_for_zero(t4ddg_1(1, 0, 2, 2)
1562  - (t4ddg_2(1, 2, 0, 2) + t4ddg_3(1, 2, 0, 2)),
1563  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(1,2,0,2)");
1564  test_for_zero(t4ddg_1(1, 1, 2, 0)
1565  - (t4ddg_2(1, 2, 1, 0) + t4ddg_3(1, 0, 1, 2)),
1566  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(1,2,1,0)");
1567  test_for_zero(t4ddg_1(1, 1, 2, 1)
1568  - (t4ddg_2(1, 2, 1, 1) + t4ddg_3(1, 1, 1, 2)),
1569  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(1,2,1,1)");
1570  test_for_zero(t4ddg_1(1, 1, 2, 2)
1571  - (t4ddg_2(1, 2, 1, 2) + t4ddg_3(1, 2, 1, 2)),
1572  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(1,2,1,2)");
1573  test_for_zero(t4ddg_1(1, 2, 2, 0)
1574  - (t4ddg_2(1, 2, 2, 0) + t4ddg_3(1, 0, 2, 2)),
1575  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(1,2,2,0)");
1576  test_for_zero(t4ddg_1(1, 2, 2, 1)
1577  - (t4ddg_2(1, 2, 2, 1) + t4ddg_3(1, 1, 2, 2)),
1578  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(1,2,2,1)");
1579  test_for_zero(t4ddg_1(1, 2, 2, 2)
1580  - (t4ddg_2(1, 2, 2, 2) + t4ddg_3(1, 2, 2, 2)),
1581  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(1,2,2,2)");
1582  test_for_zero(t4ddg_1(2, 0, 0, 0)
1583  - (t4ddg_2(2, 0, 0, 0) + t4ddg_3(2, 0, 0, 0)),
1584  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(2,0,0,0)");
1585  test_for_zero(t4ddg_1(2, 0, 0, 1)
1586  - (t4ddg_2(2, 0, 0, 1) + t4ddg_3(2, 1, 0, 0)),
1587  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(2,0,0,1)");
1588  test_for_zero(t4ddg_1(2, 0, 0, 2)
1589  - (t4ddg_2(2, 0, 0, 2) + t4ddg_3(2, 2, 0, 0)),
1590  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(2,0,0,2)");
1591  test_for_zero(t4ddg_1(2, 1, 0, 0)
1592  - (t4ddg_2(2, 0, 1, 0) + t4ddg_3(2, 0, 1, 0)),
1593  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(2,0,1,0)");
1594  test_for_zero(t4ddg_1(2, 1, 0, 1)
1595  - (t4ddg_2(2, 0, 1, 1) + t4ddg_3(2, 1, 1, 0)),
1596  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(2,0,1,1)");
1597  test_for_zero(t4ddg_1(2, 1, 0, 2)
1598  - (t4ddg_2(2, 0, 1, 2) + t4ddg_3(2, 2, 1, 0)),
1599  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(2,0,1,2)");
1600  test_for_zero(t4ddg_1(2, 2, 0, 0)
1601  - (t4ddg_2(2, 0, 2, 0) + t4ddg_3(2, 0, 2, 0)),
1602  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(2,0,2,0)");
1603  test_for_zero(t4ddg_1(2, 2, 0, 1)
1604  - (t4ddg_2(2, 0, 2, 1) + t4ddg_3(2, 1, 2, 0)),
1605  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(2,0,2,1)");
1606  test_for_zero(t4ddg_1(2, 2, 0, 2)
1607  - (t4ddg_2(2, 0, 2, 2) + t4ddg_3(2, 2, 2, 0)),
1608  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(2,0,2,2)");
1609  test_for_zero(t4ddg_1(2, 0, 1, 0)
1610  - (t4ddg_2(2, 1, 0, 0) + t4ddg_3(2, 0, 0, 1)),
1611  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(2,1,0,0)");
1612  test_for_zero(t4ddg_1(2, 0, 1, 1)
1613  - (t4ddg_2(2, 1, 0, 1) + t4ddg_3(2, 1, 0, 1)),
1614  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(2,1,0,1)");
1615  test_for_zero(t4ddg_1(2, 0, 1, 2)
1616  - (t4ddg_2(2, 1, 0, 2) + t4ddg_3(2, 2, 0, 1)),
1617  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(2,1,0,2)");
1618  test_for_zero(t4ddg_1(2, 1, 1, 0)
1619  - (t4ddg_2(2, 1, 1, 0) + t4ddg_3(2, 0, 1, 1)),
1620  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(2,1,1,0)");
1621  test_for_zero(t4ddg_1(2, 1, 1, 1)
1622  - (t4ddg_2(2, 1, 1, 1) + t4ddg_3(2, 1, 1, 1)),
1623  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(2,1,1,1)");
1624  test_for_zero(t4ddg_1(2, 1, 1, 2)
1625  - (t4ddg_2(2, 1, 1, 2) + t4ddg_3(2, 2, 1, 1)),
1626  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(2,1,1,2)");
1627  test_for_zero(t4ddg_1(2, 2, 1, 0)
1628  - (t4ddg_2(2, 1, 2, 0) + t4ddg_3(2, 0, 2, 1)),
1629  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(2,1,2,0)");
1630  test_for_zero(t4ddg_1(2, 2, 1, 1)
1631  - (t4ddg_2(2, 1, 2, 1) + t4ddg_3(2, 1, 2, 1)),
1632  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(2,1,2,1)");
1633  test_for_zero(t4ddg_1(2, 2, 1, 2)
1634  - (t4ddg_2(2, 1, 2, 2) + t4ddg_3(2, 2, 2, 1)),
1635  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(2,1,2,2)");
1636  test_for_zero(t4ddg_1(2, 0, 2, 0)
1637  - (t4ddg_2(2, 2, 0, 0) + t4ddg_3(2, 0, 0, 2)),
1638  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(2,2,0,0)");
1639  test_for_zero(t4ddg_1(2, 0, 2, 1)
1640  - (t4ddg_2(2, 2, 0, 1) + t4ddg_3(2, 1, 0, 2)),
1641  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(2,2,0,1)");
1642  test_for_zero(t4ddg_1(2, 0, 2, 2)
1643  - (t4ddg_2(2, 2, 0, 2) + t4ddg_3(2, 2, 0, 2)),
1644  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(2,2,0,2)");
1645  test_for_zero(t4ddg_1(2, 1, 2, 0)
1646  - (t4ddg_2(2, 2, 1, 0) + t4ddg_3(2, 0, 1, 2)),
1647  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(2,2,1,0)");
1648  test_for_zero(t4ddg_1(2, 1, 2, 1)
1649  - (t4ddg_2(2, 2, 1, 1) + t4ddg_3(2, 1, 1, 2)),
1650  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(2,2,1,1)");
1651  test_for_zero(t4ddg_1(2, 1, 2, 2)
1652  - (t4ddg_2(2, 2, 1, 2) + t4ddg_3(2, 2, 1, 2)),
1653  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(2,2,1,2)");
1654  test_for_zero(t4ddg_1(2, 2, 2, 0)
1655  - (t4ddg_2(2, 2, 2, 0) + t4ddg_3(2, 0, 2, 2)),
1656  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(2,2,2,0)");
1657  test_for_zero(t4ddg_1(2, 2, 2, 1)
1658  - (t4ddg_2(2, 2, 2, 1) + t4ddg_3(2, 1, 2, 2)),
1659  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(2,2,2,1)");
1660  test_for_zero(t4ddg_1(2, 2, 2, 2)
1661  - (t4ddg_2(2, 2, 2, 2) + t4ddg_3(2, 2, 2, 2)),
1662  "T4ddg(i,k,j,l)||T4ddg(i,l,j,k)(2,2,2,2)");
1663 }
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
i
FTensor::Index< 'i', SPACE_DIM > i
Definition: hcurl_divergence_operator_2d.cpp:27
FTensor::Index< 'i', 3 >
FTensor::Dg
Definition: Dg_value.hpp:9
std
Definition: enable_if.hpp:5
j
FTensor::Index< 'j', 3 > j
Definition: matrix_function.cpp:19
FTensor::Ddg< double, 3, 3 >
test_T4ddg_01
void test_T4ddg_01(const Tensor2_symmetric< double, 3 > &t2s_2, const Tensor2_symmetric< double, 3 > &t2s_3, const Dg< double, 3, 3 > &t3dg_2, const Dg< double, 3, 3 > &t3dg_3)
Definition: test_T4ddg_01.cpp:7
m
FTensor::Index< 'm', 3 > m
Definition: shallow_wave.cpp:80
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