v0.14.0
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>
4 using namespace FTensor;
5 using namespace std;
6 
7 void 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
JSON compatible output.
Definition: Christof_constructor.hpp:6
FTensor::Tensor1< double, 3 >
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
FTensor::Tensor2_symmetric
Definition: Tensor2_symmetric_value.hpp:13
FTensor::Tensor2< double, 3, 3 >
FTensor::Number< 0 >
test_for_zero
void test_for_zero(const T &t, const std::string &s)
Definition: test_for_zero.hpp:7
i
FTensor::Index< 'i', SPACE_DIM > i
Definition: hcurl_divergence_operator_2d.cpp:27
FTensor::Index< 'i', 3 >
FTensor::Dg
Definition: Dg_value.hpp:9
std
Definition: enable_if.hpp:5
j
FTensor::Index< 'j', 3 > j
Definition: matrix_function.cpp:19
FTensor::Ddg< double, 3, 3 >
k
FTensor::Index< 'k', 3 > k
Definition: matrix_function.cpp:20
l
FTensor::Index< 'l', 3 > l
Definition: matrix_function.cpp:21
Generated by Doxygen 1.8.17 and hosted at University of Glasgow