v0.15.0
Loading...
Searching...
No Matches
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>
4using namespace FTensor;
5using namespace std;
6
7void 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::Tensor2_symmetric
Definition Tensor2_symmetric_value.hpp:14
i
FTensor::Index< 'i', SPACE_DIM > i
Definition hcurl_divergence_operator_2d.cpp:27
l
FTensor::Index< 'l', 3 > l
Definition matrix_function.cpp:21
j
FTensor::Index< 'j', 3 > j
Definition matrix_function.cpp:19
k
FTensor::Index< 'k', 3 > k
Definition matrix_function.cpp:20
FTensor
Tensors class implemented by Walter Landry.
Definition FTensor.hpp:51
test_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
test_for_zero
void test_for_zero(const T &t, const std::string &s)
Definition test_for_zero.hpp:7
Generated by Doxygen 1.9.8 and hosted at University of Glasgow