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