v0.14.0
test_T4ddg_02.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_02(const Tensor2_symmetric<double, 3> &t2s_2)
8 {
9  Index<'i', 3> i;
10  Index<'j', 3> j;
11  Index<'k', 3> k;
12  Index<'l', 3> l;
13  Index<'m', 3> m;
14  Index<'n', 3> n;
15 
16  Number<0> N0;
17  Number<1> N1;
18  Number<2> N2;
19 
20  Ddg<double, 3, 3> t4ddg_1, t4ddg_2, t4ddg_3;
21 
22  t4ddg_2(i, k, j, l) = t2s_2(i, k) * t2s_2(j, l);
23  t4ddg_3(i, l, j, k) = t2s_2(i, l) * t2s_2(j, k);
24  t4ddg_1(i, k, m, n) = (t4ddg_2(i, j, k, l) ^ t4ddg_3(j, l, m, n));
25  test_for_zero(t4ddg_1(0, 0, 0, 0)
26  - (t4ddg_2(0, 0, 0, 0) * t4ddg_3(0, 0, 0, 0)
27  + t4ddg_2(0, 0, 0, 1) * t4ddg_3(0, 1, 0, 0)
28  + t4ddg_2(0, 0, 0, 2) * t4ddg_3(0, 2, 0, 0)
29  + t4ddg_2(0, 1, 0, 0) * t4ddg_3(1, 0, 0, 0)
30  + t4ddg_2(0, 1, 0, 1) * t4ddg_3(1, 1, 0, 0)
31  + t4ddg_2(0, 1, 0, 2) * t4ddg_3(1, 2, 0, 0)
32  + t4ddg_2(0, 2, 0, 0) * t4ddg_3(2, 0, 0, 0)
33  + t4ddg_2(0, 2, 0, 1) * t4ddg_3(2, 1, 0, 0)
34  + t4ddg_2(0, 2, 0, 2) * t4ddg_3(2, 2, 0, 0)),
35  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,0,0,0)");
36  test_for_zero(t4ddg_1(0, 0, 0, 1)
37  - (t4ddg_2(0, 0, 0, 0) * t4ddg_3(0, 0, 0, 1)
38  + t4ddg_2(0, 0, 0, 1) * t4ddg_3(0, 1, 0, 1)
39  + t4ddg_2(0, 0, 0, 2) * t4ddg_3(0, 2, 0, 1)
40  + t4ddg_2(0, 1, 0, 0) * t4ddg_3(1, 0, 0, 1)
41  + t4ddg_2(0, 1, 0, 1) * t4ddg_3(1, 1, 0, 1)
42  + t4ddg_2(0, 1, 0, 2) * t4ddg_3(1, 2, 0, 1)
43  + t4ddg_2(0, 2, 0, 0) * t4ddg_3(2, 0, 0, 1)
44  + t4ddg_2(0, 2, 0, 1) * t4ddg_3(2, 1, 0, 1)
45  + t4ddg_2(0, 2, 0, 2) * t4ddg_3(2, 2, 0, 1)),
46  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,0,0,1)");
47  test_for_zero(t4ddg_1(0, 0, 0, 2)
48  - (t4ddg_2(0, 0, 0, 0) * t4ddg_3(0, 0, 0, 2)
49  + t4ddg_2(0, 0, 0, 1) * t4ddg_3(0, 1, 0, 2)
50  + t4ddg_2(0, 0, 0, 2) * t4ddg_3(0, 2, 0, 2)
51  + t4ddg_2(0, 1, 0, 0) * t4ddg_3(1, 0, 0, 2)
52  + t4ddg_2(0, 1, 0, 1) * t4ddg_3(1, 1, 0, 2)
53  + t4ddg_2(0, 1, 0, 2) * t4ddg_3(1, 2, 0, 2)
54  + t4ddg_2(0, 2, 0, 0) * t4ddg_3(2, 0, 0, 2)
55  + t4ddg_2(0, 2, 0, 1) * t4ddg_3(2, 1, 0, 2)
56  + t4ddg_2(0, 2, 0, 2) * t4ddg_3(2, 2, 0, 2)),
57  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,0,0,2)");
58  test_for_zero(t4ddg_1(0, 1, 0, 0)
59  - (t4ddg_2(0, 0, 1, 0) * t4ddg_3(0, 0, 0, 0)
60  + t4ddg_2(0, 0, 1, 1) * t4ddg_3(0, 1, 0, 0)
61  + t4ddg_2(0, 0, 1, 2) * t4ddg_3(0, 2, 0, 0)
62  + t4ddg_2(0, 1, 1, 0) * t4ddg_3(1, 0, 0, 0)
63  + t4ddg_2(0, 1, 1, 1) * t4ddg_3(1, 1, 0, 0)
64  + t4ddg_2(0, 1, 1, 2) * t4ddg_3(1, 2, 0, 0)
65  + t4ddg_2(0, 2, 1, 0) * t4ddg_3(2, 0, 0, 0)
66  + t4ddg_2(0, 2, 1, 1) * t4ddg_3(2, 1, 0, 0)
67  + t4ddg_2(0, 2, 1, 2) * t4ddg_3(2, 2, 0, 0)),
68  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,0,1,0)");
69  test_for_zero(t4ddg_1(0, 1, 0, 1)
70  - (t4ddg_2(0, 0, 1, 0) * t4ddg_3(0, 0, 0, 1)
71  + t4ddg_2(0, 0, 1, 1) * t4ddg_3(0, 1, 0, 1)
72  + t4ddg_2(0, 0, 1, 2) * t4ddg_3(0, 2, 0, 1)
73  + t4ddg_2(0, 1, 1, 0) * t4ddg_3(1, 0, 0, 1)
74  + t4ddg_2(0, 1, 1, 1) * t4ddg_3(1, 1, 0, 1)
75  + t4ddg_2(0, 1, 1, 2) * t4ddg_3(1, 2, 0, 1)
76  + t4ddg_2(0, 2, 1, 0) * t4ddg_3(2, 0, 0, 1)
77  + t4ddg_2(0, 2, 1, 1) * t4ddg_3(2, 1, 0, 1)
78  + t4ddg_2(0, 2, 1, 2) * t4ddg_3(2, 2, 0, 1)),
79  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,0,1,1)");
80  test_for_zero(t4ddg_1(0, 1, 0, 2)
81  - (t4ddg_2(0, 0, 1, 0) * t4ddg_3(0, 0, 0, 2)
82  + t4ddg_2(0, 0, 1, 1) * t4ddg_3(0, 1, 0, 2)
83  + t4ddg_2(0, 0, 1, 2) * t4ddg_3(0, 2, 0, 2)
84  + t4ddg_2(0, 1, 1, 0) * t4ddg_3(1, 0, 0, 2)
85  + t4ddg_2(0, 1, 1, 1) * t4ddg_3(1, 1, 0, 2)
86  + t4ddg_2(0, 1, 1, 2) * t4ddg_3(1, 2, 0, 2)
87  + t4ddg_2(0, 2, 1, 0) * t4ddg_3(2, 0, 0, 2)
88  + t4ddg_2(0, 2, 1, 1) * t4ddg_3(2, 1, 0, 2)
89  + t4ddg_2(0, 2, 1, 2) * t4ddg_3(2, 2, 0, 2)),
90  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,0,1,2)");
91  test_for_zero(t4ddg_1(0, 2, 0, 0)
92  - (t4ddg_2(0, 0, 2, 0) * t4ddg_3(0, 0, 0, 0)
93  + t4ddg_2(0, 0, 2, 1) * t4ddg_3(0, 1, 0, 0)
94  + t4ddg_2(0, 0, 2, 2) * t4ddg_3(0, 2, 0, 0)
95  + t4ddg_2(0, 1, 2, 0) * t4ddg_3(1, 0, 0, 0)
96  + t4ddg_2(0, 1, 2, 1) * t4ddg_3(1, 1, 0, 0)
97  + t4ddg_2(0, 1, 2, 2) * t4ddg_3(1, 2, 0, 0)
98  + t4ddg_2(0, 2, 2, 0) * t4ddg_3(2, 0, 0, 0)
99  + t4ddg_2(0, 2, 2, 1) * t4ddg_3(2, 1, 0, 0)
100  + t4ddg_2(0, 2, 2, 2) * t4ddg_3(2, 2, 0, 0)),
101  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,0,2,0)");
102  test_for_zero(t4ddg_1(0, 2, 0, 1)
103  - (t4ddg_2(0, 0, 2, 0) * t4ddg_3(0, 0, 0, 1)
104  + t4ddg_2(0, 0, 2, 1) * t4ddg_3(0, 1, 0, 1)
105  + t4ddg_2(0, 0, 2, 2) * t4ddg_3(0, 2, 0, 1)
106  + t4ddg_2(0, 1, 2, 0) * t4ddg_3(1, 0, 0, 1)
107  + t4ddg_2(0, 1, 2, 1) * t4ddg_3(1, 1, 0, 1)
108  + t4ddg_2(0, 1, 2, 2) * t4ddg_3(1, 2, 0, 1)
109  + t4ddg_2(0, 2, 2, 0) * t4ddg_3(2, 0, 0, 1)
110  + t4ddg_2(0, 2, 2, 1) * t4ddg_3(2, 1, 0, 1)
111  + t4ddg_2(0, 2, 2, 2) * t4ddg_3(2, 2, 0, 1)),
112  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,0,2,1)");
113  test_for_zero(t4ddg_1(0, 2, 0, 2)
114  - (t4ddg_2(0, 0, 2, 0) * t4ddg_3(0, 0, 0, 2)
115  + t4ddg_2(0, 0, 2, 1) * t4ddg_3(0, 1, 0, 2)
116  + t4ddg_2(0, 0, 2, 2) * t4ddg_3(0, 2, 0, 2)
117  + t4ddg_2(0, 1, 2, 0) * t4ddg_3(1, 0, 0, 2)
118  + t4ddg_2(0, 1, 2, 1) * t4ddg_3(1, 1, 0, 2)
119  + t4ddg_2(0, 1, 2, 2) * t4ddg_3(1, 2, 0, 2)
120  + t4ddg_2(0, 2, 2, 0) * t4ddg_3(2, 0, 0, 2)
121  + t4ddg_2(0, 2, 2, 1) * t4ddg_3(2, 1, 0, 2)
122  + t4ddg_2(0, 2, 2, 2) * t4ddg_3(2, 2, 0, 2)),
123  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,0,2,2)");
124  test_for_zero(t4ddg_1(0, 0, 1, 0)
125  - (t4ddg_2(0, 0, 0, 0) * t4ddg_3(0, 0, 1, 0)
126  + t4ddg_2(0, 0, 0, 1) * t4ddg_3(0, 1, 1, 0)
127  + t4ddg_2(0, 0, 0, 2) * t4ddg_3(0, 2, 1, 0)
128  + t4ddg_2(0, 1, 0, 0) * t4ddg_3(1, 0, 1, 0)
129  + t4ddg_2(0, 1, 0, 1) * t4ddg_3(1, 1, 1, 0)
130  + t4ddg_2(0, 1, 0, 2) * t4ddg_3(1, 2, 1, 0)
131  + t4ddg_2(0, 2, 0, 0) * t4ddg_3(2, 0, 1, 0)
132  + t4ddg_2(0, 2, 0, 1) * t4ddg_3(2, 1, 1, 0)
133  + t4ddg_2(0, 2, 0, 2) * t4ddg_3(2, 2, 1, 0)),
134  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,1,0,0)");
135  test_for_zero(t4ddg_1(0, 0, 1, 1)
136  - (t4ddg_2(0, 0, 0, 0) * t4ddg_3(0, 0, 1, 1)
137  + t4ddg_2(0, 0, 0, 1) * t4ddg_3(0, 1, 1, 1)
138  + t4ddg_2(0, 0, 0, 2) * t4ddg_3(0, 2, 1, 1)
139  + t4ddg_2(0, 1, 0, 0) * t4ddg_3(1, 0, 1, 1)
140  + t4ddg_2(0, 1, 0, 1) * t4ddg_3(1, 1, 1, 1)
141  + t4ddg_2(0, 1, 0, 2) * t4ddg_3(1, 2, 1, 1)
142  + t4ddg_2(0, 2, 0, 0) * t4ddg_3(2, 0, 1, 1)
143  + t4ddg_2(0, 2, 0, 1) * t4ddg_3(2, 1, 1, 1)
144  + t4ddg_2(0, 2, 0, 2) * t4ddg_3(2, 2, 1, 1)),
145  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,1,0,1)");
146  test_for_zero(t4ddg_1(0, 0, 1, 2)
147  - (t4ddg_2(0, 0, 0, 0) * t4ddg_3(0, 0, 1, 2)
148  + t4ddg_2(0, 0, 0, 1) * t4ddg_3(0, 1, 1, 2)
149  + t4ddg_2(0, 0, 0, 2) * t4ddg_3(0, 2, 1, 2)
150  + t4ddg_2(0, 1, 0, 0) * t4ddg_3(1, 0, 1, 2)
151  + t4ddg_2(0, 1, 0, 1) * t4ddg_3(1, 1, 1, 2)
152  + t4ddg_2(0, 1, 0, 2) * t4ddg_3(1, 2, 1, 2)
153  + t4ddg_2(0, 2, 0, 0) * t4ddg_3(2, 0, 1, 2)
154  + t4ddg_2(0, 2, 0, 1) * t4ddg_3(2, 1, 1, 2)
155  + t4ddg_2(0, 2, 0, 2) * t4ddg_3(2, 2, 1, 2)),
156  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,1,0,2)");
157  test_for_zero(t4ddg_1(0, 1, 1, 0)
158  - (t4ddg_2(0, 0, 1, 0) * t4ddg_3(0, 0, 1, 0)
159  + t4ddg_2(0, 0, 1, 1) * t4ddg_3(0, 1, 1, 0)
160  + t4ddg_2(0, 0, 1, 2) * t4ddg_3(0, 2, 1, 0)
161  + t4ddg_2(0, 1, 1, 0) * t4ddg_3(1, 0, 1, 0)
162  + t4ddg_2(0, 1, 1, 1) * t4ddg_3(1, 1, 1, 0)
163  + t4ddg_2(0, 1, 1, 2) * t4ddg_3(1, 2, 1, 0)
164  + t4ddg_2(0, 2, 1, 0) * t4ddg_3(2, 0, 1, 0)
165  + t4ddg_2(0, 2, 1, 1) * t4ddg_3(2, 1, 1, 0)
166  + t4ddg_2(0, 2, 1, 2) * t4ddg_3(2, 2, 1, 0)),
167  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,1,1,0)");
168  test_for_zero(t4ddg_1(0, 1, 1, 1)
169  - (t4ddg_2(0, 0, 1, 0) * t4ddg_3(0, 0, 1, 1)
170  + t4ddg_2(0, 0, 1, 1) * t4ddg_3(0, 1, 1, 1)
171  + t4ddg_2(0, 0, 1, 2) * t4ddg_3(0, 2, 1, 1)
172  + t4ddg_2(0, 1, 1, 0) * t4ddg_3(1, 0, 1, 1)
173  + t4ddg_2(0, 1, 1, 1) * t4ddg_3(1, 1, 1, 1)
174  + t4ddg_2(0, 1, 1, 2) * t4ddg_3(1, 2, 1, 1)
175  + t4ddg_2(0, 2, 1, 0) * t4ddg_3(2, 0, 1, 1)
176  + t4ddg_2(0, 2, 1, 1) * t4ddg_3(2, 1, 1, 1)
177  + t4ddg_2(0, 2, 1, 2) * t4ddg_3(2, 2, 1, 1)),
178  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,1,1,1)");
179  test_for_zero(t4ddg_1(0, 1, 1, 2)
180  - (t4ddg_2(0, 0, 1, 0) * t4ddg_3(0, 0, 1, 2)
181  + t4ddg_2(0, 0, 1, 1) * t4ddg_3(0, 1, 1, 2)
182  + t4ddg_2(0, 0, 1, 2) * t4ddg_3(0, 2, 1, 2)
183  + t4ddg_2(0, 1, 1, 0) * t4ddg_3(1, 0, 1, 2)
184  + t4ddg_2(0, 1, 1, 1) * t4ddg_3(1, 1, 1, 2)
185  + t4ddg_2(0, 1, 1, 2) * t4ddg_3(1, 2, 1, 2)
186  + t4ddg_2(0, 2, 1, 0) * t4ddg_3(2, 0, 1, 2)
187  + t4ddg_2(0, 2, 1, 1) * t4ddg_3(2, 1, 1, 2)
188  + t4ddg_2(0, 2, 1, 2) * t4ddg_3(2, 2, 1, 2)),
189  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,1,1,2)");
190  test_for_zero(t4ddg_1(0, 2, 1, 0)
191  - (t4ddg_2(0, 0, 2, 0) * t4ddg_3(0, 0, 1, 0)
192  + t4ddg_2(0, 0, 2, 1) * t4ddg_3(0, 1, 1, 0)
193  + t4ddg_2(0, 0, 2, 2) * t4ddg_3(0, 2, 1, 0)
194  + t4ddg_2(0, 1, 2, 0) * t4ddg_3(1, 0, 1, 0)
195  + t4ddg_2(0, 1, 2, 1) * t4ddg_3(1, 1, 1, 0)
196  + t4ddg_2(0, 1, 2, 2) * t4ddg_3(1, 2, 1, 0)
197  + t4ddg_2(0, 2, 2, 0) * t4ddg_3(2, 0, 1, 0)
198  + t4ddg_2(0, 2, 2, 1) * t4ddg_3(2, 1, 1, 0)
199  + t4ddg_2(0, 2, 2, 2) * t4ddg_3(2, 2, 1, 0)),
200  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,1,2,0)");
201  test_for_zero(t4ddg_1(0, 2, 1, 1)
202  - (t4ddg_2(0, 0, 2, 0) * t4ddg_3(0, 0, 1, 1)
203  + t4ddg_2(0, 0, 2, 1) * t4ddg_3(0, 1, 1, 1)
204  + t4ddg_2(0, 0, 2, 2) * t4ddg_3(0, 2, 1, 1)
205  + t4ddg_2(0, 1, 2, 0) * t4ddg_3(1, 0, 1, 1)
206  + t4ddg_2(0, 1, 2, 1) * t4ddg_3(1, 1, 1, 1)
207  + t4ddg_2(0, 1, 2, 2) * t4ddg_3(1, 2, 1, 1)
208  + t4ddg_2(0, 2, 2, 0) * t4ddg_3(2, 0, 1, 1)
209  + t4ddg_2(0, 2, 2, 1) * t4ddg_3(2, 1, 1, 1)
210  + t4ddg_2(0, 2, 2, 2) * t4ddg_3(2, 2, 1, 1)),
211  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,1,2,1)");
212  test_for_zero(t4ddg_1(0, 2, 1, 2)
213  - (t4ddg_2(0, 0, 2, 0) * t4ddg_3(0, 0, 1, 2)
214  + t4ddg_2(0, 0, 2, 1) * t4ddg_3(0, 1, 1, 2)
215  + t4ddg_2(0, 0, 2, 2) * t4ddg_3(0, 2, 1, 2)
216  + t4ddg_2(0, 1, 2, 0) * t4ddg_3(1, 0, 1, 2)
217  + t4ddg_2(0, 1, 2, 1) * t4ddg_3(1, 1, 1, 2)
218  + t4ddg_2(0, 1, 2, 2) * t4ddg_3(1, 2, 1, 2)
219  + t4ddg_2(0, 2, 2, 0) * t4ddg_3(2, 0, 1, 2)
220  + t4ddg_2(0, 2, 2, 1) * t4ddg_3(2, 1, 1, 2)
221  + t4ddg_2(0, 2, 2, 2) * t4ddg_3(2, 2, 1, 2)),
222  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,1,2,2)");
223  test_for_zero(t4ddg_1(0, 0, 2, 0)
224  - (t4ddg_2(0, 0, 0, 0) * t4ddg_3(0, 0, 2, 0)
225  + t4ddg_2(0, 0, 0, 1) * t4ddg_3(0, 1, 2, 0)
226  + t4ddg_2(0, 0, 0, 2) * t4ddg_3(0, 2, 2, 0)
227  + t4ddg_2(0, 1, 0, 0) * t4ddg_3(1, 0, 2, 0)
228  + t4ddg_2(0, 1, 0, 1) * t4ddg_3(1, 1, 2, 0)
229  + t4ddg_2(0, 1, 0, 2) * t4ddg_3(1, 2, 2, 0)
230  + t4ddg_2(0, 2, 0, 0) * t4ddg_3(2, 0, 2, 0)
231  + t4ddg_2(0, 2, 0, 1) * t4ddg_3(2, 1, 2, 0)
232  + t4ddg_2(0, 2, 0, 2) * t4ddg_3(2, 2, 2, 0)),
233  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,2,0,0)");
234  test_for_zero(t4ddg_1(0, 0, 2, 1)
235  - (t4ddg_2(0, 0, 0, 0) * t4ddg_3(0, 0, 2, 1)
236  + t4ddg_2(0, 0, 0, 1) * t4ddg_3(0, 1, 2, 1)
237  + t4ddg_2(0, 0, 0, 2) * t4ddg_3(0, 2, 2, 1)
238  + t4ddg_2(0, 1, 0, 0) * t4ddg_3(1, 0, 2, 1)
239  + t4ddg_2(0, 1, 0, 1) * t4ddg_3(1, 1, 2, 1)
240  + t4ddg_2(0, 1, 0, 2) * t4ddg_3(1, 2, 2, 1)
241  + t4ddg_2(0, 2, 0, 0) * t4ddg_3(2, 0, 2, 1)
242  + t4ddg_2(0, 2, 0, 1) * t4ddg_3(2, 1, 2, 1)
243  + t4ddg_2(0, 2, 0, 2) * t4ddg_3(2, 2, 2, 1)),
244  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,2,0,1)");
245  test_for_zero(t4ddg_1(0, 0, 2, 2)
246  - (t4ddg_2(0, 0, 0, 0) * t4ddg_3(0, 0, 2, 2)
247  + t4ddg_2(0, 0, 0, 1) * t4ddg_3(0, 1, 2, 2)
248  + t4ddg_2(0, 0, 0, 2) * t4ddg_3(0, 2, 2, 2)
249  + t4ddg_2(0, 1, 0, 0) * t4ddg_3(1, 0, 2, 2)
250  + t4ddg_2(0, 1, 0, 1) * t4ddg_3(1, 1, 2, 2)
251  + t4ddg_2(0, 1, 0, 2) * t4ddg_3(1, 2, 2, 2)
252  + t4ddg_2(0, 2, 0, 0) * t4ddg_3(2, 0, 2, 2)
253  + t4ddg_2(0, 2, 0, 1) * t4ddg_3(2, 1, 2, 2)
254  + t4ddg_2(0, 2, 0, 2) * t4ddg_3(2, 2, 2, 2)),
255  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,2,0,2)");
256  test_for_zero(t4ddg_1(0, 1, 2, 0)
257  - (t4ddg_2(0, 0, 1, 0) * t4ddg_3(0, 0, 2, 0)
258  + t4ddg_2(0, 0, 1, 1) * t4ddg_3(0, 1, 2, 0)
259  + t4ddg_2(0, 0, 1, 2) * t4ddg_3(0, 2, 2, 0)
260  + t4ddg_2(0, 1, 1, 0) * t4ddg_3(1, 0, 2, 0)
261  + t4ddg_2(0, 1, 1, 1) * t4ddg_3(1, 1, 2, 0)
262  + t4ddg_2(0, 1, 1, 2) * t4ddg_3(1, 2, 2, 0)
263  + t4ddg_2(0, 2, 1, 0) * t4ddg_3(2, 0, 2, 0)
264  + t4ddg_2(0, 2, 1, 1) * t4ddg_3(2, 1, 2, 0)
265  + t4ddg_2(0, 2, 1, 2) * t4ddg_3(2, 2, 2, 0)),
266  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,2,1,0)");
267  test_for_zero(t4ddg_1(0, 1, 2, 1)
268  - (t4ddg_2(0, 0, 1, 0) * t4ddg_3(0, 0, 2, 1)
269  + t4ddg_2(0, 0, 1, 1) * t4ddg_3(0, 1, 2, 1)
270  + t4ddg_2(0, 0, 1, 2) * t4ddg_3(0, 2, 2, 1)
271  + t4ddg_2(0, 1, 1, 0) * t4ddg_3(1, 0, 2, 1)
272  + t4ddg_2(0, 1, 1, 1) * t4ddg_3(1, 1, 2, 1)
273  + t4ddg_2(0, 1, 1, 2) * t4ddg_3(1, 2, 2, 1)
274  + t4ddg_2(0, 2, 1, 0) * t4ddg_3(2, 0, 2, 1)
275  + t4ddg_2(0, 2, 1, 1) * t4ddg_3(2, 1, 2, 1)
276  + t4ddg_2(0, 2, 1, 2) * t4ddg_3(2, 2, 2, 1)),
277  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,2,1,1)");
278  test_for_zero(t4ddg_1(0, 1, 2, 2)
279  - (t4ddg_2(0, 0, 1, 0) * t4ddg_3(0, 0, 2, 2)
280  + t4ddg_2(0, 0, 1, 1) * t4ddg_3(0, 1, 2, 2)
281  + t4ddg_2(0, 0, 1, 2) * t4ddg_3(0, 2, 2, 2)
282  + t4ddg_2(0, 1, 1, 0) * t4ddg_3(1, 0, 2, 2)
283  + t4ddg_2(0, 1, 1, 1) * t4ddg_3(1, 1, 2, 2)
284  + t4ddg_2(0, 1, 1, 2) * t4ddg_3(1, 2, 2, 2)
285  + t4ddg_2(0, 2, 1, 0) * t4ddg_3(2, 0, 2, 2)
286  + t4ddg_2(0, 2, 1, 1) * t4ddg_3(2, 1, 2, 2)
287  + t4ddg_2(0, 2, 1, 2) * t4ddg_3(2, 2, 2, 2)),
288  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,2,1,2)");
289  test_for_zero(t4ddg_1(0, 2, 2, 0)
290  - (t4ddg_2(0, 0, 2, 0) * t4ddg_3(0, 0, 2, 0)
291  + t4ddg_2(0, 0, 2, 1) * t4ddg_3(0, 1, 2, 0)
292  + t4ddg_2(0, 0, 2, 2) * t4ddg_3(0, 2, 2, 0)
293  + t4ddg_2(0, 1, 2, 0) * t4ddg_3(1, 0, 2, 0)
294  + t4ddg_2(0, 1, 2, 1) * t4ddg_3(1, 1, 2, 0)
295  + t4ddg_2(0, 1, 2, 2) * t4ddg_3(1, 2, 2, 0)
296  + t4ddg_2(0, 2, 2, 0) * t4ddg_3(2, 0, 2, 0)
297  + t4ddg_2(0, 2, 2, 1) * t4ddg_3(2, 1, 2, 0)
298  + t4ddg_2(0, 2, 2, 2) * t4ddg_3(2, 2, 2, 0)),
299  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,2,2,0)");
300  test_for_zero(t4ddg_1(0, 2, 2, 1)
301  - (t4ddg_2(0, 0, 2, 0) * t4ddg_3(0, 0, 2, 1)
302  + t4ddg_2(0, 0, 2, 1) * t4ddg_3(0, 1, 2, 1)
303  + t4ddg_2(0, 0, 2, 2) * t4ddg_3(0, 2, 2, 1)
304  + t4ddg_2(0, 1, 2, 0) * t4ddg_3(1, 0, 2, 1)
305  + t4ddg_2(0, 1, 2, 1) * t4ddg_3(1, 1, 2, 1)
306  + t4ddg_2(0, 1, 2, 2) * t4ddg_3(1, 2, 2, 1)
307  + t4ddg_2(0, 2, 2, 0) * t4ddg_3(2, 0, 2, 1)
308  + t4ddg_2(0, 2, 2, 1) * t4ddg_3(2, 1, 2, 1)
309  + t4ddg_2(0, 2, 2, 2) * t4ddg_3(2, 2, 2, 1)),
310  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,2,2,1)");
311  test_for_zero(t4ddg_1(0, 2, 2, 2)
312  - (t4ddg_2(0, 0, 2, 0) * t4ddg_3(0, 0, 2, 2)
313  + t4ddg_2(0, 0, 2, 1) * t4ddg_3(0, 1, 2, 2)
314  + t4ddg_2(0, 0, 2, 2) * t4ddg_3(0, 2, 2, 2)
315  + t4ddg_2(0, 1, 2, 0) * t4ddg_3(1, 0, 2, 2)
316  + t4ddg_2(0, 1, 2, 1) * t4ddg_3(1, 1, 2, 2)
317  + t4ddg_2(0, 1, 2, 2) * t4ddg_3(1, 2, 2, 2)
318  + t4ddg_2(0, 2, 2, 0) * t4ddg_3(2, 0, 2, 2)
319  + t4ddg_2(0, 2, 2, 1) * t4ddg_3(2, 1, 2, 2)
320  + t4ddg_2(0, 2, 2, 2) * t4ddg_3(2, 2, 2, 2)),
321  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,2,2,2)");
322  test_for_zero(t4ddg_1(1, 0, 0, 0)
323  - (t4ddg_2(1, 0, 0, 0) * t4ddg_3(0, 0, 0, 0)
324  + t4ddg_2(1, 0, 0, 1) * t4ddg_3(0, 1, 0, 0)
325  + t4ddg_2(1, 0, 0, 2) * t4ddg_3(0, 2, 0, 0)
326  + t4ddg_2(1, 1, 0, 0) * t4ddg_3(1, 0, 0, 0)
327  + t4ddg_2(1, 1, 0, 1) * t4ddg_3(1, 1, 0, 0)
328  + t4ddg_2(1, 1, 0, 2) * t4ddg_3(1, 2, 0, 0)
329  + t4ddg_2(1, 2, 0, 0) * t4ddg_3(2, 0, 0, 0)
330  + t4ddg_2(1, 2, 0, 1) * t4ddg_3(2, 1, 0, 0)
331  + t4ddg_2(1, 2, 0, 2) * t4ddg_3(2, 2, 0, 0)),
332  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,0,0,0)");
333  test_for_zero(t4ddg_1(1, 0, 0, 1)
334  - (t4ddg_2(1, 0, 0, 0) * t4ddg_3(0, 0, 0, 1)
335  + t4ddg_2(1, 0, 0, 1) * t4ddg_3(0, 1, 0, 1)
336  + t4ddg_2(1, 0, 0, 2) * t4ddg_3(0, 2, 0, 1)
337  + t4ddg_2(1, 1, 0, 0) * t4ddg_3(1, 0, 0, 1)
338  + t4ddg_2(1, 1, 0, 1) * t4ddg_3(1, 1, 0, 1)
339  + t4ddg_2(1, 1, 0, 2) * t4ddg_3(1, 2, 0, 1)
340  + t4ddg_2(1, 2, 0, 0) * t4ddg_3(2, 0, 0, 1)
341  + t4ddg_2(1, 2, 0, 1) * t4ddg_3(2, 1, 0, 1)
342  + t4ddg_2(1, 2, 0, 2) * t4ddg_3(2, 2, 0, 1)),
343  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,0,0,1)");
344  test_for_zero(t4ddg_1(1, 0, 0, 2)
345  - (t4ddg_2(1, 0, 0, 0) * t4ddg_3(0, 0, 0, 2)
346  + t4ddg_2(1, 0, 0, 1) * t4ddg_3(0, 1, 0, 2)
347  + t4ddg_2(1, 0, 0, 2) * t4ddg_3(0, 2, 0, 2)
348  + t4ddg_2(1, 1, 0, 0) * t4ddg_3(1, 0, 0, 2)
349  + t4ddg_2(1, 1, 0, 1) * t4ddg_3(1, 1, 0, 2)
350  + t4ddg_2(1, 1, 0, 2) * t4ddg_3(1, 2, 0, 2)
351  + t4ddg_2(1, 2, 0, 0) * t4ddg_3(2, 0, 0, 2)
352  + t4ddg_2(1, 2, 0, 1) * t4ddg_3(2, 1, 0, 2)
353  + t4ddg_2(1, 2, 0, 2) * t4ddg_3(2, 2, 0, 2)),
354  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,0,0,2)");
355  test_for_zero(t4ddg_1(1, 1, 0, 0)
356  - (t4ddg_2(1, 0, 1, 0) * t4ddg_3(0, 0, 0, 0)
357  + t4ddg_2(1, 0, 1, 1) * t4ddg_3(0, 1, 0, 0)
358  + t4ddg_2(1, 0, 1, 2) * t4ddg_3(0, 2, 0, 0)
359  + t4ddg_2(1, 1, 1, 0) * t4ddg_3(1, 0, 0, 0)
360  + t4ddg_2(1, 1, 1, 1) * t4ddg_3(1, 1, 0, 0)
361  + t4ddg_2(1, 1, 1, 2) * t4ddg_3(1, 2, 0, 0)
362  + t4ddg_2(1, 2, 1, 0) * t4ddg_3(2, 0, 0, 0)
363  + t4ddg_2(1, 2, 1, 1) * t4ddg_3(2, 1, 0, 0)
364  + t4ddg_2(1, 2, 1, 2) * t4ddg_3(2, 2, 0, 0)),
365  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,0,1,0)");
366  test_for_zero(t4ddg_1(1, 1, 0, 1)
367  - (t4ddg_2(1, 0, 1, 0) * t4ddg_3(0, 0, 0, 1)
368  + t4ddg_2(1, 0, 1, 1) * t4ddg_3(0, 1, 0, 1)
369  + t4ddg_2(1, 0, 1, 2) * t4ddg_3(0, 2, 0, 1)
370  + t4ddg_2(1, 1, 1, 0) * t4ddg_3(1, 0, 0, 1)
371  + t4ddg_2(1, 1, 1, 1) * t4ddg_3(1, 1, 0, 1)
372  + t4ddg_2(1, 1, 1, 2) * t4ddg_3(1, 2, 0, 1)
373  + t4ddg_2(1, 2, 1, 0) * t4ddg_3(2, 0, 0, 1)
374  + t4ddg_2(1, 2, 1, 1) * t4ddg_3(2, 1, 0, 1)
375  + t4ddg_2(1, 2, 1, 2) * t4ddg_3(2, 2, 0, 1)),
376  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,0,1,1)");
377  test_for_zero(t4ddg_1(1, 1, 0, 2)
378  - (t4ddg_2(1, 0, 1, 0) * t4ddg_3(0, 0, 0, 2)
379  + t4ddg_2(1, 0, 1, 1) * t4ddg_3(0, 1, 0, 2)
380  + t4ddg_2(1, 0, 1, 2) * t4ddg_3(0, 2, 0, 2)
381  + t4ddg_2(1, 1, 1, 0) * t4ddg_3(1, 0, 0, 2)
382  + t4ddg_2(1, 1, 1, 1) * t4ddg_3(1, 1, 0, 2)
383  + t4ddg_2(1, 1, 1, 2) * t4ddg_3(1, 2, 0, 2)
384  + t4ddg_2(1, 2, 1, 0) * t4ddg_3(2, 0, 0, 2)
385  + t4ddg_2(1, 2, 1, 1) * t4ddg_3(2, 1, 0, 2)
386  + t4ddg_2(1, 2, 1, 2) * t4ddg_3(2, 2, 0, 2)),
387  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,0,1,2)");
388  test_for_zero(t4ddg_1(1, 2, 0, 0)
389  - (t4ddg_2(1, 0, 2, 0) * t4ddg_3(0, 0, 0, 0)
390  + t4ddg_2(1, 0, 2, 1) * t4ddg_3(0, 1, 0, 0)
391  + t4ddg_2(1, 0, 2, 2) * t4ddg_3(0, 2, 0, 0)
392  + t4ddg_2(1, 1, 2, 0) * t4ddg_3(1, 0, 0, 0)
393  + t4ddg_2(1, 1, 2, 1) * t4ddg_3(1, 1, 0, 0)
394  + t4ddg_2(1, 1, 2, 2) * t4ddg_3(1, 2, 0, 0)
395  + t4ddg_2(1, 2, 2, 0) * t4ddg_3(2, 0, 0, 0)
396  + t4ddg_2(1, 2, 2, 1) * t4ddg_3(2, 1, 0, 0)
397  + t4ddg_2(1, 2, 2, 2) * t4ddg_3(2, 2, 0, 0)),
398  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,0,2,0)");
399  test_for_zero(t4ddg_1(1, 2, 0, 1)
400  - (t4ddg_2(1, 0, 2, 0) * t4ddg_3(0, 0, 0, 1)
401  + t4ddg_2(1, 0, 2, 1) * t4ddg_3(0, 1, 0, 1)
402  + t4ddg_2(1, 0, 2, 2) * t4ddg_3(0, 2, 0, 1)
403  + t4ddg_2(1, 1, 2, 0) * t4ddg_3(1, 0, 0, 1)
404  + t4ddg_2(1, 1, 2, 1) * t4ddg_3(1, 1, 0, 1)
405  + t4ddg_2(1, 1, 2, 2) * t4ddg_3(1, 2, 0, 1)
406  + t4ddg_2(1, 2, 2, 0) * t4ddg_3(2, 0, 0, 1)
407  + t4ddg_2(1, 2, 2, 1) * t4ddg_3(2, 1, 0, 1)
408  + t4ddg_2(1, 2, 2, 2) * t4ddg_3(2, 2, 0, 1)),
409  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,0,2,1)");
410  test_for_zero(t4ddg_1(1, 2, 0, 2)
411  - (t4ddg_2(1, 0, 2, 0) * t4ddg_3(0, 0, 0, 2)
412  + t4ddg_2(1, 0, 2, 1) * t4ddg_3(0, 1, 0, 2)
413  + t4ddg_2(1, 0, 2, 2) * t4ddg_3(0, 2, 0, 2)
414  + t4ddg_2(1, 1, 2, 0) * t4ddg_3(1, 0, 0, 2)
415  + t4ddg_2(1, 1, 2, 1) * t4ddg_3(1, 1, 0, 2)
416  + t4ddg_2(1, 1, 2, 2) * t4ddg_3(1, 2, 0, 2)
417  + t4ddg_2(1, 2, 2, 0) * t4ddg_3(2, 0, 0, 2)
418  + t4ddg_2(1, 2, 2, 1) * t4ddg_3(2, 1, 0, 2)
419  + t4ddg_2(1, 2, 2, 2) * t4ddg_3(2, 2, 0, 2)),
420  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,0,2,2)");
421  test_for_zero(t4ddg_1(1, 0, 1, 0)
422  - (t4ddg_2(1, 0, 0, 0) * t4ddg_3(0, 0, 1, 0)
423  + t4ddg_2(1, 0, 0, 1) * t4ddg_3(0, 1, 1, 0)
424  + t4ddg_2(1, 0, 0, 2) * t4ddg_3(0, 2, 1, 0)
425  + t4ddg_2(1, 1, 0, 0) * t4ddg_3(1, 0, 1, 0)
426  + t4ddg_2(1, 1, 0, 1) * t4ddg_3(1, 1, 1, 0)
427  + t4ddg_2(1, 1, 0, 2) * t4ddg_3(1, 2, 1, 0)
428  + t4ddg_2(1, 2, 0, 0) * t4ddg_3(2, 0, 1, 0)
429  + t4ddg_2(1, 2, 0, 1) * t4ddg_3(2, 1, 1, 0)
430  + t4ddg_2(1, 2, 0, 2) * t4ddg_3(2, 2, 1, 0)),
431  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,1,0,0)");
432  test_for_zero(t4ddg_1(1, 0, 1, 1)
433  - (t4ddg_2(1, 0, 0, 0) * t4ddg_3(0, 0, 1, 1)
434  + t4ddg_2(1, 0, 0, 1) * t4ddg_3(0, 1, 1, 1)
435  + t4ddg_2(1, 0, 0, 2) * t4ddg_3(0, 2, 1, 1)
436  + t4ddg_2(1, 1, 0, 0) * t4ddg_3(1, 0, 1, 1)
437  + t4ddg_2(1, 1, 0, 1) * t4ddg_3(1, 1, 1, 1)
438  + t4ddg_2(1, 1, 0, 2) * t4ddg_3(1, 2, 1, 1)
439  + t4ddg_2(1, 2, 0, 0) * t4ddg_3(2, 0, 1, 1)
440  + t4ddg_2(1, 2, 0, 1) * t4ddg_3(2, 1, 1, 1)
441  + t4ddg_2(1, 2, 0, 2) * t4ddg_3(2, 2, 1, 1)),
442  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,1,0,1)");
443  test_for_zero(t4ddg_1(1, 0, 1, 2)
444  - (t4ddg_2(1, 0, 0, 0) * t4ddg_3(0, 0, 1, 2)
445  + t4ddg_2(1, 0, 0, 1) * t4ddg_3(0, 1, 1, 2)
446  + t4ddg_2(1, 0, 0, 2) * t4ddg_3(0, 2, 1, 2)
447  + t4ddg_2(1, 1, 0, 0) * t4ddg_3(1, 0, 1, 2)
448  + t4ddg_2(1, 1, 0, 1) * t4ddg_3(1, 1, 1, 2)
449  + t4ddg_2(1, 1, 0, 2) * t4ddg_3(1, 2, 1, 2)
450  + t4ddg_2(1, 2, 0, 0) * t4ddg_3(2, 0, 1, 2)
451  + t4ddg_2(1, 2, 0, 1) * t4ddg_3(2, 1, 1, 2)
452  + t4ddg_2(1, 2, 0, 2) * t4ddg_3(2, 2, 1, 2)),
453  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,1,0,2)");
454  test_for_zero(t4ddg_1(1, 1, 1, 0)
455  - (t4ddg_2(1, 0, 1, 0) * t4ddg_3(0, 0, 1, 0)
456  + t4ddg_2(1, 0, 1, 1) * t4ddg_3(0, 1, 1, 0)
457  + t4ddg_2(1, 0, 1, 2) * t4ddg_3(0, 2, 1, 0)
458  + t4ddg_2(1, 1, 1, 0) * t4ddg_3(1, 0, 1, 0)
459  + t4ddg_2(1, 1, 1, 1) * t4ddg_3(1, 1, 1, 0)
460  + t4ddg_2(1, 1, 1, 2) * t4ddg_3(1, 2, 1, 0)
461  + t4ddg_2(1, 2, 1, 0) * t4ddg_3(2, 0, 1, 0)
462  + t4ddg_2(1, 2, 1, 1) * t4ddg_3(2, 1, 1, 0)
463  + t4ddg_2(1, 2, 1, 2) * t4ddg_3(2, 2, 1, 0)),
464  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,1,1,0)");
465  test_for_zero(t4ddg_1(1, 1, 1, 1)
466  - (t4ddg_2(1, 0, 1, 0) * t4ddg_3(0, 0, 1, 1)
467  + t4ddg_2(1, 0, 1, 1) * t4ddg_3(0, 1, 1, 1)
468  + t4ddg_2(1, 0, 1, 2) * t4ddg_3(0, 2, 1, 1)
469  + t4ddg_2(1, 1, 1, 0) * t4ddg_3(1, 0, 1, 1)
470  + t4ddg_2(1, 1, 1, 1) * t4ddg_3(1, 1, 1, 1)
471  + t4ddg_2(1, 1, 1, 2) * t4ddg_3(1, 2, 1, 1)
472  + t4ddg_2(1, 2, 1, 0) * t4ddg_3(2, 0, 1, 1)
473  + t4ddg_2(1, 2, 1, 1) * t4ddg_3(2, 1, 1, 1)
474  + t4ddg_2(1, 2, 1, 2) * t4ddg_3(2, 2, 1, 1)),
475  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,1,1,1)");
476  test_for_zero(t4ddg_1(1, 1, 1, 2)
477  - (t4ddg_2(1, 0, 1, 0) * t4ddg_3(0, 0, 1, 2)
478  + t4ddg_2(1, 0, 1, 1) * t4ddg_3(0, 1, 1, 2)
479  + t4ddg_2(1, 0, 1, 2) * t4ddg_3(0, 2, 1, 2)
480  + t4ddg_2(1, 1, 1, 0) * t4ddg_3(1, 0, 1, 2)
481  + t4ddg_2(1, 1, 1, 1) * t4ddg_3(1, 1, 1, 2)
482  + t4ddg_2(1, 1, 1, 2) * t4ddg_3(1, 2, 1, 2)
483  + t4ddg_2(1, 2, 1, 0) * t4ddg_3(2, 0, 1, 2)
484  + t4ddg_2(1, 2, 1, 1) * t4ddg_3(2, 1, 1, 2)
485  + t4ddg_2(1, 2, 1, 2) * t4ddg_3(2, 2, 1, 2)),
486  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,1,1,2)");
487  test_for_zero(t4ddg_1(1, 2, 1, 0)
488  - (t4ddg_2(1, 0, 2, 0) * t4ddg_3(0, 0, 1, 0)
489  + t4ddg_2(1, 0, 2, 1) * t4ddg_3(0, 1, 1, 0)
490  + t4ddg_2(1, 0, 2, 2) * t4ddg_3(0, 2, 1, 0)
491  + t4ddg_2(1, 1, 2, 0) * t4ddg_3(1, 0, 1, 0)
492  + t4ddg_2(1, 1, 2, 1) * t4ddg_3(1, 1, 1, 0)
493  + t4ddg_2(1, 1, 2, 2) * t4ddg_3(1, 2, 1, 0)
494  + t4ddg_2(1, 2, 2, 0) * t4ddg_3(2, 0, 1, 0)
495  + t4ddg_2(1, 2, 2, 1) * t4ddg_3(2, 1, 1, 0)
496  + t4ddg_2(1, 2, 2, 2) * t4ddg_3(2, 2, 1, 0)),
497  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,1,2,0)");
498  test_for_zero(t4ddg_1(1, 2, 1, 1)
499  - (t4ddg_2(1, 0, 2, 0) * t4ddg_3(0, 0, 1, 1)
500  + t4ddg_2(1, 0, 2, 1) * t4ddg_3(0, 1, 1, 1)
501  + t4ddg_2(1, 0, 2, 2) * t4ddg_3(0, 2, 1, 1)
502  + t4ddg_2(1, 1, 2, 0) * t4ddg_3(1, 0, 1, 1)
503  + t4ddg_2(1, 1, 2, 1) * t4ddg_3(1, 1, 1, 1)
504  + t4ddg_2(1, 1, 2, 2) * t4ddg_3(1, 2, 1, 1)
505  + t4ddg_2(1, 2, 2, 0) * t4ddg_3(2, 0, 1, 1)
506  + t4ddg_2(1, 2, 2, 1) * t4ddg_3(2, 1, 1, 1)
507  + t4ddg_2(1, 2, 2, 2) * t4ddg_3(2, 2, 1, 1)),
508  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,1,2,1)");
509  test_for_zero(t4ddg_1(1, 2, 1, 2)
510  - (t4ddg_2(1, 0, 2, 0) * t4ddg_3(0, 0, 1, 2)
511  + t4ddg_2(1, 0, 2, 1) * t4ddg_3(0, 1, 1, 2)
512  + t4ddg_2(1, 0, 2, 2) * t4ddg_3(0, 2, 1, 2)
513  + t4ddg_2(1, 1, 2, 0) * t4ddg_3(1, 0, 1, 2)
514  + t4ddg_2(1, 1, 2, 1) * t4ddg_3(1, 1, 1, 2)
515  + t4ddg_2(1, 1, 2, 2) * t4ddg_3(1, 2, 1, 2)
516  + t4ddg_2(1, 2, 2, 0) * t4ddg_3(2, 0, 1, 2)
517  + t4ddg_2(1, 2, 2, 1) * t4ddg_3(2, 1, 1, 2)
518  + t4ddg_2(1, 2, 2, 2) * t4ddg_3(2, 2, 1, 2)),
519  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,1,2,2)");
520  test_for_zero(t4ddg_1(1, 0, 2, 0)
521  - (t4ddg_2(1, 0, 0, 0) * t4ddg_3(0, 0, 2, 0)
522  + t4ddg_2(1, 0, 0, 1) * t4ddg_3(0, 1, 2, 0)
523  + t4ddg_2(1, 0, 0, 2) * t4ddg_3(0, 2, 2, 0)
524  + t4ddg_2(1, 1, 0, 0) * t4ddg_3(1, 0, 2, 0)
525  + t4ddg_2(1, 1, 0, 1) * t4ddg_3(1, 1, 2, 0)
526  + t4ddg_2(1, 1, 0, 2) * t4ddg_3(1, 2, 2, 0)
527  + t4ddg_2(1, 2, 0, 0) * t4ddg_3(2, 0, 2, 0)
528  + t4ddg_2(1, 2, 0, 1) * t4ddg_3(2, 1, 2, 0)
529  + t4ddg_2(1, 2, 0, 2) * t4ddg_3(2, 2, 2, 0)),
530  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,2,0,0)");
531  test_for_zero(t4ddg_1(1, 0, 2, 1)
532  - (t4ddg_2(1, 0, 0, 0) * t4ddg_3(0, 0, 2, 1)
533  + t4ddg_2(1, 0, 0, 1) * t4ddg_3(0, 1, 2, 1)
534  + t4ddg_2(1, 0, 0, 2) * t4ddg_3(0, 2, 2, 1)
535  + t4ddg_2(1, 1, 0, 0) * t4ddg_3(1, 0, 2, 1)
536  + t4ddg_2(1, 1, 0, 1) * t4ddg_3(1, 1, 2, 1)
537  + t4ddg_2(1, 1, 0, 2) * t4ddg_3(1, 2, 2, 1)
538  + t4ddg_2(1, 2, 0, 0) * t4ddg_3(2, 0, 2, 1)
539  + t4ddg_2(1, 2, 0, 1) * t4ddg_3(2, 1, 2, 1)
540  + t4ddg_2(1, 2, 0, 2) * t4ddg_3(2, 2, 2, 1)),
541  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,2,0,1)");
542  test_for_zero(t4ddg_1(1, 0, 2, 2)
543  - (t4ddg_2(1, 0, 0, 0) * t4ddg_3(0, 0, 2, 2)
544  + t4ddg_2(1, 0, 0, 1) * t4ddg_3(0, 1, 2, 2)
545  + t4ddg_2(1, 0, 0, 2) * t4ddg_3(0, 2, 2, 2)
546  + t4ddg_2(1, 1, 0, 0) * t4ddg_3(1, 0, 2, 2)
547  + t4ddg_2(1, 1, 0, 1) * t4ddg_3(1, 1, 2, 2)
548  + t4ddg_2(1, 1, 0, 2) * t4ddg_3(1, 2, 2, 2)
549  + t4ddg_2(1, 2, 0, 0) * t4ddg_3(2, 0, 2, 2)
550  + t4ddg_2(1, 2, 0, 1) * t4ddg_3(2, 1, 2, 2)
551  + t4ddg_2(1, 2, 0, 2) * t4ddg_3(2, 2, 2, 2)),
552  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,2,0,2)");
553  test_for_zero(t4ddg_1(1, 1, 2, 0)
554  - (t4ddg_2(1, 0, 1, 0) * t4ddg_3(0, 0, 2, 0)
555  + t4ddg_2(1, 0, 1, 1) * t4ddg_3(0, 1, 2, 0)
556  + t4ddg_2(1, 0, 1, 2) * t4ddg_3(0, 2, 2, 0)
557  + t4ddg_2(1, 1, 1, 0) * t4ddg_3(1, 0, 2, 0)
558  + t4ddg_2(1, 1, 1, 1) * t4ddg_3(1, 1, 2, 0)
559  + t4ddg_2(1, 1, 1, 2) * t4ddg_3(1, 2, 2, 0)
560  + t4ddg_2(1, 2, 1, 0) * t4ddg_3(2, 0, 2, 0)
561  + t4ddg_2(1, 2, 1, 1) * t4ddg_3(2, 1, 2, 0)
562  + t4ddg_2(1, 2, 1, 2) * t4ddg_3(2, 2, 2, 0)),
563  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,2,1,0)");
564  test_for_zero(t4ddg_1(1, 1, 2, 1)
565  - (t4ddg_2(1, 0, 1, 0) * t4ddg_3(0, 0, 2, 1)
566  + t4ddg_2(1, 0, 1, 1) * t4ddg_3(0, 1, 2, 1)
567  + t4ddg_2(1, 0, 1, 2) * t4ddg_3(0, 2, 2, 1)
568  + t4ddg_2(1, 1, 1, 0) * t4ddg_3(1, 0, 2, 1)
569  + t4ddg_2(1, 1, 1, 1) * t4ddg_3(1, 1, 2, 1)
570  + t4ddg_2(1, 1, 1, 2) * t4ddg_3(1, 2, 2, 1)
571  + t4ddg_2(1, 2, 1, 0) * t4ddg_3(2, 0, 2, 1)
572  + t4ddg_2(1, 2, 1, 1) * t4ddg_3(2, 1, 2, 1)
573  + t4ddg_2(1, 2, 1, 2) * t4ddg_3(2, 2, 2, 1)),
574  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,2,1,1)");
575  test_for_zero(t4ddg_1(1, 1, 2, 2)
576  - (t4ddg_2(1, 0, 1, 0) * t4ddg_3(0, 0, 2, 2)
577  + t4ddg_2(1, 0, 1, 1) * t4ddg_3(0, 1, 2, 2)
578  + t4ddg_2(1, 0, 1, 2) * t4ddg_3(0, 2, 2, 2)
579  + t4ddg_2(1, 1, 1, 0) * t4ddg_3(1, 0, 2, 2)
580  + t4ddg_2(1, 1, 1, 1) * t4ddg_3(1, 1, 2, 2)
581  + t4ddg_2(1, 1, 1, 2) * t4ddg_3(1, 2, 2, 2)
582  + t4ddg_2(1, 2, 1, 0) * t4ddg_3(2, 0, 2, 2)
583  + t4ddg_2(1, 2, 1, 1) * t4ddg_3(2, 1, 2, 2)
584  + t4ddg_2(1, 2, 1, 2) * t4ddg_3(2, 2, 2, 2)),
585  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,2,1,2)");
586  test_for_zero(t4ddg_1(1, 2, 2, 0)
587  - (t4ddg_2(1, 0, 2, 0) * t4ddg_3(0, 0, 2, 0)
588  + t4ddg_2(1, 0, 2, 1) * t4ddg_3(0, 1, 2, 0)
589  + t4ddg_2(1, 0, 2, 2) * t4ddg_3(0, 2, 2, 0)
590  + t4ddg_2(1, 1, 2, 0) * t4ddg_3(1, 0, 2, 0)
591  + t4ddg_2(1, 1, 2, 1) * t4ddg_3(1, 1, 2, 0)
592  + t4ddg_2(1, 1, 2, 2) * t4ddg_3(1, 2, 2, 0)
593  + t4ddg_2(1, 2, 2, 0) * t4ddg_3(2, 0, 2, 0)
594  + t4ddg_2(1, 2, 2, 1) * t4ddg_3(2, 1, 2, 0)
595  + t4ddg_2(1, 2, 2, 2) * t4ddg_3(2, 2, 2, 0)),
596  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,2,2,0)");
597  test_for_zero(t4ddg_1(1, 2, 2, 1)
598  - (t4ddg_2(1, 0, 2, 0) * t4ddg_3(0, 0, 2, 1)
599  + t4ddg_2(1, 0, 2, 1) * t4ddg_3(0, 1, 2, 1)
600  + t4ddg_2(1, 0, 2, 2) * t4ddg_3(0, 2, 2, 1)
601  + t4ddg_2(1, 1, 2, 0) * t4ddg_3(1, 0, 2, 1)
602  + t4ddg_2(1, 1, 2, 1) * t4ddg_3(1, 1, 2, 1)
603  + t4ddg_2(1, 1, 2, 2) * t4ddg_3(1, 2, 2, 1)
604  + t4ddg_2(1, 2, 2, 0) * t4ddg_3(2, 0, 2, 1)
605  + t4ddg_2(1, 2, 2, 1) * t4ddg_3(2, 1, 2, 1)
606  + t4ddg_2(1, 2, 2, 2) * t4ddg_3(2, 2, 2, 1)),
607  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,2,2,1)");
608  test_for_zero(t4ddg_1(1, 2, 2, 2)
609  - (t4ddg_2(1, 0, 2, 0) * t4ddg_3(0, 0, 2, 2)
610  + t4ddg_2(1, 0, 2, 1) * t4ddg_3(0, 1, 2, 2)
611  + t4ddg_2(1, 0, 2, 2) * t4ddg_3(0, 2, 2, 2)
612  + t4ddg_2(1, 1, 2, 0) * t4ddg_3(1, 0, 2, 2)
613  + t4ddg_2(1, 1, 2, 1) * t4ddg_3(1, 1, 2, 2)
614  + t4ddg_2(1, 1, 2, 2) * t4ddg_3(1, 2, 2, 2)
615  + t4ddg_2(1, 2, 2, 0) * t4ddg_3(2, 0, 2, 2)
616  + t4ddg_2(1, 2, 2, 1) * t4ddg_3(2, 1, 2, 2)
617  + t4ddg_2(1, 2, 2, 2) * t4ddg_3(2, 2, 2, 2)),
618  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,2,2,2)");
619  test_for_zero(t4ddg_1(2, 0, 0, 0)
620  - (t4ddg_2(2, 0, 0, 0) * t4ddg_3(0, 0, 0, 0)
621  + t4ddg_2(2, 0, 0, 1) * t4ddg_3(0, 1, 0, 0)
622  + t4ddg_2(2, 0, 0, 2) * t4ddg_3(0, 2, 0, 0)
623  + t4ddg_2(2, 1, 0, 0) * t4ddg_3(1, 0, 0, 0)
624  + t4ddg_2(2, 1, 0, 1) * t4ddg_3(1, 1, 0, 0)
625  + t4ddg_2(2, 1, 0, 2) * t4ddg_3(1, 2, 0, 0)
626  + t4ddg_2(2, 2, 0, 0) * t4ddg_3(2, 0, 0, 0)
627  + t4ddg_2(2, 2, 0, 1) * t4ddg_3(2, 1, 0, 0)
628  + t4ddg_2(2, 2, 0, 2) * t4ddg_3(2, 2, 0, 0)),
629  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,0,0,0)");
630  test_for_zero(t4ddg_1(2, 0, 0, 1)
631  - (t4ddg_2(2, 0, 0, 0) * t4ddg_3(0, 0, 0, 1)
632  + t4ddg_2(2, 0, 0, 1) * t4ddg_3(0, 1, 0, 1)
633  + t4ddg_2(2, 0, 0, 2) * t4ddg_3(0, 2, 0, 1)
634  + t4ddg_2(2, 1, 0, 0) * t4ddg_3(1, 0, 0, 1)
635  + t4ddg_2(2, 1, 0, 1) * t4ddg_3(1, 1, 0, 1)
636  + t4ddg_2(2, 1, 0, 2) * t4ddg_3(1, 2, 0, 1)
637  + t4ddg_2(2, 2, 0, 0) * t4ddg_3(2, 0, 0, 1)
638  + t4ddg_2(2, 2, 0, 1) * t4ddg_3(2, 1, 0, 1)
639  + t4ddg_2(2, 2, 0, 2) * t4ddg_3(2, 2, 0, 1)),
640  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,0,0,1)");
641  test_for_zero(t4ddg_1(2, 0, 0, 2)
642  - (t4ddg_2(2, 0, 0, 0) * t4ddg_3(0, 0, 0, 2)
643  + t4ddg_2(2, 0, 0, 1) * t4ddg_3(0, 1, 0, 2)
644  + t4ddg_2(2, 0, 0, 2) * t4ddg_3(0, 2, 0, 2)
645  + t4ddg_2(2, 1, 0, 0) * t4ddg_3(1, 0, 0, 2)
646  + t4ddg_2(2, 1, 0, 1) * t4ddg_3(1, 1, 0, 2)
647  + t4ddg_2(2, 1, 0, 2) * t4ddg_3(1, 2, 0, 2)
648  + t4ddg_2(2, 2, 0, 0) * t4ddg_3(2, 0, 0, 2)
649  + t4ddg_2(2, 2, 0, 1) * t4ddg_3(2, 1, 0, 2)
650  + t4ddg_2(2, 2, 0, 2) * t4ddg_3(2, 2, 0, 2)),
651  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,0,0,2)");
652  test_for_zero(t4ddg_1(2, 1, 0, 0)
653  - (t4ddg_2(2, 0, 1, 0) * t4ddg_3(0, 0, 0, 0)
654  + t4ddg_2(2, 0, 1, 1) * t4ddg_3(0, 1, 0, 0)
655  + t4ddg_2(2, 0, 1, 2) * t4ddg_3(0, 2, 0, 0)
656  + t4ddg_2(2, 1, 1, 0) * t4ddg_3(1, 0, 0, 0)
657  + t4ddg_2(2, 1, 1, 1) * t4ddg_3(1, 1, 0, 0)
658  + t4ddg_2(2, 1, 1, 2) * t4ddg_3(1, 2, 0, 0)
659  + t4ddg_2(2, 2, 1, 0) * t4ddg_3(2, 0, 0, 0)
660  + t4ddg_2(2, 2, 1, 1) * t4ddg_3(2, 1, 0, 0)
661  + t4ddg_2(2, 2, 1, 2) * t4ddg_3(2, 2, 0, 0)),
662  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,0,1,0)");
663  test_for_zero(t4ddg_1(2, 1, 0, 1)
664  - (t4ddg_2(2, 0, 1, 0) * t4ddg_3(0, 0, 0, 1)
665  + t4ddg_2(2, 0, 1, 1) * t4ddg_3(0, 1, 0, 1)
666  + t4ddg_2(2, 0, 1, 2) * t4ddg_3(0, 2, 0, 1)
667  + t4ddg_2(2, 1, 1, 0) * t4ddg_3(1, 0, 0, 1)
668  + t4ddg_2(2, 1, 1, 1) * t4ddg_3(1, 1, 0, 1)
669  + t4ddg_2(2, 1, 1, 2) * t4ddg_3(1, 2, 0, 1)
670  + t4ddg_2(2, 2, 1, 0) * t4ddg_3(2, 0, 0, 1)
671  + t4ddg_2(2, 2, 1, 1) * t4ddg_3(2, 1, 0, 1)
672  + t4ddg_2(2, 2, 1, 2) * t4ddg_3(2, 2, 0, 1)),
673  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,0,1,1)");
674  test_for_zero(t4ddg_1(2, 1, 0, 2)
675  - (t4ddg_2(2, 0, 1, 0) * t4ddg_3(0, 0, 0, 2)
676  + t4ddg_2(2, 0, 1, 1) * t4ddg_3(0, 1, 0, 2)
677  + t4ddg_2(2, 0, 1, 2) * t4ddg_3(0, 2, 0, 2)
678  + t4ddg_2(2, 1, 1, 0) * t4ddg_3(1, 0, 0, 2)
679  + t4ddg_2(2, 1, 1, 1) * t4ddg_3(1, 1, 0, 2)
680  + t4ddg_2(2, 1, 1, 2) * t4ddg_3(1, 2, 0, 2)
681  + t4ddg_2(2, 2, 1, 0) * t4ddg_3(2, 0, 0, 2)
682  + t4ddg_2(2, 2, 1, 1) * t4ddg_3(2, 1, 0, 2)
683  + t4ddg_2(2, 2, 1, 2) * t4ddg_3(2, 2, 0, 2)),
684  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,0,1,2)");
685  test_for_zero(t4ddg_1(2, 2, 0, 0)
686  - (t4ddg_2(2, 0, 2, 0) * t4ddg_3(0, 0, 0, 0)
687  + t4ddg_2(2, 0, 2, 1) * t4ddg_3(0, 1, 0, 0)
688  + t4ddg_2(2, 0, 2, 2) * t4ddg_3(0, 2, 0, 0)
689  + t4ddg_2(2, 1, 2, 0) * t4ddg_3(1, 0, 0, 0)
690  + t4ddg_2(2, 1, 2, 1) * t4ddg_3(1, 1, 0, 0)
691  + t4ddg_2(2, 1, 2, 2) * t4ddg_3(1, 2, 0, 0)
692  + t4ddg_2(2, 2, 2, 0) * t4ddg_3(2, 0, 0, 0)
693  + t4ddg_2(2, 2, 2, 1) * t4ddg_3(2, 1, 0, 0)
694  + t4ddg_2(2, 2, 2, 2) * t4ddg_3(2, 2, 0, 0)),
695  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,0,2,0)");
696  test_for_zero(t4ddg_1(2, 2, 0, 1)
697  - (t4ddg_2(2, 0, 2, 0) * t4ddg_3(0, 0, 0, 1)
698  + t4ddg_2(2, 0, 2, 1) * t4ddg_3(0, 1, 0, 1)
699  + t4ddg_2(2, 0, 2, 2) * t4ddg_3(0, 2, 0, 1)
700  + t4ddg_2(2, 1, 2, 0) * t4ddg_3(1, 0, 0, 1)
701  + t4ddg_2(2, 1, 2, 1) * t4ddg_3(1, 1, 0, 1)
702  + t4ddg_2(2, 1, 2, 2) * t4ddg_3(1, 2, 0, 1)
703  + t4ddg_2(2, 2, 2, 0) * t4ddg_3(2, 0, 0, 1)
704  + t4ddg_2(2, 2, 2, 1) * t4ddg_3(2, 1, 0, 1)
705  + t4ddg_2(2, 2, 2, 2) * t4ddg_3(2, 2, 0, 1)),
706  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,0,2,1)");
707  test_for_zero(t4ddg_1(2, 2, 0, 2)
708  - (t4ddg_2(2, 0, 2, 0) * t4ddg_3(0, 0, 0, 2)
709  + t4ddg_2(2, 0, 2, 1) * t4ddg_3(0, 1, 0, 2)
710  + t4ddg_2(2, 0, 2, 2) * t4ddg_3(0, 2, 0, 2)
711  + t4ddg_2(2, 1, 2, 0) * t4ddg_3(1, 0, 0, 2)
712  + t4ddg_2(2, 1, 2, 1) * t4ddg_3(1, 1, 0, 2)
713  + t4ddg_2(2, 1, 2, 2) * t4ddg_3(1, 2, 0, 2)
714  + t4ddg_2(2, 2, 2, 0) * t4ddg_3(2, 0, 0, 2)
715  + t4ddg_2(2, 2, 2, 1) * t4ddg_3(2, 1, 0, 2)
716  + t4ddg_2(2, 2, 2, 2) * t4ddg_3(2, 2, 0, 2)),
717  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,0,2,2)");
718  test_for_zero(t4ddg_1(2, 0, 1, 0)
719  - (t4ddg_2(2, 0, 0, 0) * t4ddg_3(0, 0, 1, 0)
720  + t4ddg_2(2, 0, 0, 1) * t4ddg_3(0, 1, 1, 0)
721  + t4ddg_2(2, 0, 0, 2) * t4ddg_3(0, 2, 1, 0)
722  + t4ddg_2(2, 1, 0, 0) * t4ddg_3(1, 0, 1, 0)
723  + t4ddg_2(2, 1, 0, 1) * t4ddg_3(1, 1, 1, 0)
724  + t4ddg_2(2, 1, 0, 2) * t4ddg_3(1, 2, 1, 0)
725  + t4ddg_2(2, 2, 0, 0) * t4ddg_3(2, 0, 1, 0)
726  + t4ddg_2(2, 2, 0, 1) * t4ddg_3(2, 1, 1, 0)
727  + t4ddg_2(2, 2, 0, 2) * t4ddg_3(2, 2, 1, 0)),
728  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,1,0,0)");
729  test_for_zero(t4ddg_1(2, 0, 1, 1)
730  - (t4ddg_2(2, 0, 0, 0) * t4ddg_3(0, 0, 1, 1)
731  + t4ddg_2(2, 0, 0, 1) * t4ddg_3(0, 1, 1, 1)
732  + t4ddg_2(2, 0, 0, 2) * t4ddg_3(0, 2, 1, 1)
733  + t4ddg_2(2, 1, 0, 0) * t4ddg_3(1, 0, 1, 1)
734  + t4ddg_2(2, 1, 0, 1) * t4ddg_3(1, 1, 1, 1)
735  + t4ddg_2(2, 1, 0, 2) * t4ddg_3(1, 2, 1, 1)
736  + t4ddg_2(2, 2, 0, 0) * t4ddg_3(2, 0, 1, 1)
737  + t4ddg_2(2, 2, 0, 1) * t4ddg_3(2, 1, 1, 1)
738  + t4ddg_2(2, 2, 0, 2) * t4ddg_3(2, 2, 1, 1)),
739  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,1,0,1)");
740  test_for_zero(t4ddg_1(2, 0, 1, 2)
741  - (t4ddg_2(2, 0, 0, 0) * t4ddg_3(0, 0, 1, 2)
742  + t4ddg_2(2, 0, 0, 1) * t4ddg_3(0, 1, 1, 2)
743  + t4ddg_2(2, 0, 0, 2) * t4ddg_3(0, 2, 1, 2)
744  + t4ddg_2(2, 1, 0, 0) * t4ddg_3(1, 0, 1, 2)
745  + t4ddg_2(2, 1, 0, 1) * t4ddg_3(1, 1, 1, 2)
746  + t4ddg_2(2, 1, 0, 2) * t4ddg_3(1, 2, 1, 2)
747  + t4ddg_2(2, 2, 0, 0) * t4ddg_3(2, 0, 1, 2)
748  + t4ddg_2(2, 2, 0, 1) * t4ddg_3(2, 1, 1, 2)
749  + t4ddg_2(2, 2, 0, 2) * t4ddg_3(2, 2, 1, 2)),
750  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,1,0,2)");
751  test_for_zero(t4ddg_1(2, 1, 1, 0)
752  - (t4ddg_2(2, 0, 1, 0) * t4ddg_3(0, 0, 1, 0)
753  + t4ddg_2(2, 0, 1, 1) * t4ddg_3(0, 1, 1, 0)
754  + t4ddg_2(2, 0, 1, 2) * t4ddg_3(0, 2, 1, 0)
755  + t4ddg_2(2, 1, 1, 0) * t4ddg_3(1, 0, 1, 0)
756  + t4ddg_2(2, 1, 1, 1) * t4ddg_3(1, 1, 1, 0)
757  + t4ddg_2(2, 1, 1, 2) * t4ddg_3(1, 2, 1, 0)
758  + t4ddg_2(2, 2, 1, 0) * t4ddg_3(2, 0, 1, 0)
759  + t4ddg_2(2, 2, 1, 1) * t4ddg_3(2, 1, 1, 0)
760  + t4ddg_2(2, 2, 1, 2) * t4ddg_3(2, 2, 1, 0)),
761  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,1,1,0)");
762  test_for_zero(t4ddg_1(2, 1, 1, 1)
763  - (t4ddg_2(2, 0, 1, 0) * t4ddg_3(0, 0, 1, 1)
764  + t4ddg_2(2, 0, 1, 1) * t4ddg_3(0, 1, 1, 1)
765  + t4ddg_2(2, 0, 1, 2) * t4ddg_3(0, 2, 1, 1)
766  + t4ddg_2(2, 1, 1, 0) * t4ddg_3(1, 0, 1, 1)
767  + t4ddg_2(2, 1, 1, 1) * t4ddg_3(1, 1, 1, 1)
768  + t4ddg_2(2, 1, 1, 2) * t4ddg_3(1, 2, 1, 1)
769  + t4ddg_2(2, 2, 1, 0) * t4ddg_3(2, 0, 1, 1)
770  + t4ddg_2(2, 2, 1, 1) * t4ddg_3(2, 1, 1, 1)
771  + t4ddg_2(2, 2, 1, 2) * t4ddg_3(2, 2, 1, 1)),
772  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,1,1,1)");
773  test_for_zero(t4ddg_1(2, 1, 1, 2)
774  - (t4ddg_2(2, 0, 1, 0) * t4ddg_3(0, 0, 1, 2)
775  + t4ddg_2(2, 0, 1, 1) * t4ddg_3(0, 1, 1, 2)
776  + t4ddg_2(2, 0, 1, 2) * t4ddg_3(0, 2, 1, 2)
777  + t4ddg_2(2, 1, 1, 0) * t4ddg_3(1, 0, 1, 2)
778  + t4ddg_2(2, 1, 1, 1) * t4ddg_3(1, 1, 1, 2)
779  + t4ddg_2(2, 1, 1, 2) * t4ddg_3(1, 2, 1, 2)
780  + t4ddg_2(2, 2, 1, 0) * t4ddg_3(2, 0, 1, 2)
781  + t4ddg_2(2, 2, 1, 1) * t4ddg_3(2, 1, 1, 2)
782  + t4ddg_2(2, 2, 1, 2) * t4ddg_3(2, 2, 1, 2)),
783  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,1,1,2)");
784  test_for_zero(t4ddg_1(2, 2, 1, 0)
785  - (t4ddg_2(2, 0, 2, 0) * t4ddg_3(0, 0, 1, 0)
786  + t4ddg_2(2, 0, 2, 1) * t4ddg_3(0, 1, 1, 0)
787  + t4ddg_2(2, 0, 2, 2) * t4ddg_3(0, 2, 1, 0)
788  + t4ddg_2(2, 1, 2, 0) * t4ddg_3(1, 0, 1, 0)
789  + t4ddg_2(2, 1, 2, 1) * t4ddg_3(1, 1, 1, 0)
790  + t4ddg_2(2, 1, 2, 2) * t4ddg_3(1, 2, 1, 0)
791  + t4ddg_2(2, 2, 2, 0) * t4ddg_3(2, 0, 1, 0)
792  + t4ddg_2(2, 2, 2, 1) * t4ddg_3(2, 1, 1, 0)
793  + t4ddg_2(2, 2, 2, 2) * t4ddg_3(2, 2, 1, 0)),
794  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,1,2,0)");
795  test_for_zero(t4ddg_1(2, 2, 1, 1)
796  - (t4ddg_2(2, 0, 2, 0) * t4ddg_3(0, 0, 1, 1)
797  + t4ddg_2(2, 0, 2, 1) * t4ddg_3(0, 1, 1, 1)
798  + t4ddg_2(2, 0, 2, 2) * t4ddg_3(0, 2, 1, 1)
799  + t4ddg_2(2, 1, 2, 0) * t4ddg_3(1, 0, 1, 1)
800  + t4ddg_2(2, 1, 2, 1) * t4ddg_3(1, 1, 1, 1)
801  + t4ddg_2(2, 1, 2, 2) * t4ddg_3(1, 2, 1, 1)
802  + t4ddg_2(2, 2, 2, 0) * t4ddg_3(2, 0, 1, 1)
803  + t4ddg_2(2, 2, 2, 1) * t4ddg_3(2, 1, 1, 1)
804  + t4ddg_2(2, 2, 2, 2) * t4ddg_3(2, 2, 1, 1)),
805  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,1,2,1)");
806  test_for_zero(t4ddg_1(2, 2, 1, 2)
807  - (t4ddg_2(2, 0, 2, 0) * t4ddg_3(0, 0, 1, 2)
808  + t4ddg_2(2, 0, 2, 1) * t4ddg_3(0, 1, 1, 2)
809  + t4ddg_2(2, 0, 2, 2) * t4ddg_3(0, 2, 1, 2)
810  + t4ddg_2(2, 1, 2, 0) * t4ddg_3(1, 0, 1, 2)
811  + t4ddg_2(2, 1, 2, 1) * t4ddg_3(1, 1, 1, 2)
812  + t4ddg_2(2, 1, 2, 2) * t4ddg_3(1, 2, 1, 2)
813  + t4ddg_2(2, 2, 2, 0) * t4ddg_3(2, 0, 1, 2)
814  + t4ddg_2(2, 2, 2, 1) * t4ddg_3(2, 1, 1, 2)
815  + t4ddg_2(2, 2, 2, 2) * t4ddg_3(2, 2, 1, 2)),
816  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,1,2,2)");
817  test_for_zero(t4ddg_1(2, 0, 2, 0)
818  - (t4ddg_2(2, 0, 0, 0) * t4ddg_3(0, 0, 2, 0)
819  + t4ddg_2(2, 0, 0, 1) * t4ddg_3(0, 1, 2, 0)
820  + t4ddg_2(2, 0, 0, 2) * t4ddg_3(0, 2, 2, 0)
821  + t4ddg_2(2, 1, 0, 0) * t4ddg_3(1, 0, 2, 0)
822  + t4ddg_2(2, 1, 0, 1) * t4ddg_3(1, 1, 2, 0)
823  + t4ddg_2(2, 1, 0, 2) * t4ddg_3(1, 2, 2, 0)
824  + t4ddg_2(2, 2, 0, 0) * t4ddg_3(2, 0, 2, 0)
825  + t4ddg_2(2, 2, 0, 1) * t4ddg_3(2, 1, 2, 0)
826  + t4ddg_2(2, 2, 0, 2) * t4ddg_3(2, 2, 2, 0)),
827  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,2,0,0)");
828  test_for_zero(t4ddg_1(2, 0, 2, 1)
829  - (t4ddg_2(2, 0, 0, 0) * t4ddg_3(0, 0, 2, 1)
830  + t4ddg_2(2, 0, 0, 1) * t4ddg_3(0, 1, 2, 1)
831  + t4ddg_2(2, 0, 0, 2) * t4ddg_3(0, 2, 2, 1)
832  + t4ddg_2(2, 1, 0, 0) * t4ddg_3(1, 0, 2, 1)
833  + t4ddg_2(2, 1, 0, 1) * t4ddg_3(1, 1, 2, 1)
834  + t4ddg_2(2, 1, 0, 2) * t4ddg_3(1, 2, 2, 1)
835  + t4ddg_2(2, 2, 0, 0) * t4ddg_3(2, 0, 2, 1)
836  + t4ddg_2(2, 2, 0, 1) * t4ddg_3(2, 1, 2, 1)
837  + t4ddg_2(2, 2, 0, 2) * t4ddg_3(2, 2, 2, 1)),
838  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,2,0,1)");
839  test_for_zero(t4ddg_1(2, 0, 2, 2)
840  - (t4ddg_2(2, 0, 0, 0) * t4ddg_3(0, 0, 2, 2)
841  + t4ddg_2(2, 0, 0, 1) * t4ddg_3(0, 1, 2, 2)
842  + t4ddg_2(2, 0, 0, 2) * t4ddg_3(0, 2, 2, 2)
843  + t4ddg_2(2, 1, 0, 0) * t4ddg_3(1, 0, 2, 2)
844  + t4ddg_2(2, 1, 0, 1) * t4ddg_3(1, 1, 2, 2)
845  + t4ddg_2(2, 1, 0, 2) * t4ddg_3(1, 2, 2, 2)
846  + t4ddg_2(2, 2, 0, 0) * t4ddg_3(2, 0, 2, 2)
847  + t4ddg_2(2, 2, 0, 1) * t4ddg_3(2, 1, 2, 2)
848  + t4ddg_2(2, 2, 0, 2) * t4ddg_3(2, 2, 2, 2)),
849  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,2,0,2)");
850  test_for_zero(t4ddg_1(2, 1, 2, 0)
851  - (t4ddg_2(2, 0, 1, 0) * t4ddg_3(0, 0, 2, 0)
852  + t4ddg_2(2, 0, 1, 1) * t4ddg_3(0, 1, 2, 0)
853  + t4ddg_2(2, 0, 1, 2) * t4ddg_3(0, 2, 2, 0)
854  + t4ddg_2(2, 1, 1, 0) * t4ddg_3(1, 0, 2, 0)
855  + t4ddg_2(2, 1, 1, 1) * t4ddg_3(1, 1, 2, 0)
856  + t4ddg_2(2, 1, 1, 2) * t4ddg_3(1, 2, 2, 0)
857  + t4ddg_2(2, 2, 1, 0) * t4ddg_3(2, 0, 2, 0)
858  + t4ddg_2(2, 2, 1, 1) * t4ddg_3(2, 1, 2, 0)
859  + t4ddg_2(2, 2, 1, 2) * t4ddg_3(2, 2, 2, 0)),
860  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,2,1,0)");
861  test_for_zero(t4ddg_1(2, 1, 2, 1)
862  - (t4ddg_2(2, 0, 1, 0) * t4ddg_3(0, 0, 2, 1)
863  + t4ddg_2(2, 0, 1, 1) * t4ddg_3(0, 1, 2, 1)
864  + t4ddg_2(2, 0, 1, 2) * t4ddg_3(0, 2, 2, 1)
865  + t4ddg_2(2, 1, 1, 0) * t4ddg_3(1, 0, 2, 1)
866  + t4ddg_2(2, 1, 1, 1) * t4ddg_3(1, 1, 2, 1)
867  + t4ddg_2(2, 1, 1, 2) * t4ddg_3(1, 2, 2, 1)
868  + t4ddg_2(2, 2, 1, 0) * t4ddg_3(2, 0, 2, 1)
869  + t4ddg_2(2, 2, 1, 1) * t4ddg_3(2, 1, 2, 1)
870  + t4ddg_2(2, 2, 1, 2) * t4ddg_3(2, 2, 2, 1)),
871  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,2,1,1)");
872  test_for_zero(t4ddg_1(2, 1, 2, 2)
873  - (t4ddg_2(2, 0, 1, 0) * t4ddg_3(0, 0, 2, 2)
874  + t4ddg_2(2, 0, 1, 1) * t4ddg_3(0, 1, 2, 2)
875  + t4ddg_2(2, 0, 1, 2) * t4ddg_3(0, 2, 2, 2)
876  + t4ddg_2(2, 1, 1, 0) * t4ddg_3(1, 0, 2, 2)
877  + t4ddg_2(2, 1, 1, 1) * t4ddg_3(1, 1, 2, 2)
878  + t4ddg_2(2, 1, 1, 2) * t4ddg_3(1, 2, 2, 2)
879  + t4ddg_2(2, 2, 1, 0) * t4ddg_3(2, 0, 2, 2)
880  + t4ddg_2(2, 2, 1, 1) * t4ddg_3(2, 1, 2, 2)
881  + t4ddg_2(2, 2, 1, 2) * t4ddg_3(2, 2, 2, 2)),
882  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,2,1,2)");
883  test_for_zero(t4ddg_1(2, 2, 2, 0)
884  - (t4ddg_2(2, 0, 2, 0) * t4ddg_3(0, 0, 2, 0)
885  + t4ddg_2(2, 0, 2, 1) * t4ddg_3(0, 1, 2, 0)
886  + t4ddg_2(2, 0, 2, 2) * t4ddg_3(0, 2, 2, 0)
887  + t4ddg_2(2, 1, 2, 0) * t4ddg_3(1, 0, 2, 0)
888  + t4ddg_2(2, 1, 2, 1) * t4ddg_3(1, 1, 2, 0)
889  + t4ddg_2(2, 1, 2, 2) * t4ddg_3(1, 2, 2, 0)
890  + t4ddg_2(2, 2, 2, 0) * t4ddg_3(2, 0, 2, 0)
891  + t4ddg_2(2, 2, 2, 1) * t4ddg_3(2, 1, 2, 0)
892  + t4ddg_2(2, 2, 2, 2) * t4ddg_3(2, 2, 2, 0)),
893  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,2,2,0)");
894  test_for_zero(t4ddg_1(2, 2, 2, 1)
895  - (t4ddg_2(2, 0, 2, 0) * t4ddg_3(0, 0, 2, 1)
896  + t4ddg_2(2, 0, 2, 1) * t4ddg_3(0, 1, 2, 1)
897  + t4ddg_2(2, 0, 2, 2) * t4ddg_3(0, 2, 2, 1)
898  + t4ddg_2(2, 1, 2, 0) * t4ddg_3(1, 0, 2, 1)
899  + t4ddg_2(2, 1, 2, 1) * t4ddg_3(1, 1, 2, 1)
900  + t4ddg_2(2, 1, 2, 2) * t4ddg_3(1, 2, 2, 1)
901  + t4ddg_2(2, 2, 2, 0) * t4ddg_3(2, 0, 2, 1)
902  + t4ddg_2(2, 2, 2, 1) * t4ddg_3(2, 1, 2, 1)
903  + t4ddg_2(2, 2, 2, 2) * t4ddg_3(2, 2, 2, 1)),
904  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,2,2,1)");
905  test_for_zero(t4ddg_1(2, 2, 2, 2)
906  - (t4ddg_2(2, 0, 2, 0) * t4ddg_3(0, 0, 2, 2)
907  + t4ddg_2(2, 0, 2, 1) * t4ddg_3(0, 1, 2, 2)
908  + t4ddg_2(2, 0, 2, 2) * t4ddg_3(0, 2, 2, 2)
909  + t4ddg_2(2, 1, 2, 0) * t4ddg_3(1, 0, 2, 2)
910  + t4ddg_2(2, 1, 2, 1) * t4ddg_3(1, 1, 2, 2)
911  + t4ddg_2(2, 1, 2, 2) * t4ddg_3(1, 2, 2, 2)
912  + t4ddg_2(2, 2, 2, 0) * t4ddg_3(2, 0, 2, 2)
913  + t4ddg_2(2, 2, 2, 1) * t4ddg_3(2, 1, 2, 2)
914  + t4ddg_2(2, 2, 2, 2) * t4ddg_3(2, 2, 2, 2)),
915  "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,2,2,2)");
916 }
FTensor
JSON compatible output.
Definition: Christof_constructor.hpp:6
test_T4ddg_02
void test_T4ddg_02(const Tensor2_symmetric< double, 3 > &t2s_2)
Definition: test_T4ddg_02.cpp:7
FTensor::Tensor2_symmetric
Definition: Tensor2_symmetric_value.hpp:13
FTensor::Number< 0 >
test_for_zero
void test_for_zero(const T &t, const std::string &s)
Definition: test_for_zero.hpp:7
i
FTensor::Index< 'i', SPACE_DIM > i
Definition: hcurl_divergence_operator_2d.cpp:27
FTensor::Index< 'i', 3 >
convert.n
n
Definition: convert.py:82
std
Definition: enable_if.hpp:5
j
FTensor::Index< 'j', 3 > j
Definition: matrix_function.cpp:19
FTensor::Ddg< double, 3, 3 >
m
FTensor::Index< 'm', 3 > m
Definition: shallow_wave.cpp:80
k
FTensor::Index< 'k', 3 > k
Definition: matrix_function.cpp:20
l
FTensor::Index< 'l', 3 > l
Definition: matrix_function.cpp:21
Generated by Doxygen 1.8.17 and hosted at University of Glasgow