v0.14.0
test_T4ddg_03.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_03(Tensor2<double, 3, 3> &t2_1,
8  Tensor2_symmetric<double, 3> &t2s_1,
9  const Tensor2_symmetric<double, 3> &t2s_2,
10  const Tensor2_symmetric<double, 3> &t2s_3)
11 {
12  Index<'i', 3> i;
13  Index<'j', 3> j;
14  Index<'k', 3> k;
15  Index<'l', 3> l;
16 
17  Number<0> N0;
18  Number<1> N1;
19  Number<2> N2;
20 
21  Ddg<double, 3, 3> t4ddg_1, t4ddg_2;
22 
23  t4ddg_2(i, k, j, l) = t2s_2(i, k) * t2s_3(j, l);
24  t2_1(i, k) = t4ddg_2(i, j, k, l) * t2s_2(j, l);
25  test_for_zero(t2_1(0, 0)
26  - (t4ddg_2(0, 0, 0, 0) * t2s_2(0, 0)
27  + t4ddg_2(0, 0, 0, 1) * t2s_2(0, 1)
28  + t4ddg_2(0, 0, 0, 2) * t2s_2(0, 2)
29  + t4ddg_2(0, 1, 0, 0) * t2s_2(1, 0)
30  + t4ddg_2(0, 1, 0, 1) * t2s_2(1, 1)
31  + t4ddg_2(0, 1, 0, 2) * t2s_2(1, 2)
32  + t4ddg_2(0, 2, 0, 0) * t2s_2(2, 0)
33  + t4ddg_2(0, 2, 0, 1) * t2s_2(2, 1)
34  + t4ddg_2(0, 2, 0, 2) * t2s_2(2, 2)),
35  "T4ddg(i,j,k,l)*T2s(j,l)(0,0)");
36  test_for_zero(t2_1(0, 1)
37  - (t4ddg_2(0, 0, 1, 0) * t2s_2(0, 0)
38  + t4ddg_2(0, 0, 1, 1) * t2s_2(0, 1)
39  + t4ddg_2(0, 0, 1, 2) * t2s_2(0, 2)
40  + t4ddg_2(0, 1, 1, 0) * t2s_2(1, 0)
41  + t4ddg_2(0, 1, 1, 1) * t2s_2(1, 1)
42  + t4ddg_2(0, 1, 1, 2) * t2s_2(1, 2)
43  + t4ddg_2(0, 2, 1, 0) * t2s_2(2, 0)
44  + t4ddg_2(0, 2, 1, 1) * t2s_2(2, 1)
45  + t4ddg_2(0, 2, 1, 2) * t2s_2(2, 2)),
46  "T4ddg(i,j,k,l)*T2s(j,l)(0,1)");
47  test_for_zero(t2_1(0, 2)
48  - (t4ddg_2(0, 0, 2, 0) * t2s_2(0, 0)
49  + t4ddg_2(0, 0, 2, 1) * t2s_2(0, 1)
50  + t4ddg_2(0, 0, 2, 2) * t2s_2(0, 2)
51  + t4ddg_2(0, 1, 2, 0) * t2s_2(1, 0)
52  + t4ddg_2(0, 1, 2, 1) * t2s_2(1, 1)
53  + t4ddg_2(0, 1, 2, 2) * t2s_2(1, 2)
54  + t4ddg_2(0, 2, 2, 0) * t2s_2(2, 0)
55  + t4ddg_2(0, 2, 2, 1) * t2s_2(2, 1)
56  + t4ddg_2(0, 2, 2, 2) * t2s_2(2, 2)),
57  "T4ddg(i,j,k,l)*T2s(j,l)(0,2)");
58  test_for_zero(t2_1(1, 0)
59  - (t4ddg_2(1, 0, 0, 0) * t2s_2(0, 0)
60  + t4ddg_2(1, 0, 0, 1) * t2s_2(0, 1)
61  + t4ddg_2(1, 0, 0, 2) * t2s_2(0, 2)
62  + t4ddg_2(1, 1, 0, 0) * t2s_2(1, 0)
63  + t4ddg_2(1, 1, 0, 1) * t2s_2(1, 1)
64  + t4ddg_2(1, 1, 0, 2) * t2s_2(1, 2)
65  + t4ddg_2(1, 2, 0, 0) * t2s_2(2, 0)
66  + t4ddg_2(1, 2, 0, 1) * t2s_2(2, 1)
67  + t4ddg_2(1, 2, 0, 2) * t2s_2(2, 2)),
68  "T4ddg(i,j,k,l)*T2s(j,l)(1,0)");
69  test_for_zero(t2_1(1, 1)
70  - (t4ddg_2(1, 0, 1, 0) * t2s_2(0, 0)
71  + t4ddg_2(1, 0, 1, 1) * t2s_2(0, 1)
72  + t4ddg_2(1, 0, 1, 2) * t2s_2(0, 2)
73  + t4ddg_2(1, 1, 1, 0) * t2s_2(1, 0)
74  + t4ddg_2(1, 1, 1, 1) * t2s_2(1, 1)
75  + t4ddg_2(1, 1, 1, 2) * t2s_2(1, 2)
76  + t4ddg_2(1, 2, 1, 0) * t2s_2(2, 0)
77  + t4ddg_2(1, 2, 1, 1) * t2s_2(2, 1)
78  + t4ddg_2(1, 2, 1, 2) * t2s_2(2, 2)),
79  "T4ddg(i,j,k,l)*T2s(j,l)(1,1)");
80  test_for_zero(t2_1(1, 2)
81  - (t4ddg_2(1, 0, 2, 0) * t2s_2(0, 0)
82  + t4ddg_2(1, 0, 2, 1) * t2s_2(0, 1)
83  + t4ddg_2(1, 0, 2, 2) * t2s_2(0, 2)
84  + t4ddg_2(1, 1, 2, 0) * t2s_2(1, 0)
85  + t4ddg_2(1, 1, 2, 1) * t2s_2(1, 1)
86  + t4ddg_2(1, 1, 2, 2) * t2s_2(1, 2)
87  + t4ddg_2(1, 2, 2, 0) * t2s_2(2, 0)
88  + t4ddg_2(1, 2, 2, 1) * t2s_2(2, 1)
89  + t4ddg_2(1, 2, 2, 2) * t2s_2(2, 2)),
90  "T4ddg(i,j,k,l)*T2s(j,l)(1,2)");
91  test_for_zero(t2_1(2, 0)
92  - (t4ddg_2(2, 0, 0, 0) * t2s_2(0, 0)
93  + t4ddg_2(2, 0, 0, 1) * t2s_2(0, 1)
94  + t4ddg_2(2, 0, 0, 2) * t2s_2(0, 2)
95  + t4ddg_2(2, 1, 0, 0) * t2s_2(1, 0)
96  + t4ddg_2(2, 1, 0, 1) * t2s_2(1, 1)
97  + t4ddg_2(2, 1, 0, 2) * t2s_2(1, 2)
98  + t4ddg_2(2, 2, 0, 0) * t2s_2(2, 0)
99  + t4ddg_2(2, 2, 0, 1) * t2s_2(2, 1)
100  + t4ddg_2(2, 2, 0, 2) * t2s_2(2, 2)),
101  "T4ddg(i,j,k,l)*T2s(j,l)(2,0)");
102  test_for_zero(t2_1(2, 1)
103  - (t4ddg_2(2, 0, 1, 0) * t2s_2(0, 0)
104  + t4ddg_2(2, 0, 1, 1) * t2s_2(0, 1)
105  + t4ddg_2(2, 0, 1, 2) * t2s_2(0, 2)
106  + t4ddg_2(2, 1, 1, 0) * t2s_2(1, 0)
107  + t4ddg_2(2, 1, 1, 1) * t2s_2(1, 1)
108  + t4ddg_2(2, 1, 1, 2) * t2s_2(1, 2)
109  + t4ddg_2(2, 2, 1, 0) * t2s_2(2, 0)
110  + t4ddg_2(2, 2, 1, 1) * t2s_2(2, 1)
111  + t4ddg_2(2, 2, 1, 2) * t2s_2(2, 2)),
112  "T4ddg(i,j,k,l)*T2s(j,l)(2,1)");
113  test_for_zero(t2_1(2, 2)
114  - (t4ddg_2(2, 0, 2, 0) * t2s_2(0, 0)
115  + t4ddg_2(2, 0, 2, 1) * t2s_2(0, 1)
116  + t4ddg_2(2, 0, 2, 2) * t2s_2(0, 2)
117  + t4ddg_2(2, 1, 2, 0) * t2s_2(1, 0)
118  + t4ddg_2(2, 1, 2, 1) * t2s_2(1, 1)
119  + t4ddg_2(2, 1, 2, 2) * t2s_2(1, 2)
120  + t4ddg_2(2, 2, 2, 0) * t2s_2(2, 0)
121  + t4ddg_2(2, 2, 2, 1) * t2s_2(2, 1)
122  + t4ddg_2(2, 2, 2, 2) * t2s_2(2, 2)),
123  "T4ddg(i,j,k,l)*T2s(j,l)(2,2)");
124 
125  t2_1(i, k) = t2s_3(j, l) * t4ddg_2(i, j, k, l);
126  test_for_zero(t2_1(0, 0)
127  - (t4ddg_2(0, 0, 0, 0) * t2s_3(0, 0)
128  + t4ddg_2(0, 0, 0, 1) * t2s_3(0, 1)
129  + t4ddg_2(0, 0, 0, 2) * t2s_3(0, 2)
130  + t4ddg_2(0, 1, 0, 0) * t2s_3(1, 0)
131  + t4ddg_2(0, 1, 0, 1) * t2s_3(1, 1)
132  + t4ddg_2(0, 1, 0, 2) * t2s_3(1, 2)
133  + t4ddg_2(0, 2, 0, 0) * t2s_3(2, 0)
134  + t4ddg_2(0, 2, 0, 1) * t2s_3(2, 1)
135  + t4ddg_2(0, 2, 0, 2) * t2s_3(2, 2)),
136  "T2s(j,l)*T4ddg(i,j,k,l)(0,0)");
137  test_for_zero(t2_1(0, 1)
138  - (t4ddg_2(0, 0, 1, 0) * t2s_3(0, 0)
139  + t4ddg_2(0, 0, 1, 1) * t2s_3(0, 1)
140  + t4ddg_2(0, 0, 1, 2) * t2s_3(0, 2)
141  + t4ddg_2(0, 1, 1, 0) * t2s_3(1, 0)
142  + t4ddg_2(0, 1, 1, 1) * t2s_3(1, 1)
143  + t4ddg_2(0, 1, 1, 2) * t2s_3(1, 2)
144  + t4ddg_2(0, 2, 1, 0) * t2s_3(2, 0)
145  + t4ddg_2(0, 2, 1, 1) * t2s_3(2, 1)
146  + t4ddg_2(0, 2, 1, 2) * t2s_3(2, 2)),
147  "T2s(j,l)*T4ddg(i,j,k,l)(0,1)");
148  test_for_zero(t2_1(0, 2)
149  - (t4ddg_2(0, 0, 2, 0) * t2s_3(0, 0)
150  + t4ddg_2(0, 0, 2, 1) * t2s_3(0, 1)
151  + t4ddg_2(0, 0, 2, 2) * t2s_3(0, 2)
152  + t4ddg_2(0, 1, 2, 0) * t2s_3(1, 0)
153  + t4ddg_2(0, 1, 2, 1) * t2s_3(1, 1)
154  + t4ddg_2(0, 1, 2, 2) * t2s_3(1, 2)
155  + t4ddg_2(0, 2, 2, 0) * t2s_3(2, 0)
156  + t4ddg_2(0, 2, 2, 1) * t2s_3(2, 1)
157  + t4ddg_2(0, 2, 2, 2) * t2s_3(2, 2)),
158  "T2s(j,l)*T4ddg(i,j,k,l)(0,2)");
159  test_for_zero(t2_1(1, 0)
160  - (t4ddg_2(1, 0, 0, 0) * t2s_3(0, 0)
161  + t4ddg_2(1, 0, 0, 1) * t2s_3(0, 1)
162  + t4ddg_2(1, 0, 0, 2) * t2s_3(0, 2)
163  + t4ddg_2(1, 1, 0, 0) * t2s_3(1, 0)
164  + t4ddg_2(1, 1, 0, 1) * t2s_3(1, 1)
165  + t4ddg_2(1, 1, 0, 2) * t2s_3(1, 2)
166  + t4ddg_2(1, 2, 0, 0) * t2s_3(2, 0)
167  + t4ddg_2(1, 2, 0, 1) * t2s_3(2, 1)
168  + t4ddg_2(1, 2, 0, 2) * t2s_3(2, 2)),
169  "T2s(j,l)*T4ddg(i,j,k,l)(1,0)");
170  test_for_zero(t2_1(1, 1)
171  - (t4ddg_2(1, 0, 1, 0) * t2s_3(0, 0)
172  + t4ddg_2(1, 0, 1, 1) * t2s_3(0, 1)
173  + t4ddg_2(1, 0, 1, 2) * t2s_3(0, 2)
174  + t4ddg_2(1, 1, 1, 0) * t2s_3(1, 0)
175  + t4ddg_2(1, 1, 1, 1) * t2s_3(1, 1)
176  + t4ddg_2(1, 1, 1, 2) * t2s_3(1, 2)
177  + t4ddg_2(1, 2, 1, 0) * t2s_3(2, 0)
178  + t4ddg_2(1, 2, 1, 1) * t2s_3(2, 1)
179  + t4ddg_2(1, 2, 1, 2) * t2s_3(2, 2)),
180  "T2s(j,l)*T4ddg(i,j,k,l)(1,1)");
181  test_for_zero(t2_1(1, 2)
182  - (t4ddg_2(1, 0, 2, 0) * t2s_3(0, 0)
183  + t4ddg_2(1, 0, 2, 1) * t2s_3(0, 1)
184  + t4ddg_2(1, 0, 2, 2) * t2s_3(0, 2)
185  + t4ddg_2(1, 1, 2, 0) * t2s_3(1, 0)
186  + t4ddg_2(1, 1, 2, 1) * t2s_3(1, 1)
187  + t4ddg_2(1, 1, 2, 2) * t2s_3(1, 2)
188  + t4ddg_2(1, 2, 2, 0) * t2s_3(2, 0)
189  + t4ddg_2(1, 2, 2, 1) * t2s_3(2, 1)
190  + t4ddg_2(1, 2, 2, 2) * t2s_3(2, 2)),
191  "T2s(j,l)*T4ddg(i,j,k,l)(1,2)");
192  test_for_zero(t2_1(2, 0)
193  - (t4ddg_2(2, 0, 0, 0) * t2s_3(0, 0)
194  + t4ddg_2(2, 0, 0, 1) * t2s_3(0, 1)
195  + t4ddg_2(2, 0, 0, 2) * t2s_3(0, 2)
196  + t4ddg_2(2, 1, 0, 0) * t2s_3(1, 0)
197  + t4ddg_2(2, 1, 0, 1) * t2s_3(1, 1)
198  + t4ddg_2(2, 1, 0, 2) * t2s_3(1, 2)
199  + t4ddg_2(2, 2, 0, 0) * t2s_3(2, 0)
200  + t4ddg_2(2, 2, 0, 1) * t2s_3(2, 1)
201  + t4ddg_2(2, 2, 0, 2) * t2s_3(2, 2)),
202  "T2s(j,l)*T4ddg(i,j,k,l)(2,0)");
203  test_for_zero(t2_1(2, 1)
204  - (t4ddg_2(2, 0, 1, 0) * t2s_3(0, 0)
205  + t4ddg_2(2, 0, 1, 1) * t2s_3(0, 1)
206  + t4ddg_2(2, 0, 1, 2) * t2s_3(0, 2)
207  + t4ddg_2(2, 1, 1, 0) * t2s_3(1, 0)
208  + t4ddg_2(2, 1, 1, 1) * t2s_3(1, 1)
209  + t4ddg_2(2, 1, 1, 2) * t2s_3(1, 2)
210  + t4ddg_2(2, 2, 1, 0) * t2s_3(2, 0)
211  + t4ddg_2(2, 2, 1, 1) * t2s_3(2, 1)
212  + t4ddg_2(2, 2, 1, 2) * t2s_3(2, 2)),
213  "T2s(j,l)*T4ddg(i,j,k,l)(2,1)");
214  test_for_zero(t2_1(2, 2)
215  - (t4ddg_2(2, 0, 2, 0) * t2s_3(0, 0)
216  + t4ddg_2(2, 0, 2, 1) * t2s_3(0, 1)
217  + t4ddg_2(2, 0, 2, 2) * t2s_3(0, 2)
218  + t4ddg_2(2, 1, 2, 0) * t2s_3(1, 0)
219  + t4ddg_2(2, 1, 2, 1) * t2s_3(1, 1)
220  + t4ddg_2(2, 1, 2, 2) * t2s_3(1, 2)
221  + t4ddg_2(2, 2, 2, 0) * t2s_3(2, 0)
222  + t4ddg_2(2, 2, 2, 1) * t2s_3(2, 1)
223  + t4ddg_2(2, 2, 2, 2) * t2s_3(2, 2)),
224  "T2s(j,l)*T4ddg(i,j,k,l)(2,2)");
225 
226  t2s_1(j, l) = t4ddg_2(i, k, j, l) * t2s_2(i, k);
227  test_for_zero(t2s_1(0, 0)
228  - (t4ddg_2(0, 0, 0, 0) * t2s_2(0, 0)
229  + t4ddg_2(0, 1, 0, 0) * t2s_2(0, 1)
230  + t4ddg_2(0, 2, 0, 0) * t2s_2(0, 2)
231  + t4ddg_2(1, 0, 0, 0) * t2s_2(1, 0)
232  + t4ddg_2(1, 1, 0, 0) * t2s_2(1, 1)
233  + t4ddg_2(1, 2, 0, 0) * t2s_2(1, 2)
234  + t4ddg_2(2, 0, 0, 0) * t2s_2(2, 0)
235  + t4ddg_2(2, 1, 0, 0) * t2s_2(2, 1)
236  + t4ddg_2(2, 2, 0, 0) * t2s_2(2, 2)),
237  "T4ddg(i,j,k,l)*T2s(i,k)(0,0)");
238  test_for_zero(t2s_1(0, 1)
239  - (t4ddg_2(0, 0, 0, 1) * t2s_2(0, 0)
240  + t4ddg_2(0, 1, 0, 1) * t2s_2(0, 1)
241  + t4ddg_2(0, 2, 0, 1) * t2s_2(0, 2)
242  + t4ddg_2(1, 0, 0, 1) * t2s_2(1, 0)
243  + t4ddg_2(1, 1, 0, 1) * t2s_2(1, 1)
244  + t4ddg_2(1, 2, 0, 1) * t2s_2(1, 2)
245  + t4ddg_2(2, 0, 0, 1) * t2s_2(2, 0)
246  + t4ddg_2(2, 1, 0, 1) * t2s_2(2, 1)
247  + t4ddg_2(2, 2, 0, 1) * t2s_2(2, 2)),
248  "T4ddg(i,j,k,l)*T2s(i,k)(0,1)");
249  test_for_zero(t2s_1(0, 2)
250  - (t4ddg_2(0, 0, 0, 2) * t2s_2(0, 0)
251  + t4ddg_2(0, 1, 0, 2) * t2s_2(0, 1)
252  + t4ddg_2(0, 2, 0, 2) * t2s_2(0, 2)
253  + t4ddg_2(1, 0, 0, 2) * t2s_2(1, 0)
254  + t4ddg_2(1, 1, 0, 2) * t2s_2(1, 1)
255  + t4ddg_2(1, 2, 0, 2) * t2s_2(1, 2)
256  + t4ddg_2(2, 0, 0, 2) * t2s_2(2, 0)
257  + t4ddg_2(2, 1, 0, 2) * t2s_2(2, 1)
258  + t4ddg_2(2, 2, 0, 2) * t2s_2(2, 2)),
259  "T4ddg(i,j,k,l)*T2s(i,k)(0,2)");
260  test_for_zero(t2s_1(1, 0)
261  - (t4ddg_2(0, 0, 1, 0) * t2s_2(0, 0)
262  + t4ddg_2(0, 1, 1, 0) * t2s_2(0, 1)
263  + t4ddg_2(0, 2, 1, 0) * t2s_2(0, 2)
264  + t4ddg_2(1, 0, 1, 0) * t2s_2(1, 0)
265  + t4ddg_2(1, 1, 1, 0) * t2s_2(1, 1)
266  + t4ddg_2(1, 2, 1, 0) * t2s_2(1, 2)
267  + t4ddg_2(2, 0, 1, 0) * t2s_2(2, 0)
268  + t4ddg_2(2, 1, 1, 0) * t2s_2(2, 1)
269  + t4ddg_2(2, 2, 1, 0) * t2s_2(2, 2)),
270  "T4ddg(i,j,k,l)*T2s(i,k)(1,0)");
271  test_for_zero(t2s_1(1, 1)
272  - (t4ddg_2(0, 0, 1, 1) * t2s_2(0, 0)
273  + t4ddg_2(0, 1, 1, 1) * t2s_2(0, 1)
274  + t4ddg_2(0, 2, 1, 1) * t2s_2(0, 2)
275  + t4ddg_2(1, 0, 1, 1) * t2s_2(1, 0)
276  + t4ddg_2(1, 1, 1, 1) * t2s_2(1, 1)
277  + t4ddg_2(1, 2, 1, 1) * t2s_2(1, 2)
278  + t4ddg_2(2, 0, 1, 1) * t2s_2(2, 0)
279  + t4ddg_2(2, 1, 1, 1) * t2s_2(2, 1)
280  + t4ddg_2(2, 2, 1, 1) * t2s_2(2, 2)),
281  "T4ddg(i,j,k,l)*T2s(i,k)(1,1)");
282  test_for_zero(t2s_1(1, 2)
283  - (t4ddg_2(0, 0, 1, 2) * t2s_2(0, 0)
284  + t4ddg_2(0, 1, 1, 2) * t2s_2(0, 1)
285  + t4ddg_2(0, 2, 1, 2) * t2s_2(0, 2)
286  + t4ddg_2(1, 0, 1, 2) * t2s_2(1, 0)
287  + t4ddg_2(1, 1, 1, 2) * t2s_2(1, 1)
288  + t4ddg_2(1, 2, 1, 2) * t2s_2(1, 2)
289  + t4ddg_2(2, 0, 1, 2) * t2s_2(2, 0)
290  + t4ddg_2(2, 1, 1, 2) * t2s_2(2, 1)
291  + t4ddg_2(2, 2, 1, 2) * t2s_2(2, 2)),
292  "T4ddg(i,j,k,l)*T2s(i,k)(1,2)");
293  test_for_zero(t2s_1(2, 0)
294  - (t4ddg_2(0, 0, 2, 0) * t2s_2(0, 0)
295  + t4ddg_2(0, 1, 2, 0) * t2s_2(0, 1)
296  + t4ddg_2(0, 2, 2, 0) * t2s_2(0, 2)
297  + t4ddg_2(1, 0, 2, 0) * t2s_2(1, 0)
298  + t4ddg_2(1, 1, 2, 0) * t2s_2(1, 1)
299  + t4ddg_2(1, 2, 2, 0) * t2s_2(1, 2)
300  + t4ddg_2(2, 0, 2, 0) * t2s_2(2, 0)
301  + t4ddg_2(2, 1, 2, 0) * t2s_2(2, 1)
302  + t4ddg_2(2, 2, 2, 0) * t2s_2(2, 2)),
303  "T4ddg(i,j,k,l)*T2s(i,k)(2,0)");
304  test_for_zero(t2s_1(2, 1)
305  - (t4ddg_2(0, 0, 2, 1) * t2s_2(0, 0)
306  + t4ddg_2(0, 1, 2, 1) * t2s_2(0, 1)
307  + t4ddg_2(0, 2, 2, 1) * t2s_2(0, 2)
308  + t4ddg_2(1, 0, 2, 1) * t2s_2(1, 0)
309  + t4ddg_2(1, 1, 2, 1) * t2s_2(1, 1)
310  + t4ddg_2(1, 2, 2, 1) * t2s_2(1, 2)
311  + t4ddg_2(2, 0, 2, 1) * t2s_2(2, 0)
312  + t4ddg_2(2, 1, 2, 1) * t2s_2(2, 1)
313  + t4ddg_2(2, 2, 2, 1) * t2s_2(2, 2)),
314  "T4ddg(i,j,k,l)*T2s(i,k)(2,1)");
315  test_for_zero(t2s_1(2, 2)
316  - (t4ddg_2(0, 0, 2, 2) * t2s_2(0, 0)
317  + t4ddg_2(0, 1, 2, 2) * t2s_2(0, 1)
318  + t4ddg_2(0, 2, 2, 2) * t2s_2(0, 2)
319  + t4ddg_2(1, 0, 2, 2) * t2s_2(1, 0)
320  + t4ddg_2(1, 1, 2, 2) * t2s_2(1, 1)
321  + t4ddg_2(1, 2, 2, 2) * t2s_2(1, 2)
322  + t4ddg_2(2, 0, 2, 2) * t2s_2(2, 0)
323  + t4ddg_2(2, 1, 2, 2) * t2s_2(2, 1)
324  + t4ddg_2(2, 2, 2, 2) * t2s_2(2, 2)),
325  "T4ddg(i,j,k,l)*T2s(i,k)(2,2)");
326 
327  t2s_1(j, l) = t2s_3(i, k) * t4ddg_2(i, k, j, l);
328  test_for_zero(t2s_1(0, 0)
329  - (t4ddg_2(0, 0, 0, 0) * t2s_3(0, 0)
330  + t4ddg_2(0, 1, 0, 0) * t2s_3(0, 1)
331  + t4ddg_2(0, 2, 0, 0) * t2s_3(0, 2)
332  + t4ddg_2(1, 0, 0, 0) * t2s_3(1, 0)
333  + t4ddg_2(1, 1, 0, 0) * t2s_3(1, 1)
334  + t4ddg_2(1, 2, 0, 0) * t2s_3(1, 2)
335  + t4ddg_2(2, 0, 0, 0) * t2s_3(2, 0)
336  + t4ddg_2(2, 1, 0, 0) * t2s_3(2, 1)
337  + t4ddg_2(2, 2, 0, 0) * t2s_3(2, 2)),
338  "T2s(i,k)*T4ddg(i,j,k,l)(0,0)");
339  test_for_zero(t2s_1(0, 1)
340  - (t4ddg_2(0, 0, 0, 1) * t2s_3(0, 0)
341  + t4ddg_2(0, 1, 0, 1) * t2s_3(0, 1)
342  + t4ddg_2(0, 2, 0, 1) * t2s_3(0, 2)
343  + t4ddg_2(1, 0, 0, 1) * t2s_3(1, 0)
344  + t4ddg_2(1, 1, 0, 1) * t2s_3(1, 1)
345  + t4ddg_2(1, 2, 0, 1) * t2s_3(1, 2)
346  + t4ddg_2(2, 0, 0, 1) * t2s_3(2, 0)
347  + t4ddg_2(2, 1, 0, 1) * t2s_3(2, 1)
348  + t4ddg_2(2, 2, 0, 1) * t2s_3(2, 2)),
349  "T2s(i,k)*T4ddg(i,j,k,l)(0,1)");
350  test_for_zero(t2s_1(0, 2)
351  - (t4ddg_2(0, 0, 0, 2) * t2s_3(0, 0)
352  + t4ddg_2(0, 1, 0, 2) * t2s_3(0, 1)
353  + t4ddg_2(0, 2, 0, 2) * t2s_3(0, 2)
354  + t4ddg_2(1, 0, 0, 2) * t2s_3(1, 0)
355  + t4ddg_2(1, 1, 0, 2) * t2s_3(1, 1)
356  + t4ddg_2(1, 2, 0, 2) * t2s_3(1, 2)
357  + t4ddg_2(2, 0, 0, 2) * t2s_3(2, 0)
358  + t4ddg_2(2, 1, 0, 2) * t2s_3(2, 1)
359  + t4ddg_2(2, 2, 0, 2) * t2s_3(2, 2)),
360  "T2s(i,k)*T4ddg(i,j,k,l)(0,2)");
361  test_for_zero(t2s_1(1, 0)
362  - (t4ddg_2(0, 0, 1, 0) * t2s_3(0, 0)
363  + t4ddg_2(0, 1, 1, 0) * t2s_3(0, 1)
364  + t4ddg_2(0, 2, 1, 0) * t2s_3(0, 2)
365  + t4ddg_2(1, 0, 1, 0) * t2s_3(1, 0)
366  + t4ddg_2(1, 1, 1, 0) * t2s_3(1, 1)
367  + t4ddg_2(1, 2, 1, 0) * t2s_3(1, 2)
368  + t4ddg_2(2, 0, 1, 0) * t2s_3(2, 0)
369  + t4ddg_2(2, 1, 1, 0) * t2s_3(2, 1)
370  + t4ddg_2(2, 2, 1, 0) * t2s_3(2, 2)),
371  "T2s(i,k)*T4ddg(i,j,k,l)(1,0)");
372  test_for_zero(t2s_1(1, 1)
373  - (t4ddg_2(0, 0, 1, 1) * t2s_3(0, 0)
374  + t4ddg_2(0, 1, 1, 1) * t2s_3(0, 1)
375  + t4ddg_2(0, 2, 1, 1) * t2s_3(0, 2)
376  + t4ddg_2(1, 0, 1, 1) * t2s_3(1, 0)
377  + t4ddg_2(1, 1, 1, 1) * t2s_3(1, 1)
378  + t4ddg_2(1, 2, 1, 1) * t2s_3(1, 2)
379  + t4ddg_2(2, 0, 1, 1) * t2s_3(2, 0)
380  + t4ddg_2(2, 1, 1, 1) * t2s_3(2, 1)
381  + t4ddg_2(2, 2, 1, 1) * t2s_3(2, 2)),
382  "T2s(i,k)*T4ddg(i,j,k,l)(1,1)");
383  test_for_zero(t2s_1(1, 2)
384  - (t4ddg_2(0, 0, 1, 2) * t2s_3(0, 0)
385  + t4ddg_2(0, 1, 1, 2) * t2s_3(0, 1)
386  + t4ddg_2(0, 2, 1, 2) * t2s_3(0, 2)
387  + t4ddg_2(1, 0, 1, 2) * t2s_3(1, 0)
388  + t4ddg_2(1, 1, 1, 2) * t2s_3(1, 1)
389  + t4ddg_2(1, 2, 1, 2) * t2s_3(1, 2)
390  + t4ddg_2(2, 0, 1, 2) * t2s_3(2, 0)
391  + t4ddg_2(2, 1, 1, 2) * t2s_3(2, 1)
392  + t4ddg_2(2, 2, 1, 2) * t2s_3(2, 2)),
393  "T2s(i,k)*T4ddg(i,j,k,l)(1,2)");
394  test_for_zero(t2s_1(2, 0)
395  - (t4ddg_2(0, 0, 2, 0) * t2s_3(0, 0)
396  + t4ddg_2(0, 1, 2, 0) * t2s_3(0, 1)
397  + t4ddg_2(0, 2, 2, 0) * t2s_3(0, 2)
398  + t4ddg_2(1, 0, 2, 0) * t2s_3(1, 0)
399  + t4ddg_2(1, 1, 2, 0) * t2s_3(1, 1)
400  + t4ddg_2(1, 2, 2, 0) * t2s_3(1, 2)
401  + t4ddg_2(2, 0, 2, 0) * t2s_3(2, 0)
402  + t4ddg_2(2, 1, 2, 0) * t2s_3(2, 1)
403  + t4ddg_2(2, 2, 2, 0) * t2s_3(2, 2)),
404  "T2s(i,k)*T4ddg(i,j,k,l)(2,0)");
405  test_for_zero(t2s_1(2, 1)
406  - (t4ddg_2(0, 0, 2, 1) * t2s_3(0, 0)
407  + t4ddg_2(0, 1, 2, 1) * t2s_3(0, 1)
408  + t4ddg_2(0, 2, 2, 1) * t2s_3(0, 2)
409  + t4ddg_2(1, 0, 2, 1) * t2s_3(1, 0)
410  + t4ddg_2(1, 1, 2, 1) * t2s_3(1, 1)
411  + t4ddg_2(1, 2, 2, 1) * t2s_3(1, 2)
412  + t4ddg_2(2, 0, 2, 1) * t2s_3(2, 0)
413  + t4ddg_2(2, 1, 2, 1) * t2s_3(2, 1)
414  + t4ddg_2(2, 2, 2, 1) * t2s_3(2, 2)),
415  "T2s(i,k)*T4ddg(i,j,k,l)(2,1)");
416  test_for_zero(t2s_1(2, 2)
417  - (t4ddg_2(0, 0, 2, 2) * t2s_3(0, 0)
418  + t4ddg_2(0, 1, 2, 2) * t2s_3(0, 1)
419  + t4ddg_2(0, 2, 2, 2) * t2s_3(0, 2)
420  + t4ddg_2(1, 0, 2, 2) * t2s_3(1, 0)
421  + t4ddg_2(1, 1, 2, 2) * t2s_3(1, 1)
422  + t4ddg_2(1, 2, 2, 2) * t2s_3(1, 2)
423  + t4ddg_2(2, 0, 2, 2) * t2s_3(2, 0)
424  + t4ddg_2(2, 1, 2, 2) * t2s_3(2, 1)
425  + t4ddg_2(2, 2, 2, 2) * t2s_3(2, 2)),
426  "T2s(i,k)*T4ddg(i,j,k,l)(2,2)");
427 
428  t4ddg_1(i, j, k, l) = (t4ddg_2(i, j, k, l) & t2s_3(i, j));
429  test_for_zero(t4ddg_1(0, 0, 0, 0) - (t4ddg_2(0, 0, 0, 0) * t2s_3(0, 0)),
430  "T4ddg(i,j,k,l)&T2s(i,j)(0,0,0,0)");
431  test_for_zero(t4ddg_1(0, 0, 0, 1) - (t4ddg_2(0, 0, 0, 1) * t2s_3(0, 0)),
432  "T4ddg(i,j,k,l)&T2s(i,j)(0,0,0,1)");
433  test_for_zero(t4ddg_1(0, 0, 0, 2) - (t4ddg_2(0, 0, 0, 2) * t2s_3(0, 0)),
434  "T4ddg(i,j,k,l)&T2s(i,j)(0,0,0,2)");
435  test_for_zero(t4ddg_1(0, 0, 1, 0) - (t4ddg_2(0, 0, 1, 0) * t2s_3(0, 0)),
436  "T4ddg(i,j,k,l)&T2s(i,j)(0,0,1,0)");
437  test_for_zero(t4ddg_1(0, 0, 1, 1) - (t4ddg_2(0, 0, 1, 1) * t2s_3(0, 0)),
438  "T4ddg(i,j,k,l)&T2s(i,j)(0,0,1,1)");
439  test_for_zero(t4ddg_1(0, 0, 1, 2) - (t4ddg_2(0, 0, 1, 2) * t2s_3(0, 0)),
440  "T4ddg(i,j,k,l)&T2s(i,j)(0,0,1,2)");
441  test_for_zero(t4ddg_1(0, 0, 2, 0) - (t4ddg_2(0, 0, 2, 0) * t2s_3(0, 0)),
442  "T4ddg(i,j,k,l)&T2s(i,j)(0,0,2,0)");
443  test_for_zero(t4ddg_1(0, 0, 2, 1) - (t4ddg_2(0, 0, 2, 1) * t2s_3(0, 0)),
444  "T4ddg(i,j,k,l)&T2s(i,j)(0,0,2,1)");
445  test_for_zero(t4ddg_1(0, 0, 2, 2) - (t4ddg_2(0, 0, 2, 2) * t2s_3(0, 0)),
446  "T4ddg(i,j,k,l)&T2s(i,j)(0,0,2,2)");
447  test_for_zero(t4ddg_1(0, 1, 0, 0) - (t4ddg_2(0, 1, 0, 0) * t2s_3(0, 1)),
448  "T4ddg(i,j,k,l)&T2s(i,j)(0,1,0,0)");
449  test_for_zero(t4ddg_1(0, 1, 0, 1) - (t4ddg_2(0, 1, 0, 1) * t2s_3(0, 1)),
450  "T4ddg(i,j,k,l)&T2s(i,j)(0,1,0,1)");
451  test_for_zero(t4ddg_1(0, 1, 0, 2) - (t4ddg_2(0, 1, 0, 2) * t2s_3(0, 1)),
452  "T4ddg(i,j,k,l)&T2s(i,j)(0,1,0,2)");
453  test_for_zero(t4ddg_1(0, 1, 1, 0) - (t4ddg_2(0, 1, 1, 0) * t2s_3(0, 1)),
454  "T4ddg(i,j,k,l)&T2s(i,j)(0,1,1,0)");
455  test_for_zero(t4ddg_1(0, 1, 1, 1) - (t4ddg_2(0, 1, 1, 1) * t2s_3(0, 1)),
456  "T4ddg(i,j,k,l)&T2s(i,j)(0,1,1,1)");
457  test_for_zero(t4ddg_1(0, 1, 1, 2) - (t4ddg_2(0, 1, 1, 2) * t2s_3(0, 1)),
458  "T4ddg(i,j,k,l)&T2s(i,j)(0,1,1,2)");
459  test_for_zero(t4ddg_1(0, 1, 2, 0) - (t4ddg_2(0, 1, 2, 0) * t2s_3(0, 1)),
460  "T4ddg(i,j,k,l)&T2s(i,j)(0,1,2,0)");
461  test_for_zero(t4ddg_1(0, 1, 2, 1) - (t4ddg_2(0, 1, 2, 1) * t2s_3(0, 1)),
462  "T4ddg(i,j,k,l)&T2s(i,j)(0,1,2,1)");
463  test_for_zero(t4ddg_1(0, 1, 2, 2) - (t4ddg_2(0, 1, 2, 2) * t2s_3(0, 1)),
464  "T4ddg(i,j,k,l)&T2s(i,j)(0,1,2,2)");
465  test_for_zero(t4ddg_1(0, 2, 0, 0) - (t4ddg_2(0, 2, 0, 0) * t2s_3(0, 2)),
466  "T4ddg(i,j,k,l)&T2s(i,j)(0,2,0,0)");
467  test_for_zero(t4ddg_1(0, 2, 0, 1) - (t4ddg_2(0, 2, 0, 1) * t2s_3(0, 2)),
468  "T4ddg(i,j,k,l)&T2s(i,j)(0,2,0,1)");
469  test_for_zero(t4ddg_1(0, 2, 0, 2) - (t4ddg_2(0, 2, 0, 2) * t2s_3(0, 2)),
470  "T4ddg(i,j,k,l)&T2s(i,j)(0,2,0,2)");
471  test_for_zero(t4ddg_1(0, 2, 1, 0) - (t4ddg_2(0, 2, 1, 0) * t2s_3(0, 2)),
472  "T4ddg(i,j,k,l)&T2s(i,j)(0,2,1,0)");
473  test_for_zero(t4ddg_1(0, 2, 1, 1) - (t4ddg_2(0, 2, 1, 1) * t2s_3(0, 2)),
474  "T4ddg(i,j,k,l)&T2s(i,j)(0,2,1,1)");
475  test_for_zero(t4ddg_1(0, 2, 1, 2) - (t4ddg_2(0, 2, 1, 2) * t2s_3(0, 2)),
476  "T4ddg(i,j,k,l)&T2s(i,j)(0,2,1,2)");
477  test_for_zero(t4ddg_1(0, 2, 2, 0) - (t4ddg_2(0, 2, 2, 0) * t2s_3(0, 2)),
478  "T4ddg(i,j,k,l)&T2s(i,j)(0,2,2,0)");
479  test_for_zero(t4ddg_1(0, 2, 2, 1) - (t4ddg_2(0, 2, 2, 1) * t2s_3(0, 2)),
480  "T4ddg(i,j,k,l)&T2s(i,j)(0,2,2,1)");
481  test_for_zero(t4ddg_1(0, 2, 2, 2) - (t4ddg_2(0, 2, 2, 2) * t2s_3(0, 2)),
482  "T4ddg(i,j,k,l)&T2s(i,j)(0,2,2,2)");
483  test_for_zero(t4ddg_1(1, 0, 0, 0) - (t4ddg_2(1, 0, 0, 0) * t2s_3(1, 0)),
484  "T4ddg(i,j,k,l)&T2s(i,j)(1,0,0,0)");
485  test_for_zero(t4ddg_1(1, 0, 0, 1) - (t4ddg_2(1, 0, 0, 1) * t2s_3(1, 0)),
486  "T4ddg(i,j,k,l)&T2s(i,j)(1,0,0,1)");
487  test_for_zero(t4ddg_1(1, 0, 0, 2) - (t4ddg_2(1, 0, 0, 2) * t2s_3(1, 0)),
488  "T4ddg(i,j,k,l)&T2s(i,j)(1,0,0,2)");
489  test_for_zero(t4ddg_1(1, 0, 1, 0) - (t4ddg_2(1, 0, 1, 0) * t2s_3(1, 0)),
490  "T4ddg(i,j,k,l)&T2s(i,j)(1,0,1,0)");
491  test_for_zero(t4ddg_1(1, 0, 1, 1) - (t4ddg_2(1, 0, 1, 1) * t2s_3(1, 0)),
492  "T4ddg(i,j,k,l)&T2s(i,j)(1,0,1,1)");
493  test_for_zero(t4ddg_1(1, 0, 1, 2) - (t4ddg_2(1, 0, 1, 2) * t2s_3(1, 0)),
494  "T4ddg(i,j,k,l)&T2s(i,j)(1,0,1,2)");
495  test_for_zero(t4ddg_1(1, 0, 2, 0) - (t4ddg_2(1, 0, 2, 0) * t2s_3(1, 0)),
496  "T4ddg(i,j,k,l)&T2s(i,j)(1,0,2,0)");
497  test_for_zero(t4ddg_1(1, 0, 2, 1) - (t4ddg_2(1, 0, 2, 1) * t2s_3(1, 0)),
498  "T4ddg(i,j,k,l)&T2s(i,j)(1,0,2,1)");
499  test_for_zero(t4ddg_1(1, 0, 2, 2) - (t4ddg_2(1, 0, 2, 2) * t2s_3(1, 0)),
500  "T4ddg(i,j,k,l)&T2s(i,j)(1,0,2,2)");
501  test_for_zero(t4ddg_1(1, 1, 0, 0) - (t4ddg_2(1, 1, 0, 0) * t2s_3(1, 1)),
502  "T4ddg(i,j,k,l)&T2s(i,j)(1,1,0,0)");
503  test_for_zero(t4ddg_1(1, 1, 0, 1) - (t4ddg_2(1, 1, 0, 1) * t2s_3(1, 1)),
504  "T4ddg(i,j,k,l)&T2s(i,j)(1,1,0,1)");
505  test_for_zero(t4ddg_1(1, 1, 0, 2) - (t4ddg_2(1, 1, 0, 2) * t2s_3(1, 1)),
506  "T4ddg(i,j,k,l)&T2s(i,j)(1,1,0,2)");
507  test_for_zero(t4ddg_1(1, 1, 1, 0) - (t4ddg_2(1, 1, 1, 0) * t2s_3(1, 1)),
508  "T4ddg(i,j,k,l)&T2s(i,j)(1,1,1,0)");
509  test_for_zero(t4ddg_1(1, 1, 1, 1) - (t4ddg_2(1, 1, 1, 1) * t2s_3(1, 1)),
510  "T4ddg(i,j,k,l)&T2s(i,j)(1,1,1,1)");
511  test_for_zero(t4ddg_1(1, 1, 1, 2) - (t4ddg_2(1, 1, 1, 2) * t2s_3(1, 1)),
512  "T4ddg(i,j,k,l)&T2s(i,j)(1,1,1,2)");
513  test_for_zero(t4ddg_1(1, 1, 2, 0) - (t4ddg_2(1, 1, 2, 0) * t2s_3(1, 1)),
514  "T4ddg(i,j,k,l)&T2s(i,j)(1,1,2,0)");
515  test_for_zero(t4ddg_1(1, 1, 2, 1) - (t4ddg_2(1, 1, 2, 1) * t2s_3(1, 1)),
516  "T4ddg(i,j,k,l)&T2s(i,j)(1,1,2,1)");
517  test_for_zero(t4ddg_1(1, 1, 2, 2) - (t4ddg_2(1, 1, 2, 2) * t2s_3(1, 1)),
518  "T4ddg(i,j,k,l)&T2s(i,j)(1,1,2,2)");
519  test_for_zero(t4ddg_1(1, 2, 0, 0) - (t4ddg_2(1, 2, 0, 0) * t2s_3(1, 2)),
520  "T4ddg(i,j,k,l)&T2s(i,j)(1,2,0,0)");
521  test_for_zero(t4ddg_1(1, 2, 0, 1) - (t4ddg_2(1, 2, 0, 1) * t2s_3(1, 2)),
522  "T4ddg(i,j,k,l)&T2s(i,j)(1,2,0,1)");
523  test_for_zero(t4ddg_1(1, 2, 0, 2) - (t4ddg_2(1, 2, 0, 2) * t2s_3(1, 2)),
524  "T4ddg(i,j,k,l)&T2s(i,j)(1,2,0,2)");
525  test_for_zero(t4ddg_1(1, 2, 1, 0) - (t4ddg_2(1, 2, 1, 0) * t2s_3(1, 2)),
526  "T4ddg(i,j,k,l)&T2s(i,j)(1,2,1,0)");
527  test_for_zero(t4ddg_1(1, 2, 1, 1) - (t4ddg_2(1, 2, 1, 1) * t2s_3(1, 2)),
528  "T4ddg(i,j,k,l)&T2s(i,j)(1,2,1,1)");
529  test_for_zero(t4ddg_1(1, 2, 1, 2) - (t4ddg_2(1, 2, 1, 2) * t2s_3(1, 2)),
530  "T4ddg(i,j,k,l)&T2s(i,j)(1,2,1,2)");
531  test_for_zero(t4ddg_1(1, 2, 2, 0) - (t4ddg_2(1, 2, 2, 0) * t2s_3(1, 2)),
532  "T4ddg(i,j,k,l)&T2s(i,j)(1,2,2,0)");
533  test_for_zero(t4ddg_1(1, 2, 2, 1) - (t4ddg_2(1, 2, 2, 1) * t2s_3(1, 2)),
534  "T4ddg(i,j,k,l)&T2s(i,j)(1,2,2,1)");
535  test_for_zero(t4ddg_1(1, 2, 2, 2) - (t4ddg_2(1, 2, 2, 2) * t2s_3(1, 2)),
536  "T4ddg(i,j,k,l)&T2s(i,j)(1,2,2,2)");
537  test_for_zero(t4ddg_1(2, 0, 0, 0) - (t4ddg_2(2, 0, 0, 0) * t2s_3(2, 0)),
538  "T4ddg(i,j,k,l)&T2s(i,j)(2,0,0,0)");
539  test_for_zero(t4ddg_1(2, 0, 0, 1) - (t4ddg_2(2, 0, 0, 1) * t2s_3(2, 0)),
540  "T4ddg(i,j,k,l)&T2s(i,j)(2,0,0,1)");
541  test_for_zero(t4ddg_1(2, 0, 0, 2) - (t4ddg_2(2, 0, 0, 2) * t2s_3(2, 0)),
542  "T4ddg(i,j,k,l)&T2s(i,j)(2,0,0,2)");
543  test_for_zero(t4ddg_1(2, 0, 1, 0) - (t4ddg_2(2, 0, 1, 0) * t2s_3(2, 0)),
544  "T4ddg(i,j,k,l)&T2s(i,j)(2,0,1,0)");
545  test_for_zero(t4ddg_1(2, 0, 1, 1) - (t4ddg_2(2, 0, 1, 1) * t2s_3(2, 0)),
546  "T4ddg(i,j,k,l)&T2s(i,j)(2,0,1,1)");
547  test_for_zero(t4ddg_1(2, 0, 1, 2) - (t4ddg_2(2, 0, 1, 2) * t2s_3(2, 0)),
548  "T4ddg(i,j,k,l)&T2s(i,j)(2,0,1,2)");
549  test_for_zero(t4ddg_1(2, 0, 2, 0) - (t4ddg_2(2, 0, 2, 0) * t2s_3(2, 0)),
550  "T4ddg(i,j,k,l)&T2s(i,j)(2,0,2,0)");
551  test_for_zero(t4ddg_1(2, 0, 2, 1) - (t4ddg_2(2, 0, 2, 1) * t2s_3(2, 0)),
552  "T4ddg(i,j,k,l)&T2s(i,j)(2,0,2,1)");
553  test_for_zero(t4ddg_1(2, 0, 2, 2) - (t4ddg_2(2, 0, 2, 2) * t2s_3(2, 0)),
554  "T4ddg(i,j,k,l)&T2s(i,j)(2,0,2,2)");
555  test_for_zero(t4ddg_1(2, 1, 0, 0) - (t4ddg_2(2, 1, 0, 0) * t2s_3(2, 1)),
556  "T4ddg(i,j,k,l)&T2s(i,j)(2,1,0,0)");
557  test_for_zero(t4ddg_1(2, 1, 0, 1) - (t4ddg_2(2, 1, 0, 1) * t2s_3(2, 1)),
558  "T4ddg(i,j,k,l)&T2s(i,j)(2,1,0,1)");
559  test_for_zero(t4ddg_1(2, 1, 0, 2) - (t4ddg_2(2, 1, 0, 2) * t2s_3(2, 1)),
560  "T4ddg(i,j,k,l)&T2s(i,j)(2,1,0,2)");
561  test_for_zero(t4ddg_1(2, 1, 1, 0) - (t4ddg_2(2, 1, 1, 0) * t2s_3(2, 1)),
562  "T4ddg(i,j,k,l)&T2s(i,j)(2,1,1,0)");
563  test_for_zero(t4ddg_1(2, 1, 1, 1) - (t4ddg_2(2, 1, 1, 1) * t2s_3(2, 1)),
564  "T4ddg(i,j,k,l)&T2s(i,j)(2,1,1,1)");
565  test_for_zero(t4ddg_1(2, 1, 1, 2) - (t4ddg_2(2, 1, 1, 2) * t2s_3(2, 1)),
566  "T4ddg(i,j,k,l)&T2s(i,j)(2,1,1,2)");
567  test_for_zero(t4ddg_1(2, 1, 2, 0) - (t4ddg_2(2, 1, 2, 0) * t2s_3(2, 1)),
568  "T4ddg(i,j,k,l)&T2s(i,j)(2,1,2,0)");
569  test_for_zero(t4ddg_1(2, 1, 2, 1) - (t4ddg_2(2, 1, 2, 1) * t2s_3(2, 1)),
570  "T4ddg(i,j,k,l)&T2s(i,j)(2,1,2,1)");
571  test_for_zero(t4ddg_1(2, 1, 2, 2) - (t4ddg_2(2, 1, 2, 2) * t2s_3(2, 1)),
572  "T4ddg(i,j,k,l)&T2s(i,j)(2,1,2,2)");
573  test_for_zero(t4ddg_1(2, 2, 0, 0) - (t4ddg_2(2, 2, 0, 0) * t2s_3(2, 2)),
574  "T4ddg(i,j,k,l)&T2s(i,j)(2,2,0,0)");
575  test_for_zero(t4ddg_1(2, 2, 0, 1) - (t4ddg_2(2, 2, 0, 1) * t2s_3(2, 2)),
576  "T4ddg(i,j,k,l)&T2s(i,j)(2,2,0,1)");
577  test_for_zero(t4ddg_1(2, 2, 0, 2) - (t4ddg_2(2, 2, 0, 2) * t2s_3(2, 2)),
578  "T4ddg(i,j,k,l)&T2s(i,j)(2,2,0,2)");
579  test_for_zero(t4ddg_1(2, 2, 1, 0) - (t4ddg_2(2, 2, 1, 0) * t2s_3(2, 2)),
580  "T4ddg(i,j,k,l)&T2s(i,j)(2,2,1,0)");
581  test_for_zero(t4ddg_1(2, 2, 1, 1) - (t4ddg_2(2, 2, 1, 1) * t2s_3(2, 2)),
582  "T4ddg(i,j,k,l)&T2s(i,j)(2,2,1,1)");
583  test_for_zero(t4ddg_1(2, 2, 1, 2) - (t4ddg_2(2, 2, 1, 2) * t2s_3(2, 2)),
584  "T4ddg(i,j,k,l)&T2s(i,j)(2,2,1,2)");
585  test_for_zero(t4ddg_1(2, 2, 2, 0) - (t4ddg_2(2, 2, 2, 0) * t2s_3(2, 2)),
586  "T4ddg(i,j,k,l)&T2s(i,j)(2,2,2,0)");
587  test_for_zero(t4ddg_1(2, 2, 2, 1) - (t4ddg_2(2, 2, 2, 1) * t2s_3(2, 2)),
588  "T4ddg(i,j,k,l)&T2s(i,j)(2,2,2,1)");
589  test_for_zero(t4ddg_1(2, 2, 2, 2) - (t4ddg_2(2, 2, 2, 2) * t2s_3(2, 2)),
590  "T4ddg(i,j,k,l)&T2s(i,j)(2,2,2,2)");
591 
592  t4ddg_1(i, j, k, l) = (t2s_2(i, j) & t4ddg_2(i, j, k, l));
593  test_for_zero(t4ddg_1(0, 0, 0, 0) - (t4ddg_2(0, 0, 0, 0) * t2s_2(0, 0)),
594  "T2s(i,j)&T4ddg(i,j,k,l)(0,0,0,0)");
595  test_for_zero(t4ddg_1(0, 0, 0, 1) - (t4ddg_2(0, 0, 0, 1) * t2s_2(0, 0)),
596  "T2s(i,j)&T4ddg(i,j,k,l)(0,0,0,1)");
597  test_for_zero(t4ddg_1(0, 0, 0, 2) - (t4ddg_2(0, 0, 0, 2) * t2s_2(0, 0)),
598  "T2s(i,j)&T4ddg(i,j,k,l)(0,0,0,2)");
599  test_for_zero(t4ddg_1(0, 0, 1, 0) - (t4ddg_2(0, 0, 1, 0) * t2s_2(0, 0)),
600  "T2s(i,j)&T4ddg(i,j,k,l)(0,0,1,0)");
601  test_for_zero(t4ddg_1(0, 0, 1, 1) - (t4ddg_2(0, 0, 1, 1) * t2s_2(0, 0)),
602  "T2s(i,j)&T4ddg(i,j,k,l)(0,0,1,1)");
603  test_for_zero(t4ddg_1(0, 0, 1, 2) - (t4ddg_2(0, 0, 1, 2) * t2s_2(0, 0)),
604  "T2s(i,j)&T4ddg(i,j,k,l)(0,0,1,2)");
605  test_for_zero(t4ddg_1(0, 0, 2, 0) - (t4ddg_2(0, 0, 2, 0) * t2s_2(0, 0)),
606  "T2s(i,j)&T4ddg(i,j,k,l)(0,0,2,0)");
607  test_for_zero(t4ddg_1(0, 0, 2, 1) - (t4ddg_2(0, 0, 2, 1) * t2s_2(0, 0)),
608  "T2s(i,j)&T4ddg(i,j,k,l)(0,0,2,1)");
609  test_for_zero(t4ddg_1(0, 0, 2, 2) - (t4ddg_2(0, 0, 2, 2) * t2s_2(0, 0)),
610  "T2s(i,j)&T4ddg(i,j,k,l)(0,0,2,2)");
611  test_for_zero(t4ddg_1(0, 1, 0, 0) - (t4ddg_2(0, 1, 0, 0) * t2s_2(0, 1)),
612  "T2s(i,j)&T4ddg(i,j,k,l)(0,1,0,0)");
613  test_for_zero(t4ddg_1(0, 1, 0, 1) - (t4ddg_2(0, 1, 0, 1) * t2s_2(0, 1)),
614  "T2s(i,j)&T4ddg(i,j,k,l)(0,1,0,1)");
615  test_for_zero(t4ddg_1(0, 1, 0, 2) - (t4ddg_2(0, 1, 0, 2) * t2s_2(0, 1)),
616  "T2s(i,j)&T4ddg(i,j,k,l)(0,1,0,2)");
617  test_for_zero(t4ddg_1(0, 1, 1, 0) - (t4ddg_2(0, 1, 1, 0) * t2s_2(0, 1)),
618  "T2s(i,j)&T4ddg(i,j,k,l)(0,1,1,0)");
619  test_for_zero(t4ddg_1(0, 1, 1, 1) - (t4ddg_2(0, 1, 1, 1) * t2s_2(0, 1)),
620  "T2s(i,j)&T4ddg(i,j,k,l)(0,1,1,1)");
621  test_for_zero(t4ddg_1(0, 1, 1, 2) - (t4ddg_2(0, 1, 1, 2) * t2s_2(0, 1)),
622  "T2s(i,j)&T4ddg(i,j,k,l)(0,1,1,2)");
623  test_for_zero(t4ddg_1(0, 1, 2, 0) - (t4ddg_2(0, 1, 2, 0) * t2s_2(0, 1)),
624  "T2s(i,j)&T4ddg(i,j,k,l)(0,1,2,0)");
625  test_for_zero(t4ddg_1(0, 1, 2, 1) - (t4ddg_2(0, 1, 2, 1) * t2s_2(0, 1)),
626  "T2s(i,j)&T4ddg(i,j,k,l)(0,1,2,1)");
627  test_for_zero(t4ddg_1(0, 1, 2, 2) - (t4ddg_2(0, 1, 2, 2) * t2s_2(0, 1)),
628  "T2s(i,j)&T4ddg(i,j,k,l)(0,1,2,2)");
629  test_for_zero(t4ddg_1(0, 2, 0, 0) - (t4ddg_2(0, 2, 0, 0) * t2s_2(0, 2)),
630  "T2s(i,j)&T4ddg(i,j,k,l)(0,2,0,0)");
631  test_for_zero(t4ddg_1(0, 2, 0, 1) - (t4ddg_2(0, 2, 0, 1) * t2s_2(0, 2)),
632  "T2s(i,j)&T4ddg(i,j,k,l)(0,2,0,1)");
633  test_for_zero(t4ddg_1(0, 2, 0, 2) - (t4ddg_2(0, 2, 0, 2) * t2s_2(0, 2)),
634  "T2s(i,j)&T4ddg(i,j,k,l)(0,2,0,2)");
635  test_for_zero(t4ddg_1(0, 2, 1, 0) - (t4ddg_2(0, 2, 1, 0) * t2s_2(0, 2)),
636  "T2s(i,j)&T4ddg(i,j,k,l)(0,2,1,0)");
637  test_for_zero(t4ddg_1(0, 2, 1, 1) - (t4ddg_2(0, 2, 1, 1) * t2s_2(0, 2)),
638  "T2s(i,j)&T4ddg(i,j,k,l)(0,2,1,1)");
639  test_for_zero(t4ddg_1(0, 2, 1, 2) - (t4ddg_2(0, 2, 1, 2) * t2s_2(0, 2)),
640  "T2s(i,j)&T4ddg(i,j,k,l)(0,2,1,2)");
641  test_for_zero(t4ddg_1(0, 2, 2, 0) - (t4ddg_2(0, 2, 2, 0) * t2s_2(0, 2)),
642  "T2s(i,j)&T4ddg(i,j,k,l)(0,2,2,0)");
643  test_for_zero(t4ddg_1(0, 2, 2, 1) - (t4ddg_2(0, 2, 2, 1) * t2s_2(0, 2)),
644  "T2s(i,j)&T4ddg(i,j,k,l)(0,2,2,1)");
645  test_for_zero(t4ddg_1(0, 2, 2, 2) - (t4ddg_2(0, 2, 2, 2) * t2s_2(0, 2)),
646  "T2s(i,j)&T4ddg(i,j,k,l)(0,2,2,2)");
647  test_for_zero(t4ddg_1(1, 0, 0, 0) - (t4ddg_2(1, 0, 0, 0) * t2s_2(1, 0)),
648  "T2s(i,j)&T4ddg(i,j,k,l)(1,0,0,0)");
649  test_for_zero(t4ddg_1(1, 0, 0, 1) - (t4ddg_2(1, 0, 0, 1) * t2s_2(1, 0)),
650  "T2s(i,j)&T4ddg(i,j,k,l)(1,0,0,1)");
651  test_for_zero(t4ddg_1(1, 0, 0, 2) - (t4ddg_2(1, 0, 0, 2) * t2s_2(1, 0)),
652  "T2s(i,j)&T4ddg(i,j,k,l)(1,0,0,2)");
653  test_for_zero(t4ddg_1(1, 0, 1, 0) - (t4ddg_2(1, 0, 1, 0) * t2s_2(1, 0)),
654  "T2s(i,j)&T4ddg(i,j,k,l)(1,0,1,0)");
655  test_for_zero(t4ddg_1(1, 0, 1, 1) - (t4ddg_2(1, 0, 1, 1) * t2s_2(1, 0)),
656  "T2s(i,j)&T4ddg(i,j,k,l)(1,0,1,1)");
657  test_for_zero(t4ddg_1(1, 0, 1, 2) - (t4ddg_2(1, 0, 1, 2) * t2s_2(1, 0)),
658  "T2s(i,j)&T4ddg(i,j,k,l)(1,0,1,2)");
659  test_for_zero(t4ddg_1(1, 0, 2, 0) - (t4ddg_2(1, 0, 2, 0) * t2s_2(1, 0)),
660  "T2s(i,j)&T4ddg(i,j,k,l)(1,0,2,0)");
661  test_for_zero(t4ddg_1(1, 0, 2, 1) - (t4ddg_2(1, 0, 2, 1) * t2s_2(1, 0)),
662  "T2s(i,j)&T4ddg(i,j,k,l)(1,0,2,1)");
663  test_for_zero(t4ddg_1(1, 0, 2, 2) - (t4ddg_2(1, 0, 2, 2) * t2s_2(1, 0)),
664  "T2s(i,j)&T4ddg(i,j,k,l)(1,0,2,2)");
665  test_for_zero(t4ddg_1(1, 1, 0, 0) - (t4ddg_2(1, 1, 0, 0) * t2s_2(1, 1)),
666  "T2s(i,j)&T4ddg(i,j,k,l)(1,1,0,0)");
667  test_for_zero(t4ddg_1(1, 1, 0, 1) - (t4ddg_2(1, 1, 0, 1) * t2s_2(1, 1)),
668  "T2s(i,j)&T4ddg(i,j,k,l)(1,1,0,1)");
669  test_for_zero(t4ddg_1(1, 1, 0, 2) - (t4ddg_2(1, 1, 0, 2) * t2s_2(1, 1)),
670  "T2s(i,j)&T4ddg(i,j,k,l)(1,1,0,2)");
671  test_for_zero(t4ddg_1(1, 1, 1, 0) - (t4ddg_2(1, 1, 1, 0) * t2s_2(1, 1)),
672  "T2s(i,j)&T4ddg(i,j,k,l)(1,1,1,0)");
673  test_for_zero(t4ddg_1(1, 1, 1, 1) - (t4ddg_2(1, 1, 1, 1) * t2s_2(1, 1)),
674  "T2s(i,j)&T4ddg(i,j,k,l)(1,1,1,1)");
675  test_for_zero(t4ddg_1(1, 1, 1, 2) - (t4ddg_2(1, 1, 1, 2) * t2s_2(1, 1)),
676  "T2s(i,j)&T4ddg(i,j,k,l)(1,1,1,2)");
677  test_for_zero(t4ddg_1(1, 1, 2, 0) - (t4ddg_2(1, 1, 2, 0) * t2s_2(1, 1)),
678  "T2s(i,j)&T4ddg(i,j,k,l)(1,1,2,0)");
679  test_for_zero(t4ddg_1(1, 1, 2, 1) - (t4ddg_2(1, 1, 2, 1) * t2s_2(1, 1)),
680  "T2s(i,j)&T4ddg(i,j,k,l)(1,1,2,1)");
681  test_for_zero(t4ddg_1(1, 1, 2, 2) - (t4ddg_2(1, 1, 2, 2) * t2s_2(1, 1)),
682  "T2s(i,j)&T4ddg(i,j,k,l)(1,1,2,2)");
683  test_for_zero(t4ddg_1(1, 2, 0, 0) - (t4ddg_2(1, 2, 0, 0) * t2s_2(1, 2)),
684  "T2s(i,j)&T4ddg(i,j,k,l)(1,2,0,0)");
685  test_for_zero(t4ddg_1(1, 2, 0, 1) - (t4ddg_2(1, 2, 0, 1) * t2s_2(1, 2)),
686  "T2s(i,j)&T4ddg(i,j,k,l)(1,2,0,1)");
687  test_for_zero(t4ddg_1(1, 2, 0, 2) - (t4ddg_2(1, 2, 0, 2) * t2s_2(1, 2)),
688  "T2s(i,j)&T4ddg(i,j,k,l)(1,2,0,2)");
689  test_for_zero(t4ddg_1(1, 2, 1, 0) - (t4ddg_2(1, 2, 1, 0) * t2s_2(1, 2)),
690  "T2s(i,j)&T4ddg(i,j,k,l)(1,2,1,0)");
691  test_for_zero(t4ddg_1(1, 2, 1, 1) - (t4ddg_2(1, 2, 1, 1) * t2s_2(1, 2)),
692  "T2s(i,j)&T4ddg(i,j,k,l)(1,2,1,1)");
693  test_for_zero(t4ddg_1(1, 2, 1, 2) - (t4ddg_2(1, 2, 1, 2) * t2s_2(1, 2)),
694  "T2s(i,j)&T4ddg(i,j,k,l)(1,2,1,2)");
695  test_for_zero(t4ddg_1(1, 2, 2, 0) - (t4ddg_2(1, 2, 2, 0) * t2s_2(1, 2)),
696  "T2s(i,j)&T4ddg(i,j,k,l)(1,2,2,0)");
697  test_for_zero(t4ddg_1(1, 2, 2, 1) - (t4ddg_2(1, 2, 2, 1) * t2s_2(1, 2)),
698  "T2s(i,j)&T4ddg(i,j,k,l)(1,2,2,1)");
699  test_for_zero(t4ddg_1(1, 2, 2, 2) - (t4ddg_2(1, 2, 2, 2) * t2s_2(1, 2)),
700  "T2s(i,j)&T4ddg(i,j,k,l)(1,2,2,2)");
701  test_for_zero(t4ddg_1(2, 0, 0, 0) - (t4ddg_2(2, 0, 0, 0) * t2s_2(2, 0)),
702  "T2s(i,j)&T4ddg(i,j,k,l)(2,0,0,0)");
703  test_for_zero(t4ddg_1(2, 0, 0, 1) - (t4ddg_2(2, 0, 0, 1) * t2s_2(2, 0)),
704  "T2s(i,j)&T4ddg(i,j,k,l)(2,0,0,1)");
705  test_for_zero(t4ddg_1(2, 0, 0, 2) - (t4ddg_2(2, 0, 0, 2) * t2s_2(2, 0)),
706  "T2s(i,j)&T4ddg(i,j,k,l)(2,0,0,2)");
707  test_for_zero(t4ddg_1(2, 0, 1, 0) - (t4ddg_2(2, 0, 1, 0) * t2s_2(2, 0)),
708  "T2s(i,j)&T4ddg(i,j,k,l)(2,0,1,0)");
709  test_for_zero(t4ddg_1(2, 0, 1, 1) - (t4ddg_2(2, 0, 1, 1) * t2s_2(2, 0)),
710  "T2s(i,j)&T4ddg(i,j,k,l)(2,0,1,1)");
711  test_for_zero(t4ddg_1(2, 0, 1, 2) - (t4ddg_2(2, 0, 1, 2) * t2s_2(2, 0)),
712  "T2s(i,j)&T4ddg(i,j,k,l)(2,0,1,2)");
713  test_for_zero(t4ddg_1(2, 0, 2, 0) - (t4ddg_2(2, 0, 2, 0) * t2s_2(2, 0)),
714  "T2s(i,j)&T4ddg(i,j,k,l)(2,0,2,0)");
715  test_for_zero(t4ddg_1(2, 0, 2, 1) - (t4ddg_2(2, 0, 2, 1) * t2s_2(2, 0)),
716  "T2s(i,j)&T4ddg(i,j,k,l)(2,0,2,1)");
717  test_for_zero(t4ddg_1(2, 0, 2, 2) - (t4ddg_2(2, 0, 2, 2) * t2s_2(2, 0)),
718  "T2s(i,j)&T4ddg(i,j,k,l)(2,0,2,2)");
719  test_for_zero(t4ddg_1(2, 1, 0, 0) - (t4ddg_2(2, 1, 0, 0) * t2s_2(2, 1)),
720  "T2s(i,j)&T4ddg(i,j,k,l)(2,1,0,0)");
721  test_for_zero(t4ddg_1(2, 1, 0, 1) - (t4ddg_2(2, 1, 0, 1) * t2s_2(2, 1)),
722  "T2s(i,j)&T4ddg(i,j,k,l)(2,1,0,1)");
723  test_for_zero(t4ddg_1(2, 1, 0, 2) - (t4ddg_2(2, 1, 0, 2) * t2s_2(2, 1)),
724  "T2s(i,j)&T4ddg(i,j,k,l)(2,1,0,2)");
725  test_for_zero(t4ddg_1(2, 1, 1, 0) - (t4ddg_2(2, 1, 1, 0) * t2s_2(2, 1)),
726  "T2s(i,j)&T4ddg(i,j,k,l)(2,1,1,0)");
727  test_for_zero(t4ddg_1(2, 1, 1, 1) - (t4ddg_2(2, 1, 1, 1) * t2s_2(2, 1)),
728  "T2s(i,j)&T4ddg(i,j,k,l)(2,1,1,1)");
729  test_for_zero(t4ddg_1(2, 1, 1, 2) - (t4ddg_2(2, 1, 1, 2) * t2s_2(2, 1)),
730  "T2s(i,j)&T4ddg(i,j,k,l)(2,1,1,2)");
731  test_for_zero(t4ddg_1(2, 1, 2, 0) - (t4ddg_2(2, 1, 2, 0) * t2s_2(2, 1)),
732  "T2s(i,j)&T4ddg(i,j,k,l)(2,1,2,0)");
733  test_for_zero(t4ddg_1(2, 1, 2, 1) - (t4ddg_2(2, 1, 2, 1) * t2s_2(2, 1)),
734  "T2s(i,j)&T4ddg(i,j,k,l)(2,1,2,1)");
735  test_for_zero(t4ddg_1(2, 1, 2, 2) - (t4ddg_2(2, 1, 2, 2) * t2s_2(2, 1)),
736  "T2s(i,j)&T4ddg(i,j,k,l)(2,1,2,2)");
737  test_for_zero(t4ddg_1(2, 2, 0, 0) - (t4ddg_2(2, 2, 0, 0) * t2s_2(2, 2)),
738  "T2s(i,j)&T4ddg(i,j,k,l)(2,2,0,0)");
739  test_for_zero(t4ddg_1(2, 2, 0, 1) - (t4ddg_2(2, 2, 0, 1) * t2s_2(2, 2)),
740  "T2s(i,j)&T4ddg(i,j,k,l)(2,2,0,1)");
741  test_for_zero(t4ddg_1(2, 2, 0, 2) - (t4ddg_2(2, 2, 0, 2) * t2s_2(2, 2)),
742  "T2s(i,j)&T4ddg(i,j,k,l)(2,2,0,2)");
743  test_for_zero(t4ddg_1(2, 2, 1, 0) - (t4ddg_2(2, 2, 1, 0) * t2s_2(2, 2)),
744  "T2s(i,j)&T4ddg(i,j,k,l)(2,2,1,0)");
745  test_for_zero(t4ddg_1(2, 2, 1, 1) - (t4ddg_2(2, 2, 1, 1) * t2s_2(2, 2)),
746  "T2s(i,j)&T4ddg(i,j,k,l)(2,2,1,1)");
747  test_for_zero(t4ddg_1(2, 2, 1, 2) - (t4ddg_2(2, 2, 1, 2) * t2s_2(2, 2)),
748  "T2s(i,j)&T4ddg(i,j,k,l)(2,2,1,2)");
749  test_for_zero(t4ddg_1(2, 2, 2, 0) - (t4ddg_2(2, 2, 2, 0) * t2s_2(2, 2)),
750  "T2s(i,j)&T4ddg(i,j,k,l)(2,2,2,0)");
751  test_for_zero(t4ddg_1(2, 2, 2, 1) - (t4ddg_2(2, 2, 2, 1) * t2s_2(2, 2)),
752  "T2s(i,j)&T4ddg(i,j,k,l)(2,2,2,1)");
753  test_for_zero(t4ddg_1(2, 2, 2, 2) - (t4ddg_2(2, 2, 2, 2) * t2s_2(2, 2)),
754  "T2s(i,j)&T4ddg(i,j,k,l)(2,2,2,2)");
755 
756  t4ddg_1(i, j, k, l) = (t4ddg_2(i, j, k, l) & t2s_3(k, l));
757  test_for_zero(t4ddg_1(0, 0, 0, 0) - (t4ddg_2(0, 0, 0, 0) * t2s_3(0, 0)),
758  "T4ddg(i,j,k,l)&T2s(k,l)(0,0,0,0)");
759  test_for_zero(t4ddg_1(0, 0, 0, 1) - (t4ddg_2(0, 0, 0, 1) * t2s_3(0, 1)),
760  "T4ddg(i,j,k,l)&T2s(k,l)(0,0,0,1)");
761  test_for_zero(t4ddg_1(0, 0, 0, 2) - (t4ddg_2(0, 0, 0, 2) * t2s_3(0, 2)),
762  "T4ddg(i,j,k,l)&T2s(k,l)(0,0,0,2)");
763  test_for_zero(t4ddg_1(0, 0, 1, 0) - (t4ddg_2(0, 0, 1, 0) * t2s_3(1, 0)),
764  "T4ddg(i,j,k,l)&T2s(k,l)(0,0,1,0)");
765  test_for_zero(t4ddg_1(0, 0, 1, 1) - (t4ddg_2(0, 0, 1, 1) * t2s_3(1, 1)),
766  "T4ddg(i,j,k,l)&T2s(k,l)(0,0,1,1)");
767  test_for_zero(t4ddg_1(0, 0, 1, 2) - (t4ddg_2(0, 0, 1, 2) * t2s_3(1, 2)),
768  "T4ddg(i,j,k,l)&T2s(k,l)(0,0,1,2)");
769  test_for_zero(t4ddg_1(0, 0, 2, 0) - (t4ddg_2(0, 0, 2, 0) * t2s_3(2, 0)),
770  "T4ddg(i,j,k,l)&T2s(k,l)(0,0,2,0)");
771  test_for_zero(t4ddg_1(0, 0, 2, 1) - (t4ddg_2(0, 0, 2, 1) * t2s_3(2, 1)),
772  "T4ddg(i,j,k,l)&T2s(k,l)(0,0,2,1)");
773  test_for_zero(t4ddg_1(0, 0, 2, 2) - (t4ddg_2(0, 0, 2, 2) * t2s_3(2, 2)),
774  "T4ddg(i,j,k,l)&T2s(k,l)(0,0,2,2)");
775  test_for_zero(t4ddg_1(0, 1, 0, 0) - (t4ddg_2(0, 1, 0, 0) * t2s_3(0, 0)),
776  "T4ddg(i,j,k,l)&T2s(k,l)(0,1,0,0)");
777  test_for_zero(t4ddg_1(0, 1, 0, 1) - (t4ddg_2(0, 1, 0, 1) * t2s_3(0, 1)),
778  "T4ddg(i,j,k,l)&T2s(k,l)(0,1,0,1)");
779  test_for_zero(t4ddg_1(0, 1, 0, 2) - (t4ddg_2(0, 1, 0, 2) * t2s_3(0, 2)),
780  "T4ddg(i,j,k,l)&T2s(k,l)(0,1,0,2)");
781  test_for_zero(t4ddg_1(0, 1, 1, 0) - (t4ddg_2(0, 1, 1, 0) * t2s_3(1, 0)),
782  "T4ddg(i,j,k,l)&T2s(k,l)(0,1,1,0)");
783  test_for_zero(t4ddg_1(0, 1, 1, 1) - (t4ddg_2(0, 1, 1, 1) * t2s_3(1, 1)),
784  "T4ddg(i,j,k,l)&T2s(k,l)(0,1,1,1)");
785  test_for_zero(t4ddg_1(0, 1, 1, 2) - (t4ddg_2(0, 1, 1, 2) * t2s_3(1, 2)),
786  "T4ddg(i,j,k,l)&T2s(k,l)(0,1,1,2)");
787  test_for_zero(t4ddg_1(0, 1, 2, 0) - (t4ddg_2(0, 1, 2, 0) * t2s_3(2, 0)),
788  "T4ddg(i,j,k,l)&T2s(k,l)(0,1,2,0)");
789  test_for_zero(t4ddg_1(0, 1, 2, 1) - (t4ddg_2(0, 1, 2, 1) * t2s_3(2, 1)),
790  "T4ddg(i,j,k,l)&T2s(k,l)(0,1,2,1)");
791  test_for_zero(t4ddg_1(0, 1, 2, 2) - (t4ddg_2(0, 1, 2, 2) * t2s_3(2, 2)),
792  "T4ddg(i,j,k,l)&T2s(k,l)(0,1,2,2)");
793  test_for_zero(t4ddg_1(0, 2, 0, 0) - (t4ddg_2(0, 2, 0, 0) * t2s_3(0, 0)),
794  "T4ddg(i,j,k,l)&T2s(k,l)(0,2,0,0)");
795  test_for_zero(t4ddg_1(0, 2, 0, 1) - (t4ddg_2(0, 2, 0, 1) * t2s_3(0, 1)),
796  "T4ddg(i,j,k,l)&T2s(k,l)(0,2,0,1)");
797  test_for_zero(t4ddg_1(0, 2, 0, 2) - (t4ddg_2(0, 2, 0, 2) * t2s_3(0, 2)),
798  "T4ddg(i,j,k,l)&T2s(k,l)(0,2,0,2)");
799  test_for_zero(t4ddg_1(0, 2, 1, 0) - (t4ddg_2(0, 2, 1, 0) * t2s_3(1, 0)),
800  "T4ddg(i,j,k,l)&T2s(k,l)(0,2,1,0)");
801  test_for_zero(t4ddg_1(0, 2, 1, 1) - (t4ddg_2(0, 2, 1, 1) * t2s_3(1, 1)),
802  "T4ddg(i,j,k,l)&T2s(k,l)(0,2,1,1)");
803  test_for_zero(t4ddg_1(0, 2, 1, 2) - (t4ddg_2(0, 2, 1, 2) * t2s_3(1, 2)),
804  "T4ddg(i,j,k,l)&T2s(k,l)(0,2,1,2)");
805  test_for_zero(t4ddg_1(0, 2, 2, 0) - (t4ddg_2(0, 2, 2, 0) * t2s_3(2, 0)),
806  "T4ddg(i,j,k,l)&T2s(k,l)(0,2,2,0)");
807  test_for_zero(t4ddg_1(0, 2, 2, 1) - (t4ddg_2(0, 2, 2, 1) * t2s_3(2, 1)),
808  "T4ddg(i,j,k,l)&T2s(k,l)(0,2,2,1)");
809  test_for_zero(t4ddg_1(0, 2, 2, 2) - (t4ddg_2(0, 2, 2, 2) * t2s_3(2, 2)),
810  "T4ddg(i,j,k,l)&T2s(k,l)(0,2,2,2)");
811  test_for_zero(t4ddg_1(1, 0, 0, 0) - (t4ddg_2(1, 0, 0, 0) * t2s_3(0, 0)),
812  "T4ddg(i,j,k,l)&T2s(k,l)(1,0,0,0)");
813  test_for_zero(t4ddg_1(1, 0, 0, 1) - (t4ddg_2(1, 0, 0, 1) * t2s_3(0, 1)),
814  "T4ddg(i,j,k,l)&T2s(k,l)(1,0,0,1)");
815  test_for_zero(t4ddg_1(1, 0, 0, 2) - (t4ddg_2(1, 0, 0, 2) * t2s_3(0, 2)),
816  "T4ddg(i,j,k,l)&T2s(k,l)(1,0,0,2)");
817  test_for_zero(t4ddg_1(1, 0, 1, 0) - (t4ddg_2(1, 0, 1, 0) * t2s_3(1, 0)),
818  "T4ddg(i,j,k,l)&T2s(k,l)(1,0,1,0)");
819  test_for_zero(t4ddg_1(1, 0, 1, 1) - (t4ddg_2(1, 0, 1, 1) * t2s_3(1, 1)),
820  "T4ddg(i,j,k,l)&T2s(k,l)(1,0,1,1)");
821  test_for_zero(t4ddg_1(1, 0, 1, 2) - (t4ddg_2(1, 0, 1, 2) * t2s_3(1, 2)),
822  "T4ddg(i,j,k,l)&T2s(k,l)(1,0,1,2)");
823  test_for_zero(t4ddg_1(1, 0, 2, 0) - (t4ddg_2(1, 0, 2, 0) * t2s_3(2, 0)),
824  "T4ddg(i,j,k,l)&T2s(k,l)(1,0,2,0)");
825  test_for_zero(t4ddg_1(1, 0, 2, 1) - (t4ddg_2(1, 0, 2, 1) * t2s_3(2, 1)),
826  "T4ddg(i,j,k,l)&T2s(k,l)(1,0,2,1)");
827  test_for_zero(t4ddg_1(1, 0, 2, 2) - (t4ddg_2(1, 0, 2, 2) * t2s_3(2, 2)),
828  "T4ddg(i,j,k,l)&T2s(k,l)(1,0,2,2)");
829  test_for_zero(t4ddg_1(1, 1, 0, 0) - (t4ddg_2(1, 1, 0, 0) * t2s_3(0, 0)),
830  "T4ddg(i,j,k,l)&T2s(k,l)(1,1,0,0)");
831  test_for_zero(t4ddg_1(1, 1, 0, 1) - (t4ddg_2(1, 1, 0, 1) * t2s_3(0, 1)),
832  "T4ddg(i,j,k,l)&T2s(k,l)(1,1,0,1)");
833  test_for_zero(t4ddg_1(1, 1, 0, 2) - (t4ddg_2(1, 1, 0, 2) * t2s_3(0, 2)),
834  "T4ddg(i,j,k,l)&T2s(k,l)(1,1,0,2)");
835  test_for_zero(t4ddg_1(1, 1, 1, 0) - (t4ddg_2(1, 1, 1, 0) * t2s_3(1, 0)),
836  "T4ddg(i,j,k,l)&T2s(k,l)(1,1,1,0)");
837  test_for_zero(t4ddg_1(1, 1, 1, 1) - (t4ddg_2(1, 1, 1, 1) * t2s_3(1, 1)),
838  "T4ddg(i,j,k,l)&T2s(k,l)(1,1,1,1)");
839  test_for_zero(t4ddg_1(1, 1, 1, 2) - (t4ddg_2(1, 1, 1, 2) * t2s_3(1, 2)),
840  "T4ddg(i,j,k,l)&T2s(k,l)(1,1,1,2)");
841  test_for_zero(t4ddg_1(1, 1, 2, 0) - (t4ddg_2(1, 1, 2, 0) * t2s_3(2, 0)),
842  "T4ddg(i,j,k,l)&T2s(k,l)(1,1,2,0)");
843  test_for_zero(t4ddg_1(1, 1, 2, 1) - (t4ddg_2(1, 1, 2, 1) * t2s_3(2, 1)),
844  "T4ddg(i,j,k,l)&T2s(k,l)(1,1,2,1)");
845  test_for_zero(t4ddg_1(1, 1, 2, 2) - (t4ddg_2(1, 1, 2, 2) * t2s_3(2, 2)),
846  "T4ddg(i,j,k,l)&T2s(k,l)(1,1,2,2)");
847  test_for_zero(t4ddg_1(1, 2, 0, 0) - (t4ddg_2(1, 2, 0, 0) * t2s_3(0, 0)),
848  "T4ddg(i,j,k,l)&T2s(k,l)(1,2,0,0)");
849  test_for_zero(t4ddg_1(1, 2, 0, 1) - (t4ddg_2(1, 2, 0, 1) * t2s_3(0, 1)),
850  "T4ddg(i,j,k,l)&T2s(k,l)(1,2,0,1)");
851  test_for_zero(t4ddg_1(1, 2, 0, 2) - (t4ddg_2(1, 2, 0, 2) * t2s_3(0, 2)),
852  "T4ddg(i,j,k,l)&T2s(k,l)(1,2,0,2)");
853  test_for_zero(t4ddg_1(1, 2, 1, 0) - (t4ddg_2(1, 2, 1, 0) * t2s_3(1, 0)),
854  "T4ddg(i,j,k,l)&T2s(k,l)(1,2,1,0)");
855  test_for_zero(t4ddg_1(1, 2, 1, 1) - (t4ddg_2(1, 2, 1, 1) * t2s_3(1, 1)),
856  "T4ddg(i,j,k,l)&T2s(k,l)(1,2,1,1)");
857  test_for_zero(t4ddg_1(1, 2, 1, 2) - (t4ddg_2(1, 2, 1, 2) * t2s_3(1, 2)),
858  "T4ddg(i,j,k,l)&T2s(k,l)(1,2,1,2)");
859  test_for_zero(t4ddg_1(1, 2, 2, 0) - (t4ddg_2(1, 2, 2, 0) * t2s_3(2, 0)),
860  "T4ddg(i,j,k,l)&T2s(k,l)(1,2,2,0)");
861  test_for_zero(t4ddg_1(1, 2, 2, 1) - (t4ddg_2(1, 2, 2, 1) * t2s_3(2, 1)),
862  "T4ddg(i,j,k,l)&T2s(k,l)(1,2,2,1)");
863  test_for_zero(t4ddg_1(1, 2, 2, 2) - (t4ddg_2(1, 2, 2, 2) * t2s_3(2, 2)),
864  "T4ddg(i,j,k,l)&T2s(k,l)(1,2,2,2)");
865  test_for_zero(t4ddg_1(2, 0, 0, 0) - (t4ddg_2(2, 0, 0, 0) * t2s_3(0, 0)),
866  "T4ddg(i,j,k,l)&T2s(k,l)(2,0,0,0)");
867  test_for_zero(t4ddg_1(2, 0, 0, 1) - (t4ddg_2(2, 0, 0, 1) * t2s_3(0, 1)),
868  "T4ddg(i,j,k,l)&T2s(k,l)(2,0,0,1)");
869  test_for_zero(t4ddg_1(2, 0, 0, 2) - (t4ddg_2(2, 0, 0, 2) * t2s_3(0, 2)),
870  "T4ddg(i,j,k,l)&T2s(k,l)(2,0,0,2)");
871  test_for_zero(t4ddg_1(2, 0, 1, 0) - (t4ddg_2(2, 0, 1, 0) * t2s_3(1, 0)),
872  "T4ddg(i,j,k,l)&T2s(k,l)(2,0,1,0)");
873  test_for_zero(t4ddg_1(2, 0, 1, 1) - (t4ddg_2(2, 0, 1, 1) * t2s_3(1, 1)),
874  "T4ddg(i,j,k,l)&T2s(k,l)(2,0,1,1)");
875  test_for_zero(t4ddg_1(2, 0, 1, 2) - (t4ddg_2(2, 0, 1, 2) * t2s_3(1, 2)),
876  "T4ddg(i,j,k,l)&T2s(k,l)(2,0,1,2)");
877  test_for_zero(t4ddg_1(2, 0, 2, 0) - (t4ddg_2(2, 0, 2, 0) * t2s_3(2, 0)),
878  "T4ddg(i,j,k,l)&T2s(k,l)(2,0,2,0)");
879  test_for_zero(t4ddg_1(2, 0, 2, 1) - (t4ddg_2(2, 0, 2, 1) * t2s_3(2, 1)),
880  "T4ddg(i,j,k,l)&T2s(k,l)(2,0,2,1)");
881  test_for_zero(t4ddg_1(2, 0, 2, 2) - (t4ddg_2(2, 0, 2, 2) * t2s_3(2, 2)),
882  "T4ddg(i,j,k,l)&T2s(k,l)(2,0,2,2)");
883  test_for_zero(t4ddg_1(2, 1, 0, 0) - (t4ddg_2(2, 1, 0, 0) * t2s_3(0, 0)),
884  "T4ddg(i,j,k,l)&T2s(k,l)(2,1,0,0)");
885  test_for_zero(t4ddg_1(2, 1, 0, 1) - (t4ddg_2(2, 1, 0, 1) * t2s_3(0, 1)),
886  "T4ddg(i,j,k,l)&T2s(k,l)(2,1,0,1)");
887  test_for_zero(t4ddg_1(2, 1, 0, 2) - (t4ddg_2(2, 1, 0, 2) * t2s_3(0, 2)),
888  "T4ddg(i,j,k,l)&T2s(k,l)(2,1,0,2)");
889  test_for_zero(t4ddg_1(2, 1, 1, 0) - (t4ddg_2(2, 1, 1, 0) * t2s_3(1, 0)),
890  "T4ddg(i,j,k,l)&T2s(k,l)(2,1,1,0)");
891  test_for_zero(t4ddg_1(2, 1, 1, 1) - (t4ddg_2(2, 1, 1, 1) * t2s_3(1, 1)),
892  "T4ddg(i,j,k,l)&T2s(k,l)(2,1,1,1)");
893  test_for_zero(t4ddg_1(2, 1, 1, 2) - (t4ddg_2(2, 1, 1, 2) * t2s_3(1, 2)),
894  "T4ddg(i,j,k,l)&T2s(k,l)(2,1,1,2)");
895  test_for_zero(t4ddg_1(2, 1, 2, 0) - (t4ddg_2(2, 1, 2, 0) * t2s_3(2, 0)),
896  "T4ddg(i,j,k,l)&T2s(k,l)(2,1,2,0)");
897  test_for_zero(t4ddg_1(2, 1, 2, 1) - (t4ddg_2(2, 1, 2, 1) * t2s_3(2, 1)),
898  "T4ddg(i,j,k,l)&T2s(k,l)(2,1,2,1)");
899  test_for_zero(t4ddg_1(2, 1, 2, 2) - (t4ddg_2(2, 1, 2, 2) * t2s_3(2, 2)),
900  "T4ddg(i,j,k,l)&T2s(k,l)(2,1,2,2)");
901  test_for_zero(t4ddg_1(2, 2, 0, 0) - (t4ddg_2(2, 2, 0, 0) * t2s_3(0, 0)),
902  "T4ddg(i,j,k,l)&T2s(k,l)(2,2,0,0)");
903  test_for_zero(t4ddg_1(2, 2, 0, 1) - (t4ddg_2(2, 2, 0, 1) * t2s_3(0, 1)),
904  "T4ddg(i,j,k,l)&T2s(k,l)(2,2,0,1)");
905  test_for_zero(t4ddg_1(2, 2, 0, 2) - (t4ddg_2(2, 2, 0, 2) * t2s_3(0, 2)),
906  "T4ddg(i,j,k,l)&T2s(k,l)(2,2,0,2)");
907  test_for_zero(t4ddg_1(2, 2, 1, 0) - (t4ddg_2(2, 2, 1, 0) * t2s_3(1, 0)),
908  "T4ddg(i,j,k,l)&T2s(k,l)(2,2,1,0)");
909  test_for_zero(t4ddg_1(2, 2, 1, 1) - (t4ddg_2(2, 2, 1, 1) * t2s_3(1, 1)),
910  "T4ddg(i,j,k,l)&T2s(k,l)(2,2,1,1)");
911  test_for_zero(t4ddg_1(2, 2, 1, 2) - (t4ddg_2(2, 2, 1, 2) * t2s_3(1, 2)),
912  "T4ddg(i,j,k,l)&T2s(k,l)(2,2,1,2)");
913  test_for_zero(t4ddg_1(2, 2, 2, 0) - (t4ddg_2(2, 2, 2, 0) * t2s_3(2, 0)),
914  "T4ddg(i,j,k,l)&T2s(k,l)(2,2,2,0)");
915  test_for_zero(t4ddg_1(2, 2, 2, 1) - (t4ddg_2(2, 2, 2, 1) * t2s_3(2, 1)),
916  "T4ddg(i,j,k,l)&T2s(k,l)(2,2,2,1)");
917  test_for_zero(t4ddg_1(2, 2, 2, 2) - (t4ddg_2(2, 2, 2, 2) * t2s_3(2, 2)),
918  "T4ddg(i,j,k,l)&T2s(k,l)(2,2,2,2)");
919 
920  t4ddg_1(i, j, k, l) = (t2s_2(k, l) & t4ddg_2(i, j, k, l));
921  test_for_zero(t4ddg_1(0, 0, 0, 0) - (t4ddg_2(0, 0, 0, 0) * t2s_2(0, 0)),
922  "T2s(k,l)&T4ddg(i,j,k,l)(0,0,0,0)");
923  test_for_zero(t4ddg_1(0, 0, 0, 1) - (t4ddg_2(0, 0, 0, 1) * t2s_2(0, 1)),
924  "T2s(k,l)&T4ddg(i,j,k,l)(0,0,0,1)");
925  test_for_zero(t4ddg_1(0, 0, 0, 2) - (t4ddg_2(0, 0, 0, 2) * t2s_2(0, 2)),
926  "T2s(k,l)&T4ddg(i,j,k,l)(0,0,0,2)");
927  test_for_zero(t4ddg_1(0, 0, 1, 0) - (t4ddg_2(0, 0, 1, 0) * t2s_2(1, 0)),
928  "T2s(k,l)&T4ddg(i,j,k,l)(0,0,1,0)");
929  test_for_zero(t4ddg_1(0, 0, 1, 1) - (t4ddg_2(0, 0, 1, 1) * t2s_2(1, 1)),
930  "T2s(k,l)&T4ddg(i,j,k,l)(0,0,1,1)");
931  test_for_zero(t4ddg_1(0, 0, 1, 2) - (t4ddg_2(0, 0, 1, 2) * t2s_2(1, 2)),
932  "T2s(k,l)&T4ddg(i,j,k,l)(0,0,1,2)");
933  test_for_zero(t4ddg_1(0, 0, 2, 0) - (t4ddg_2(0, 0, 2, 0) * t2s_2(2, 0)),
934  "T2s(k,l)&T4ddg(i,j,k,l)(0,0,2,0)");
935  test_for_zero(t4ddg_1(0, 0, 2, 1) - (t4ddg_2(0, 0, 2, 1) * t2s_2(2, 1)),
936  "T2s(k,l)&T4ddg(i,j,k,l)(0,0,2,1)");
937  test_for_zero(t4ddg_1(0, 0, 2, 2) - (t4ddg_2(0, 0, 2, 2) * t2s_2(2, 2)),
938  "T2s(k,l)&T4ddg(i,j,k,l)(0,0,2,2)");
939  test_for_zero(t4ddg_1(0, 1, 0, 0) - (t4ddg_2(0, 1, 0, 0) * t2s_2(0, 0)),
940  "T2s(k,l)&T4ddg(i,j,k,l)(0,1,0,0)");
941  test_for_zero(t4ddg_1(0, 1, 0, 1) - (t4ddg_2(0, 1, 0, 1) * t2s_2(0, 1)),
942  "T2s(k,l)&T4ddg(i,j,k,l)(0,1,0,1)");
943  test_for_zero(t4ddg_1(0, 1, 0, 2) - (t4ddg_2(0, 1, 0, 2) * t2s_2(0, 2)),
944  "T2s(k,l)&T4ddg(i,j,k,l)(0,1,0,2)");
945  test_for_zero(t4ddg_1(0, 1, 1, 0) - (t4ddg_2(0, 1, 1, 0) * t2s_2(1, 0)),
946  "T2s(k,l)&T4ddg(i,j,k,l)(0,1,1,0)");
947  test_for_zero(t4ddg_1(0, 1, 1, 1) - (t4ddg_2(0, 1, 1, 1) * t2s_2(1, 1)),
948  "T2s(k,l)&T4ddg(i,j,k,l)(0,1,1,1)");
949  test_for_zero(t4ddg_1(0, 1, 1, 2) - (t4ddg_2(0, 1, 1, 2) * t2s_2(1, 2)),
950  "T2s(k,l)&T4ddg(i,j,k,l)(0,1,1,2)");
951  test_for_zero(t4ddg_1(0, 1, 2, 0) - (t4ddg_2(0, 1, 2, 0) * t2s_2(2, 0)),
952  "T2s(k,l)&T4ddg(i,j,k,l)(0,1,2,0)");
953  test_for_zero(t4ddg_1(0, 1, 2, 1) - (t4ddg_2(0, 1, 2, 1) * t2s_2(2, 1)),
954  "T2s(k,l)&T4ddg(i,j,k,l)(0,1,2,1)");
955  test_for_zero(t4ddg_1(0, 1, 2, 2) - (t4ddg_2(0, 1, 2, 2) * t2s_2(2, 2)),
956  "T2s(k,l)&T4ddg(i,j,k,l)(0,1,2,2)");
957  test_for_zero(t4ddg_1(0, 2, 0, 0) - (t4ddg_2(0, 2, 0, 0) * t2s_2(0, 0)),
958  "T2s(k,l)&T4ddg(i,j,k,l)(0,2,0,0)");
959  test_for_zero(t4ddg_1(0, 2, 0, 1) - (t4ddg_2(0, 2, 0, 1) * t2s_2(0, 1)),
960  "T2s(k,l)&T4ddg(i,j,k,l)(0,2,0,1)");
961  test_for_zero(t4ddg_1(0, 2, 0, 2) - (t4ddg_2(0, 2, 0, 2) * t2s_2(0, 2)),
962  "T2s(k,l)&T4ddg(i,j,k,l)(0,2,0,2)");
963  test_for_zero(t4ddg_1(0, 2, 1, 0) - (t4ddg_2(0, 2, 1, 0) * t2s_2(1, 0)),
964  "T2s(k,l)&T4ddg(i,j,k,l)(0,2,1,0)");
965  test_for_zero(t4ddg_1(0, 2, 1, 1) - (t4ddg_2(0, 2, 1, 1) * t2s_2(1, 1)),
966  "T2s(k,l)&T4ddg(i,j,k,l)(0,2,1,1)");
967  test_for_zero(t4ddg_1(0, 2, 1, 2) - (t4ddg_2(0, 2, 1, 2) * t2s_2(1, 2)),
968  "T2s(k,l)&T4ddg(i,j,k,l)(0,2,1,2)");
969  test_for_zero(t4ddg_1(0, 2, 2, 0) - (t4ddg_2(0, 2, 2, 0) * t2s_2(2, 0)),
970  "T2s(k,l)&T4ddg(i,j,k,l)(0,2,2,0)");
971  test_for_zero(t4ddg_1(0, 2, 2, 1) - (t4ddg_2(0, 2, 2, 1) * t2s_2(2, 1)),
972  "T2s(k,l)&T4ddg(i,j,k,l)(0,2,2,1)");
973  test_for_zero(t4ddg_1(0, 2, 2, 2) - (t4ddg_2(0, 2, 2, 2) * t2s_2(2, 2)),
974  "T2s(k,l)&T4ddg(i,j,k,l)(0,2,2,2)");
975  test_for_zero(t4ddg_1(1, 0, 0, 0) - (t4ddg_2(1, 0, 0, 0) * t2s_2(0, 0)),
976  "T2s(k,l)&T4ddg(i,j,k,l)(1,0,0,0)");
977  test_for_zero(t4ddg_1(1, 0, 0, 1) - (t4ddg_2(1, 0, 0, 1) * t2s_2(0, 1)),
978  "T2s(k,l)&T4ddg(i,j,k,l)(1,0,0,1)");
979  test_for_zero(t4ddg_1(1, 0, 0, 2) - (t4ddg_2(1, 0, 0, 2) * t2s_2(0, 2)),
980  "T2s(k,l)&T4ddg(i,j,k,l)(1,0,0,2)");
981  test_for_zero(t4ddg_1(1, 0, 1, 0) - (t4ddg_2(1, 0, 1, 0) * t2s_2(1, 0)),
982  "T2s(k,l)&T4ddg(i,j,k,l)(1,0,1,0)");
983  test_for_zero(t4ddg_1(1, 0, 1, 1) - (t4ddg_2(1, 0, 1, 1) * t2s_2(1, 1)),
984  "T2s(k,l)&T4ddg(i,j,k,l)(1,0,1,1)");
985  test_for_zero(t4ddg_1(1, 0, 1, 2) - (t4ddg_2(1, 0, 1, 2) * t2s_2(1, 2)),
986  "T2s(k,l)&T4ddg(i,j,k,l)(1,0,1,2)");
987  test_for_zero(t4ddg_1(1, 0, 2, 0) - (t4ddg_2(1, 0, 2, 0) * t2s_2(2, 0)),
988  "T2s(k,l)&T4ddg(i,j,k,l)(1,0,2,0)");
989  test_for_zero(t4ddg_1(1, 0, 2, 1) - (t4ddg_2(1, 0, 2, 1) * t2s_2(2, 1)),
990  "T2s(k,l)&T4ddg(i,j,k,l)(1,0,2,1)");
991  test_for_zero(t4ddg_1(1, 0, 2, 2) - (t4ddg_2(1, 0, 2, 2) * t2s_2(2, 2)),
992  "T2s(k,l)&T4ddg(i,j,k,l)(1,0,2,2)");
993  test_for_zero(t4ddg_1(1, 1, 0, 0) - (t4ddg_2(1, 1, 0, 0) * t2s_2(0, 0)),
994  "T2s(k,l)&T4ddg(i,j,k,l)(1,1,0,0)");
995  test_for_zero(t4ddg_1(1, 1, 0, 1) - (t4ddg_2(1, 1, 0, 1) * t2s_2(0, 1)),
996  "T2s(k,l)&T4ddg(i,j,k,l)(1,1,0,1)");
997  test_for_zero(t4ddg_1(1, 1, 0, 2) - (t4ddg_2(1, 1, 0, 2) * t2s_2(0, 2)),
998  "T2s(k,l)&T4ddg(i,j,k,l)(1,1,0,2)");
999  test_for_zero(t4ddg_1(1, 1, 1, 0) - (t4ddg_2(1, 1, 1, 0) * t2s_2(1, 0)),
1000  "T2s(k,l)&T4ddg(i,j,k,l)(1,1,1,0)");
1001  test_for_zero(t4ddg_1(1, 1, 1, 1) - (t4ddg_2(1, 1, 1, 1) * t2s_2(1, 1)),
1002  "T2s(k,l)&T4ddg(i,j,k,l)(1,1,1,1)");
1003  test_for_zero(t4ddg_1(1, 1, 1, 2) - (t4ddg_2(1, 1, 1, 2) * t2s_2(1, 2)),
1004  "T2s(k,l)&T4ddg(i,j,k,l)(1,1,1,2)");
1005  test_for_zero(t4ddg_1(1, 1, 2, 0) - (t4ddg_2(1, 1, 2, 0) * t2s_2(2, 0)),
1006  "T2s(k,l)&T4ddg(i,j,k,l)(1,1,2,0)");
1007  test_for_zero(t4ddg_1(1, 1, 2, 1) - (t4ddg_2(1, 1, 2, 1) * t2s_2(2, 1)),
1008  "T2s(k,l)&T4ddg(i,j,k,l)(1,1,2,1)");
1009  test_for_zero(t4ddg_1(1, 1, 2, 2) - (t4ddg_2(1, 1, 2, 2) * t2s_2(2, 2)),
1010  "T2s(k,l)&T4ddg(i,j,k,l)(1,1,2,2)");
1011  test_for_zero(t4ddg_1(1, 2, 0, 0) - (t4ddg_2(1, 2, 0, 0) * t2s_2(0, 0)),
1012  "T2s(k,l)&T4ddg(i,j,k,l)(1,2,0,0)");
1013  test_for_zero(t4ddg_1(1, 2, 0, 1) - (t4ddg_2(1, 2, 0, 1) * t2s_2(0, 1)),
1014  "T2s(k,l)&T4ddg(i,j,k,l)(1,2,0,1)");
1015  test_for_zero(t4ddg_1(1, 2, 0, 2) - (t4ddg_2(1, 2, 0, 2) * t2s_2(0, 2)),
1016  "T2s(k,l)&T4ddg(i,j,k,l)(1,2,0,2)");
1017  test_for_zero(t4ddg_1(1, 2, 1, 0) - (t4ddg_2(1, 2, 1, 0) * t2s_2(1, 0)),
1018  "T2s(k,l)&T4ddg(i,j,k,l)(1,2,1,0)");
1019  test_for_zero(t4ddg_1(1, 2, 1, 1) - (t4ddg_2(1, 2, 1, 1) * t2s_2(1, 1)),
1020  "T2s(k,l)&T4ddg(i,j,k,l)(1,2,1,1)");
1021  test_for_zero(t4ddg_1(1, 2, 1, 2) - (t4ddg_2(1, 2, 1, 2) * t2s_2(1, 2)),
1022  "T2s(k,l)&T4ddg(i,j,k,l)(1,2,1,2)");
1023  test_for_zero(t4ddg_1(1, 2, 2, 0) - (t4ddg_2(1, 2, 2, 0) * t2s_2(2, 0)),
1024  "T2s(k,l)&T4ddg(i,j,k,l)(1,2,2,0)");
1025  test_for_zero(t4ddg_1(1, 2, 2, 1) - (t4ddg_2(1, 2, 2, 1) * t2s_2(2, 1)),
1026  "T2s(k,l)&T4ddg(i,j,k,l)(1,2,2,1)");
1027  test_for_zero(t4ddg_1(1, 2, 2, 2) - (t4ddg_2(1, 2, 2, 2) * t2s_2(2, 2)),
1028  "T2s(k,l)&T4ddg(i,j,k,l)(1,2,2,2)");
1029  test_for_zero(t4ddg_1(2, 0, 0, 0) - (t4ddg_2(2, 0, 0, 0) * t2s_2(0, 0)),
1030  "T2s(k,l)&T4ddg(i,j,k,l)(2,0,0,0)");
1031  test_for_zero(t4ddg_1(2, 0, 0, 1) - (t4ddg_2(2, 0, 0, 1) * t2s_2(0, 1)),
1032  "T2s(k,l)&T4ddg(i,j,k,l)(2,0,0,1)");
1033  test_for_zero(t4ddg_1(2, 0, 0, 2) - (t4ddg_2(2, 0, 0, 2) * t2s_2(0, 2)),
1034  "T2s(k,l)&T4ddg(i,j,k,l)(2,0,0,2)");
1035  test_for_zero(t4ddg_1(2, 0, 1, 0) - (t4ddg_2(2, 0, 1, 0) * t2s_2(1, 0)),
1036  "T2s(k,l)&T4ddg(i,j,k,l)(2,0,1,0)");
1037  test_for_zero(t4ddg_1(2, 0, 1, 1) - (t4ddg_2(2, 0, 1, 1) * t2s_2(1, 1)),
1038  "T2s(k,l)&T4ddg(i,j,k,l)(2,0,1,1)");
1039  test_for_zero(t4ddg_1(2, 0, 1, 2) - (t4ddg_2(2, 0, 1, 2) * t2s_2(1, 2)),
1040  "T2s(k,l)&T4ddg(i,j,k,l)(2,0,1,2)");
1041  test_for_zero(t4ddg_1(2, 0, 2, 0) - (t4ddg_2(2, 0, 2, 0) * t2s_2(2, 0)),
1042  "T2s(k,l)&T4ddg(i,j,k,l)(2,0,2,0)");
1043  test_for_zero(t4ddg_1(2, 0, 2, 1) - (t4ddg_2(2, 0, 2, 1) * t2s_2(2, 1)),
1044  "T2s(k,l)&T4ddg(i,j,k,l)(2,0,2,1)");
1045  test_for_zero(t4ddg_1(2, 0, 2, 2) - (t4ddg_2(2, 0, 2, 2) * t2s_2(2, 2)),
1046  "T2s(k,l)&T4ddg(i,j,k,l)(2,0,2,2)");
1047  test_for_zero(t4ddg_1(2, 1, 0, 0) - (t4ddg_2(2, 1, 0, 0) * t2s_2(0, 0)),
1048  "T2s(k,l)&T4ddg(i,j,k,l)(2,1,0,0)");
1049  test_for_zero(t4ddg_1(2, 1, 0, 1) - (t4ddg_2(2, 1, 0, 1) * t2s_2(0, 1)),
1050  "T2s(k,l)&T4ddg(i,j,k,l)(2,1,0,1)");
1051  test_for_zero(t4ddg_1(2, 1, 0, 2) - (t4ddg_2(2, 1, 0, 2) * t2s_2(0, 2)),
1052  "T2s(k,l)&T4ddg(i,j,k,l)(2,1,0,2)");
1053  test_for_zero(t4ddg_1(2, 1, 1, 0) - (t4ddg_2(2, 1, 1, 0) * t2s_2(1, 0)),
1054  "T2s(k,l)&T4ddg(i,j,k,l)(2,1,1,0)");
1055  test_for_zero(t4ddg_1(2, 1, 1, 1) - (t4ddg_2(2, 1, 1, 1) * t2s_2(1, 1)),
1056  "T2s(k,l)&T4ddg(i,j,k,l)(2,1,1,1)");
1057  test_for_zero(t4ddg_1(2, 1, 1, 2) - (t4ddg_2(2, 1, 1, 2) * t2s_2(1, 2)),
1058  "T2s(k,l)&T4ddg(i,j,k,l)(2,1,1,2)");
1059  test_for_zero(t4ddg_1(2, 1, 2, 0) - (t4ddg_2(2, 1, 2, 0) * t2s_2(2, 0)),
1060  "T2s(k,l)&T4ddg(i,j,k,l)(2,1,2,0)");
1061  test_for_zero(t4ddg_1(2, 1, 2, 1) - (t4ddg_2(2, 1, 2, 1) * t2s_2(2, 1)),
1062  "T2s(k,l)&T4ddg(i,j,k,l)(2,1,2,1)");
1063  test_for_zero(t4ddg_1(2, 1, 2, 2) - (t4ddg_2(2, 1, 2, 2) * t2s_2(2, 2)),
1064  "T2s(k,l)&T4ddg(i,j,k,l)(2,1,2,2)");
1065  test_for_zero(t4ddg_1(2, 2, 0, 0) - (t4ddg_2(2, 2, 0, 0) * t2s_2(0, 0)),
1066  "T2s(k,l)&T4ddg(i,j,k,l)(2,2,0,0)");
1067  test_for_zero(t4ddg_1(2, 2, 0, 1) - (t4ddg_2(2, 2, 0, 1) * t2s_2(0, 1)),
1068  "T2s(k,l)&T4ddg(i,j,k,l)(2,2,0,1)");
1069  test_for_zero(t4ddg_1(2, 2, 0, 2) - (t4ddg_2(2, 2, 0, 2) * t2s_2(0, 2)),
1070  "T2s(k,l)&T4ddg(i,j,k,l)(2,2,0,2)");
1071  test_for_zero(t4ddg_1(2, 2, 1, 0) - (t4ddg_2(2, 2, 1, 0) * t2s_2(1, 0)),
1072  "T2s(k,l)&T4ddg(i,j,k,l)(2,2,1,0)");
1073  test_for_zero(t4ddg_1(2, 2, 1, 1) - (t4ddg_2(2, 2, 1, 1) * t2s_2(1, 1)),
1074  "T2s(k,l)&T4ddg(i,j,k,l)(2,2,1,1)");
1075  test_for_zero(t4ddg_1(2, 2, 1, 2) - (t4ddg_2(2, 2, 1, 2) * t2s_2(1, 2)),
1076  "T2s(k,l)&T4ddg(i,j,k,l)(2,2,1,2)");
1077  test_for_zero(t4ddg_1(2, 2, 2, 0) - (t4ddg_2(2, 2, 2, 0) * t2s_2(2, 0)),
1078  "T2s(k,l)&T4ddg(i,j,k,l)(2,2,2,0)");
1079  test_for_zero(t4ddg_1(2, 2, 2, 1) - (t4ddg_2(2, 2, 2, 1) * t2s_2(2, 1)),
1080  "T2s(k,l)&T4ddg(i,j,k,l)(2,2,2,1)");
1081  test_for_zero(t4ddg_1(2, 2, 2, 2) - (t4ddg_2(2, 2, 2, 2) * t2s_2(2, 2)),
1082  "T2s(k,l)&T4ddg(i,j,k,l)(2,2,2,2)");
1083 
1084  t4ddg_1(i, j, k, l) = (t4ddg_2(i, j, k, l) % t2s_3(i, j));
1085  test_for_zero(t4ddg_1(0, 0, 0, 0) - (t4ddg_2(0, 0, 0, 0) / t2s_3(0, 0)),
1086  "T4ddg(i,j,k,l)%T2s(i,j)(0,0,0,0)");
1087  test_for_zero(t4ddg_1(0, 0, 0, 1) - (t4ddg_2(0, 0, 0, 1) / t2s_3(0, 0)),
1088  "T4ddg(i,j,k,l)%T2s(i,j)(0,0,0,1)");
1089  test_for_zero(t4ddg_1(0, 0, 0, 2) - (t4ddg_2(0, 0, 0, 2) / t2s_3(0, 0)),
1090  "T4ddg(i,j,k,l)%T2s(i,j)(0,0,0,2)");
1091  test_for_zero(t4ddg_1(0, 0, 1, 0) - (t4ddg_2(0, 0, 1, 0) / t2s_3(0, 0)),
1092  "T4ddg(i,j,k,l)%T2s(i,j)(0,0,1,0)");
1093  test_for_zero(t4ddg_1(0, 0, 1, 1) - (t4ddg_2(0, 0, 1, 1) / t2s_3(0, 0)),
1094  "T4ddg(i,j,k,l)%T2s(i,j)(0,0,1,1)");
1095  test_for_zero(t4ddg_1(0, 0, 1, 2) - (t4ddg_2(0, 0, 1, 2) / t2s_3(0, 0)),
1096  "T4ddg(i,j,k,l)%T2s(i,j)(0,0,1,2)");
1097  test_for_zero(t4ddg_1(0, 0, 2, 0) - (t4ddg_2(0, 0, 2, 0) / t2s_3(0, 0)),
1098  "T4ddg(i,j,k,l)%T2s(i,j)(0,0,2,0)");
1099  test_for_zero(t4ddg_1(0, 0, 2, 1) - (t4ddg_2(0, 0, 2, 1) / t2s_3(0, 0)),
1100  "T4ddg(i,j,k,l)%T2s(i,j)(0,0,2,1)");
1101  test_for_zero(t4ddg_1(0, 0, 2, 2) - (t4ddg_2(0, 0, 2, 2) / t2s_3(0, 0)),
1102  "T4ddg(i,j,k,l)%T2s(i,j)(0,0,2,2)");
1103  test_for_zero(t4ddg_1(0, 1, 0, 0) - (t4ddg_2(0, 1, 0, 0) / t2s_3(0, 1)),
1104  "T4ddg(i,j,k,l)%T2s(i,j)(0,1,0,0)");
1105  test_for_zero(t4ddg_1(0, 1, 0, 1) - (t4ddg_2(0, 1, 0, 1) / t2s_3(0, 1)),
1106  "T4ddg(i,j,k,l)%T2s(i,j)(0,1,0,1)");
1107  test_for_zero(t4ddg_1(0, 1, 0, 2) - (t4ddg_2(0, 1, 0, 2) / t2s_3(0, 1)),
1108  "T4ddg(i,j,k,l)%T2s(i,j)(0,1,0,2)");
1109  test_for_zero(t4ddg_1(0, 1, 1, 0) - (t4ddg_2(0, 1, 1, 0) / t2s_3(0, 1)),
1110  "T4ddg(i,j,k,l)%T2s(i,j)(0,1,1,0)");
1111  test_for_zero(t4ddg_1(0, 1, 1, 1) - (t4ddg_2(0, 1, 1, 1) / t2s_3(0, 1)),
1112  "T4ddg(i,j,k,l)%T2s(i,j)(0,1,1,1)");
1113  test_for_zero(t4ddg_1(0, 1, 1, 2) - (t4ddg_2(0, 1, 1, 2) / t2s_3(0, 1)),
1114  "T4ddg(i,j,k,l)%T2s(i,j)(0,1,1,2)");
1115  test_for_zero(t4ddg_1(0, 1, 2, 0) - (t4ddg_2(0, 1, 2, 0) / t2s_3(0, 1)),
1116  "T4ddg(i,j,k,l)%T2s(i,j)(0,1,2,0)");
1117  test_for_zero(t4ddg_1(0, 1, 2, 1) - (t4ddg_2(0, 1, 2, 1) / t2s_3(0, 1)),
1118  "T4ddg(i,j,k,l)%T2s(i,j)(0,1,2,1)");
1119  test_for_zero(t4ddg_1(0, 1, 2, 2) - (t4ddg_2(0, 1, 2, 2) / t2s_3(0, 1)),
1120  "T4ddg(i,j,k,l)%T2s(i,j)(0,1,2,2)");
1121  test_for_zero(t4ddg_1(0, 2, 0, 0) - (t4ddg_2(0, 2, 0, 0) / t2s_3(0, 2)),
1122  "T4ddg(i,j,k,l)%T2s(i,j)(0,2,0,0)");
1123  test_for_zero(t4ddg_1(0, 2, 0, 1) - (t4ddg_2(0, 2, 0, 1) / t2s_3(0, 2)),
1124  "T4ddg(i,j,k,l)%T2s(i,j)(0,2,0,1)");
1125  test_for_zero(t4ddg_1(0, 2, 0, 2) - (t4ddg_2(0, 2, 0, 2) / t2s_3(0, 2)),
1126  "T4ddg(i,j,k,l)%T2s(i,j)(0,2,0,2)");
1127  test_for_zero(t4ddg_1(0, 2, 1, 0) - (t4ddg_2(0, 2, 1, 0) / t2s_3(0, 2)),
1128  "T4ddg(i,j,k,l)%T2s(i,j)(0,2,1,0)");
1129  test_for_zero(t4ddg_1(0, 2, 1, 1) - (t4ddg_2(0, 2, 1, 1) / t2s_3(0, 2)),
1130  "T4ddg(i,j,k,l)%T2s(i,j)(0,2,1,1)");
1131  test_for_zero(t4ddg_1(0, 2, 1, 2) - (t4ddg_2(0, 2, 1, 2) / t2s_3(0, 2)),
1132  "T4ddg(i,j,k,l)%T2s(i,j)(0,2,1,2)");
1133  test_for_zero(t4ddg_1(0, 2, 2, 0) - (t4ddg_2(0, 2, 2, 0) / t2s_3(0, 2)),
1134  "T4ddg(i,j,k,l)%T2s(i,j)(0,2,2,0)");
1135  test_for_zero(t4ddg_1(0, 2, 2, 1) - (t4ddg_2(0, 2, 2, 1) / t2s_3(0, 2)),
1136  "T4ddg(i,j,k,l)%T2s(i,j)(0,2,2,1)");
1137  test_for_zero(t4ddg_1(0, 2, 2, 2) - (t4ddg_2(0, 2, 2, 2) / t2s_3(0, 2)),
1138  "T4ddg(i,j,k,l)%T2s(i,j)(0,2,2,2)");
1139  test_for_zero(t4ddg_1(1, 0, 0, 0) - (t4ddg_2(1, 0, 0, 0) / t2s_3(1, 0)),
1140  "T4ddg(i,j,k,l)%T2s(i,j)(1,0,0,0)");
1141  test_for_zero(t4ddg_1(1, 0, 0, 1) - (t4ddg_2(1, 0, 0, 1) / t2s_3(1, 0)),
1142  "T4ddg(i,j,k,l)%T2s(i,j)(1,0,0,1)");
1143  test_for_zero(t4ddg_1(1, 0, 0, 2) - (t4ddg_2(1, 0, 0, 2) / t2s_3(1, 0)),
1144  "T4ddg(i,j,k,l)%T2s(i,j)(1,0,0,2)");
1145  test_for_zero(t4ddg_1(1, 0, 1, 0) - (t4ddg_2(1, 0, 1, 0) / t2s_3(1, 0)),
1146  "T4ddg(i,j,k,l)%T2s(i,j)(1,0,1,0)");
1147  test_for_zero(t4ddg_1(1, 0, 1, 1) - (t4ddg_2(1, 0, 1, 1) / t2s_3(1, 0)),
1148  "T4ddg(i,j,k,l)%T2s(i,j)(1,0,1,1)");
1149  test_for_zero(t4ddg_1(1, 0, 1, 2) - (t4ddg_2(1, 0, 1, 2) / t2s_3(1, 0)),
1150  "T4ddg(i,j,k,l)%T2s(i,j)(1,0,1,2)");
1151  test_for_zero(t4ddg_1(1, 0, 2, 0) - (t4ddg_2(1, 0, 2, 0) / t2s_3(1, 0)),
1152  "T4ddg(i,j,k,l)%T2s(i,j)(1,0,2,0)");
1153  test_for_zero(t4ddg_1(1, 0, 2, 1) - (t4ddg_2(1, 0, 2, 1) / t2s_3(1, 0)),
1154  "T4ddg(i,j,k,l)%T2s(i,j)(1,0,2,1)");
1155  test_for_zero(t4ddg_1(1, 0, 2, 2) - (t4ddg_2(1, 0, 2, 2) / t2s_3(1, 0)),
1156  "T4ddg(i,j,k,l)%T2s(i,j)(1,0,2,2)");
1157  test_for_zero(t4ddg_1(1, 1, 0, 0) - (t4ddg_2(1, 1, 0, 0) / t2s_3(1, 1)),
1158  "T4ddg(i,j,k,l)%T2s(i,j)(1,1,0,0)");
1159  test_for_zero(t4ddg_1(1, 1, 0, 1) - (t4ddg_2(1, 1, 0, 1) / t2s_3(1, 1)),
1160  "T4ddg(i,j,k,l)%T2s(i,j)(1,1,0,1)");
1161  test_for_zero(t4ddg_1(1, 1, 0, 2) - (t4ddg_2(1, 1, 0, 2) / t2s_3(1, 1)),
1162  "T4ddg(i,j,k,l)%T2s(i,j)(1,1,0,2)");
1163  test_for_zero(t4ddg_1(1, 1, 1, 0) - (t4ddg_2(1, 1, 1, 0) / t2s_3(1, 1)),
1164  "T4ddg(i,j,k,l)%T2s(i,j)(1,1,1,0)");
1165  test_for_zero(t4ddg_1(1, 1, 1, 1) - (t4ddg_2(1, 1, 1, 1) / t2s_3(1, 1)),
1166  "T4ddg(i,j,k,l)%T2s(i,j)(1,1,1,1)");
1167  test_for_zero(t4ddg_1(1, 1, 1, 2) - (t4ddg_2(1, 1, 1, 2) / t2s_3(1, 1)),
1168  "T4ddg(i,j,k,l)%T2s(i,j)(1,1,1,2)");
1169  test_for_zero(t4ddg_1(1, 1, 2, 0) - (t4ddg_2(1, 1, 2, 0) / t2s_3(1, 1)),
1170  "T4ddg(i,j,k,l)%T2s(i,j)(1,1,2,0)");
1171  test_for_zero(t4ddg_1(1, 1, 2, 1) - (t4ddg_2(1, 1, 2, 1) / t2s_3(1, 1)),
1172  "T4ddg(i,j,k,l)%T2s(i,j)(1,1,2,1)");
1173  test_for_zero(t4ddg_1(1, 1, 2, 2) - (t4ddg_2(1, 1, 2, 2) / t2s_3(1, 1)),
1174  "T4ddg(i,j,k,l)%T2s(i,j)(1,1,2,2)");
1175  test_for_zero(t4ddg_1(1, 2, 0, 0) - (t4ddg_2(1, 2, 0, 0) / t2s_3(1, 2)),
1176  "T4ddg(i,j,k,l)%T2s(i,j)(1,2,0,0)");
1177  test_for_zero(t4ddg_1(1, 2, 0, 1) - (t4ddg_2(1, 2, 0, 1) / t2s_3(1, 2)),
1178  "T4ddg(i,j,k,l)%T2s(i,j)(1,2,0,1)");
1179  test_for_zero(t4ddg_1(1, 2, 0, 2) - (t4ddg_2(1, 2, 0, 2) / t2s_3(1, 2)),
1180  "T4ddg(i,j,k,l)%T2s(i,j)(1,2,0,2)");
1181  test_for_zero(t4ddg_1(1, 2, 1, 0) - (t4ddg_2(1, 2, 1, 0) / t2s_3(1, 2)),
1182  "T4ddg(i,j,k,l)%T2s(i,j)(1,2,1,0)");
1183  test_for_zero(t4ddg_1(1, 2, 1, 1) - (t4ddg_2(1, 2, 1, 1) / t2s_3(1, 2)),
1184  "T4ddg(i,j,k,l)%T2s(i,j)(1,2,1,1)");
1185  test_for_zero(t4ddg_1(1, 2, 1, 2) - (t4ddg_2(1, 2, 1, 2) / t2s_3(1, 2)),
1186  "T4ddg(i,j,k,l)%T2s(i,j)(1,2,1,2)");
1187  test_for_zero(t4ddg_1(1, 2, 2, 0) - (t4ddg_2(1, 2, 2, 0) / t2s_3(1, 2)),
1188  "T4ddg(i,j,k,l)%T2s(i,j)(1,2,2,0)");
1189  test_for_zero(t4ddg_1(1, 2, 2, 1) - (t4ddg_2(1, 2, 2, 1) / t2s_3(1, 2)),
1190  "T4ddg(i,j,k,l)%T2s(i,j)(1,2,2,1)");
1191  test_for_zero(t4ddg_1(1, 2, 2, 2) - (t4ddg_2(1, 2, 2, 2) / t2s_3(1, 2)),
1192  "T4ddg(i,j,k,l)%T2s(i,j)(1,2,2,2)");
1193  test_for_zero(t4ddg_1(2, 0, 0, 0) - (t4ddg_2(2, 0, 0, 0) / t2s_3(2, 0)),
1194  "T4ddg(i,j,k,l)%T2s(i,j)(2,0,0,0)");
1195  test_for_zero(t4ddg_1(2, 0, 0, 1) - (t4ddg_2(2, 0, 0, 1) / t2s_3(2, 0)),
1196  "T4ddg(i,j,k,l)%T2s(i,j)(2,0,0,1)");
1197  test_for_zero(t4ddg_1(2, 0, 0, 2) - (t4ddg_2(2, 0, 0, 2) / t2s_3(2, 0)),
1198  "T4ddg(i,j,k,l)%T2s(i,j)(2,0,0,2)");
1199  test_for_zero(t4ddg_1(2, 0, 1, 0) - (t4ddg_2(2, 0, 1, 0) / t2s_3(2, 0)),
1200  "T4ddg(i,j,k,l)%T2s(i,j)(2,0,1,0)");
1201  test_for_zero(t4ddg_1(2, 0, 1, 1) - (t4ddg_2(2, 0, 1, 1) / t2s_3(2, 0)),
1202  "T4ddg(i,j,k,l)%T2s(i,j)(2,0,1,1)");
1203  test_for_zero(t4ddg_1(2, 0, 1, 2) - (t4ddg_2(2, 0, 1, 2) / t2s_3(2, 0)),
1204  "T4ddg(i,j,k,l)%T2s(i,j)(2,0,1,2)");
1205  test_for_zero(t4ddg_1(2, 0, 2, 0) - (t4ddg_2(2, 0, 2, 0) / t2s_3(2, 0)),
1206  "T4ddg(i,j,k,l)%T2s(i,j)(2,0,2,0)");
1207  test_for_zero(t4ddg_1(2, 0, 2, 1) - (t4ddg_2(2, 0, 2, 1) / t2s_3(2, 0)),
1208  "T4ddg(i,j,k,l)%T2s(i,j)(2,0,2,1)");
1209  test_for_zero(t4ddg_1(2, 0, 2, 2) - (t4ddg_2(2, 0, 2, 2) / t2s_3(2, 0)),
1210  "T4ddg(i,j,k,l)%T2s(i,j)(2,0,2,2)");
1211  test_for_zero(t4ddg_1(2, 1, 0, 0) - (t4ddg_2(2, 1, 0, 0) / t2s_3(2, 1)),
1212  "T4ddg(i,j,k,l)%T2s(i,j)(2,1,0,0)");
1213  test_for_zero(t4ddg_1(2, 1, 0, 1) - (t4ddg_2(2, 1, 0, 1) / t2s_3(2, 1)),
1214  "T4ddg(i,j,k,l)%T2s(i,j)(2,1,0,1)");
1215  test_for_zero(t4ddg_1(2, 1, 0, 2) - (t4ddg_2(2, 1, 0, 2) / t2s_3(2, 1)),
1216  "T4ddg(i,j,k,l)%T2s(i,j)(2,1,0,2)");
1217  test_for_zero(t4ddg_1(2, 1, 1, 0) - (t4ddg_2(2, 1, 1, 0) / t2s_3(2, 1)),
1218  "T4ddg(i,j,k,l)%T2s(i,j)(2,1,1,0)");
1219  test_for_zero(t4ddg_1(2, 1, 1, 1) - (t4ddg_2(2, 1, 1, 1) / t2s_3(2, 1)),
1220  "T4ddg(i,j,k,l)%T2s(i,j)(2,1,1,1)");
1221  test_for_zero(t4ddg_1(2, 1, 1, 2) - (t4ddg_2(2, 1, 1, 2) / t2s_3(2, 1)),
1222  "T4ddg(i,j,k,l)%T2s(i,j)(2,1,1,2)");
1223  test_for_zero(t4ddg_1(2, 1, 2, 0) - (t4ddg_2(2, 1, 2, 0) / t2s_3(2, 1)),
1224  "T4ddg(i,j,k,l)%T2s(i,j)(2,1,2,0)");
1225  test_for_zero(t4ddg_1(2, 1, 2, 1) - (t4ddg_2(2, 1, 2, 1) / t2s_3(2, 1)),
1226  "T4ddg(i,j,k,l)%T2s(i,j)(2,1,2,1)");
1227  test_for_zero(t4ddg_1(2, 1, 2, 2) - (t4ddg_2(2, 1, 2, 2) / t2s_3(2, 1)),
1228  "T4ddg(i,j,k,l)%T2s(i,j)(2,1,2,2)");
1229  test_for_zero(t4ddg_1(2, 2, 0, 0) - (t4ddg_2(2, 2, 0, 0) / t2s_3(2, 2)),
1230  "T4ddg(i,j,k,l)%T2s(i,j)(2,2,0,0)");
1231  test_for_zero(t4ddg_1(2, 2, 0, 1) - (t4ddg_2(2, 2, 0, 1) / t2s_3(2, 2)),
1232  "T4ddg(i,j,k,l)%T2s(i,j)(2,2,0,1)");
1233  test_for_zero(t4ddg_1(2, 2, 0, 2) - (t4ddg_2(2, 2, 0, 2) / t2s_3(2, 2)),
1234  "T4ddg(i,j,k,l)%T2s(i,j)(2,2,0,2)");
1235  test_for_zero(t4ddg_1(2, 2, 1, 0) - (t4ddg_2(2, 2, 1, 0) / t2s_3(2, 2)),
1236  "T4ddg(i,j,k,l)%T2s(i,j)(2,2,1,0)");
1237  test_for_zero(t4ddg_1(2, 2, 1, 1) - (t4ddg_2(2, 2, 1, 1) / t2s_3(2, 2)),
1238  "T4ddg(i,j,k,l)%T2s(i,j)(2,2,1,1)");
1239  test_for_zero(t4ddg_1(2, 2, 1, 2) - (t4ddg_2(2, 2, 1, 2) / t2s_3(2, 2)),
1240  "T4ddg(i,j,k,l)%T2s(i,j)(2,2,1,2)");
1241  test_for_zero(t4ddg_1(2, 2, 2, 0) - (t4ddg_2(2, 2, 2, 0) / t2s_3(2, 2)),
1242  "T4ddg(i,j,k,l)%T2s(i,j)(2,2,2,0)");
1243  test_for_zero(t4ddg_1(2, 2, 2, 1) - (t4ddg_2(2, 2, 2, 1) / t2s_3(2, 2)),
1244  "T4ddg(i,j,k,l)%T2s(i,j)(2,2,2,1)");
1245  test_for_zero(t4ddg_1(2, 2, 2, 2) - (t4ddg_2(2, 2, 2, 2) / t2s_3(2, 2)),
1246  "T4ddg(i,j,k,l)%T2s(i,j)(2,2,2,2)");
1247 
1248  t4ddg_1(i, j, k, l) = (t2s_2(i, j) % t4ddg_2(i, j, k, l));
1249  test_for_zero(t4ddg_1(0, 0, 0, 0) - (t4ddg_2(0, 0, 0, 0) / t2s_2(0, 0)),
1250  "T2s(i,j)%T4ddg(i,j,k,l)(0,0,0,0)");
1251  test_for_zero(t4ddg_1(0, 0, 0, 1) - (t4ddg_2(0, 0, 0, 1) / t2s_2(0, 0)),
1252  "T2s(i,j)%T4ddg(i,j,k,l)(0,0,0,1)");
1253  test_for_zero(t4ddg_1(0, 0, 0, 2) - (t4ddg_2(0, 0, 0, 2) / t2s_2(0, 0)),
1254  "T2s(i,j)%T4ddg(i,j,k,l)(0,0,0,2)");
1255  test_for_zero(t4ddg_1(0, 0, 1, 0) - (t4ddg_2(0, 0, 1, 0) / t2s_2(0, 0)),
1256  "T2s(i,j)%T4ddg(i,j,k,l)(0,0,1,0)");
1257  test_for_zero(t4ddg_1(0, 0, 1, 1) - (t4ddg_2(0, 0, 1, 1) / t2s_2(0, 0)),
1258  "T2s(i,j)%T4ddg(i,j,k,l)(0,0,1,1)");
1259  test_for_zero(t4ddg_1(0, 0, 1, 2) - (t4ddg_2(0, 0, 1, 2) / t2s_2(0, 0)),
1260  "T2s(i,j)%T4ddg(i,j,k,l)(0,0,1,2)");
1261  test_for_zero(t4ddg_1(0, 0, 2, 0) - (t4ddg_2(0, 0, 2, 0) / t2s_2(0, 0)),
1262  "T2s(i,j)%T4ddg(i,j,k,l)(0,0,2,0)");
1263  test_for_zero(t4ddg_1(0, 0, 2, 1) - (t4ddg_2(0, 0, 2, 1) / t2s_2(0, 0)),
1264  "T2s(i,j)%T4ddg(i,j,k,l)(0,0,2,1)");
1265  test_for_zero(t4ddg_1(0, 0, 2, 2) - (t4ddg_2(0, 0, 2, 2) / t2s_2(0, 0)),
1266  "T2s(i,j)%T4ddg(i,j,k,l)(0,0,2,2)");
1267  test_for_zero(t4ddg_1(0, 1, 0, 0) - (t4ddg_2(0, 1, 0, 0) / t2s_2(0, 1)),
1268  "T2s(i,j)%T4ddg(i,j,k,l)(0,1,0,0)");
1269  test_for_zero(t4ddg_1(0, 1, 0, 1) - (t4ddg_2(0, 1, 0, 1) / t2s_2(0, 1)),
1270  "T2s(i,j)%T4ddg(i,j,k,l)(0,1,0,1)");
1271  test_for_zero(t4ddg_1(0, 1, 0, 2) - (t4ddg_2(0, 1, 0, 2) / t2s_2(0, 1)),
1272  "T2s(i,j)%T4ddg(i,j,k,l)(0,1,0,2)");
1273  test_for_zero(t4ddg_1(0, 1, 1, 0) - (t4ddg_2(0, 1, 1, 0) / t2s_2(0, 1)),
1274  "T2s(i,j)%T4ddg(i,j,k,l)(0,1,1,0)");
1275  test_for_zero(t4ddg_1(0, 1, 1, 1) - (t4ddg_2(0, 1, 1, 1) / t2s_2(0, 1)),
1276  "T2s(i,j)%T4ddg(i,j,k,l)(0,1,1,1)");
1277  test_for_zero(t4ddg_1(0, 1, 1, 2) - (t4ddg_2(0, 1, 1, 2) / t2s_2(0, 1)),
1278  "T2s(i,j)%T4ddg(i,j,k,l)(0,1,1,2)");
1279  test_for_zero(t4ddg_1(0, 1, 2, 0) - (t4ddg_2(0, 1, 2, 0) / t2s_2(0, 1)),
1280  "T2s(i,j)%T4ddg(i,j,k,l)(0,1,2,0)");
1281  test_for_zero(t4ddg_1(0, 1, 2, 1) - (t4ddg_2(0, 1, 2, 1) / t2s_2(0, 1)),
1282  "T2s(i,j)%T4ddg(i,j,k,l)(0,1,2,1)");
1283  test_for_zero(t4ddg_1(0, 1, 2, 2) - (t4ddg_2(0, 1, 2, 2) / t2s_2(0, 1)),
1284  "T2s(i,j)%T4ddg(i,j,k,l)(0,1,2,2)");
1285  test_for_zero(t4ddg_1(0, 2, 0, 0) - (t4ddg_2(0, 2, 0, 0) / t2s_2(0, 2)),
1286  "T2s(i,j)%T4ddg(i,j,k,l)(0,2,0,0)");
1287  test_for_zero(t4ddg_1(0, 2, 0, 1) - (t4ddg_2(0, 2, 0, 1) / t2s_2(0, 2)),
1288  "T2s(i,j)%T4ddg(i,j,k,l)(0,2,0,1)");
1289  test_for_zero(t4ddg_1(0, 2, 0, 2) - (t4ddg_2(0, 2, 0, 2) / t2s_2(0, 2)),
1290  "T2s(i,j)%T4ddg(i,j,k,l)(0,2,0,2)");
1291  test_for_zero(t4ddg_1(0, 2, 1, 0) - (t4ddg_2(0, 2, 1, 0) / t2s_2(0, 2)),
1292  "T2s(i,j)%T4ddg(i,j,k,l)(0,2,1,0)");
1293  test_for_zero(t4ddg_1(0, 2, 1, 1) - (t4ddg_2(0, 2, 1, 1) / t2s_2(0, 2)),
1294  "T2s(i,j)%T4ddg(i,j,k,l)(0,2,1,1)");
1295  test_for_zero(t4ddg_1(0, 2, 1, 2) - (t4ddg_2(0, 2, 1, 2) / t2s_2(0, 2)),
1296  "T2s(i,j)%T4ddg(i,j,k,l)(0,2,1,2)");
1297  test_for_zero(t4ddg_1(0, 2, 2, 0) - (t4ddg_2(0, 2, 2, 0) / t2s_2(0, 2)),
1298  "T2s(i,j)%T4ddg(i,j,k,l)(0,2,2,0)");
1299  test_for_zero(t4ddg_1(0, 2, 2, 1) - (t4ddg_2(0, 2, 2, 1) / t2s_2(0, 2)),
1300  "T2s(i,j)%T4ddg(i,j,k,l)(0,2,2,1)");
1301  test_for_zero(t4ddg_1(0, 2, 2, 2) - (t4ddg_2(0, 2, 2, 2) / t2s_2(0, 2)),
1302  "T2s(i,j)%T4ddg(i,j,k,l)(0,2,2,2)");
1303  test_for_zero(t4ddg_1(1, 0, 0, 0) - (t4ddg_2(1, 0, 0, 0) / t2s_2(1, 0)),
1304  "T2s(i,j)%T4ddg(i,j,k,l)(1,0,0,0)");
1305  test_for_zero(t4ddg_1(1, 0, 0, 1) - (t4ddg_2(1, 0, 0, 1) / t2s_2(1, 0)),
1306  "T2s(i,j)%T4ddg(i,j,k,l)(1,0,0,1)");
1307  test_for_zero(t4ddg_1(1, 0, 0, 2) - (t4ddg_2(1, 0, 0, 2) / t2s_2(1, 0)),
1308  "T2s(i,j)%T4ddg(i,j,k,l)(1,0,0,2)");
1309  test_for_zero(t4ddg_1(1, 0, 1, 0) - (t4ddg_2(1, 0, 1, 0) / t2s_2(1, 0)),
1310  "T2s(i,j)%T4ddg(i,j,k,l)(1,0,1,0)");
1311  test_for_zero(t4ddg_1(1, 0, 1, 1) - (t4ddg_2(1, 0, 1, 1) / t2s_2(1, 0)),
1312  "T2s(i,j)%T4ddg(i,j,k,l)(1,0,1,1)");
1313  test_for_zero(t4ddg_1(1, 0, 1, 2) - (t4ddg_2(1, 0, 1, 2) / t2s_2(1, 0)),
1314  "T2s(i,j)%T4ddg(i,j,k,l)(1,0,1,2)");
1315  test_for_zero(t4ddg_1(1, 0, 2, 0) - (t4ddg_2(1, 0, 2, 0) / t2s_2(1, 0)),
1316  "T2s(i,j)%T4ddg(i,j,k,l)(1,0,2,0)");
1317  test_for_zero(t4ddg_1(1, 0, 2, 1) - (t4ddg_2(1, 0, 2, 1) / t2s_2(1, 0)),
1318  "T2s(i,j)%T4ddg(i,j,k,l)(1,0,2,1)");
1319  test_for_zero(t4ddg_1(1, 0, 2, 2) - (t4ddg_2(1, 0, 2, 2) / t2s_2(1, 0)),
1320  "T2s(i,j)%T4ddg(i,j,k,l)(1,0,2,2)");
1321  test_for_zero(t4ddg_1(1, 1, 0, 0) - (t4ddg_2(1, 1, 0, 0) / t2s_2(1, 1)),
1322  "T2s(i,j)%T4ddg(i,j,k,l)(1,1,0,0)");
1323  test_for_zero(t4ddg_1(1, 1, 0, 1) - (t4ddg_2(1, 1, 0, 1) / t2s_2(1, 1)),
1324  "T2s(i,j)%T4ddg(i,j,k,l)(1,1,0,1)");
1325  test_for_zero(t4ddg_1(1, 1, 0, 2) - (t4ddg_2(1, 1, 0, 2) / t2s_2(1, 1)),
1326  "T2s(i,j)%T4ddg(i,j,k,l)(1,1,0,2)");
1327  test_for_zero(t4ddg_1(1, 1, 1, 0) - (t4ddg_2(1, 1, 1, 0) / t2s_2(1, 1)),
1328  "T2s(i,j)%T4ddg(i,j,k,l)(1,1,1,0)");
1329  test_for_zero(t4ddg_1(1, 1, 1, 1) - (t4ddg_2(1, 1, 1, 1) / t2s_2(1, 1)),
1330  "T2s(i,j)%T4ddg(i,j,k,l)(1,1,1,1)");
1331  test_for_zero(t4ddg_1(1, 1, 1, 2) - (t4ddg_2(1, 1, 1, 2) / t2s_2(1, 1)),
1332  "T2s(i,j)%T4ddg(i,j,k,l)(1,1,1,2)");
1333  test_for_zero(t4ddg_1(1, 1, 2, 0) - (t4ddg_2(1, 1, 2, 0) / t2s_2(1, 1)),
1334  "T2s(i,j)%T4ddg(i,j,k,l)(1,1,2,0)");
1335  test_for_zero(t4ddg_1(1, 1, 2, 1) - (t4ddg_2(1, 1, 2, 1) / t2s_2(1, 1)),
1336  "T2s(i,j)%T4ddg(i,j,k,l)(1,1,2,1)");
1337  test_for_zero(t4ddg_1(1, 1, 2, 2) - (t4ddg_2(1, 1, 2, 2) / t2s_2(1, 1)),
1338  "T2s(i,j)%T4ddg(i,j,k,l)(1,1,2,2)");
1339  test_for_zero(t4ddg_1(1, 2, 0, 0) - (t4ddg_2(1, 2, 0, 0) / t2s_2(1, 2)),
1340  "T2s(i,j)%T4ddg(i,j,k,l)(1,2,0,0)");
1341  test_for_zero(t4ddg_1(1, 2, 0, 1) - (t4ddg_2(1, 2, 0, 1) / t2s_2(1, 2)),
1342  "T2s(i,j)%T4ddg(i,j,k,l)(1,2,0,1)");
1343  test_for_zero(t4ddg_1(1, 2, 0, 2) - (t4ddg_2(1, 2, 0, 2) / t2s_2(1, 2)),
1344  "T2s(i,j)%T4ddg(i,j,k,l)(1,2,0,2)");
1345  test_for_zero(t4ddg_1(1, 2, 1, 0) - (t4ddg_2(1, 2, 1, 0) / t2s_2(1, 2)),
1346  "T2s(i,j)%T4ddg(i,j,k,l)(1,2,1,0)");
1347  test_for_zero(t4ddg_1(1, 2, 1, 1) - (t4ddg_2(1, 2, 1, 1) / t2s_2(1, 2)),
1348  "T2s(i,j)%T4ddg(i,j,k,l)(1,2,1,1)");
1349  test_for_zero(t4ddg_1(1, 2, 1, 2) - (t4ddg_2(1, 2, 1, 2) / t2s_2(1, 2)),
1350  "T2s(i,j)%T4ddg(i,j,k,l)(1,2,1,2)");
1351  test_for_zero(t4ddg_1(1, 2, 2, 0) - (t4ddg_2(1, 2, 2, 0) / t2s_2(1, 2)),
1352  "T2s(i,j)%T4ddg(i,j,k,l)(1,2,2,0)");
1353  test_for_zero(t4ddg_1(1, 2, 2, 1) - (t4ddg_2(1, 2, 2, 1) / t2s_2(1, 2)),
1354  "T2s(i,j)%T4ddg(i,j,k,l)(1,2,2,1)");
1355  test_for_zero(t4ddg_1(1, 2, 2, 2) - (t4ddg_2(1, 2, 2, 2) / t2s_2(1, 2)),
1356  "T2s(i,j)%T4ddg(i,j,k,l)(1,2,2,2)");
1357  test_for_zero(t4ddg_1(2, 0, 0, 0) - (t4ddg_2(2, 0, 0, 0) / t2s_2(2, 0)),
1358  "T2s(i,j)%T4ddg(i,j,k,l)(2,0,0,0)");
1359  test_for_zero(t4ddg_1(2, 0, 0, 1) - (t4ddg_2(2, 0, 0, 1) / t2s_2(2, 0)),
1360  "T2s(i,j)%T4ddg(i,j,k,l)(2,0,0,1)");
1361  test_for_zero(t4ddg_1(2, 0, 0, 2) - (t4ddg_2(2, 0, 0, 2) / t2s_2(2, 0)),
1362  "T2s(i,j)%T4ddg(i,j,k,l)(2,0,0,2)");
1363  test_for_zero(t4ddg_1(2, 0, 1, 0) - (t4ddg_2(2, 0, 1, 0) / t2s_2(2, 0)),
1364  "T2s(i,j)%T4ddg(i,j,k,l)(2,0,1,0)");
1365  test_for_zero(t4ddg_1(2, 0, 1, 1) - (t4ddg_2(2, 0, 1, 1) / t2s_2(2, 0)),
1366  "T2s(i,j)%T4ddg(i,j,k,l)(2,0,1,1)");
1367  test_for_zero(t4ddg_1(2, 0, 1, 2) - (t4ddg_2(2, 0, 1, 2) / t2s_2(2, 0)),
1368  "T2s(i,j)%T4ddg(i,j,k,l)(2,0,1,2)");
1369  test_for_zero(t4ddg_1(2, 0, 2, 0) - (t4ddg_2(2, 0, 2, 0) / t2s_2(2, 0)),
1370  "T2s(i,j)%T4ddg(i,j,k,l)(2,0,2,0)");
1371  test_for_zero(t4ddg_1(2, 0, 2, 1) - (t4ddg_2(2, 0, 2, 1) / t2s_2(2, 0)),
1372  "T2s(i,j)%T4ddg(i,j,k,l)(2,0,2,1)");
1373  test_for_zero(t4ddg_1(2, 0, 2, 2) - (t4ddg_2(2, 0, 2, 2) / t2s_2(2, 0)),
1374  "T2s(i,j)%T4ddg(i,j,k,l)(2,0,2,2)");
1375  test_for_zero(t4ddg_1(2, 1, 0, 0) - (t4ddg_2(2, 1, 0, 0) / t2s_2(2, 1)),
1376  "T2s(i,j)%T4ddg(i,j,k,l)(2,1,0,0)");
1377  test_for_zero(t4ddg_1(2, 1, 0, 1) - (t4ddg_2(2, 1, 0, 1) / t2s_2(2, 1)),
1378  "T2s(i,j)%T4ddg(i,j,k,l)(2,1,0,1)");
1379  test_for_zero(t4ddg_1(2, 1, 0, 2) - (t4ddg_2(2, 1, 0, 2) / t2s_2(2, 1)),
1380  "T2s(i,j)%T4ddg(i,j,k,l)(2,1,0,2)");
1381  test_for_zero(t4ddg_1(2, 1, 1, 0) - (t4ddg_2(2, 1, 1, 0) / t2s_2(2, 1)),
1382  "T2s(i,j)%T4ddg(i,j,k,l)(2,1,1,0)");
1383  test_for_zero(t4ddg_1(2, 1, 1, 1) - (t4ddg_2(2, 1, 1, 1) / t2s_2(2, 1)),
1384  "T2s(i,j)%T4ddg(i,j,k,l)(2,1,1,1)");
1385  test_for_zero(t4ddg_1(2, 1, 1, 2) - (t4ddg_2(2, 1, 1, 2) / t2s_2(2, 1)),
1386  "T2s(i,j)%T4ddg(i,j,k,l)(2,1,1,2)");
1387  test_for_zero(t4ddg_1(2, 1, 2, 0) - (t4ddg_2(2, 1, 2, 0) / t2s_2(2, 1)),
1388  "T2s(i,j)%T4ddg(i,j,k,l)(2,1,2,0)");
1389  test_for_zero(t4ddg_1(2, 1, 2, 1) - (t4ddg_2(2, 1, 2, 1) / t2s_2(2, 1)),
1390  "T2s(i,j)%T4ddg(i,j,k,l)(2,1,2,1)");
1391  test_for_zero(t4ddg_1(2, 1, 2, 2) - (t4ddg_2(2, 1, 2, 2) / t2s_2(2, 1)),
1392  "T2s(i,j)%T4ddg(i,j,k,l)(2,1,2,2)");
1393  test_for_zero(t4ddg_1(2, 2, 0, 0) - (t4ddg_2(2, 2, 0, 0) / t2s_2(2, 2)),
1394  "T2s(i,j)%T4ddg(i,j,k,l)(2,2,0,0)");
1395  test_for_zero(t4ddg_1(2, 2, 0, 1) - (t4ddg_2(2, 2, 0, 1) / t2s_2(2, 2)),
1396  "T2s(i,j)%T4ddg(i,j,k,l)(2,2,0,1)");
1397  test_for_zero(t4ddg_1(2, 2, 0, 2) - (t4ddg_2(2, 2, 0, 2) / t2s_2(2, 2)),
1398  "T2s(i,j)%T4ddg(i,j,k,l)(2,2,0,2)");
1399  test_for_zero(t4ddg_1(2, 2, 1, 0) - (t4ddg_2(2, 2, 1, 0) / t2s_2(2, 2)),
1400  "T2s(i,j)%T4ddg(i,j,k,l)(2,2,1,0)");
1401  test_for_zero(t4ddg_1(2, 2, 1, 1) - (t4ddg_2(2, 2, 1, 1) / t2s_2(2, 2)),
1402  "T2s(i,j)%T4ddg(i,j,k,l)(2,2,1,1)");
1403  test_for_zero(t4ddg_1(2, 2, 1, 2) - (t4ddg_2(2, 2, 1, 2) / t2s_2(2, 2)),
1404  "T2s(i,j)%T4ddg(i,j,k,l)(2,2,1,2)");
1405  test_for_zero(t4ddg_1(2, 2, 2, 0) - (t4ddg_2(2, 2, 2, 0) / t2s_2(2, 2)),
1406  "T2s(i,j)%T4ddg(i,j,k,l)(2,2,2,0)");
1407  test_for_zero(t4ddg_1(2, 2, 2, 1) - (t4ddg_2(2, 2, 2, 1) / t2s_2(2, 2)),
1408  "T2s(i,j)%T4ddg(i,j,k,l)(2,2,2,1)");
1409  test_for_zero(t4ddg_1(2, 2, 2, 2) - (t4ddg_2(2, 2, 2, 2) / t2s_2(2, 2)),
1410  "T2s(i,j)%T4ddg(i,j,k,l)(2,2,2,2)");
1411 
1412  t4ddg_1(i, j, k, l) = (t4ddg_2(i, j, k, l) % t2s_3(k, l));
1413  test_for_zero(t4ddg_1(0, 0, 0, 0) - (t4ddg_2(0, 0, 0, 0) / t2s_3(0, 0)),
1414  "T4ddg(i,j,k,l)%T2s(k,l)(0,0,0,0)");
1415  test_for_zero(t4ddg_1(0, 0, 0, 1) - (t4ddg_2(0, 0, 0, 1) / t2s_3(0, 1)),
1416  "T4ddg(i,j,k,l)%T2s(k,l)(0,0,0,1)");
1417  test_for_zero(t4ddg_1(0, 0, 0, 2) - (t4ddg_2(0, 0, 0, 2) / t2s_3(0, 2)),
1418  "T4ddg(i,j,k,l)%T2s(k,l)(0,0,0,2)");
1419  test_for_zero(t4ddg_1(0, 0, 1, 0) - (t4ddg_2(0, 0, 1, 0) / t2s_3(1, 0)),
1420  "T4ddg(i,j,k,l)%T2s(k,l)(0,0,1,0)");
1421  test_for_zero(t4ddg_1(0, 0, 1, 1) - (t4ddg_2(0, 0, 1, 1) / t2s_3(1, 1)),
1422  "T4ddg(i,j,k,l)%T2s(k,l)(0,0,1,1)");
1423  test_for_zero(t4ddg_1(0, 0, 1, 2) - (t4ddg_2(0, 0, 1, 2) / t2s_3(1, 2)),
1424  "T4ddg(i,j,k,l)%T2s(k,l)(0,0,1,2)");
1425  test_for_zero(t4ddg_1(0, 0, 2, 0) - (t4ddg_2(0, 0, 2, 0) / t2s_3(2, 0)),
1426  "T4ddg(i,j,k,l)%T2s(k,l)(0,0,2,0)");
1427  test_for_zero(t4ddg_1(0, 0, 2, 1) - (t4ddg_2(0, 0, 2, 1) / t2s_3(2, 1)),
1428  "T4ddg(i,j,k,l)%T2s(k,l)(0,0,2,1)");
1429  test_for_zero(t4ddg_1(0, 0, 2, 2) - (t4ddg_2(0, 0, 2, 2) / t2s_3(2, 2)),
1430  "T4ddg(i,j,k,l)%T2s(k,l)(0,0,2,2)");
1431  test_for_zero(t4ddg_1(0, 1, 0, 0) - (t4ddg_2(0, 1, 0, 0) / t2s_3(0, 0)),
1432  "T4ddg(i,j,k,l)%T2s(k,l)(0,1,0,0)");
1433  test_for_zero(t4ddg_1(0, 1, 0, 1) - (t4ddg_2(0, 1, 0, 1) / t2s_3(0, 1)),
1434  "T4ddg(i,j,k,l)%T2s(k,l)(0,1,0,1)");
1435  test_for_zero(t4ddg_1(0, 1, 0, 2) - (t4ddg_2(0, 1, 0, 2) / t2s_3(0, 2)),
1436  "T4ddg(i,j,k,l)%T2s(k,l)(0,1,0,2)");
1437  test_for_zero(t4ddg_1(0, 1, 1, 0) - (t4ddg_2(0, 1, 1, 0) / t2s_3(1, 0)),
1438  "T4ddg(i,j,k,l)%T2s(k,l)(0,1,1,0)");
1439  test_for_zero(t4ddg_1(0, 1, 1, 1) - (t4ddg_2(0, 1, 1, 1) / t2s_3(1, 1)),
1440  "T4ddg(i,j,k,l)%T2s(k,l)(0,1,1,1)");
1441  test_for_zero(t4ddg_1(0, 1, 1, 2) - (t4ddg_2(0, 1, 1, 2) / t2s_3(1, 2)),
1442  "T4ddg(i,j,k,l)%T2s(k,l)(0,1,1,2)");
1443  test_for_zero(t4ddg_1(0, 1, 2, 0) - (t4ddg_2(0, 1, 2, 0) / t2s_3(2, 0)),
1444  "T4ddg(i,j,k,l)%T2s(k,l)(0,1,2,0)");
1445  test_for_zero(t4ddg_1(0, 1, 2, 1) - (t4ddg_2(0, 1, 2, 1) / t2s_3(2, 1)),
1446  "T4ddg(i,j,k,l)%T2s(k,l)(0,1,2,1)");
1447  test_for_zero(t4ddg_1(0, 1, 2, 2) - (t4ddg_2(0, 1, 2, 2) / t2s_3(2, 2)),
1448  "T4ddg(i,j,k,l)%T2s(k,l)(0,1,2,2)");
1449  test_for_zero(t4ddg_1(0, 2, 0, 0) - (t4ddg_2(0, 2, 0, 0) / t2s_3(0, 0)),
1450  "T4ddg(i,j,k,l)%T2s(k,l)(0,2,0,0)");
1451  test_for_zero(t4ddg_1(0, 2, 0, 1) - (t4ddg_2(0, 2, 0, 1) / t2s_3(0, 1)),
1452  "T4ddg(i,j,k,l)%T2s(k,l)(0,2,0,1)");
1453  test_for_zero(t4ddg_1(0, 2, 0, 2) - (t4ddg_2(0, 2, 0, 2) / t2s_3(0, 2)),
1454  "T4ddg(i,j,k,l)%T2s(k,l)(0,2,0,2)");
1455  test_for_zero(t4ddg_1(0, 2, 1, 0) - (t4ddg_2(0, 2, 1, 0) / t2s_3(1, 0)),
1456  "T4ddg(i,j,k,l)%T2s(k,l)(0,2,1,0)");
1457  test_for_zero(t4ddg_1(0, 2, 1, 1) - (t4ddg_2(0, 2, 1, 1) / t2s_3(1, 1)),
1458  "T4ddg(i,j,k,l)%T2s(k,l)(0,2,1,1)");
1459  test_for_zero(t4ddg_1(0, 2, 1, 2) - (t4ddg_2(0, 2, 1, 2) / t2s_3(1, 2)),
1460  "T4ddg(i,j,k,l)%T2s(k,l)(0,2,1,2)");
1461  test_for_zero(t4ddg_1(0, 2, 2, 0) - (t4ddg_2(0, 2, 2, 0) / t2s_3(2, 0)),
1462  "T4ddg(i,j,k,l)%T2s(k,l)(0,2,2,0)");
1463  test_for_zero(t4ddg_1(0, 2, 2, 1) - (t4ddg_2(0, 2, 2, 1) / t2s_3(2, 1)),
1464  "T4ddg(i,j,k,l)%T2s(k,l)(0,2,2,1)");
1465  test_for_zero(t4ddg_1(0, 2, 2, 2) - (t4ddg_2(0, 2, 2, 2) / t2s_3(2, 2)),
1466  "T4ddg(i,j,k,l)%T2s(k,l)(0,2,2,2)");
1467  test_for_zero(t4ddg_1(1, 0, 0, 0) - (t4ddg_2(1, 0, 0, 0) / t2s_3(0, 0)),
1468  "T4ddg(i,j,k,l)%T2s(k,l)(1,0,0,0)");
1469  test_for_zero(t4ddg_1(1, 0, 0, 1) - (t4ddg_2(1, 0, 0, 1) / t2s_3(0, 1)),
1470  "T4ddg(i,j,k,l)%T2s(k,l)(1,0,0,1)");
1471  test_for_zero(t4ddg_1(1, 0, 0, 2) - (t4ddg_2(1, 0, 0, 2) / t2s_3(0, 2)),
1472  "T4ddg(i,j,k,l)%T2s(k,l)(1,0,0,2)");
1473  test_for_zero(t4ddg_1(1, 0, 1, 0) - (t4ddg_2(1, 0, 1, 0) / t2s_3(1, 0)),
1474  "T4ddg(i,j,k,l)%T2s(k,l)(1,0,1,0)");
1475  test_for_zero(t4ddg_1(1, 0, 1, 1) - (t4ddg_2(1, 0, 1, 1) / t2s_3(1, 1)),
1476  "T4ddg(i,j,k,l)%T2s(k,l)(1,0,1,1)");
1477  test_for_zero(t4ddg_1(1, 0, 1, 2) - (t4ddg_2(1, 0, 1, 2) / t2s_3(1, 2)),
1478  "T4ddg(i,j,k,l)%T2s(k,l)(1,0,1,2)");
1479  test_for_zero(t4ddg_1(1, 0, 2, 0) - (t4ddg_2(1, 0, 2, 0) / t2s_3(2, 0)),
1480  "T4ddg(i,j,k,l)%T2s(k,l)(1,0,2,0)");
1481  test_for_zero(t4ddg_1(1, 0, 2, 1) - (t4ddg_2(1, 0, 2, 1) / t2s_3(2, 1)),
1482  "T4ddg(i,j,k,l)%T2s(k,l)(1,0,2,1)");
1483  test_for_zero(t4ddg_1(1, 0, 2, 2) - (t4ddg_2(1, 0, 2, 2) / t2s_3(2, 2)),
1484  "T4ddg(i,j,k,l)%T2s(k,l)(1,0,2,2)");
1485  test_for_zero(t4ddg_1(1, 1, 0, 0) - (t4ddg_2(1, 1, 0, 0) / t2s_3(0, 0)),
1486  "T4ddg(i,j,k,l)%T2s(k,l)(1,1,0,0)");
1487  test_for_zero(t4ddg_1(1, 1, 0, 1) - (t4ddg_2(1, 1, 0, 1) / t2s_3(0, 1)),
1488  "T4ddg(i,j,k,l)%T2s(k,l)(1,1,0,1)");
1489  test_for_zero(t4ddg_1(1, 1, 0, 2) - (t4ddg_2(1, 1, 0, 2) / t2s_3(0, 2)),
1490  "T4ddg(i,j,k,l)%T2s(k,l)(1,1,0,2)");
1491  test_for_zero(t4ddg_1(1, 1, 1, 0) - (t4ddg_2(1, 1, 1, 0) / t2s_3(1, 0)),
1492  "T4ddg(i,j,k,l)%T2s(k,l)(1,1,1,0)");
1493  test_for_zero(t4ddg_1(1, 1, 1, 1) - (t4ddg_2(1, 1, 1, 1) / t2s_3(1, 1)),
1494  "T4ddg(i,j,k,l)%T2s(k,l)(1,1,1,1)");
1495  test_for_zero(t4ddg_1(1, 1, 1, 2) - (t4ddg_2(1, 1, 1, 2) / t2s_3(1, 2)),
1496  "T4ddg(i,j,k,l)%T2s(k,l)(1,1,1,2)");
1497  test_for_zero(t4ddg_1(1, 1, 2, 0) - (t4ddg_2(1, 1, 2, 0) / t2s_3(2, 0)),
1498  "T4ddg(i,j,k,l)%T2s(k,l)(1,1,2,0)");
1499  test_for_zero(t4ddg_1(1, 1, 2, 1) - (t4ddg_2(1, 1, 2, 1) / t2s_3(2, 1)),
1500  "T4ddg(i,j,k,l)%T2s(k,l)(1,1,2,1)");
1501  test_for_zero(t4ddg_1(1, 1, 2, 2) - (t4ddg_2(1, 1, 2, 2) / t2s_3(2, 2)),
1502  "T4ddg(i,j,k,l)%T2s(k,l)(1,1,2,2)");
1503  test_for_zero(t4ddg_1(1, 2, 0, 0) - (t4ddg_2(1, 2, 0, 0) / t2s_3(0, 0)),
1504  "T4ddg(i,j,k,l)%T2s(k,l)(1,2,0,0)");
1505  test_for_zero(t4ddg_1(1, 2, 0, 1) - (t4ddg_2(1, 2, 0, 1) / t2s_3(0, 1)),
1506  "T4ddg(i,j,k,l)%T2s(k,l)(1,2,0,1)");
1507  test_for_zero(t4ddg_1(1, 2, 0, 2) - (t4ddg_2(1, 2, 0, 2) / t2s_3(0, 2)),
1508  "T4ddg(i,j,k,l)%T2s(k,l)(1,2,0,2)");
1509  test_for_zero(t4ddg_1(1, 2, 1, 0) - (t4ddg_2(1, 2, 1, 0) / t2s_3(1, 0)),
1510  "T4ddg(i,j,k,l)%T2s(k,l)(1,2,1,0)");
1511  test_for_zero(t4ddg_1(1, 2, 1, 1) - (t4ddg_2(1, 2, 1, 1) / t2s_3(1, 1)),
1512  "T4ddg(i,j,k,l)%T2s(k,l)(1,2,1,1)");
1513  test_for_zero(t4ddg_1(1, 2, 1, 2) - (t4ddg_2(1, 2, 1, 2) / t2s_3(1, 2)),
1514  "T4ddg(i,j,k,l)%T2s(k,l)(1,2,1,2)");
1515  test_for_zero(t4ddg_1(1, 2, 2, 0) - (t4ddg_2(1, 2, 2, 0) / t2s_3(2, 0)),
1516  "T4ddg(i,j,k,l)%T2s(k,l)(1,2,2,0)");
1517  test_for_zero(t4ddg_1(1, 2, 2, 1) - (t4ddg_2(1, 2, 2, 1) / t2s_3(2, 1)),
1518  "T4ddg(i,j,k,l)%T2s(k,l)(1,2,2,1)");
1519  test_for_zero(t4ddg_1(1, 2, 2, 2) - (t4ddg_2(1, 2, 2, 2) / t2s_3(2, 2)),
1520  "T4ddg(i,j,k,l)%T2s(k,l)(1,2,2,2)");
1521  test_for_zero(t4ddg_1(2, 0, 0, 0) - (t4ddg_2(2, 0, 0, 0) / t2s_3(0, 0)),
1522  "T4ddg(i,j,k,l)%T2s(k,l)(2,0,0,0)");
1523  test_for_zero(t4ddg_1(2, 0, 0, 1) - (t4ddg_2(2, 0, 0, 1) / t2s_3(0, 1)),
1524  "T4ddg(i,j,k,l)%T2s(k,l)(2,0,0,1)");
1525  test_for_zero(t4ddg_1(2, 0, 0, 2) - (t4ddg_2(2, 0, 0, 2) / t2s_3(0, 2)),
1526  "T4ddg(i,j,k,l)%T2s(k,l)(2,0,0,2)");
1527  test_for_zero(t4ddg_1(2, 0, 1, 0) - (t4ddg_2(2, 0, 1, 0) / t2s_3(1, 0)),
1528  "T4ddg(i,j,k,l)%T2s(k,l)(2,0,1,0)");
1529  test_for_zero(t4ddg_1(2, 0, 1, 1) - (t4ddg_2(2, 0, 1, 1) / t2s_3(1, 1)),
1530  "T4ddg(i,j,k,l)%T2s(k,l)(2,0,1,1)");
1531  test_for_zero(t4ddg_1(2, 0, 1, 2) - (t4ddg_2(2, 0, 1, 2) / t2s_3(1, 2)),
1532  "T4ddg(i,j,k,l)%T2s(k,l)(2,0,1,2)");
1533  test_for_zero(t4ddg_1(2, 0, 2, 0) - (t4ddg_2(2, 0, 2, 0) / t2s_3(2, 0)),
1534  "T4ddg(i,j,k,l)%T2s(k,l)(2,0,2,0)");
1535  test_for_zero(t4ddg_1(2, 0, 2, 1) - (t4ddg_2(2, 0, 2, 1) / t2s_3(2, 1)),
1536  "T4ddg(i,j,k,l)%T2s(k,l)(2,0,2,1)");
1537  test_for_zero(t4ddg_1(2, 0, 2, 2) - (t4ddg_2(2, 0, 2, 2) / t2s_3(2, 2)),
1538  "T4ddg(i,j,k,l)%T2s(k,l)(2,0,2,2)");
1539  test_for_zero(t4ddg_1(2, 1, 0, 0) - (t4ddg_2(2, 1, 0, 0) / t2s_3(0, 0)),
1540  "T4ddg(i,j,k,l)%T2s(k,l)(2,1,0,0)");
1541  test_for_zero(t4ddg_1(2, 1, 0, 1) - (t4ddg_2(2, 1, 0, 1) / t2s_3(0, 1)),
1542  "T4ddg(i,j,k,l)%T2s(k,l)(2,1,0,1)");
1543  test_for_zero(t4ddg_1(2, 1, 0, 2) - (t4ddg_2(2, 1, 0, 2) / t2s_3(0, 2)),
1544  "T4ddg(i,j,k,l)%T2s(k,l)(2,1,0,2)");
1545  test_for_zero(t4ddg_1(2, 1, 1, 0) - (t4ddg_2(2, 1, 1, 0) / t2s_3(1, 0)),
1546  "T4ddg(i,j,k,l)%T2s(k,l)(2,1,1,0)");
1547  test_for_zero(t4ddg_1(2, 1, 1, 1) - (t4ddg_2(2, 1, 1, 1) / t2s_3(1, 1)),
1548  "T4ddg(i,j,k,l)%T2s(k,l)(2,1,1,1)");
1549  test_for_zero(t4ddg_1(2, 1, 1, 2) - (t4ddg_2(2, 1, 1, 2) / t2s_3(1, 2)),
1550  "T4ddg(i,j,k,l)%T2s(k,l)(2,1,1,2)");
1551  test_for_zero(t4ddg_1(2, 1, 2, 0) - (t4ddg_2(2, 1, 2, 0) / t2s_3(2, 0)),
1552  "T4ddg(i,j,k,l)%T2s(k,l)(2,1,2,0)");
1553  test_for_zero(t4ddg_1(2, 1, 2, 1) - (t4ddg_2(2, 1, 2, 1) / t2s_3(2, 1)),
1554  "T4ddg(i,j,k,l)%T2s(k,l)(2,1,2,1)");
1555  test_for_zero(t4ddg_1(2, 1, 2, 2) - (t4ddg_2(2, 1, 2, 2) / t2s_3(2, 2)),
1556  "T4ddg(i,j,k,l)%T2s(k,l)(2,1,2,2)");
1557  test_for_zero(t4ddg_1(2, 2, 0, 0) - (t4ddg_2(2, 2, 0, 0) / t2s_3(0, 0)),
1558  "T4ddg(i,j,k,l)%T2s(k,l)(2,2,0,0)");
1559  test_for_zero(t4ddg_1(2, 2, 0, 1) - (t4ddg_2(2, 2, 0, 1) / t2s_3(0, 1)),
1560  "T4ddg(i,j,k,l)%T2s(k,l)(2,2,0,1)");
1561  test_for_zero(t4ddg_1(2, 2, 0, 2) - (t4ddg_2(2, 2, 0, 2) / t2s_3(0, 2)),
1562  "T4ddg(i,j,k,l)%T2s(k,l)(2,2,0,2)");
1563  test_for_zero(t4ddg_1(2, 2, 1, 0) - (t4ddg_2(2, 2, 1, 0) / t2s_3(1, 0)),
1564  "T4ddg(i,j,k,l)%T2s(k,l)(2,2,1,0)");
1565  test_for_zero(t4ddg_1(2, 2, 1, 1) - (t4ddg_2(2, 2, 1, 1) / t2s_3(1, 1)),
1566  "T4ddg(i,j,k,l)%T2s(k,l)(2,2,1,1)");
1567  test_for_zero(t4ddg_1(2, 2, 1, 2) - (t4ddg_2(2, 2, 1, 2) / t2s_3(1, 2)),
1568  "T4ddg(i,j,k,l)%T2s(k,l)(2,2,1,2)");
1569  test_for_zero(t4ddg_1(2, 2, 2, 0) - (t4ddg_2(2, 2, 2, 0) / t2s_3(2, 0)),
1570  "T4ddg(i,j,k,l)%T2s(k,l)(2,2,2,0)");
1571  test_for_zero(t4ddg_1(2, 2, 2, 1) - (t4ddg_2(2, 2, 2, 1) / t2s_3(2, 1)),
1572  "T4ddg(i,j,k,l)%T2s(k,l)(2,2,2,1)");
1573  test_for_zero(t4ddg_1(2, 2, 2, 2) - (t4ddg_2(2, 2, 2, 2) / t2s_3(2, 2)),
1574  "T4ddg(i,j,k,l)%T2s(k,l)(2,2,2,2)");
1575 
1576  t4ddg_1(i, j, k, l) = (t2s_2(k, l) % t4ddg_2(i, j, k, l));
1577  test_for_zero(t4ddg_1(0, 0, 0, 0) - (t4ddg_2(0, 0, 0, 0) / t2s_2(0, 0)),
1578  "T2s(k,l)%T4ddg(i,j,k,l)(0,0,0,0)");
1579  test_for_zero(t4ddg_1(0, 0, 0, 1) - (t4ddg_2(0, 0, 0, 1) / t2s_2(0, 1)),
1580  "T2s(k,l)%T4ddg(i,j,k,l)(0,0,0,1)");
1581  test_for_zero(t4ddg_1(0, 0, 0, 2) - (t4ddg_2(0, 0, 0, 2) / t2s_2(0, 2)),
1582  "T2s(k,l)%T4ddg(i,j,k,l)(0,0,0,2)");
1583  test_for_zero(t4ddg_1(0, 0, 1, 0) - (t4ddg_2(0, 0, 1, 0) / t2s_2(1, 0)),
1584  "T2s(k,l)%T4ddg(i,j,k,l)(0,0,1,0)");
1585  test_for_zero(t4ddg_1(0, 0, 1, 1) - (t4ddg_2(0, 0, 1, 1) / t2s_2(1, 1)),
1586  "T2s(k,l)%T4ddg(i,j,k,l)(0,0,1,1)");
1587  test_for_zero(t4ddg_1(0, 0, 1, 2) - (t4ddg_2(0, 0, 1, 2) / t2s_2(1, 2)),
1588  "T2s(k,l)%T4ddg(i,j,k,l)(0,0,1,2)");
1589  test_for_zero(t4ddg_1(0, 0, 2, 0) - (t4ddg_2(0, 0, 2, 0) / t2s_2(2, 0)),
1590  "T2s(k,l)%T4ddg(i,j,k,l)(0,0,2,0)");
1591  test_for_zero(t4ddg_1(0, 0, 2, 1) - (t4ddg_2(0, 0, 2, 1) / t2s_2(2, 1)),
1592  "T2s(k,l)%T4ddg(i,j,k,l)(0,0,2,1)");
1593  test_for_zero(t4ddg_1(0, 0, 2, 2) - (t4ddg_2(0, 0, 2, 2) / t2s_2(2, 2)),
1594  "T2s(k,l)%T4ddg(i,j,k,l)(0,0,2,2)");
1595  test_for_zero(t4ddg_1(0, 1, 0, 0) - (t4ddg_2(0, 1, 0, 0) / t2s_2(0, 0)),
1596  "T2s(k,l)%T4ddg(i,j,k,l)(0,1,0,0)");
1597  test_for_zero(t4ddg_1(0, 1, 0, 1) - (t4ddg_2(0, 1, 0, 1) / t2s_2(0, 1)),
1598  "T2s(k,l)%T4ddg(i,j,k,l)(0,1,0,1)");
1599  test_for_zero(t4ddg_1(0, 1, 0, 2) - (t4ddg_2(0, 1, 0, 2) / t2s_2(0, 2)),
1600  "T2s(k,l)%T4ddg(i,j,k,l)(0,1,0,2)");
1601  test_for_zero(t4ddg_1(0, 1, 1, 0) - (t4ddg_2(0, 1, 1, 0) / t2s_2(1, 0)),
1602  "T2s(k,l)%T4ddg(i,j,k,l)(0,1,1,0)");
1603  test_for_zero(t4ddg_1(0, 1, 1, 1) - (t4ddg_2(0, 1, 1, 1) / t2s_2(1, 1)),
1604  "T2s(k,l)%T4ddg(i,j,k,l)(0,1,1,1)");
1605  test_for_zero(t4ddg_1(0, 1, 1, 2) - (t4ddg_2(0, 1, 1, 2) / t2s_2(1, 2)),
1606  "T2s(k,l)%T4ddg(i,j,k,l)(0,1,1,2)");
1607  test_for_zero(t4ddg_1(0, 1, 2, 0) - (t4ddg_2(0, 1, 2, 0) / t2s_2(2, 0)),
1608  "T2s(k,l)%T4ddg(i,j,k,l)(0,1,2,0)");
1609  test_for_zero(t4ddg_1(0, 1, 2, 1) - (t4ddg_2(0, 1, 2, 1) / t2s_2(2, 1)),
1610  "T2s(k,l)%T4ddg(i,j,k,l)(0,1,2,1)");
1611  test_for_zero(t4ddg_1(0, 1, 2, 2) - (t4ddg_2(0, 1, 2, 2) / t2s_2(2, 2)),
1612  "T2s(k,l)%T4ddg(i,j,k,l)(0,1,2,2)");
1613  test_for_zero(t4ddg_1(0, 2, 0, 0) - (t4ddg_2(0, 2, 0, 0) / t2s_2(0, 0)),
1614  "T2s(k,l)%T4ddg(i,j,k,l)(0,2,0,0)");
1615  test_for_zero(t4ddg_1(0, 2, 0, 1) - (t4ddg_2(0, 2, 0, 1) / t2s_2(0, 1)),
1616  "T2s(k,l)%T4ddg(i,j,k,l)(0,2,0,1)");
1617  test_for_zero(t4ddg_1(0, 2, 0, 2) - (t4ddg_2(0, 2, 0, 2) / t2s_2(0, 2)),
1618  "T2s(k,l)%T4ddg(i,j,k,l)(0,2,0,2)");
1619  test_for_zero(t4ddg_1(0, 2, 1, 0) - (t4ddg_2(0, 2, 1, 0) / t2s_2(1, 0)),
1620  "T2s(k,l)%T4ddg(i,j,k,l)(0,2,1,0)");
1621  test_for_zero(t4ddg_1(0, 2, 1, 1) - (t4ddg_2(0, 2, 1, 1) / t2s_2(1, 1)),
1622  "T2s(k,l)%T4ddg(i,j,k,l)(0,2,1,1)");
1623  test_for_zero(t4ddg_1(0, 2, 1, 2) - (t4ddg_2(0, 2, 1, 2) / t2s_2(1, 2)),
1624  "T2s(k,l)%T4ddg(i,j,k,l)(0,2,1,2)");
1625  test_for_zero(t4ddg_1(0, 2, 2, 0) - (t4ddg_2(0, 2, 2, 0) / t2s_2(2, 0)),
1626  "T2s(k,l)%T4ddg(i,j,k,l)(0,2,2,0)");
1627  test_for_zero(t4ddg_1(0, 2, 2, 1) - (t4ddg_2(0, 2, 2, 1) / t2s_2(2, 1)),
1628  "T2s(k,l)%T4ddg(i,j,k,l)(0,2,2,1)");
1629  test_for_zero(t4ddg_1(0, 2, 2, 2) - (t4ddg_2(0, 2, 2, 2) / t2s_2(2, 2)),
1630  "T2s(k,l)%T4ddg(i,j,k,l)(0,2,2,2)");
1631  test_for_zero(t4ddg_1(1, 0, 0, 0) - (t4ddg_2(1, 0, 0, 0) / t2s_2(0, 0)),
1632  "T2s(k,l)%T4ddg(i,j,k,l)(1,0,0,0)");
1633  test_for_zero(t4ddg_1(1, 0, 0, 1) - (t4ddg_2(1, 0, 0, 1) / t2s_2(0, 1)),
1634  "T2s(k,l)%T4ddg(i,j,k,l)(1,0,0,1)");
1635  test_for_zero(t4ddg_1(1, 0, 0, 2) - (t4ddg_2(1, 0, 0, 2) / t2s_2(0, 2)),
1636  "T2s(k,l)%T4ddg(i,j,k,l)(1,0,0,2)");
1637  test_for_zero(t4ddg_1(1, 0, 1, 0) - (t4ddg_2(1, 0, 1, 0) / t2s_2(1, 0)),
1638  "T2s(k,l)%T4ddg(i,j,k,l)(1,0,1,0)");
1639  test_for_zero(t4ddg_1(1, 0, 1, 1) - (t4ddg_2(1, 0, 1, 1) / t2s_2(1, 1)),
1640  "T2s(k,l)%T4ddg(i,j,k,l)(1,0,1,1)");
1641  test_for_zero(t4ddg_1(1, 0, 1, 2) - (t4ddg_2(1, 0, 1, 2) / t2s_2(1, 2)),
1642  "T2s(k,l)%T4ddg(i,j,k,l)(1,0,1,2)");
1643  test_for_zero(t4ddg_1(1, 0, 2, 0) - (t4ddg_2(1, 0, 2, 0) / t2s_2(2, 0)),
1644  "T2s(k,l)%T4ddg(i,j,k,l)(1,0,2,0)");
1645  test_for_zero(t4ddg_1(1, 0, 2, 1) - (t4ddg_2(1, 0, 2, 1) / t2s_2(2, 1)),
1646  "T2s(k,l)%T4ddg(i,j,k,l)(1,0,2,1)");
1647  test_for_zero(t4ddg_1(1, 0, 2, 2) - (t4ddg_2(1, 0, 2, 2) / t2s_2(2, 2)),
1648  "T2s(k,l)%T4ddg(i,j,k,l)(1,0,2,2)");
1649  test_for_zero(t4ddg_1(1, 1, 0, 0) - (t4ddg_2(1, 1, 0, 0) / t2s_2(0, 0)),
1650  "T2s(k,l)%T4ddg(i,j,k,l)(1,1,0,0)");
1651  test_for_zero(t4ddg_1(1, 1, 0, 1) - (t4ddg_2(1, 1, 0, 1) / t2s_2(0, 1)),
1652  "T2s(k,l)%T4ddg(i,j,k,l)(1,1,0,1)");
1653  test_for_zero(t4ddg_1(1, 1, 0, 2) - (t4ddg_2(1, 1, 0, 2) / t2s_2(0, 2)),
1654  "T2s(k,l)%T4ddg(i,j,k,l)(1,1,0,2)");
1655  test_for_zero(t4ddg_1(1, 1, 1, 0) - (t4ddg_2(1, 1, 1, 0) / t2s_2(1, 0)),
1656  "T2s(k,l)%T4ddg(i,j,k,l)(1,1,1,0)");
1657  test_for_zero(t4ddg_1(1, 1, 1, 1) - (t4ddg_2(1, 1, 1, 1) / t2s_2(1, 1)),
1658  "T2s(k,l)%T4ddg(i,j,k,l)(1,1,1,1)");
1659  test_for_zero(t4ddg_1(1, 1, 1, 2) - (t4ddg_2(1, 1, 1, 2) / t2s_2(1, 2)),
1660  "T2s(k,l)%T4ddg(i,j,k,l)(1,1,1,2)");
1661  test_for_zero(t4ddg_1(1, 1, 2, 0) - (t4ddg_2(1, 1, 2, 0) / t2s_2(2, 0)),
1662  "T2s(k,l)%T4ddg(i,j,k,l)(1,1,2,0)");
1663  test_for_zero(t4ddg_1(1, 1, 2, 1) - (t4ddg_2(1, 1, 2, 1) / t2s_2(2, 1)),
1664  "T2s(k,l)%T4ddg(i,j,k,l)(1,1,2,1)");
1665  test_for_zero(t4ddg_1(1, 1, 2, 2) - (t4ddg_2(1, 1, 2, 2) / t2s_2(2, 2)),
1666  "T2s(k,l)%T4ddg(i,j,k,l)(1,1,2,2)");
1667  test_for_zero(t4ddg_1(1, 2, 0, 0) - (t4ddg_2(1, 2, 0, 0) / t2s_2(0, 0)),
1668  "T2s(k,l)%T4ddg(i,j,k,l)(1,2,0,0)");
1669  test_for_zero(t4ddg_1(1, 2, 0, 1) - (t4ddg_2(1, 2, 0, 1) / t2s_2(0, 1)),
1670  "T2s(k,l)%T4ddg(i,j,k,l)(1,2,0,1)");
1671  test_for_zero(t4ddg_1(1, 2, 0, 2) - (t4ddg_2(1, 2, 0, 2) / t2s_2(0, 2)),
1672  "T2s(k,l)%T4ddg(i,j,k,l)(1,2,0,2)");
1673  test_for_zero(t4ddg_1(1, 2, 1, 0) - (t4ddg_2(1, 2, 1, 0) / t2s_2(1, 0)),
1674  "T2s(k,l)%T4ddg(i,j,k,l)(1,2,1,0)");
1675  test_for_zero(t4ddg_1(1, 2, 1, 1) - (t4ddg_2(1, 2, 1, 1) / t2s_2(1, 1)),
1676  "T2s(k,l)%T4ddg(i,j,k,l)(1,2,1,1)");
1677  test_for_zero(t4ddg_1(1, 2, 1, 2) - (t4ddg_2(1, 2, 1, 2) / t2s_2(1, 2)),
1678  "T2s(k,l)%T4ddg(i,j,k,l)(1,2,1,2)");
1679  test_for_zero(t4ddg_1(1, 2, 2, 0) - (t4ddg_2(1, 2, 2, 0) / t2s_2(2, 0)),
1680  "T2s(k,l)%T4ddg(i,j,k,l)(1,2,2,0)");
1681  test_for_zero(t4ddg_1(1, 2, 2, 1) - (t4ddg_2(1, 2, 2, 1) / t2s_2(2, 1)),
1682  "T2s(k,l)%T4ddg(i,j,k,l)(1,2,2,1)");
1683  test_for_zero(t4ddg_1(1, 2, 2, 2) - (t4ddg_2(1, 2, 2, 2) / t2s_2(2, 2)),
1684  "T2s(k,l)%T4ddg(i,j,k,l)(1,2,2,2)");
1685  test_for_zero(t4ddg_1(2, 0, 0, 0) - (t4ddg_2(2, 0, 0, 0) / t2s_2(0, 0)),
1686  "T2s(k,l)%T4ddg(i,j,k,l)(2,0,0,0)");
1687  test_for_zero(t4ddg_1(2, 0, 0, 1) - (t4ddg_2(2, 0, 0, 1) / t2s_2(0, 1)),
1688  "T2s(k,l)%T4ddg(i,j,k,l)(2,0,0,1)");
1689  test_for_zero(t4ddg_1(2, 0, 0, 2) - (t4ddg_2(2, 0, 0, 2) / t2s_2(0, 2)),
1690  "T2s(k,l)%T4ddg(i,j,k,l)(2,0,0,2)");
1691  test_for_zero(t4ddg_1(2, 0, 1, 0) - (t4ddg_2(2, 0, 1, 0) / t2s_2(1, 0)),
1692  "T2s(k,l)%T4ddg(i,j,k,l)(2,0,1,0)");
1693  test_for_zero(t4ddg_1(2, 0, 1, 1) - (t4ddg_2(2, 0, 1, 1) / t2s_2(1, 1)),
1694  "T2s(k,l)%T4ddg(i,j,k,l)(2,0,1,1)");
1695  test_for_zero(t4ddg_1(2, 0, 1, 2) - (t4ddg_2(2, 0, 1, 2) / t2s_2(1, 2)),
1696  "T2s(k,l)%T4ddg(i,j,k,l)(2,0,1,2)");
1697  test_for_zero(t4ddg_1(2, 0, 2, 0) - (t4ddg_2(2, 0, 2, 0) / t2s_2(2, 0)),
1698  "T2s(k,l)%T4ddg(i,j,k,l)(2,0,2,0)");
1699  test_for_zero(t4ddg_1(2, 0, 2, 1) - (t4ddg_2(2, 0, 2, 1) / t2s_2(2, 1)),
1700  "T2s(k,l)%T4ddg(i,j,k,l)(2,0,2,1)");
1701  test_for_zero(t4ddg_1(2, 0, 2, 2) - (t4ddg_2(2, 0, 2, 2) / t2s_2(2, 2)),
1702  "T2s(k,l)%T4ddg(i,j,k,l)(2,0,2,2)");
1703  test_for_zero(t4ddg_1(2, 1, 0, 0) - (t4ddg_2(2, 1, 0, 0) / t2s_2(0, 0)),
1704  "T2s(k,l)%T4ddg(i,j,k,l)(2,1,0,0)");
1705  test_for_zero(t4ddg_1(2, 1, 0, 1) - (t4ddg_2(2, 1, 0, 1) / t2s_2(0, 1)),
1706  "T2s(k,l)%T4ddg(i,j,k,l)(2,1,0,1)");
1707  test_for_zero(t4ddg_1(2, 1, 0, 2) - (t4ddg_2(2, 1, 0, 2) / t2s_2(0, 2)),
1708  "T2s(k,l)%T4ddg(i,j,k,l)(2,1,0,2)");
1709  test_for_zero(t4ddg_1(2, 1, 1, 0) - (t4ddg_2(2, 1, 1, 0) / t2s_2(1, 0)),
1710  "T2s(k,l)%T4ddg(i,j,k,l)(2,1,1,0)");
1711  test_for_zero(t4ddg_1(2, 1, 1, 1) - (t4ddg_2(2, 1, 1, 1) / t2s_2(1, 1)),
1712  "T2s(k,l)%T4ddg(i,j,k,l)(2,1,1,1)");
1713  test_for_zero(t4ddg_1(2, 1, 1, 2) - (t4ddg_2(2, 1, 1, 2) / t2s_2(1, 2)),
1714  "T2s(k,l)%T4ddg(i,j,k,l)(2,1,1,2)");
1715  test_for_zero(t4ddg_1(2, 1, 2, 0) - (t4ddg_2(2, 1, 2, 0) / t2s_2(2, 0)),
1716  "T2s(k,l)%T4ddg(i,j,k,l)(2,1,2,0)");
1717  test_for_zero(t4ddg_1(2, 1, 2, 1) - (t4ddg_2(2, 1, 2, 1) / t2s_2(2, 1)),
1718  "T2s(k,l)%T4ddg(i,j,k,l)(2,1,2,1)");
1719  test_for_zero(t4ddg_1(2, 1, 2, 2) - (t4ddg_2(2, 1, 2, 2) / t2s_2(2, 2)),
1720  "T2s(k,l)%T4ddg(i,j,k,l)(2,1,2,2)");
1721  test_for_zero(t4ddg_1(2, 2, 0, 0) - (t4ddg_2(2, 2, 0, 0) / t2s_2(0, 0)),
1722  "T2s(k,l)%T4ddg(i,j,k,l)(2,2,0,0)");
1723  test_for_zero(t4ddg_1(2, 2, 0, 1) - (t4ddg_2(2, 2, 0, 1) / t2s_2(0, 1)),
1724  "T2s(k,l)%T4ddg(i,j,k,l)(2,2,0,1)");
1725  test_for_zero(t4ddg_1(2, 2, 0, 2) - (t4ddg_2(2, 2, 0, 2) / t2s_2(0, 2)),
1726  "T2s(k,l)%T4ddg(i,j,k,l)(2,2,0,2)");
1727  test_for_zero(t4ddg_1(2, 2, 1, 0) - (t4ddg_2(2, 2, 1, 0) / t2s_2(1, 0)),
1728  "T2s(k,l)%T4ddg(i,j,k,l)(2,2,1,0)");
1729  test_for_zero(t4ddg_1(2, 2, 1, 1) - (t4ddg_2(2, 2, 1, 1) / t2s_2(1, 1)),
1730  "T2s(k,l)%T4ddg(i,j,k,l)(2,2,1,1)");
1731  test_for_zero(t4ddg_1(2, 2, 1, 2) - (t4ddg_2(2, 2, 1, 2) / t2s_2(1, 2)),
1732  "T2s(k,l)%T4ddg(i,j,k,l)(2,2,1,2)");
1733  test_for_zero(t4ddg_1(2, 2, 2, 0) - (t4ddg_2(2, 2, 2, 0) / t2s_2(2, 0)),
1734  "T2s(k,l)%T4ddg(i,j,k,l)(2,2,2,0)");
1735  test_for_zero(t4ddg_1(2, 2, 2, 1) - (t4ddg_2(2, 2, 2, 1) / t2s_2(2, 1)),
1736  "T2s(k,l)%T4ddg(i,j,k,l)(2,2,2,1)");
1737  test_for_zero(t4ddg_1(2, 2, 2, 2) - (t4ddg_2(2, 2, 2, 2) / t2s_2(2, 2)),
1738  "T2s(k,l)%T4ddg(i,j,k,l)(2,2,2,2)");
1739 
1740  /* I originally put these declarations for unknown reasons, but they
1741  won't work because the result is not a Ddg. The
1742  multiplication messes up the symmetries. */
1743 
1744  // t4ddg_1(i,j,k,l)=(t4ddg_2(i,j,k,l)&t2s_3(j,l));
1745  // test_for_zero(t4ddg_1(0,0,0,0) - (t4ddg_2(0,0,0,0)*t2s_3(0,0))
1746  // test_for_zero(t4ddg_1(0,0,0,1) - (t4ddg_2(0,0,0,1)*t2s_3(0,1))
1747  // test_for_zero(t4ddg_1(0,0,0,2) - (t4ddg_2(0,0,0,2)*t2s_3(0,2))
1748  // test_for_zero(t4ddg_1(0,0,1,0) - (t4ddg_2(0,0,1,0)*t2s_3(0,0))
1749  // test_for_zero(t4ddg_1(0,0,1,1) - (t4ddg_2(0,0,1,1)*t2s_3(0,1))
1750  // test_for_zero(t4ddg_1(0,0,1,2) - (t4ddg_2(0,0,1,2)*t2s_3(0,2))
1751  // test_for_zero(t4ddg_1(0,0,2,0) - (t4ddg_2(0,0,2,0)*t2s_3(0,0))
1752  // test_for_zero(t4ddg_1(0,0,2,1) - (t4ddg_2(0,0,2,1)*t2s_3(0,1))
1753  // test_for_zero(t4ddg_1(0,0,2,2) - (t4ddg_2(0,0,2,2)*t2s_3(0,2))
1754  // test_for_zero(t4ddg_1(0,1,0,0) - (t4ddg_2(0,1,0,0)*t2s_3(1,0))
1755  // test_for_zero(t4ddg_1(0,1,0,1) - (t4ddg_2(0,1,0,1)*t2s_3(1,1))
1756  // test_for_zero(t4ddg_1(0,1,0,2) - (t4ddg_2(0,1,0,2)*t2s_3(1,2))
1757  // test_for_zero(t4ddg_1(0,1,1,0) - (t4ddg_2(0,1,1,0)*t2s_3(1,0))
1758  // test_for_zero(t4ddg_1(0,1,1,1) - (t4ddg_2(0,1,1,1)*t2s_3(1,1))
1759  // test_for_zero(t4ddg_1(0,1,1,2) - (t4ddg_2(0,1,1,2)*t2s_3(1,2))
1760  // test_for_zero(t4ddg_1(0,1,2,0) - (t4ddg_2(0,1,2,0)*t2s_3(1,0))
1761  // test_for_zero(t4ddg_1(0,1,2,1) - (t4ddg_2(0,1,2,1)*t2s_3(1,1))
1762  // test_for_zero(t4ddg_1(0,1,2,2) - (t4ddg_2(0,1,2,2)*t2s_3(1,2))
1763  // test_for_zero(t4ddg_1(0,2,0,0) - (t4ddg_2(0,2,0,0)*t2s_3(2,0))
1764  // test_for_zero(t4ddg_1(0,2,0,1) - (t4ddg_2(0,2,0,1)*t2s_3(2,1))
1765  // test_for_zero(t4ddg_1(0,2,0,2) - (t4ddg_2(0,2,0,2)*t2s_3(2,2))
1766  // test_for_zero(t4ddg_1(0,2,1,0) - (t4ddg_2(0,2,1,0)*t2s_3(2,0))
1767  // test_for_zero(t4ddg_1(0,2,1,1) - (t4ddg_2(0,2,1,1)*t2s_3(2,1))
1768  // test_for_zero(t4ddg_1(0,2,1,2) - (t4ddg_2(0,2,1,2)*t2s_3(2,2))
1769  // test_for_zero(t4ddg_1(0,2,2,0) - (t4ddg_2(0,2,2,0)*t2s_3(2,0))
1770  // test_for_zero(t4ddg_1(0,2,2,1) - (t4ddg_2(0,2,2,1)*t2s_3(2,1))
1771  // test_for_zero(t4ddg_1(0,2,2,2) - (t4ddg_2(0,2,2,2)*t2s_3(2,2))
1772  // test_for_zero(t4ddg_1(1,0,0,0) - (t4ddg_2(1,0,0,0)*t2s_3(0,0))
1773  // test_for_zero(t4ddg_1(1,0,0,1) - (t4ddg_2(1,0,0,1)*t2s_3(0,1))
1774  // test_for_zero(t4ddg_1(1,0,0,2) - (t4ddg_2(1,0,0,2)*t2s_3(0,2))
1775  // test_for_zero(t4ddg_1(1,0,1,0) - (t4ddg_2(1,0,1,0)*t2s_3(0,0))
1776  // test_for_zero(t4ddg_1(1,0,1,1) - (t4ddg_2(1,0,1,1)*t2s_3(0,1))
1777  // test_for_zero(t4ddg_1(1,0,1,2) - (t4ddg_2(1,0,1,2)*t2s_3(0,2))
1778  // test_for_zero(t4ddg_1(1,0,2,0) - (t4ddg_2(1,0,2,0)*t2s_3(0,0))
1779  // test_for_zero(t4ddg_1(1,0,2,1) - (t4ddg_2(1,0,2,1)*t2s_3(0,1))
1780  // test_for_zero(t4ddg_1(1,0,2,2) - (t4ddg_2(1,0,2,2)*t2s_3(0,2))
1781  // test_for_zero(t4ddg_1(1,1,0,0) - (t4ddg_2(1,1,0,0)*t2s_3(1,0))
1782  // test_for_zero(t4ddg_1(1,1,0,1) - (t4ddg_2(1,1,0,1)*t2s_3(1,1))
1783  // test_for_zero(t4ddg_1(1,1,0,2) - (t4ddg_2(1,1,0,2)*t2s_3(1,2))
1784  // test_for_zero(t4ddg_1(1,1,1,0) - (t4ddg_2(1,1,1,0)*t2s_3(1,0))
1785  // test_for_zero(t4ddg_1(1,1,1,1) - (t4ddg_2(1,1,1,1)*t2s_3(1,1))
1786  // test_for_zero(t4ddg_1(1,1,1,2) - (t4ddg_2(1,1,1,2)*t2s_3(1,2))
1787  // test_for_zero(t4ddg_1(1,1,2,0) - (t4ddg_2(1,1,2,0)*t2s_3(1,0))
1788  // test_for_zero(t4ddg_1(1,1,2,1) - (t4ddg_2(1,1,2,1)*t2s_3(1,1))
1789  // test_for_zero(t4ddg_1(1,1,2,2) - (t4ddg_2(1,1,2,2)*t2s_3(1,2))
1790  // test_for_zero(t4ddg_1(1,2,0,0) - (t4ddg_2(1,2,0,0)*t2s_3(2,0))
1791  // test_for_zero(t4ddg_1(1,2,0,1) - (t4ddg_2(1,2,0,1)*t2s_3(2,1))
1792  // test_for_zero(t4ddg_1(1,2,0,2) - (t4ddg_2(1,2,0,2)*t2s_3(2,2))
1793  // test_for_zero(t4ddg_1(1,2,1,0) - (t4ddg_2(1,2,1,0)*t2s_3(2,0))
1794  // test_for_zero(t4ddg_1(1,2,1,1) - (t4ddg_2(1,2,1,1)*t2s_3(2,1))
1795  // test_for_zero(t4ddg_1(1,2,1,2) - (t4ddg_2(1,2,1,2)*t2s_3(2,2))
1796  // test_for_zero(t4ddg_1(1,2,2,0) - (t4ddg_2(1,2,2,0)*t2s_3(2,0))
1797  // test_for_zero(t4ddg_1(1,2,2,1) - (t4ddg_2(1,2,2,1)*t2s_3(2,1))
1798  // test_for_zero(t4ddg_1(1,2,2,2) - (t4ddg_2(1,2,2,2)*t2s_3(2,2))
1799  // test_for_zero(t4ddg_1(2,0,0,0) - (t4ddg_2(2,0,0,0)*t2s_3(0,0))
1800  // test_for_zero(t4ddg_1(2,0,0,1) - (t4ddg_2(2,0,0,1)*t2s_3(0,1))
1801  // test_for_zero(t4ddg_1(2,0,0,2) - (t4ddg_2(2,0,0,2)*t2s_3(0,2))
1802  // test_for_zero(t4ddg_1(2,0,1,0) - (t4ddg_2(2,0,1,0)*t2s_3(0,0))
1803  // test_for_zero(t4ddg_1(2,0,1,1) - (t4ddg_2(2,0,1,1)*t2s_3(0,1))
1804  // test_for_zero(t4ddg_1(2,0,1,2) - (t4ddg_2(2,0,1,2)*t2s_3(0,2))
1805  // test_for_zero(t4ddg_1(2,0,2,0) - (t4ddg_2(2,0,2,0)*t2s_3(0,0))
1806  // test_for_zero(t4ddg_1(2,0,2,1) - (t4ddg_2(2,0,2,1)*t2s_3(0,1))
1807  // test_for_zero(t4ddg_1(2,0,2,2) - (t4ddg_2(2,0,2,2)*t2s_3(0,2))
1808  // test_for_zero(t4ddg_1(2,1,0,0) - (t4ddg_2(2,1,0,0)*t2s_3(1,0))
1809  // test_for_zero(t4ddg_1(2,1,0,1) - (t4ddg_2(2,1,0,1)*t2s_3(1,1))
1810  // test_for_zero(t4ddg_1(2,1,0,2) - (t4ddg_2(2,1,0,2)*t2s_3(1,2))
1811  // test_for_zero(t4ddg_1(2,1,1,0) - (t4ddg_2(2,1,1,0)*t2s_3(1,0))
1812  // test_for_zero(t4ddg_1(2,1,1,1) - (t4ddg_2(2,1,1,1)*t2s_3(1,1))
1813  // test_for_zero(t4ddg_1(2,1,1,2) - (t4ddg_2(2,1,1,2)*t2s_3(1,2))
1814  // test_for_zero(t4ddg_1(2,1,2,0) - (t4ddg_2(2,1,2,0)*t2s_3(1,0))
1815  // test_for_zero(t4ddg_1(2,1,2,1) - (t4ddg_2(2,1,2,1)*t2s_3(1,1))
1816  // test_for_zero(t4ddg_1(2,1,2,2) - (t4ddg_2(2,1,2,2)*t2s_3(1,2))
1817  // test_for_zero(t4ddg_1(2,2,0,0) - (t4ddg_2(2,2,0,0)*t2s_3(2,0))
1818  // test_for_zero(t4ddg_1(2,2,0,1) - (t4ddg_2(2,2,0,1)*t2s_3(2,1))
1819  // test_for_zero(t4ddg_1(2,2,0,2) - (t4ddg_2(2,2,0,2)*t2s_3(2,2))
1820  // test_for_zero(t4ddg_1(2,2,1,0) - (t4ddg_2(2,2,1,0)*t2s_3(2,0))
1821  // test_for_zero(t4ddg_1(2,2,1,1) - (t4ddg_2(2,2,1,1)*t2s_3(2,1))
1822  // test_for_zero(t4ddg_1(2,2,1,2) - (t4ddg_2(2,2,1,2)*t2s_3(2,2))
1823  // test_for_zero(t4ddg_1(2,2,2,0) - (t4ddg_2(2,2,2,0)*t2s_3(2,0))
1824  // test_for_zero(t4ddg_1(2,2,2,1) - (t4ddg_2(2,2,2,1)*t2s_3(2,1))
1825  // test_for_zero(t4ddg_1(2,2,2,2) - (t4ddg_2(2,2,2,2)*t2s_3(2,2)));
1826 
1827  // t4ddg_1(i,j,k,l)=(t2s_2(j,l)&t4ddg_2(i,j,k,l));
1828  // test_for_zero(t4ddg_1(0,0,0,0) - (t4ddg_2(0,0,0,0)*t2s_2(0,0))
1829  // test_for_zero(t4ddg_1(0,0,0,1) - (t4ddg_2(0,0,0,1)*t2s_2(0,1))
1830  // test_for_zero(t4ddg_1(0,0,0,2) - (t4ddg_2(0,0,0,2)*t2s_2(0,2))
1831  // test_for_zero(t4ddg_1(0,0,1,0) - (t4ddg_2(0,0,1,0)*t2s_2(0,0))
1832  // test_for_zero(t4ddg_1(0,0,1,1) - (t4ddg_2(0,0,1,1)*t2s_2(0,1))
1833  // test_for_zero(t4ddg_1(0,0,1,2) - (t4ddg_2(0,0,1,2)*t2s_2(0,2))
1834  // test_for_zero(t4ddg_1(0,0,2,0) - (t4ddg_2(0,0,2,0)*t2s_2(0,0))
1835  // test_for_zero(t4ddg_1(0,0,2,1) - (t4ddg_2(0,0,2,1)*t2s_2(0,1))
1836  // test_for_zero(t4ddg_1(0,0,2,2) - (t4ddg_2(0,0,2,2)*t2s_2(0,2))
1837  // test_for_zero(t4ddg_1(0,1,0,0) - (t4ddg_2(0,1,0,0)*t2s_2(1,0))
1838  // test_for_zero(t4ddg_1(0,1,0,1) - (t4ddg_2(0,1,0,1)*t2s_2(1,1))
1839  // test_for_zero(t4ddg_1(0,1,0,2) - (t4ddg_2(0,1,0,2)*t2s_2(1,2))
1840  // test_for_zero(t4ddg_1(0,1,1,0) - (t4ddg_2(0,1,1,0)*t2s_2(1,0))
1841  // test_for_zero(t4ddg_1(0,1,1,1) - (t4ddg_2(0,1,1,1)*t2s_2(1,1))
1842  // test_for_zero(t4ddg_1(0,1,1,2) - (t4ddg_2(0,1,1,2)*t2s_2(1,2))
1843  // test_for_zero(t4ddg_1(0,1,2,0) - (t4ddg_2(0,1,2,0)*t2s_2(1,0))
1844  // test_for_zero(t4ddg_1(0,1,2,1) - (t4ddg_2(0,1,2,1)*t2s_2(1,1))
1845  // test_for_zero(t4ddg_1(0,1,2,2) - (t4ddg_2(0,1,2,2)*t2s_2(1,2))
1846  // test_for_zero(t4ddg_1(0,2,0,0) - (t4ddg_2(0,2,0,0)*t2s_2(2,0))
1847  // test_for_zero(t4ddg_1(0,2,0,1) - (t4ddg_2(0,2,0,1)*t2s_2(2,1))
1848  // test_for_zero(t4ddg_1(0,2,0,2) - (t4ddg_2(0,2,0,2)*t2s_2(2,2))
1849  // test_for_zero(t4ddg_1(0,2,1,0) - (t4ddg_2(0,2,1,0)*t2s_2(2,0))
1850  // test_for_zero(t4ddg_1(0,2,1,1) - (t4ddg_2(0,2,1,1)*t2s_2(2,1))
1851  // test_for_zero(t4ddg_1(0,2,1,2) - (t4ddg_2(0,2,1,2)*t2s_2(2,2))
1852  // test_for_zero(t4ddg_1(0,2,2,0) - (t4ddg_2(0,2,2,0)*t2s_2(2,0))
1853  // test_for_zero(t4ddg_1(0,2,2,1) - (t4ddg_2(0,2,2,1)*t2s_2(2,1))
1854  // test_for_zero(t4ddg_1(0,2,2,2) - (t4ddg_2(0,2,2,2)*t2s_2(2,2))
1855  // test_for_zero(t4ddg_1(1,0,0,0) - (t4ddg_2(1,0,0,0)*t2s_2(0,0))
1856  // test_for_zero(t4ddg_1(1,0,0,1) - (t4ddg_2(1,0,0,1)*t2s_2(0,1))
1857  // test_for_zero(t4ddg_1(1,0,0,2) - (t4ddg_2(1,0,0,2)*t2s_2(0,2))
1858  // test_for_zero(t4ddg_1(1,0,1,0) - (t4ddg_2(1,0,1,0)*t2s_2(0,0))
1859  // test_for_zero(t4ddg_1(1,0,1,1) - (t4ddg_2(1,0,1,1)*t2s_2(0,1))
1860  // test_for_zero(t4ddg_1(1,0,1,2) - (t4ddg_2(1,0,1,2)*t2s_2(0,2))
1861  // test_for_zero(t4ddg_1(1,0,2,0) - (t4ddg_2(1,0,2,0)*t2s_2(0,0))
1862  // test_for_zero(t4ddg_1(1,0,2,1) - (t4ddg_2(1,0,2,1)*t2s_2(0,1))
1863  // test_for_zero(t4ddg_1(1,0,2,2) - (t4ddg_2(1,0,2,2)*t2s_2(0,2))
1864  // test_for_zero(t4ddg_1(1,1,0,0) - (t4ddg_2(1,1,0,0)*t2s_2(1,0))
1865  // test_for_zero(t4ddg_1(1,1,0,1) - (t4ddg_2(1,1,0,1)*t2s_2(1,1))
1866  // test_for_zero(t4ddg_1(1,1,0,2) - (t4ddg_2(1,1,0,2)*t2s_2(1,2))
1867  // test_for_zero(t4ddg_1(1,1,1,0) - (t4ddg_2(1,1,1,0)*t2s_2(1,0))
1868  // test_for_zero(t4ddg_1(1,1,1,1) - (t4ddg_2(1,1,1,1)*t2s_2(1,1))
1869  // test_for_zero(t4ddg_1(1,1,1,2) - (t4ddg_2(1,1,1,2)*t2s_2(1,2))
1870  // test_for_zero(t4ddg_1(1,1,2,0) - (t4ddg_2(1,1,2,0)*t2s_2(1,0))
1871  // test_for_zero(t4ddg_1(1,1,2,1) - (t4ddg_2(1,1,2,1)*t2s_2(1,1))
1872  // test_for_zero(t4ddg_1(1,1,2,2) - (t4ddg_2(1,1,2,2)*t2s_2(1,2))
1873  // test_for_zero(t4ddg_1(1,2,0,0) - (t4ddg_2(1,2,0,0)*t2s_2(2,0))
1874  // test_for_zero(t4ddg_1(1,2,0,1) - (t4ddg_2(1,2,0,1)*t2s_2(2,1))
1875  // test_for_zero(t4ddg_1(1,2,0,2) - (t4ddg_2(1,2,0,2)*t2s_2(2,2))
1876  // test_for_zero(t4ddg_1(1,2,1,0) - (t4ddg_2(1,2,1,0)*t2s_2(2,0))
1877  // test_for_zero(t4ddg_1(1,2,1,1) - (t4ddg_2(1,2,1,1)*t2s_2(2,1))
1878  // test_for_zero(t4ddg_1(1,2,1,2) - (t4ddg_2(1,2,1,2)*t2s_2(2,2))
1879  // test_for_zero(t4ddg_1(1,2,2,0) - (t4ddg_2(1,2,2,0)*t2s_2(2,0))
1880  // test_for_zero(t4ddg_1(1,2,2,1) - (t4ddg_2(1,2,2,1)*t2s_2(2,1))
1881  // test_for_zero(t4ddg_1(1,2,2,2) - (t4ddg_2(1,2,2,2)*t2s_2(2,2))
1882  // test_for_zero(t4ddg_1(2,0,0,0) - (t4ddg_2(2,0,0,0)*t2s_2(0,0))
1883  // test_for_zero(t4ddg_1(2,0,0,1) - (t4ddg_2(2,0,0,1)*t2s_2(0,1))
1884  // test_for_zero(t4ddg_1(2,0,0,2) - (t4ddg_2(2,0,0,2)*t2s_2(0,2))
1885  // test_for_zero(t4ddg_1(2,0,1,0) - (t4ddg_2(2,0,1,0)*t2s_2(0,0))
1886  // test_for_zero(t4ddg_1(2,0,1,1) - (t4ddg_2(2,0,1,1)*t2s_2(0,1))
1887  // test_for_zero(t4ddg_1(2,0,1,2) - (t4ddg_2(2,0,1,2)*t2s_2(0,2))
1888  // test_for_zero(t4ddg_1(2,0,2,0) - (t4ddg_2(2,0,2,0)*t2s_2(0,0))
1889  // test_for_zero(t4ddg_1(2,0,2,1) - (t4ddg_2(2,0,2,1)*t2s_2(0,1))
1890  // test_for_zero(t4ddg_1(2,0,2,2) - (t4ddg_2(2,0,2,2)*t2s_2(0,2))
1891  // test_for_zero(t4ddg_1(2,1,0,0) - (t4ddg_2(2,1,0,0)*t2s_2(1,0))
1892  // test_for_zero(t4ddg_1(2,1,0,1) - (t4ddg_2(2,1,0,1)*t2s_2(1,1))
1893  // test_for_zero(t4ddg_1(2,1,0,2) - (t4ddg_2(2,1,0,2)*t2s_2(1,2))
1894  // test_for_zero(t4ddg_1(2,1,1,0) - (t4ddg_2(2,1,1,0)*t2s_2(1,0))
1895  // test_for_zero(t4ddg_1(2,1,1,1) - (t4ddg_2(2,1,1,1)*t2s_2(1,1))
1896  // test_for_zero(t4ddg_1(2,1,1,2) - (t4ddg_2(2,1,1,2)*t2s_2(1,2))
1897  // test_for_zero(t4ddg_1(2,1,2,0) - (t4ddg_2(2,1,2,0)*t2s_2(1,0))
1898  // test_for_zero(t4ddg_1(2,1,2,1) - (t4ddg_2(2,1,2,1)*t2s_2(1,1))
1899  // test_for_zero(t4ddg_1(2,1,2,2) - (t4ddg_2(2,1,2,2)*t2s_2(1,2))
1900  // test_for_zero(t4ddg_1(2,2,0,0) - (t4ddg_2(2,2,0,0)*t2s_2(2,0))
1901  // test_for_zero(t4ddg_1(2,2,0,1) - (t4ddg_2(2,2,0,1)*t2s_2(2,1))
1902  // test_for_zero(t4ddg_1(2,2,0,2) - (t4ddg_2(2,2,0,2)*t2s_2(2,2))
1903  // test_for_zero(t4ddg_1(2,2,1,0) - (t4ddg_2(2,2,1,0)*t2s_2(2,0))
1904  // test_for_zero(t4ddg_1(2,2,1,1) - (t4ddg_2(2,2,1,1)*t2s_2(2,1))
1905  // test_for_zero(t4ddg_1(2,2,1,2) - (t4ddg_2(2,2,1,2)*t2s_2(2,2))
1906  // test_for_zero(t4ddg_1(2,2,2,0) - (t4ddg_2(2,2,2,0)*t2s_2(2,0))
1907  // test_for_zero(t4ddg_1(2,2,2,1) - (t4ddg_2(2,2,2,1)*t2s_2(2,1))
1908  // test_for_zero(t4ddg_1(2,2,2,2) - (t4ddg_2(2,2,2,2)*t2s_2(2,2)));
1909  // cout << endl;
1910 
1911  // t4ddg_1(i,j,k,l)=(t4ddg_2(i,j,k,l)&t2s_3(l,j));
1912  // test_for_zero(t4ddg_1(0,0,0,0) - (t4ddg_2(0,0,0,0)*t2s_3(0,0))
1913  // test_for_zero(t4ddg_1(0,0,0,1) - (t4ddg_2(0,0,0,1)*t2s_3(0,1))
1914  // test_for_zero(t4ddg_1(0,0,0,2) - (t4ddg_2(0,0,0,2)*t2s_3(0,2))
1915  // test_for_zero(t4ddg_1(0,0,1,0) - (t4ddg_2(0,0,1,0)*t2s_3(0,0))
1916  // test_for_zero(t4ddg_1(0,0,1,1) - (t4ddg_2(0,0,1,1)*t2s_3(0,1))
1917  // test_for_zero(t4ddg_1(0,0,1,2) - (t4ddg_2(0,0,1,2)*t2s_3(0,2))
1918  // test_for_zero(t4ddg_1(0,0,2,0) - (t4ddg_2(0,0,2,0)*t2s_3(0,0))
1919  // test_for_zero(t4ddg_1(0,0,2,1) - (t4ddg_2(0,0,2,1)*t2s_3(0,1))
1920  // test_for_zero(t4ddg_1(0,0,2,2) - (t4ddg_2(0,0,2,2)*t2s_3(0,2))
1921  // test_for_zero(t4ddg_1(0,1,0,0) - (t4ddg_2(0,1,0,0)*t2s_3(1,0))
1922  // test_for_zero(t4ddg_1(0,1,0,1) - (t4ddg_2(0,1,0,1)*t2s_3(1,1))
1923  // test_for_zero(t4ddg_1(0,1,0,2) - (t4ddg_2(0,1,0,2)*t2s_3(1,2))
1924  // test_for_zero(t4ddg_1(0,1,1,0) - (t4ddg_2(0,1,1,0)*t2s_3(1,0))
1925  // test_for_zero(t4ddg_1(0,1,1,1) - (t4ddg_2(0,1,1,1)*t2s_3(1,1))
1926  // test_for_zero(t4ddg_1(0,1,1,2) - (t4ddg_2(0,1,1,2)*t2s_3(1,2))
1927  // test_for_zero(t4ddg_1(0,1,2,0) - (t4ddg_2(0,1,2,0)*t2s_3(1,0))
1928  // test_for_zero(t4ddg_1(0,1,2,1) - (t4ddg_2(0,1,2,1)*t2s_3(1,1))
1929  // test_for_zero(t4ddg_1(0,1,2,2) - (t4ddg_2(0,1,2,2)*t2s_3(1,2))
1930  // test_for_zero(t4ddg_1(0,2,0,0) - (t4ddg_2(0,2,0,0)*t2s_3(2,0))
1931  // test_for_zero(t4ddg_1(0,2,0,1) - (t4ddg_2(0,2,0,1)*t2s_3(2,1))
1932  // test_for_zero(t4ddg_1(0,2,0,2) - (t4ddg_2(0,2,0,2)*t2s_3(2,2))
1933  // test_for_zero(t4ddg_1(0,2,1,0) - (t4ddg_2(0,2,1,0)*t2s_3(2,0))
1934  // test_for_zero(t4ddg_1(0,2,1,1) - (t4ddg_2(0,2,1,1)*t2s_3(2,1))
1935  // test_for_zero(t4ddg_1(0,2,1,2) - (t4ddg_2(0,2,1,2)*t2s_3(2,2))
1936  // test_for_zero(t4ddg_1(0,2,2,0) - (t4ddg_2(0,2,2,0)*t2s_3(2,0))
1937  // test_for_zero(t4ddg_1(0,2,2,1) - (t4ddg_2(0,2,2,1)*t2s_3(2,1))
1938  // test_for_zero(t4ddg_1(0,2,2,2) - (t4ddg_2(0,2,2,2)*t2s_3(2,2))
1939  // test_for_zero(t4ddg_1(1,0,0,0) - (t4ddg_2(1,0,0,0)*t2s_3(0,0))
1940  // test_for_zero(t4ddg_1(1,0,0,1) - (t4ddg_2(1,0,0,1)*t2s_3(0,1))
1941  // test_for_zero(t4ddg_1(1,0,0,2) - (t4ddg_2(1,0,0,2)*t2s_3(0,2))
1942  // test_for_zero(t4ddg_1(1,0,1,0) - (t4ddg_2(1,0,1,0)*t2s_3(0,0))
1943  // test_for_zero(t4ddg_1(1,0,1,1) - (t4ddg_2(1,0,1,1)*t2s_3(0,1))
1944  // test_for_zero(t4ddg_1(1,0,1,2) - (t4ddg_2(1,0,1,2)*t2s_3(0,2))
1945  // test_for_zero(t4ddg_1(1,0,2,0) - (t4ddg_2(1,0,2,0)*t2s_3(0,0))
1946  // test_for_zero(t4ddg_1(1,0,2,1) - (t4ddg_2(1,0,2,1)*t2s_3(0,1))
1947  // test_for_zero(t4ddg_1(1,0,2,2) - (t4ddg_2(1,0,2,2)*t2s_3(0,2))
1948  // test_for_zero(t4ddg_1(1,1,0,0) - (t4ddg_2(1,1,0,0)*t2s_3(1,0))
1949  // test_for_zero(t4ddg_1(1,1,0,1) - (t4ddg_2(1,1,0,1)*t2s_3(1,1))
1950  // test_for_zero(t4ddg_1(1,1,0,2) - (t4ddg_2(1,1,0,2)*t2s_3(1,2))
1951  // test_for_zero(t4ddg_1(1,1,1,0) - (t4ddg_2(1,1,1,0)*t2s_3(1,0))
1952  // test_for_zero(t4ddg_1(1,1,1,1) - (t4ddg_2(1,1,1,1)*t2s_3(1,1))
1953  // test_for_zero(t4ddg_1(1,1,1,2) - (t4ddg_2(1,1,1,2)*t2s_3(1,2))
1954  // test_for_zero(t4ddg_1(1,1,2,0) - (t4ddg_2(1,1,2,0)*t2s_3(1,0))
1955  // test_for_zero(t4ddg_1(1,1,2,1) - (t4ddg_2(1,1,2,1)*t2s_3(1,1))
1956  // test_for_zero(t4ddg_1(1,1,2,2) - (t4ddg_2(1,1,2,2)*t2s_3(1,2))
1957  // test_for_zero(t4ddg_1(1,2,0,0) - (t4ddg_2(1,2,0,0)*t2s_3(2,0))
1958  // test_for_zero(t4ddg_1(1,2,0,1) - (t4ddg_2(1,2,0,1)*t2s_3(2,1))
1959  // test_for_zero(t4ddg_1(1,2,0,2) - (t4ddg_2(1,2,0,2)*t2s_3(2,2))
1960  // test_for_zero(t4ddg_1(1,2,1,0) - (t4ddg_2(1,2,1,0)*t2s_3(2,0))
1961  // test_for_zero(t4ddg_1(1,2,1,1) - (t4ddg_2(1,2,1,1)*t2s_3(2,1))
1962  // test_for_zero(t4ddg_1(1,2,1,2) - (t4ddg_2(1,2,1,2)*t2s_3(2,2))
1963  // test_for_zero(t4ddg_1(1,2,2,0) - (t4ddg_2(1,2,2,0)*t2s_3(2,0))
1964  // test_for_zero(t4ddg_1(1,2,2,1) - (t4ddg_2(1,2,2,1)*t2s_3(2,1))
1965  // test_for_zero(t4ddg_1(1,2,2,2) - (t4ddg_2(1,2,2,2)*t2s_3(2,2))
1966  // test_for_zero(t4ddg_1(2,0,0,0) - (t4ddg_2(2,0,0,0)*t2s_3(0,0))
1967  // test_for_zero(t4ddg_1(2,0,0,1) - (t4ddg_2(2,0,0,1)*t2s_3(0,1))
1968  // test_for_zero(t4ddg_1(2,0,0,2) - (t4ddg_2(2,0,0,2)*t2s_3(0,2))
1969  // test_for_zero(t4ddg_1(2,0,1,0) - (t4ddg_2(2,0,1,0)*t2s_3(0,0))
1970  // test_for_zero(t4ddg_1(2,0,1,1) - (t4ddg_2(2,0,1,1)*t2s_3(0,1))
1971  // test_for_zero(t4ddg_1(2,0,1,2) - (t4ddg_2(2,0,1,2)*t2s_3(0,2))
1972  // test_for_zero(t4ddg_1(2,0,2,0) - (t4ddg_2(2,0,2,0)*t2s_3(0,0))
1973  // test_for_zero(t4ddg_1(2,0,2,1) - (t4ddg_2(2,0,2,1)*t2s_3(0,1))
1974  // test_for_zero(t4ddg_1(2,0,2,2) - (t4ddg_2(2,0,2,2)*t2s_3(0,2))
1975  // test_for_zero(t4ddg_1(2,1,0,0) - (t4ddg_2(2,1,0,0)*t2s_3(1,0))
1976  // test_for_zero(t4ddg_1(2,1,0,1) - (t4ddg_2(2,1,0,1)*t2s_3(1,1))
1977  // test_for_zero(t4ddg_1(2,1,0,2) - (t4ddg_2(2,1,0,2)*t2s_3(1,2))
1978  // test_for_zero(t4ddg_1(2,1,1,0) - (t4ddg_2(2,1,1,0)*t2s_3(1,0))
1979  // test_for_zero(t4ddg_1(2,1,1,1) - (t4ddg_2(2,1,1,1)*t2s_3(1,1))
1980  // test_for_zero(t4ddg_1(2,1,1,2) - (t4ddg_2(2,1,1,2)*t2s_3(1,2))
1981  // test_for_zero(t4ddg_1(2,1,2,0) - (t4ddg_2(2,1,2,0)*t2s_3(1,0))
1982  // test_for_zero(t4ddg_1(2,1,2,1) - (t4ddg_2(2,1,2,1)*t2s_3(1,1))
1983  // test_for_zero(t4ddg_1(2,1,2,2) - (t4ddg_2(2,1,2,2)*t2s_3(1,2))
1984  // test_for_zero(t4ddg_1(2,2,0,0) - (t4ddg_2(2,2,0,0)*t2s_3(2,0))
1985  // test_for_zero(t4ddg_1(2,2,0,1) - (t4ddg_2(2,2,0,1)*t2s_3(2,1))
1986  // test_for_zero(t4ddg_1(2,2,0,2) - (t4ddg_2(2,2,0,2)*t2s_3(2,2))
1987  // test_for_zero(t4ddg_1(2,2,1,0) - (t4ddg_2(2,2,1,0)*t2s_3(2,0))
1988  // test_for_zero(t4ddg_1(2,2,1,1) - (t4ddg_2(2,2,1,1)*t2s_3(2,1))
1989  // test_for_zero(t4ddg_1(2,2,1,2) - (t4ddg_2(2,2,1,2)*t2s_3(2,2))
1990  // test_for_zero(t4ddg_1(2,2,2,0) - (t4ddg_2(2,2,2,0)*t2s_3(2,0))
1991  // test_for_zero(t4ddg_1(2,2,2,1) - (t4ddg_2(2,2,2,1)*t2s_3(2,1))
1992  // test_for_zero(t4ddg_1(2,2,2,2) - (t4ddg_2(2,2,2,2)*t2s_3(2,2)));
1993  // cout << endl;
1994 
1995  // t4ddg_1(i,j,k,l)=(t2s_2(l,j)&t4ddg_2(i,j,k,l));
1996  // test_for_zero(t4ddg_1(0,0,0,0) - (t4ddg_2(0,0,0,0)*t2s_2(0,0))
1997  // test_for_zero(t4ddg_1(0,0,0,1) - (t4ddg_2(0,0,0,1)*t2s_2(0,1))
1998  // test_for_zero(t4ddg_1(0,0,0,2) - (t4ddg_2(0,0,0,2)*t2s_2(0,2))
1999  // test_for_zero(t4ddg_1(0,0,1,0) - (t4ddg_2(0,0,1,0)*t2s_2(0,0))
2000  // test_for_zero(t4ddg_1(0,0,1,1) - (t4ddg_2(0,0,1,1)*t2s_2(0,1))
2001  // test_for_zero(t4ddg_1(0,0,1,2) - (t4ddg_2(0,0,1,2)*t2s_2(0,2))
2002  // test_for_zero(t4ddg_1(0,0,2,0) - (t4ddg_2(0,0,2,0)*t2s_2(0,0))
2003  // test_for_zero(t4ddg_1(0,0,2,1) - (t4ddg_2(0,0,2,1)*t2s_2(0,1))
2004  // test_for_zero(t4ddg_1(0,0,2,2) - (t4ddg_2(0,0,2,2)*t2s_2(0,2))
2005  // test_for_zero(t4ddg_1(0,1,0,0) - (t4ddg_2(0,1,0,0)*t2s_2(1,0))
2006  // test_for_zero(t4ddg_1(0,1,0,1) - (t4ddg_2(0,1,0,1)*t2s_2(1,1))
2007  // test_for_zero(t4ddg_1(0,1,0,2) - (t4ddg_2(0,1,0,2)*t2s_2(1,2))
2008  // test_for_zero(t4ddg_1(0,1,1,0) - (t4ddg_2(0,1,1,0)*t2s_2(1,0))
2009  // test_for_zero(t4ddg_1(0,1,1,1) - (t4ddg_2(0,1,1,1)*t2s_2(1,1))
2010  // test_for_zero(t4ddg_1(0,1,1,2) - (t4ddg_2(0,1,1,2)*t2s_2(1,2))
2011  // test_for_zero(t4ddg_1(0,1,2,0) - (t4ddg_2(0,1,2,0)*t2s_2(1,0))
2012  // test_for_zero(t4ddg_1(0,1,2,1) - (t4ddg_2(0,1,2,1)*t2s_2(1,1))
2013  // test_for_zero(t4ddg_1(0,1,2,2) - (t4ddg_2(0,1,2,2)*t2s_2(1,2))
2014  // test_for_zero(t4ddg_1(0,2,0,0) - (t4ddg_2(0,2,0,0)*t2s_2(2,0))
2015  // test_for_zero(t4ddg_1(0,2,0,1) - (t4ddg_2(0,2,0,1)*t2s_2(2,1))
2016  // test_for_zero(t4ddg_1(0,2,0,2) - (t4ddg_2(0,2,0,2)*t2s_2(2,2))
2017  // test_for_zero(t4ddg_1(0,2,1,0) - (t4ddg_2(0,2,1,0)*t2s_2(2,0))
2018  // test_for_zero(t4ddg_1(0,2,1,1) - (t4ddg_2(0,2,1,1)*t2s_2(2,1))
2019  // test_for_zero(t4ddg_1(0,2,1,2) - (t4ddg_2(0,2,1,2)*t2s_2(2,2))
2020  // test_for_zero(t4ddg_1(0,2,2,0) - (t4ddg_2(0,2,2,0)*t2s_2(2,0))
2021  // test_for_zero(t4ddg_1(0,2,2,1) - (t4ddg_2(0,2,2,1)*t2s_2(2,1))
2022  // test_for_zero(t4ddg_1(0,2,2,2) - (t4ddg_2(0,2,2,2)*t2s_2(2,2))
2023  // test_for_zero(t4ddg_1(1,0,0,0) - (t4ddg_2(1,0,0,0)*t2s_2(0,0))
2024  // test_for_zero(t4ddg_1(1,0,0,1) - (t4ddg_2(1,0,0,1)*t2s_2(0,1))
2025  // test_for_zero(t4ddg_1(1,0,0,2) - (t4ddg_2(1,0,0,2)*t2s_2(0,2))
2026  // test_for_zero(t4ddg_1(1,0,1,0) - (t4ddg_2(1,0,1,0)*t2s_2(0,0))
2027  // test_for_zero(t4ddg_1(1,0,1,1) - (t4ddg_2(1,0,1,1)*t2s_2(0,1))
2028  // test_for_zero(t4ddg_1(1,0,1,2) - (t4ddg_2(1,0,1,2)*t2s_2(0,2))
2029  // test_for_zero(t4ddg_1(1,0,2,0) - (t4ddg_2(1,0,2,0)*t2s_2(0,0))
2030  // test_for_zero(t4ddg_1(1,0,2,1) - (t4ddg_2(1,0,2,1)*t2s_2(0,1))
2031  // test_for_zero(t4ddg_1(1,0,2,2) - (t4ddg_2(1,0,2,2)*t2s_2(0,2))
2032  // test_for_zero(t4ddg_1(1,1,0,0) - (t4ddg_2(1,1,0,0)*t2s_2(1,0))
2033  // test_for_zero(t4ddg_1(1,1,0,1) - (t4ddg_2(1,1,0,1)*t2s_2(1,1))
2034  // test_for_zero(t4ddg_1(1,1,0,2) - (t4ddg_2(1,1,0,2)*t2s_2(1,2))
2035  // test_for_zero(t4ddg_1(1,1,1,0) - (t4ddg_2(1,1,1,0)*t2s_2(1,0))
2036  // test_for_zero(t4ddg_1(1,1,1,1) - (t4ddg_2(1,1,1,1)*t2s_2(1,1))
2037  // test_for_zero(t4ddg_1(1,1,1,2) - (t4ddg_2(1,1,1,2)*t2s_2(1,2))
2038  // test_for_zero(t4ddg_1(1,1,2,0) - (t4ddg_2(1,1,2,0)*t2s_2(1,0))
2039  // test_for_zero(t4ddg_1(1,1,2,1) - (t4ddg_2(1,1,2,1)*t2s_2(1,1))
2040  // test_for_zero(t4ddg_1(1,1,2,2) - (t4ddg_2(1,1,2,2)*t2s_2(1,2))
2041  // test_for_zero(t4ddg_1(1,2,0,0) - (t4ddg_2(1,2,0,0)*t2s_2(2,0))
2042  // test_for_zero(t4ddg_1(1,2,0,1) - (t4ddg_2(1,2,0,1)*t2s_2(2,1))
2043  // test_for_zero(t4ddg_1(1,2,0,2) - (t4ddg_2(1,2,0,2)*t2s_2(2,2))
2044  // test_for_zero(t4ddg_1(1,2,1,0) - (t4ddg_2(1,2,1,0)*t2s_2(2,0))
2045  // test_for_zero(t4ddg_1(1,2,1,1) - (t4ddg_2(1,2,1,1)*t2s_2(2,1))
2046  // test_for_zero(t4ddg_1(1,2,1,2) - (t4ddg_2(1,2,1,2)*t2s_2(2,2))
2047  // test_for_zero(t4ddg_1(1,2,2,0) - (t4ddg_2(1,2,2,0)*t2s_2(2,0))
2048  // test_for_zero(t4ddg_1(1,2,2,1) - (t4ddg_2(1,2,2,1)*t2s_2(2,1))
2049  // test_for_zero(t4ddg_1(1,2,2,2) - (t4ddg_2(1,2,2,2)*t2s_2(2,2))
2050  // test_for_zero(t4ddg_1(2,0,0,0) - (t4ddg_2(2,0,0,0)*t2s_2(0,0))
2051  // test_for_zero(t4ddg_1(2,0,0,1) - (t4ddg_2(2,0,0,1)*t2s_2(0,1))
2052  // test_for_zero(t4ddg_1(2,0,0,2) - (t4ddg_2(2,0,0,2)*t2s_2(0,2))
2053  // test_for_zero(t4ddg_1(2,0,1,0) - (t4ddg_2(2,0,1,0)*t2s_2(0,0))
2054  // test_for_zero(t4ddg_1(2,0,1,1) - (t4ddg_2(2,0,1,1)*t2s_2(0,1))
2055  // test_for_zero(t4ddg_1(2,0,1,2) - (t4ddg_2(2,0,1,2)*t2s_2(0,2))
2056  // test_for_zero(t4ddg_1(2,0,2,0) - (t4ddg_2(2,0,2,0)*t2s_2(0,0))
2057  // test_for_zero(t4ddg_1(2,0,2,1) - (t4ddg_2(2,0,2,1)*t2s_2(0,1))
2058  // test_for_zero(t4ddg_1(2,0,2,2) - (t4ddg_2(2,0,2,2)*t2s_2(0,2))
2059  // test_for_zero(t4ddg_1(2,1,0,0) - (t4ddg_2(2,1,0,0)*t2s_2(1,0))
2060  // test_for_zero(t4ddg_1(2,1,0,1) - (t4ddg_2(2,1,0,1)*t2s_2(1,1))
2061  // test_for_zero(t4ddg_1(2,1,0,2) - (t4ddg_2(2,1,0,2)*t2s_2(1,2))
2062  // test_for_zero(t4ddg_1(2,1,1,0) - (t4ddg_2(2,1,1,0)*t2s_2(1,0))
2063  // test_for_zero(t4ddg_1(2,1,1,1) - (t4ddg_2(2,1,1,1)*t2s_2(1,1))
2064  // test_for_zero(t4ddg_1(2,1,1,2) - (t4ddg_2(2,1,1,2)*t2s_2(1,2))
2065  // test_for_zero(t4ddg_1(2,1,2,0) - (t4ddg_2(2,1,2,0)*t2s_2(1,0))
2066  // test_for_zero(t4ddg_1(2,1,2,1) - (t4ddg_2(2,1,2,1)*t2s_2(1,1))
2067  // test_for_zero(t4ddg_1(2,1,2,2) - (t4ddg_2(2,1,2,2)*t2s_2(1,2))
2068  // test_for_zero(t4ddg_1(2,2,0,0) - (t4ddg_2(2,2,0,0)*t2s_2(2,0))
2069  // test_for_zero(t4ddg_1(2,2,0,1) - (t4ddg_2(2,2,0,1)*t2s_2(2,1))
2070  // test_for_zero(t4ddg_1(2,2,0,2) - (t4ddg_2(2,2,0,2)*t2s_2(2,2))
2071  // test_for_zero(t4ddg_1(2,2,1,0) - (t4ddg_2(2,2,1,0)*t2s_2(2,0))
2072  // test_for_zero(t4ddg_1(2,2,1,1) - (t4ddg_2(2,2,1,1)*t2s_2(2,1))
2073  // test_for_zero(t4ddg_1(2,2,1,2) - (t4ddg_2(2,2,1,2)*t2s_2(2,2))
2074  // test_for_zero(t4ddg_1(2,2,2,0) - (t4ddg_2(2,2,2,0)*t2s_2(2,0))
2075  // test_for_zero(t4ddg_1(2,2,2,1) - (t4ddg_2(2,2,2,1)*t2s_2(2,1))
2076  // test_for_zero(t4ddg_1(2,2,2,2) - (t4ddg_2(2,2,2,2)*t2s_2(2,2)));
2077 }
FTensor
JSON compatible output.
Definition: Christof_constructor.hpp:6
FTensor::Tensor2_symmetric
Definition: Tensor2_symmetric_value.hpp:13
FTensor::Tensor2< double, 3, 3 >
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 >
std
Definition: enable_if.hpp:5
j
FTensor::Index< 'j', 3 > j
Definition: matrix_function.cpp:19
FTensor::Ddg< double, 3, 3 >
k
FTensor::Index< 'k', 3 > k
Definition: matrix_function.cpp:20
l
FTensor::Index< 'l', 3 > l
Definition: matrix_function.cpp:21
test_T4ddg_03
void test_T4ddg_03(Tensor2< double, 3, 3 > &t2_1, Tensor2_symmetric< double, 3 > &t2s_1, const Tensor2_symmetric< double, 3 > &t2s_2, const Tensor2_symmetric< double, 3 > &t2s_3)
Definition: test_T4ddg_03.cpp:7
Generated by Doxygen 1.8.17 and hosted at University of Glasgow