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