v0.14.0
test_T3dg_37.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_T3dg_37(const Tensor2<double, 3, 3> &t2_2,
8  const Tensor2_symmetric<double, 3> &t2s_2,
9  const Dg<double, 3, 3> &t3dg_2,
10  const Dg<double, 3, 3> &t3dg_3)
11 {
12  Index<'i', 3> i;
13  Index<'j', 3> j;
14  Index<'k', 3> k;
15  Index<'l', 3> l;
16 
17  Number<0> N0;
18  Number<1> N1;
19  Number<2> N2;
20 
21  test_for_zero((t3dg_2(i, j, k) + t3dg_3(k, j, i))(0, 0, 0)
22  - (t3dg_2(0, 0, 0) + t3dg_3(0, 0, 0)),
23  "T3dg(i,j,k)+T3(k,j,i)(0,0,0)");
24  test_for_zero((t3dg_2(i, j, k) + t3dg_3(k, j, i))(0, 0, 1)
25  - (t3dg_2(0, 0, 1) + t3dg_3(1, 0, 0)),
26  "T3dg(i,j,k)+T3(k,j,i)(0,0,1)");
27  test_for_zero((t3dg_2(i, j, k) + t3dg_3(k, j, i))(0, 0, 2)
28  - (t3dg_2(0, 0, 2) + t3dg_3(2, 0, 0)),
29  "T3dg(i,j,k)+T3(k,j,i)(0,0,2)");
30  test_for_zero((t3dg_2(i, j, k) + t3dg_3(k, j, i))(0, 1, 0)
31  - (t3dg_2(0, 1, 0) + t3dg_3(0, 1, 0)),
32  "T3dg(i,j,k)+T3(k,j,i)(0,1,0)");
33  test_for_zero((t3dg_2(i, j, k) + t3dg_3(k, j, i))(0, 1, 1)
34  - (t3dg_2(0, 1, 1) + t3dg_3(1, 1, 0)),
35  "T3dg(i,j,k)+T3(k,j,i)(0,1,1)");
36  test_for_zero((t3dg_2(i, j, k) + t3dg_3(k, j, i))(0, 1, 2)
37  - (t3dg_2(0, 1, 2) + t3dg_3(2, 1, 0)),
38  "T3dg(i,j,k)+T3(k,j,i)(0,1,2)");
39  test_for_zero((t3dg_2(i, j, k) + t3dg_3(k, j, i))(0, 2, 0)
40  - (t3dg_2(0, 2, 0) + t3dg_3(0, 2, 0)),
41  "T3dg(i,j,k)+T3(k,j,i)(0,2,0)");
42  test_for_zero((t3dg_2(i, j, k) + t3dg_3(k, j, i))(0, 2, 1)
43  - (t3dg_2(0, 2, 1) + t3dg_3(1, 2, 0)),
44  "T3dg(i,j,k)+T3(k,j,i)(0,2,1)");
45  test_for_zero((t3dg_2(i, j, k) + t3dg_3(k, j, i))(0, 2, 2)
46  - (t3dg_2(0, 2, 2) + t3dg_3(2, 2, 0)),
47  "T3dg(i,j,k)+T3(k,j,i)(0,2,2)");
48  test_for_zero((t3dg_2(i, j, k) + t3dg_3(k, j, i))(1, 0, 0)
49  - (t3dg_2(1, 0, 0) + t3dg_3(0, 0, 1)),
50  "T3dg(i,j,k)+T3(k,j,i)(1,0,0)");
51  test_for_zero((t3dg_2(i, j, k) + t3dg_3(k, j, i))(1, 0, 1)
52  - (t3dg_2(1, 0, 1) + t3dg_3(1, 0, 1)),
53  "T3dg(i,j,k)+T3(k,j,i)(1,0,1)");
54  test_for_zero((t3dg_2(i, j, k) + t3dg_3(k, j, i))(1, 0, 2)
55  - (t3dg_2(1, 0, 2) + t3dg_3(2, 0, 1)),
56  "T3dg(i,j,k)+T3(k,j,i)(1,0,2)");
57  test_for_zero((t3dg_2(i, j, k) + t3dg_3(k, j, i))(1, 1, 0)
58  - (t3dg_2(1, 1, 0) + t3dg_3(0, 1, 1)),
59  "T3dg(i,j,k)+T3(k,j,i)(1,1,0)");
60  test_for_zero((t3dg_2(i, j, k) + t3dg_3(k, j, i))(1, 1, 1)
61  - (t3dg_2(1, 1, 1) + t3dg_3(1, 1, 1)),
62  "T3dg(i,j,k)+T3(k,j,i)(1,1,1)");
63  test_for_zero((t3dg_2(i, j, k) + t3dg_3(k, j, i))(1, 1, 2)
64  - (t3dg_2(1, 1, 2) + t3dg_3(2, 1, 1)),
65  "T3dg(i,j,k)+T3(k,j,i)(1,1,2)");
66  test_for_zero((t3dg_2(i, j, k) + t3dg_3(k, j, i))(1, 2, 0)
67  - (t3dg_2(1, 2, 0) + t3dg_3(0, 2, 1)),
68  "T3dg(i,j,k)+T3(k,j,i)(1,2,0)");
69  test_for_zero((t3dg_2(i, j, k) + t3dg_3(k, j, i))(1, 2, 1)
70  - (t3dg_2(1, 2, 1) + t3dg_3(1, 2, 1)),
71  "T3dg(i,j,k)+T3(k,j,i)(1,2,1)");
72  test_for_zero((t3dg_2(i, j, k) + t3dg_3(k, j, i))(1, 2, 2)
73  - (t3dg_2(1, 2, 2) + t3dg_3(2, 2, 1)),
74  "T3dg(i,j,k)+T3(k,j,i)(1,2,2)");
75  test_for_zero((t3dg_2(i, j, k) + t3dg_3(k, j, i))(2, 0, 0)
76  - (t3dg_2(2, 0, 0) + t3dg_3(0, 0, 2)),
77  "T3dg(i,j,k)+T3(k,j,i)(2,0,0)");
78  test_for_zero((t3dg_2(i, j, k) + t3dg_3(k, j, i))(2, 0, 1)
79  - (t3dg_2(2, 0, 1) + t3dg_3(1, 0, 2)),
80  "T3dg(i,j,k)+T3(k,j,i)(2,0,1)");
81  test_for_zero((t3dg_2(i, j, k) + t3dg_3(k, j, i))(2, 0, 2)
82  - (t3dg_2(2, 0, 2) + t3dg_3(2, 0, 2)),
83  "T3dg(i,j,k)+T3(k,j,i)(2,0,2)");
84  test_for_zero((t3dg_2(i, j, k) + t3dg_3(k, j, i))(2, 1, 0)
85  - (t3dg_2(2, 1, 0) + t3dg_3(0, 1, 2)),
86  "T3dg(i,j,k)+T3(k,j,i)(2,1,0)");
87  test_for_zero((t3dg_2(i, j, k) + t3dg_3(k, j, i))(2, 1, 1)
88  - (t3dg_2(2, 1, 1) + t3dg_3(1, 1, 2)),
89  "T3dg(i,j,k)+T3(k,j,i)(2,1,1)");
90  test_for_zero((t3dg_2(i, j, k) + t3dg_3(k, j, i))(2, 1, 2)
91  - (t3dg_2(2, 1, 2) + t3dg_3(2, 1, 2)),
92  "T3dg(i,j,k)+T3(k,j,i)(2,1,2)");
93  test_for_zero((t3dg_2(i, j, k) + t3dg_3(k, j, i))(2, 2, 0)
94  - (t3dg_2(2, 2, 0) + t3dg_3(0, 2, 2)),
95  "T3dg(i,j,k)+T3(k,j,i)(2,2,0)");
96  test_for_zero((t3dg_2(i, j, k) + t3dg_3(k, j, i))(2, 2, 1)
97  - (t3dg_2(2, 2, 1) + t3dg_3(1, 2, 2)),
98  "T3dg(i,j,k)+T3(k,j,i)(2,2,1)");
99  test_for_zero((t3dg_2(i, j, k) + t3dg_3(k, j, i))(2, 2, 2)
100  - (t3dg_2(2, 2, 2) + t3dg_3(2, 2, 2)),
101  "T3dg(i,j,k)+T3(k,j,i)(2,2,2)");
102 
103  test_for_zero((t3dg_2(i, j, k) - t3dg_3(k, j, i))(0, 0, 0)
104  - (t3dg_2(0, 0, 0) - t3dg_3(0, 0, 0)),
105  "T3dg(i,j,k)-T3(k,j,i)(0,0,0)");
106  test_for_zero((t3dg_2(i, j, k) - t3dg_3(k, j, i))(0, 0, 1)
107  - (t3dg_2(0, 0, 1) - t3dg_3(1, 0, 0)),
108  "T3dg(i,j,k)-T3(k,j,i)(0,0,1)");
109  test_for_zero((t3dg_2(i, j, k) - t3dg_3(k, j, i))(0, 0, 2)
110  - (t3dg_2(0, 0, 2) - t3dg_3(2, 0, 0)),
111  "T3dg(i,j,k)-T3(k,j,i)(0,0,2)");
112  test_for_zero((t3dg_2(i, j, k) - t3dg_3(k, j, i))(0, 1, 0)
113  - (t3dg_2(0, 1, 0) - t3dg_3(0, 1, 0)),
114  "T3dg(i,j,k)-T3(k,j,i)(0,1,0)");
115  test_for_zero((t3dg_2(i, j, k) - t3dg_3(k, j, i))(0, 1, 1)
116  - (t3dg_2(0, 1, 1) - t3dg_3(1, 1, 0)),
117  "T3dg(i,j,k)-T3(k,j,i)(0,1,1)");
118  test_for_zero((t3dg_2(i, j, k) - t3dg_3(k, j, i))(0, 1, 2)
119  - (t3dg_2(0, 1, 2) - t3dg_3(2, 1, 0)),
120  "T3dg(i,j,k)-T3(k,j,i)(0,1,2)");
121  test_for_zero((t3dg_2(i, j, k) - t3dg_3(k, j, i))(0, 2, 0)
122  - (t3dg_2(0, 2, 0) - t3dg_3(0, 2, 0)),
123  "T3dg(i,j,k)-T3(k,j,i)(0,2,0)");
124  test_for_zero((t3dg_2(i, j, k) - t3dg_3(k, j, i))(0, 2, 1)
125  - (t3dg_2(0, 2, 1) - t3dg_3(1, 2, 0)),
126  "T3dg(i,j,k)-T3(k,j,i)(0,2,1)");
127  test_for_zero((t3dg_2(i, j, k) - t3dg_3(k, j, i))(0, 2, 2)
128  - (t3dg_2(0, 2, 2) - t3dg_3(2, 2, 0)),
129  "T3dg(i,j,k)-T3(k,j,i)(0,2,2)");
130  test_for_zero((t3dg_2(i, j, k) - t3dg_3(k, j, i))(1, 0, 0)
131  - (t3dg_2(1, 0, 0) - t3dg_3(0, 0, 1)),
132  "T3dg(i,j,k)-T3(k,j,i)(1,0,0)");
133  test_for_zero((t3dg_2(i, j, k) - t3dg_3(k, j, i))(1, 0, 1)
134  - (t3dg_2(1, 0, 1) - t3dg_3(1, 0, 1)),
135  "T3dg(i,j,k)-T3(k,j,i)(1,0,1)");
136  test_for_zero((t3dg_2(i, j, k) - t3dg_3(k, j, i))(1, 0, 2)
137  - (t3dg_2(1, 0, 2) - t3dg_3(2, 0, 1)),
138  "T3dg(i,j,k)-T3(k,j,i)(1,0,2)");
139  test_for_zero((t3dg_2(i, j, k) - t3dg_3(k, j, i))(1, 1, 0)
140  - (t3dg_2(1, 1, 0) - t3dg_3(0, 1, 1)),
141  "T3dg(i,j,k)-T3(k,j,i)(1,1,0)");
142  test_for_zero((t3dg_2(i, j, k) - t3dg_3(k, j, i))(1, 1, 1)
143  - (t3dg_2(1, 1, 1) - t3dg_3(1, 1, 1)),
144  "T3dg(i,j,k)-T3(k,j,i)(1,1,1)");
145  test_for_zero((t3dg_2(i, j, k) - t3dg_3(k, j, i))(1, 1, 2)
146  - (t3dg_2(1, 1, 2) - t3dg_3(2, 1, 1)),
147  "T3dg(i,j,k)-T3(k,j,i)(1,1,2)");
148  test_for_zero((t3dg_2(i, j, k) - t3dg_3(k, j, i))(1, 2, 0)
149  - (t3dg_2(1, 2, 0) - t3dg_3(0, 2, 1)),
150  "T3dg(i,j,k)-T3(k,j,i)(1,2,0)");
151  test_for_zero((t3dg_2(i, j, k) - t3dg_3(k, j, i))(1, 2, 1)
152  - (t3dg_2(1, 2, 1) - t3dg_3(1, 2, 1)),
153  "T3dg(i,j,k)-T3(k,j,i)(1,2,1)");
154  test_for_zero((t3dg_2(i, j, k) - t3dg_3(k, j, i))(1, 2, 2)
155  - (t3dg_2(1, 2, 2) - t3dg_3(2, 2, 1)),
156  "T3dg(i,j,k)-T3(k,j,i)(1,2,2)");
157  test_for_zero((t3dg_2(i, j, k) - t3dg_3(k, j, i))(2, 0, 0)
158  - (t3dg_2(2, 0, 0) - t3dg_3(0, 0, 2)),
159  "T3dg(i,j,k)-T3(k,j,i)(2,0,0)");
160  test_for_zero((t3dg_2(i, j, k) - t3dg_3(k, j, i))(2, 0, 1)
161  - (t3dg_2(2, 0, 1) - t3dg_3(1, 0, 2)),
162  "T3dg(i,j,k)-T3(k,j,i)(2,0,1)");
163  test_for_zero((t3dg_2(i, j, k) - t3dg_3(k, j, i))(2, 0, 2)
164  - (t3dg_2(2, 0, 2) - t3dg_3(2, 0, 2)),
165  "T3dg(i,j,k)-T3(k,j,i)(2,0,2)");
166  test_for_zero((t3dg_2(i, j, k) - t3dg_3(k, j, i))(2, 1, 0)
167  - (t3dg_2(2, 1, 0) - t3dg_3(0, 1, 2)),
168  "T3dg(i,j,k)-T3(k,j,i)(2,1,0)");
169  test_for_zero((t3dg_2(i, j, k) - t3dg_3(k, j, i))(2, 1, 1)
170  - (t3dg_2(2, 1, 1) - t3dg_3(1, 1, 2)),
171  "T3dg(i,j,k)-T3(k,j,i)(2,1,1)");
172  test_for_zero((t3dg_2(i, j, k) - t3dg_3(k, j, i))(2, 1, 2)
173  - (t3dg_2(2, 1, 2) - t3dg_3(2, 1, 2)),
174  "T3dg(i,j,k)-T3(k,j,i)(2,1,2)");
175  test_for_zero((t3dg_2(i, j, k) - t3dg_3(k, j, i))(2, 2, 0)
176  - (t3dg_2(2, 2, 0) - t3dg_3(0, 2, 2)),
177  "T3dg(i,j,k)-T3(k,j,i)(2,2,0)");
178  test_for_zero((t3dg_2(i, j, k) - t3dg_3(k, j, i))(2, 2, 1)
179  - (t3dg_2(2, 2, 1) - t3dg_3(1, 2, 2)),
180  "T3dg(i,j,k)-T3(k,j,i)(2,2,1)");
181  test_for_zero((t3dg_2(i, j, k) - t3dg_3(k, j, i))(2, 2, 2)
182  - (t3dg_2(2, 2, 2) - t3dg_3(2, 2, 2)),
183  "T3dg(i,j,k)-T3(k,j,i)(2,2,2)");
184 
185  test_for_zero((t3dg_2(i, j, k) * t2_2(j, l))(0, 0, 0)
186  - (t3dg_2(0, 0, 0) * t2_2(0, 0)
187  + t3dg_2(0, 1, 0) * t2_2(1, 0)
188  + t3dg_2(0, 2, 0) * t2_2(2, 0)),
189  "T3dg(i,j,k)*T2(j,l)(0,0,0)");
190  test_for_zero((t3dg_2(i, j, k) * t2_2(j, l))(0, 0, 1)
191  - (t3dg_2(0, 0, 0) * t2_2(0, 1)
192  + t3dg_2(0, 1, 0) * t2_2(1, 1)
193  + t3dg_2(0, 2, 0) * t2_2(2, 1)),
194  "T3dg(i,j,k)*T2(j,l)(0,0,1)");
195  test_for_zero((t3dg_2(i, j, k) * t2_2(j, l))(0, 0, 2)
196  - (t3dg_2(0, 0, 0) * t2_2(0, 2)
197  + t3dg_2(0, 1, 0) * t2_2(1, 2)
198  + t3dg_2(0, 2, 0) * t2_2(2, 2)),
199  "T3dg(i,j,k)*T2(j,l)(0,0,2)");
200  test_for_zero((t3dg_2(i, j, k) * t2_2(j, l))(0, 1, 0)
201  - (t3dg_2(0, 0, 1) * t2_2(0, 0)
202  + t3dg_2(0, 1, 1) * t2_2(1, 0)
203  + t3dg_2(0, 2, 1) * t2_2(2, 0)),
204  "T3dg(i,j,k)*T2(j,l)(0,1,0)");
205  test_for_zero((t3dg_2(i, j, k) * t2_2(j, l))(0, 1, 1)
206  - (t3dg_2(0, 0, 1) * t2_2(0, 1)
207  + t3dg_2(0, 1, 1) * t2_2(1, 1)
208  + t3dg_2(0, 2, 1) * t2_2(2, 1)),
209  "T3dg(i,j,k)*T2(j,l)(0,1,1)");
210  test_for_zero((t3dg_2(i, j, k) * t2_2(j, l))(0, 1, 2)
211  - (t3dg_2(0, 0, 1) * t2_2(0, 2)
212  + t3dg_2(0, 1, 1) * t2_2(1, 2)
213  + t3dg_2(0, 2, 1) * t2_2(2, 2)),
214  "T3dg(i,j,k)*T2(j,l)(0,1,2)");
215  test_for_zero((t3dg_2(i, j, k) * t2_2(j, l))(0, 2, 0)
216  - (t3dg_2(0, 0, 2) * t2_2(0, 0)
217  + t3dg_2(0, 1, 2) * t2_2(1, 0)
218  + t3dg_2(0, 2, 2) * t2_2(2, 0)),
219  "T3dg(i,j,k)*T2(j,l)(0,2,0)");
220  test_for_zero((t3dg_2(i, j, k) * t2_2(j, l))(0, 2, 1)
221  - (t3dg_2(0, 0, 2) * t2_2(0, 1)
222  + t3dg_2(0, 1, 2) * t2_2(1, 1)
223  + t3dg_2(0, 2, 2) * t2_2(2, 1)),
224  "T3dg(i,j,k)*T2(j,l)(0,2,1)");
225  test_for_zero((t3dg_2(i, j, k) * t2_2(j, l))(0, 2, 2)
226  - (t3dg_2(0, 0, 2) * t2_2(0, 2)
227  + t3dg_2(0, 1, 2) * t2_2(1, 2)
228  + t3dg_2(0, 2, 2) * t2_2(2, 2)),
229  "T3dg(i,j,k)*T2(j,l)(0,2,2)");
230  test_for_zero((t3dg_2(i, j, k) * t2_2(j, l))(1, 0, 0)
231  - (t3dg_2(1, 0, 0) * t2_2(0, 0)
232  + t3dg_2(1, 1, 0) * t2_2(1, 0)
233  + t3dg_2(1, 2, 0) * t2_2(2, 0)),
234  "T3dg(i,j,k)*T2(j,l)(1,0,0)");
235  test_for_zero((t3dg_2(i, j, k) * t2_2(j, l))(1, 0, 1)
236  - (t3dg_2(1, 0, 0) * t2_2(0, 1)
237  + t3dg_2(1, 1, 0) * t2_2(1, 1)
238  + t3dg_2(1, 2, 0) * t2_2(2, 1)),
239  "T3dg(i,j,k)*T2(j,l)(1,0,1)");
240  test_for_zero((t3dg_2(i, j, k) * t2_2(j, l))(1, 0, 2)
241  - (t3dg_2(1, 0, 0) * t2_2(0, 2)
242  + t3dg_2(1, 1, 0) * t2_2(1, 2)
243  + t3dg_2(1, 2, 0) * t2_2(2, 2)),
244  "T3dg(i,j,k)*T2(j,l)(1,0,2)");
245  test_for_zero((t3dg_2(i, j, k) * t2_2(j, l))(1, 1, 0)
246  - (t3dg_2(1, 0, 1) * t2_2(0, 0)
247  + t3dg_2(1, 1, 1) * t2_2(1, 0)
248  + t3dg_2(1, 2, 1) * t2_2(2, 0)),
249  "T3dg(i,j,k)*T2(j,l)(1,1,0)");
250  test_for_zero((t3dg_2(i, j, k) * t2_2(j, l))(1, 1, 1)
251  - (t3dg_2(1, 0, 1) * t2_2(0, 1)
252  + t3dg_2(1, 1, 1) * t2_2(1, 1)
253  + t3dg_2(1, 2, 1) * t2_2(2, 1)),
254  "T3dg(i,j,k)*T2(j,l)(1,1,1)");
255  test_for_zero((t3dg_2(i, j, k) * t2_2(j, l))(1, 1, 2)
256  - (t3dg_2(1, 0, 1) * t2_2(0, 2)
257  + t3dg_2(1, 1, 1) * t2_2(1, 2)
258  + t3dg_2(1, 2, 1) * t2_2(2, 2)),
259  "T3dg(i,j,k)*T2(j,l)(1,1,2)");
260  test_for_zero((t3dg_2(i, j, k) * t2_2(j, l))(1, 2, 0)
261  - (t3dg_2(1, 0, 2) * t2_2(0, 0)
262  + t3dg_2(1, 1, 2) * t2_2(1, 0)
263  + t3dg_2(1, 2, 2) * t2_2(2, 0)),
264  "T3dg(i,j,k)*T2(j,l)(1,2,0)");
265  test_for_zero((t3dg_2(i, j, k) * t2_2(j, l))(1, 2, 1)
266  - (t3dg_2(1, 0, 2) * t2_2(0, 1)
267  + t3dg_2(1, 1, 2) * t2_2(1, 1)
268  + t3dg_2(1, 2, 2) * t2_2(2, 1)),
269  "T3dg(i,j,k)*T2(j,l)(1,2,1)");
270  test_for_zero((t3dg_2(i, j, k) * t2_2(j, l))(1, 2, 2)
271  - (t3dg_2(1, 0, 2) * t2_2(0, 2)
272  + t3dg_2(1, 1, 2) * t2_2(1, 2)
273  + t3dg_2(1, 2, 2) * t2_2(2, 2)),
274  "T3dg(i,j,k)*T2(j,l)(1,2,2)");
275  test_for_zero((t3dg_2(i, j, k) * t2_2(j, l))(2, 0, 0)
276  - (t3dg_2(2, 0, 0) * t2_2(0, 0)
277  + t3dg_2(2, 1, 0) * t2_2(1, 0)
278  + t3dg_2(2, 2, 0) * t2_2(2, 0)),
279  "T3dg(i,j,k)*T2(j,l)(2,0,0)");
280  test_for_zero((t3dg_2(i, j, k) * t2_2(j, l))(2, 0, 1)
281  - (t3dg_2(2, 0, 0) * t2_2(0, 1)
282  + t3dg_2(2, 1, 0) * t2_2(1, 1)
283  + t3dg_2(2, 2, 0) * t2_2(2, 1)),
284  "T3dg(i,j,k)*T2(j,l)(2,0,1)");
285  test_for_zero((t3dg_2(i, j, k) * t2_2(j, l))(2, 0, 2)
286  - (t3dg_2(2, 0, 0) * t2_2(0, 2)
287  + t3dg_2(2, 1, 0) * t2_2(1, 2)
288  + t3dg_2(2, 2, 0) * t2_2(2, 2)),
289  "T3dg(i,j,k)*T2(j,l)(2,0,2)");
290  test_for_zero((t3dg_2(i, j, k) * t2_2(j, l))(2, 1, 0)
291  - (t3dg_2(2, 0, 1) * t2_2(0, 0)
292  + t3dg_2(2, 1, 1) * t2_2(1, 0)
293  + t3dg_2(2, 2, 1) * t2_2(2, 0)),
294  "T3dg(i,j,k)*T2(j,l)(2,1,0)");
295  test_for_zero((t3dg_2(i, j, k) * t2_2(j, l))(2, 1, 1)
296  - (t3dg_2(2, 0, 1) * t2_2(0, 1)
297  + t3dg_2(2, 1, 1) * t2_2(1, 1)
298  + t3dg_2(2, 2, 1) * t2_2(2, 1)),
299  "T3dg(i,j,k)*T2(j,l)(2,1,1)");
300  test_for_zero((t3dg_2(i, j, k) * t2_2(j, l))(2, 1, 2)
301  - (t3dg_2(2, 0, 1) * t2_2(0, 2)
302  + t3dg_2(2, 1, 1) * t2_2(1, 2)
303  + t3dg_2(2, 2, 1) * t2_2(2, 2)),
304  "T3dg(i,j,k)*T2(j,l)(2,1,2)");
305  test_for_zero((t3dg_2(i, j, k) * t2_2(j, l))(2, 2, 0)
306  - (t3dg_2(2, 0, 2) * t2_2(0, 0)
307  + t3dg_2(2, 1, 2) * t2_2(1, 0)
308  + t3dg_2(2, 2, 2) * t2_2(2, 0)),
309  "T3dg(i,j,k)*T2(j,l)(2,2,0)");
310  test_for_zero((t3dg_2(i, j, k) * t2_2(j, l))(2, 2, 1)
311  - (t3dg_2(2, 0, 2) * t2_2(0, 1)
312  + t3dg_2(2, 1, 2) * t2_2(1, 1)
313  + t3dg_2(2, 2, 2) * t2_2(2, 1)),
314  "T3dg(i,j,k)*T2(j,l)(2,2,1)");
315  test_for_zero((t3dg_2(i, j, k) * t2_2(j, l))(2, 2, 2)
316  - (t3dg_2(2, 0, 2) * t2_2(0, 2)
317  + t3dg_2(2, 1, 2) * t2_2(1, 2)
318  + t3dg_2(2, 2, 2) * t2_2(2, 2)),
319  "T3dg(i,j,k)*T2(j,l)(2,2,2)");
320 
321  test_for_zero((t3dg_3(i, j, k) * t2_2(l, j))(0, 0, 0)
322  - (t3dg_3(0, 0, 0) * t2_2(0, 0)
323  + t3dg_3(0, 1, 0) * t2_2(0, 1)
324  + t3dg_3(0, 2, 0) * t2_2(0, 2)),
325  "T3dg(i,j,k)*T2(l,j)(0,0,0)");
326  test_for_zero((t3dg_3(i, j, k) * t2_2(l, j))(0, 0, 1)
327  - (t3dg_3(0, 0, 0) * t2_2(1, 0)
328  + t3dg_3(0, 1, 0) * t2_2(1, 1)
329  + t3dg_3(0, 2, 0) * t2_2(1, 2)),
330  "T3dg(i,j,k)*T2(l,j)(0,0,1)");
331  test_for_zero((t3dg_3(i, j, k) * t2_2(l, j))(0, 0, 2)
332  - (t3dg_3(0, 0, 0) * t2_2(2, 0)
333  + t3dg_3(0, 1, 0) * t2_2(2, 1)
334  + t3dg_3(0, 2, 0) * t2_2(2, 2)),
335  "T3dg(i,j,k)*T2(l,j)(0,0,2)");
336  test_for_zero((t3dg_3(i, j, k) * t2_2(l, j))(0, 1, 0)
337  - (t3dg_3(0, 0, 1) * t2_2(0, 0)
338  + t3dg_3(0, 1, 1) * t2_2(0, 1)
339  + t3dg_3(0, 2, 1) * t2_2(0, 2)),
340  "T3dg(i,j,k)*T2(l,j)(0,1,0)");
341  test_for_zero((t3dg_3(i, j, k) * t2_2(l, j))(0, 1, 1)
342  - (t3dg_3(0, 0, 1) * t2_2(1, 0)
343  + t3dg_3(0, 1, 1) * t2_2(1, 1)
344  + t3dg_3(0, 2, 1) * t2_2(1, 2)),
345  "T3dg(i,j,k)*T2(l,j)(0,1,1)");
346  test_for_zero((t3dg_3(i, j, k) * t2_2(l, j))(0, 1, 2)
347  - (t3dg_3(0, 0, 1) * t2_2(2, 0)
348  + t3dg_3(0, 1, 1) * t2_2(2, 1)
349  + t3dg_3(0, 2, 1) * t2_2(2, 2)),
350  "T3dg(i,j,k)*T2(l,j)(0,1,2)");
351  test_for_zero((t3dg_3(i, j, k) * t2_2(l, j))(0, 2, 0)
352  - (t3dg_3(0, 0, 2) * t2_2(0, 0)
353  + t3dg_3(0, 1, 2) * t2_2(0, 1)
354  + t3dg_3(0, 2, 2) * t2_2(0, 2)),
355  "T3dg(i,j,k)*T2(l,j)(0,2,0)");
356  test_for_zero((t3dg_3(i, j, k) * t2_2(l, j))(0, 2, 1)
357  - (t3dg_3(0, 0, 2) * t2_2(1, 0)
358  + t3dg_3(0, 1, 2) * t2_2(1, 1)
359  + t3dg_3(0, 2, 2) * t2_2(1, 2)),
360  "T3dg(i,j,k)*T2(l,j)(0,2,1)");
361  test_for_zero((t3dg_3(i, j, k) * t2_2(l, j))(0, 2, 2)
362  - (t3dg_3(0, 0, 2) * t2_2(2, 0)
363  + t3dg_3(0, 1, 2) * t2_2(2, 1)
364  + t3dg_3(0, 2, 2) * t2_2(2, 2)),
365  "T3dg(i,j,k)*T2(l,j)(0,2,2)");
366  test_for_zero((t3dg_3(i, j, k) * t2_2(l, j))(1, 0, 0)
367  - (t3dg_3(1, 0, 0) * t2_2(0, 0)
368  + t3dg_3(1, 1, 0) * t2_2(0, 1)
369  + t3dg_3(1, 2, 0) * t2_2(0, 2)),
370  "T3dg(i,j,k)*T2(l,j)(1,0,0)");
371  test_for_zero((t3dg_3(i, j, k) * t2_2(l, j))(1, 0, 1)
372  - (t3dg_3(1, 0, 0) * t2_2(1, 0)
373  + t3dg_3(1, 1, 0) * t2_2(1, 1)
374  + t3dg_3(1, 2, 0) * t2_2(1, 2)),
375  "T3dg(i,j,k)*T2(l,j)(1,0,1)");
376  test_for_zero((t3dg_3(i, j, k) * t2_2(l, j))(1, 0, 2)
377  - (t3dg_3(1, 0, 0) * t2_2(2, 0)
378  + t3dg_3(1, 1, 0) * t2_2(2, 1)
379  + t3dg_3(1, 2, 0) * t2_2(2, 2)),
380  "T3dg(i,j,k)*T2(l,j)(1,0,2)");
381  test_for_zero((t3dg_3(i, j, k) * t2_2(l, j))(1, 1, 0)
382  - (t3dg_3(1, 0, 1) * t2_2(0, 0)
383  + t3dg_3(1, 1, 1) * t2_2(0, 1)
384  + t3dg_3(1, 2, 1) * t2_2(0, 2)),
385  "T3dg(i,j,k)*T2(l,j)(1,1,0)");
386  test_for_zero((t3dg_3(i, j, k) * t2_2(l, j))(1, 1, 1)
387  - (t3dg_3(1, 0, 1) * t2_2(1, 0)
388  + t3dg_3(1, 1, 1) * t2_2(1, 1)
389  + t3dg_3(1, 2, 1) * t2_2(1, 2)),
390  "T3dg(i,j,k)*T2(l,j)(1,1,1)");
391  test_for_zero((t3dg_3(i, j, k) * t2_2(l, j))(1, 1, 2)
392  - (t3dg_3(1, 0, 1) * t2_2(2, 0)
393  + t3dg_3(1, 1, 1) * t2_2(2, 1)
394  + t3dg_3(1, 2, 1) * t2_2(2, 2)),
395  "T3dg(i,j,k)*T2(l,j)(1,1,2)");
396  test_for_zero((t3dg_3(i, j, k) * t2_2(l, j))(1, 2, 0)
397  - (t3dg_3(1, 0, 2) * t2_2(0, 0)
398  + t3dg_3(1, 1, 2) * t2_2(0, 1)
399  + t3dg_3(1, 2, 2) * t2_2(0, 2)),
400  "T3dg(i,j,k)*T2(l,j)(1,2,0)");
401  test_for_zero((t3dg_3(i, j, k) * t2_2(l, j))(1, 2, 1)
402  - (t3dg_3(1, 0, 2) * t2_2(1, 0)
403  + t3dg_3(1, 1, 2) * t2_2(1, 1)
404  + t3dg_3(1, 2, 2) * t2_2(1, 2)),
405  "T3dg(i,j,k)*T2(l,j)(1,2,1)");
406  test_for_zero((t3dg_3(i, j, k) * t2_2(l, j))(1, 2, 2)
407  - (t3dg_3(1, 0, 2) * t2_2(2, 0)
408  + t3dg_3(1, 1, 2) * t2_2(2, 1)
409  + t3dg_3(1, 2, 2) * t2_2(2, 2)),
410  "T3dg(i,j,k)*T2(l,j)(1,2,2)");
411  test_for_zero((t3dg_3(i, j, k) * t2_2(l, j))(2, 0, 0)
412  - (t3dg_3(2, 0, 0) * t2_2(0, 0)
413  + t3dg_3(2, 1, 0) * t2_2(0, 1)
414  + t3dg_3(2, 2, 0) * t2_2(0, 2)),
415  "T3dg(i,j,k)*T2(l,j)(2,0,0)");
416  test_for_zero((t3dg_3(i, j, k) * t2_2(l, j))(2, 0, 1)
417  - (t3dg_3(2, 0, 0) * t2_2(1, 0)
418  + t3dg_3(2, 1, 0) * t2_2(1, 1)
419  + t3dg_3(2, 2, 0) * t2_2(1, 2)),
420  "T3dg(i,j,k)*T2(l,j)(2,0,1)");
421  test_for_zero((t3dg_3(i, j, k) * t2_2(l, j))(2, 0, 2)
422  - (t3dg_3(2, 0, 0) * t2_2(2, 0)
423  + t3dg_3(2, 1, 0) * t2_2(2, 1)
424  + t3dg_3(2, 2, 0) * t2_2(2, 2)),
425  "T3dg(i,j,k)*T2(l,j)(2,0,2)");
426  test_for_zero((t3dg_3(i, j, k) * t2_2(l, j))(2, 1, 0)
427  - (t3dg_3(2, 0, 1) * t2_2(0, 0)
428  + t3dg_3(2, 1, 1) * t2_2(0, 1)
429  + t3dg_3(2, 2, 1) * t2_2(0, 2)),
430  "T3dg(i,j,k)*T2(l,j)(2,1,0)");
431  test_for_zero((t3dg_3(i, j, k) * t2_2(l, j))(2, 1, 1)
432  - (t3dg_3(2, 0, 1) * t2_2(1, 0)
433  + t3dg_3(2, 1, 1) * t2_2(1, 1)
434  + t3dg_3(2, 2, 1) * t2_2(1, 2)),
435  "T3dg(i,j,k)*T2(l,j)(2,1,1)");
436  test_for_zero((t3dg_3(i, j, k) * t2_2(l, j))(2, 1, 2)
437  - (t3dg_3(2, 0, 1) * t2_2(2, 0)
438  + t3dg_3(2, 1, 1) * t2_2(2, 1)
439  + t3dg_3(2, 2, 1) * t2_2(2, 2)),
440  "T3dg(i,j,k)*T2(l,j)(2,1,2)");
441  test_for_zero((t3dg_3(i, j, k) * t2_2(l, j))(2, 2, 0)
442  - (t3dg_3(2, 0, 2) * t2_2(0, 0)
443  + t3dg_3(2, 1, 2) * t2_2(0, 1)
444  + t3dg_3(2, 2, 2) * t2_2(0, 2)),
445  "T3dg(i,j,k)*T2(l,j)(2,2,0)");
446  test_for_zero((t3dg_3(i, j, k) * t2_2(l, j))(2, 2, 1)
447  - (t3dg_3(2, 0, 2) * t2_2(1, 0)
448  + t3dg_3(2, 1, 2) * t2_2(1, 1)
449  + t3dg_3(2, 2, 2) * t2_2(1, 2)),
450  "T3dg(i,j,k)*T2(l,j)(2,2,1)");
451  test_for_zero((t3dg_3(i, j, k) * t2_2(l, j))(2, 2, 2)
452  - (t3dg_3(2, 0, 2) * t2_2(2, 0)
453  + t3dg_3(2, 1, 2) * t2_2(2, 1)
454  + t3dg_3(2, 2, 2) * t2_2(2, 2)),
455  "T3dg(i,j,k)*T2(l,j)(2,2,2)");
456 
457  test_for_zero((t2_2(j, l) * t3dg_3(i, j, k))(0, 0, 0)
458  - (t3dg_3(0, 0, 0) * t2_2(0, 0)
459  + t3dg_3(0, 1, 0) * t2_2(1, 0)
460  + t3dg_3(0, 2, 0) * t2_2(2, 0)),
461  "T2(j,l)*T3dg(i,j,k)(0,0,0)");
462  test_for_zero((t2_2(j, l) * t3dg_3(i, j, k))(0, 0, 1)
463  - (t3dg_3(0, 0, 0) * t2_2(0, 1)
464  + t3dg_3(0, 1, 0) * t2_2(1, 1)
465  + t3dg_3(0, 2, 0) * t2_2(2, 1)),
466  "T2(j,l)*T3dg(i,j,k)(0,0,1)");
467  test_for_zero((t2_2(j, l) * t3dg_3(i, j, k))(0, 0, 2)
468  - (t3dg_3(0, 0, 0) * t2_2(0, 2)
469  + t3dg_3(0, 1, 0) * t2_2(1, 2)
470  + t3dg_3(0, 2, 0) * t2_2(2, 2)),
471  "T2(j,l)*T3dg(i,j,k)(0,0,2)");
472  test_for_zero((t2_2(j, l) * t3dg_3(i, j, k))(0, 1, 0)
473  - (t3dg_3(0, 0, 1) * t2_2(0, 0)
474  + t3dg_3(0, 1, 1) * t2_2(1, 0)
475  + t3dg_3(0, 2, 1) * t2_2(2, 0)),
476  "T2(j,l)*T3dg(i,j,k)(0,1,0)");
477  test_for_zero((t2_2(j, l) * t3dg_3(i, j, k))(0, 1, 1)
478  - (t3dg_3(0, 0, 1) * t2_2(0, 1)
479  + t3dg_3(0, 1, 1) * t2_2(1, 1)
480  + t3dg_3(0, 2, 1) * t2_2(2, 1)),
481  "T2(j,l)*T3dg(i,j,k)(0,1,1)");
482  test_for_zero((t2_2(j, l) * t3dg_3(i, j, k))(0, 1, 2)
483  - (t3dg_3(0, 0, 1) * t2_2(0, 2)
484  + t3dg_3(0, 1, 1) * t2_2(1, 2)
485  + t3dg_3(0, 2, 1) * t2_2(2, 2)),
486  "T2(j,l)*T3dg(i,j,k)(0,1,2)");
487  test_for_zero((t2_2(j, l) * t3dg_3(i, j, k))(0, 2, 0)
488  - (t3dg_3(0, 0, 2) * t2_2(0, 0)
489  + t3dg_3(0, 1, 2) * t2_2(1, 0)
490  + t3dg_3(0, 2, 2) * t2_2(2, 0)),
491  "T2(j,l)*T3dg(i,j,k)(0,2,0)");
492  test_for_zero((t2_2(j, l) * t3dg_3(i, j, k))(0, 2, 1)
493  - (t3dg_3(0, 0, 2) * t2_2(0, 1)
494  + t3dg_3(0, 1, 2) * t2_2(1, 1)
495  + t3dg_3(0, 2, 2) * t2_2(2, 1)),
496  "T2(j,l)*T3dg(i,j,k)(0,2,1)");
497  test_for_zero((t2_2(j, l) * t3dg_3(i, j, k))(0, 2, 2)
498  - (t3dg_3(0, 0, 2) * t2_2(0, 2)
499  + t3dg_3(0, 1, 2) * t2_2(1, 2)
500  + t3dg_3(0, 2, 2) * t2_2(2, 2)),
501  "T2(j,l)*T3dg(i,j,k)(0,2,2)");
502  test_for_zero((t2_2(j, l) * t3dg_3(i, j, k))(1, 0, 0)
503  - (t3dg_3(1, 0, 0) * t2_2(0, 0)
504  + t3dg_3(1, 1, 0) * t2_2(1, 0)
505  + t3dg_3(1, 2, 0) * t2_2(2, 0)),
506  "T2(j,l)*T3dg(i,j,k)(1,0,0)");
507  test_for_zero((t2_2(j, l) * t3dg_3(i, j, k))(1, 0, 1)
508  - (t3dg_3(1, 0, 0) * t2_2(0, 1)
509  + t3dg_3(1, 1, 0) * t2_2(1, 1)
510  + t3dg_3(1, 2, 0) * t2_2(2, 1)),
511  "T2(j,l)*T3dg(i,j,k)(1,0,1)");
512  test_for_zero((t2_2(j, l) * t3dg_3(i, j, k))(1, 0, 2)
513  - (t3dg_3(1, 0, 0) * t2_2(0, 2)
514  + t3dg_3(1, 1, 0) * t2_2(1, 2)
515  + t3dg_3(1, 2, 0) * t2_2(2, 2)),
516  "T2(j,l)*T3dg(i,j,k)(1,0,2)");
517  test_for_zero((t2_2(j, l) * t3dg_3(i, j, k))(1, 1, 0)
518  - (t3dg_3(1, 0, 1) * t2_2(0, 0)
519  + t3dg_3(1, 1, 1) * t2_2(1, 0)
520  + t3dg_3(1, 2, 1) * t2_2(2, 0)),
521  "T2(j,l)*T3dg(i,j,k)(1,1,0)");
522  test_for_zero((t2_2(j, l) * t3dg_3(i, j, k))(1, 1, 1)
523  - (t3dg_3(1, 0, 1) * t2_2(0, 1)
524  + t3dg_3(1, 1, 1) * t2_2(1, 1)
525  + t3dg_3(1, 2, 1) * t2_2(2, 1)),
526  "T2(j,l)*T3dg(i,j,k)(1,1,1)");
527  test_for_zero((t2_2(j, l) * t3dg_3(i, j, k))(1, 1, 2)
528  - (t3dg_3(1, 0, 1) * t2_2(0, 2)
529  + t3dg_3(1, 1, 1) * t2_2(1, 2)
530  + t3dg_3(1, 2, 1) * t2_2(2, 2)),
531  "T2(j,l)*T3dg(i,j,k)(1,1,2)");
532  test_for_zero((t2_2(j, l) * t3dg_3(i, j, k))(1, 2, 0)
533  - (t3dg_3(1, 0, 2) * t2_2(0, 0)
534  + t3dg_3(1, 1, 2) * t2_2(1, 0)
535  + t3dg_3(1, 2, 2) * t2_2(2, 0)),
536  "T2(j,l)*T3dg(i,j,k)(1,2,0)");
537  test_for_zero((t2_2(j, l) * t3dg_3(i, j, k))(1, 2, 1)
538  - (t3dg_3(1, 0, 2) * t2_2(0, 1)
539  + t3dg_3(1, 1, 2) * t2_2(1, 1)
540  + t3dg_3(1, 2, 2) * t2_2(2, 1)),
541  "T2(j,l)*T3dg(i,j,k)(1,2,1)");
542  test_for_zero((t2_2(j, l) * t3dg_3(i, j, k))(1, 2, 2)
543  - (t3dg_3(1, 0, 2) * t2_2(0, 2)
544  + t3dg_3(1, 1, 2) * t2_2(1, 2)
545  + t3dg_3(1, 2, 2) * t2_2(2, 2)),
546  "T2(j,l)*T3dg(i,j,k)(1,2,2)");
547  test_for_zero((t2_2(j, l) * t3dg_3(i, j, k))(2, 0, 0)
548  - (t3dg_3(2, 0, 0) * t2_2(0, 0)
549  + t3dg_3(2, 1, 0) * t2_2(1, 0)
550  + t3dg_3(2, 2, 0) * t2_2(2, 0)),
551  "T2(j,l)*T3dg(i,j,k)(2,0,0)");
552  test_for_zero((t2_2(j, l) * t3dg_3(i, j, k))(2, 0, 1)
553  - (t3dg_3(2, 0, 0) * t2_2(0, 1)
554  + t3dg_3(2, 1, 0) * t2_2(1, 1)
555  + t3dg_3(2, 2, 0) * t2_2(2, 1)),
556  "T2(j,l)*T3dg(i,j,k)(2,0,1)");
557  test_for_zero((t2_2(j, l) * t3dg_3(i, j, k))(2, 0, 2)
558  - (t3dg_3(2, 0, 0) * t2_2(0, 2)
559  + t3dg_3(2, 1, 0) * t2_2(1, 2)
560  + t3dg_3(2, 2, 0) * t2_2(2, 2)),
561  "T2(j,l)*T3dg(i,j,k)(2,0,2)");
562  test_for_zero((t2_2(j, l) * t3dg_3(i, j, k))(2, 1, 0)
563  - (t3dg_3(2, 0, 1) * t2_2(0, 0)
564  + t3dg_3(2, 1, 1) * t2_2(1, 0)
565  + t3dg_3(2, 2, 1) * t2_2(2, 0)),
566  "T2(j,l)*T3dg(i,j,k)(2,1,0)");
567  test_for_zero((t2_2(j, l) * t3dg_3(i, j, k))(2, 1, 1)
568  - (t3dg_3(2, 0, 1) * t2_2(0, 1)
569  + t3dg_3(2, 1, 1) * t2_2(1, 1)
570  + t3dg_3(2, 2, 1) * t2_2(2, 1)),
571  "T2(j,l)*T3dg(i,j,k)(2,1,1)");
572  test_for_zero((t2_2(j, l) * t3dg_3(i, j, k))(2, 1, 2)
573  - (t3dg_3(2, 0, 1) * t2_2(0, 2)
574  + t3dg_3(2, 1, 1) * t2_2(1, 2)
575  + t3dg_3(2, 2, 1) * t2_2(2, 2)),
576  "T2(j,l)*T3dg(i,j,k)(2,1,2)");
577  test_for_zero((t2_2(j, l) * t3dg_3(i, j, k))(2, 2, 0)
578  - (t3dg_3(2, 0, 2) * t2_2(0, 0)
579  + t3dg_3(2, 1, 2) * t2_2(1, 0)
580  + t3dg_3(2, 2, 2) * t2_2(2, 0)),
581  "T2(j,l)*T3dg(i,j,k)(2,2,0)");
582  test_for_zero((t2_2(j, l) * t3dg_3(i, j, k))(2, 2, 1)
583  - (t3dg_3(2, 0, 2) * t2_2(0, 1)
584  + t3dg_3(2, 1, 2) * t2_2(1, 1)
585  + t3dg_3(2, 2, 2) * t2_2(2, 1)),
586  "T2(j,l)*T3dg(i,j,k)(2,2,1)");
587  test_for_zero((t2_2(j, l) * t3dg_3(i, j, k))(2, 2, 2)
588  - (t3dg_3(2, 0, 2) * t2_2(0, 2)
589  + t3dg_3(2, 1, 2) * t2_2(1, 2)
590  + t3dg_3(2, 2, 2) * t2_2(2, 2)),
591  "T2(j,l)*T3dg(i,j,k)(2,2,2)");
592 
593  test_for_zero((t2_2(l, j) * t3dg_3(i, j, k))(0, 0, 0)
594  - (t3dg_3(0, 0, 0) * t2_2(0, 0)
595  + t3dg_3(0, 1, 0) * t2_2(0, 1)
596  + t3dg_3(0, 2, 0) * t2_2(0, 2)),
597  "T2(l,j)*T3dg(i,j,k)(0,0,0)");
598  test_for_zero((t2_2(l, j) * t3dg_3(i, j, k))(0, 0, 1)
599  - (t3dg_3(0, 0, 0) * t2_2(1, 0)
600  + t3dg_3(0, 1, 0) * t2_2(1, 1)
601  + t3dg_3(0, 2, 0) * t2_2(1, 2)),
602  "T2(l,j)*T3dg(i,j,k)(0,0,1)");
603  test_for_zero((t2_2(l, j) * t3dg_3(i, j, k))(0, 0, 2)
604  - (t3dg_3(0, 0, 0) * t2_2(2, 0)
605  + t3dg_3(0, 1, 0) * t2_2(2, 1)
606  + t3dg_3(0, 2, 0) * t2_2(2, 2)),
607  "T2(l,j)*T3dg(i,j,k)(0,0,2)");
608  test_for_zero((t2_2(l, j) * t3dg_3(i, j, k))(0, 1, 0)
609  - (t3dg_3(0, 0, 1) * t2_2(0, 0)
610  + t3dg_3(0, 1, 1) * t2_2(0, 1)
611  + t3dg_3(0, 2, 1) * t2_2(0, 2)),
612  "T2(l,j)*T3dg(i,j,k)(0,1,0)");
613  test_for_zero((t2_2(l, j) * t3dg_3(i, j, k))(0, 1, 1)
614  - (t3dg_3(0, 0, 1) * t2_2(1, 0)
615  + t3dg_3(0, 1, 1) * t2_2(1, 1)
616  + t3dg_3(0, 2, 1) * t2_2(1, 2)),
617  "T2(l,j)*T3dg(i,j,k)(0,1,1)");
618  test_for_zero((t2_2(l, j) * t3dg_3(i, j, k))(0, 1, 2)
619  - (t3dg_3(0, 0, 1) * t2_2(2, 0)
620  + t3dg_3(0, 1, 1) * t2_2(2, 1)
621  + t3dg_3(0, 2, 1) * t2_2(2, 2)),
622  "T2(l,j)*T3dg(i,j,k)(0,1,2)");
623  test_for_zero((t2_2(l, j) * t3dg_3(i, j, k))(0, 2, 0)
624  - (t3dg_3(0, 0, 2) * t2_2(0, 0)
625  + t3dg_3(0, 1, 2) * t2_2(0, 1)
626  + t3dg_3(0, 2, 2) * t2_2(0, 2)),
627  "T2(l,j)*T3dg(i,j,k)(0,2,0)");
628  test_for_zero((t2_2(l, j) * t3dg_3(i, j, k))(0, 2, 1)
629  - (t3dg_3(0, 0, 2) * t2_2(1, 0)
630  + t3dg_3(0, 1, 2) * t2_2(1, 1)
631  + t3dg_3(0, 2, 2) * t2_2(1, 2)),
632  "T2(l,j)*T3dg(i,j,k)(0,2,1)");
633  test_for_zero((t2_2(l, j) * t3dg_3(i, j, k))(0, 2, 2)
634  - (t3dg_3(0, 0, 2) * t2_2(2, 0)
635  + t3dg_3(0, 1, 2) * t2_2(2, 1)
636  + t3dg_3(0, 2, 2) * t2_2(2, 2)),
637  "T2(l,j)*T3dg(i,j,k)(0,2,2)");
638  test_for_zero((t2_2(l, j) * t3dg_3(i, j, k))(1, 0, 0)
639  - (t3dg_3(1, 0, 0) * t2_2(0, 0)
640  + t3dg_3(1, 1, 0) * t2_2(0, 1)
641  + t3dg_3(1, 2, 0) * t2_2(0, 2)),
642  "T2(l,j)*T3dg(i,j,k)(1,0,0)");
643  test_for_zero((t2_2(l, j) * t3dg_3(i, j, k))(1, 0, 1)
644  - (t3dg_3(1, 0, 0) * t2_2(1, 0)
645  + t3dg_3(1, 1, 0) * t2_2(1, 1)
646  + t3dg_3(1, 2, 0) * t2_2(1, 2)),
647  "T2(l,j)*T3dg(i,j,k)(1,0,1)");
648  test_for_zero((t2_2(l, j) * t3dg_3(i, j, k))(1, 0, 2)
649  - (t3dg_3(1, 0, 0) * t2_2(2, 0)
650  + t3dg_3(1, 1, 0) * t2_2(2, 1)
651  + t3dg_3(1, 2, 0) * t2_2(2, 2)),
652  "T2(l,j)*T3dg(i,j,k)(1,0,2)");
653  test_for_zero((t2_2(l, j) * t3dg_3(i, j, k))(1, 1, 0)
654  - (t3dg_3(1, 0, 1) * t2_2(0, 0)
655  + t3dg_3(1, 1, 1) * t2_2(0, 1)
656  + t3dg_3(1, 2, 1) * t2_2(0, 2)),
657  "T2(l,j)*T3dg(i,j,k)(1,1,0)");
658  test_for_zero((t2_2(l, j) * t3dg_3(i, j, k))(1, 1, 1)
659  - (t3dg_3(1, 0, 1) * t2_2(1, 0)
660  + t3dg_3(1, 1, 1) * t2_2(1, 1)
661  + t3dg_3(1, 2, 1) * t2_2(1, 2)),
662  "T2(l,j)*T3dg(i,j,k)(1,1,1)");
663  test_for_zero((t2_2(l, j) * t3dg_3(i, j, k))(1, 1, 2)
664  - (t3dg_3(1, 0, 1) * t2_2(2, 0)
665  + t3dg_3(1, 1, 1) * t2_2(2, 1)
666  + t3dg_3(1, 2, 1) * t2_2(2, 2)),
667  "T2(l,j)*T3dg(i,j,k)(1,1,2)");
668  test_for_zero((t2_2(l, j) * t3dg_3(i, j, k))(1, 2, 0)
669  - (t3dg_3(1, 0, 2) * t2_2(0, 0)
670  + t3dg_3(1, 1, 2) * t2_2(0, 1)
671  + t3dg_3(1, 2, 2) * t2_2(0, 2)),
672  "T2(l,j)*T3dg(i,j,k)(1,2,0)");
673  test_for_zero((t2_2(l, j) * t3dg_3(i, j, k))(1, 2, 1)
674  - (t3dg_3(1, 0, 2) * t2_2(1, 0)
675  + t3dg_3(1, 1, 2) * t2_2(1, 1)
676  + t3dg_3(1, 2, 2) * t2_2(1, 2)),
677  "T2(l,j)*T3dg(i,j,k)(1,2,1)");
678  test_for_zero((t2_2(l, j) * t3dg_3(i, j, k))(1, 2, 2)
679  - (t3dg_3(1, 0, 2) * t2_2(2, 0)
680  + t3dg_3(1, 1, 2) * t2_2(2, 1)
681  + t3dg_3(1, 2, 2) * t2_2(2, 2)),
682  "T2(l,j)*T3dg(i,j,k)(1,2,2)");
683  test_for_zero((t2_2(l, j) * t3dg_3(i, j, k))(2, 0, 0)
684  - (t3dg_3(2, 0, 0) * t2_2(0, 0)
685  + t3dg_3(2, 1, 0) * t2_2(0, 1)
686  + t3dg_3(2, 2, 0) * t2_2(0, 2)),
687  "T2(l,j)*T3dg(i,j,k)(2,0,0)");
688  test_for_zero((t2_2(l, j) * t3dg_3(i, j, k))(2, 0, 1)
689  - (t3dg_3(2, 0, 0) * t2_2(1, 0)
690  + t3dg_3(2, 1, 0) * t2_2(1, 1)
691  + t3dg_3(2, 2, 0) * t2_2(1, 2)),
692  "T2(l,j)*T3dg(i,j,k)(2,0,1)");
693  test_for_zero((t2_2(l, j) * t3dg_3(i, j, k))(2, 0, 2)
694  - (t3dg_3(2, 0, 0) * t2_2(2, 0)
695  + t3dg_3(2, 1, 0) * t2_2(2, 1)
696  + t3dg_3(2, 2, 0) * t2_2(2, 2)),
697  "T2(l,j)*T3dg(i,j,k)(2,0,2)");
698  test_for_zero((t2_2(l, j) * t3dg_3(i, j, k))(2, 1, 0)
699  - (t3dg_3(2, 0, 1) * t2_2(0, 0)
700  + t3dg_3(2, 1, 1) * t2_2(0, 1)
701  + t3dg_3(2, 2, 1) * t2_2(0, 2)),
702  "T2(l,j)*T3dg(i,j,k)(2,1,0)");
703  test_for_zero((t2_2(l, j) * t3dg_3(i, j, k))(2, 1, 1)
704  - (t3dg_3(2, 0, 1) * t2_2(1, 0)
705  + t3dg_3(2, 1, 1) * t2_2(1, 1)
706  + t3dg_3(2, 2, 1) * t2_2(1, 2)),
707  "T2(l,j)*T3dg(i,j,k)(2,1,1)");
708  test_for_zero((t2_2(l, j) * t3dg_3(i, j, k))(2, 1, 2)
709  - (t3dg_3(2, 0, 1) * t2_2(2, 0)
710  + t3dg_3(2, 1, 1) * t2_2(2, 1)
711  + t3dg_3(2, 2, 1) * t2_2(2, 2)),
712  "T2(l,j)*T3dg(i,j,k)(2,1,2)");
713  test_for_zero((t2_2(l, j) * t3dg_3(i, j, k))(2, 2, 0)
714  - (t3dg_3(2, 0, 2) * t2_2(0, 0)
715  + t3dg_3(2, 1, 2) * t2_2(0, 1)
716  + t3dg_3(2, 2, 2) * t2_2(0, 2)),
717  "T2(l,j)*T3dg(i,j,k)(2,2,0)");
718  test_for_zero((t2_2(l, j) * t3dg_3(i, j, k))(2, 2, 1)
719  - (t3dg_3(2, 0, 2) * t2_2(1, 0)
720  + t3dg_3(2, 1, 2) * t2_2(1, 1)
721  + t3dg_3(2, 2, 2) * t2_2(1, 2)),
722  "T2(l,j)*T3dg(i,j,k)(2,2,1)");
723  test_for_zero((t2_2(l, j) * t3dg_3(i, j, k))(2, 2, 2)
724  - (t3dg_3(2, 0, 2) * t2_2(2, 0)
725  + t3dg_3(2, 1, 2) * t2_2(2, 1)
726  + t3dg_3(2, 2, 2) * t2_2(2, 2)),
727  "T2(l,j)*T3dg(i,j,k)(2,2,2)");
728 
729  test_for_zero((t3dg_2(i, j, k) * t2s_2(j, l))(0, 0, 0)
730  - (t3dg_2(0, 0, 0) * t2s_2(0, 0)
731  + t3dg_2(0, 1, 0) * t2s_2(1, 0)
732  + t3dg_2(0, 2, 0) * t2s_2(2, 0)),
733  "T3dg(i,j,k)*T2s(j,l)(0,0,0)");
734  test_for_zero((t3dg_2(i, j, k) * t2s_2(j, l))(0, 0, 1)
735  - (t3dg_2(0, 0, 0) * t2s_2(0, 1)
736  + t3dg_2(0, 1, 0) * t2s_2(1, 1)
737  + t3dg_2(0, 2, 0) * t2s_2(2, 1)),
738  "T3dg(i,j,k)*T2s(j,l)(0,0,1)");
739  test_for_zero((t3dg_2(i, j, k) * t2s_2(j, l))(0, 0, 2)
740  - (t3dg_2(0, 0, 0) * t2s_2(0, 2)
741  + t3dg_2(0, 1, 0) * t2s_2(1, 2)
742  + t3dg_2(0, 2, 0) * t2s_2(2, 2)),
743  "T3dg(i,j,k)*T2s(j,l)(0,0,2)");
744  test_for_zero((t3dg_2(i, j, k) * t2s_2(j, l))(0, 1, 0)
745  - (t3dg_2(0, 0, 1) * t2s_2(0, 0)
746  + t3dg_2(0, 1, 1) * t2s_2(1, 0)
747  + t3dg_2(0, 2, 1) * t2s_2(2, 0)),
748  "T3dg(i,j,k)*T2s(j,l)(0,1,0)");
749  test_for_zero((t3dg_2(i, j, k) * t2s_2(j, l))(0, 1, 1)
750  - (t3dg_2(0, 0, 1) * t2s_2(0, 1)
751  + t3dg_2(0, 1, 1) * t2s_2(1, 1)
752  + t3dg_2(0, 2, 1) * t2s_2(2, 1)),
753  "T3dg(i,j,k)*T2s(j,l)(0,1,1)");
754  test_for_zero((t3dg_2(i, j, k) * t2s_2(j, l))(0, 1, 2)
755  - (t3dg_2(0, 0, 1) * t2s_2(0, 2)
756  + t3dg_2(0, 1, 1) * t2s_2(1, 2)
757  + t3dg_2(0, 2, 1) * t2s_2(2, 2)),
758  "T3dg(i,j,k)*T2s(j,l)(0,1,2)");
759  test_for_zero((t3dg_2(i, j, k) * t2s_2(j, l))(0, 2, 0)
760  - (t3dg_2(0, 0, 2) * t2s_2(0, 0)
761  + t3dg_2(0, 1, 2) * t2s_2(1, 0)
762  + t3dg_2(0, 2, 2) * t2s_2(2, 0)),
763  "T3dg(i,j,k)*T2s(j,l)(0,2,0)");
764  test_for_zero((t3dg_2(i, j, k) * t2s_2(j, l))(0, 2, 1)
765  - (t3dg_2(0, 0, 2) * t2s_2(0, 1)
766  + t3dg_2(0, 1, 2) * t2s_2(1, 1)
767  + t3dg_2(0, 2, 2) * t2s_2(2, 1)),
768  "T3dg(i,j,k)*T2s(j,l)(0,2,1)");
769  test_for_zero((t3dg_2(i, j, k) * t2s_2(j, l))(0, 2, 2)
770  - (t3dg_2(0, 0, 2) * t2s_2(0, 2)
771  + t3dg_2(0, 1, 2) * t2s_2(1, 2)
772  + t3dg_2(0, 2, 2) * t2s_2(2, 2)),
773  "T3dg(i,j,k)*T2s(j,l)(0,2,2)");
774  test_for_zero((t3dg_2(i, j, k) * t2s_2(j, l))(1, 0, 0)
775  - (t3dg_2(1, 0, 0) * t2s_2(0, 0)
776  + t3dg_2(1, 1, 0) * t2s_2(1, 0)
777  + t3dg_2(1, 2, 0) * t2s_2(2, 0)),
778  "T3dg(i,j,k)*T2s(j,l)(1,0,0)");
779  test_for_zero((t3dg_2(i, j, k) * t2s_2(j, l))(1, 0, 1)
780  - (t3dg_2(1, 0, 0) * t2s_2(0, 1)
781  + t3dg_2(1, 1, 0) * t2s_2(1, 1)
782  + t3dg_2(1, 2, 0) * t2s_2(2, 1)),
783  "T3dg(i,j,k)*T2s(j,l)(1,0,1)");
784  test_for_zero((t3dg_2(i, j, k) * t2s_2(j, l))(1, 0, 2)
785  - (t3dg_2(1, 0, 0) * t2s_2(0, 2)
786  + t3dg_2(1, 1, 0) * t2s_2(1, 2)
787  + t3dg_2(1, 2, 0) * t2s_2(2, 2)),
788  "T3dg(i,j,k)*T2s(j,l)(1,0,2)");
789  test_for_zero((t3dg_2(i, j, k) * t2s_2(j, l))(1, 1, 0)
790  - (t3dg_2(1, 0, 1) * t2s_2(0, 0)
791  + t3dg_2(1, 1, 1) * t2s_2(1, 0)
792  + t3dg_2(1, 2, 1) * t2s_2(2, 0)),
793  "T3dg(i,j,k)*T2s(j,l)(1,1,0)");
794  test_for_zero((t3dg_2(i, j, k) * t2s_2(j, l))(1, 1, 1)
795  - (t3dg_2(1, 0, 1) * t2s_2(0, 1)
796  + t3dg_2(1, 1, 1) * t2s_2(1, 1)
797  + t3dg_2(1, 2, 1) * t2s_2(2, 1)),
798  "T3dg(i,j,k)*T2s(j,l)(1,1,1)");
799  test_for_zero((t3dg_2(i, j, k) * t2s_2(j, l))(1, 1, 2)
800  - (t3dg_2(1, 0, 1) * t2s_2(0, 2)
801  + t3dg_2(1, 1, 1) * t2s_2(1, 2)
802  + t3dg_2(1, 2, 1) * t2s_2(2, 2)),
803  "T3dg(i,j,k)*T2s(j,l)(1,1,2)");
804  test_for_zero((t3dg_2(i, j, k) * t2s_2(j, l))(1, 2, 0)
805  - (t3dg_2(1, 0, 2) * t2s_2(0, 0)
806  + t3dg_2(1, 1, 2) * t2s_2(1, 0)
807  + t3dg_2(1, 2, 2) * t2s_2(2, 0)),
808  "T3dg(i,j,k)*T2s(j,l)(1,2,0)");
809  test_for_zero((t3dg_2(i, j, k) * t2s_2(j, l))(1, 2, 1)
810  - (t3dg_2(1, 0, 2) * t2s_2(0, 1)
811  + t3dg_2(1, 1, 2) * t2s_2(1, 1)
812  + t3dg_2(1, 2, 2) * t2s_2(2, 1)),
813  "T3dg(i,j,k)*T2s(j,l)(1,2,1)");
814  test_for_zero((t3dg_2(i, j, k) * t2s_2(j, l))(1, 2, 2)
815  - (t3dg_2(1, 0, 2) * t2s_2(0, 2)
816  + t3dg_2(1, 1, 2) * t2s_2(1, 2)
817  + t3dg_2(1, 2, 2) * t2s_2(2, 2)),
818  "T3dg(i,j,k)*T2s(j,l)(1,2,2)");
819  test_for_zero((t3dg_2(i, j, k) * t2s_2(j, l))(2, 0, 0)
820  - (t3dg_2(2, 0, 0) * t2s_2(0, 0)
821  + t3dg_2(2, 1, 0) * t2s_2(1, 0)
822  + t3dg_2(2, 2, 0) * t2s_2(2, 0)),
823  "T3dg(i,j,k)*T2s(j,l)(2,0,0)");
824  test_for_zero((t3dg_2(i, j, k) * t2s_2(j, l))(2, 0, 1)
825  - (t3dg_2(2, 0, 0) * t2s_2(0, 1)
826  + t3dg_2(2, 1, 0) * t2s_2(1, 1)
827  + t3dg_2(2, 2, 0) * t2s_2(2, 1)),
828  "T3dg(i,j,k)*T2s(j,l)(2,0,1)");
829  test_for_zero((t3dg_2(i, j, k) * t2s_2(j, l))(2, 0, 2)
830  - (t3dg_2(2, 0, 0) * t2s_2(0, 2)
831  + t3dg_2(2, 1, 0) * t2s_2(1, 2)
832  + t3dg_2(2, 2, 0) * t2s_2(2, 2)),
833  "T3dg(i,j,k)*T2s(j,l)(2,0,2)");
834  test_for_zero((t3dg_2(i, j, k) * t2s_2(j, l))(2, 1, 0)
835  - (t3dg_2(2, 0, 1) * t2s_2(0, 0)
836  + t3dg_2(2, 1, 1) * t2s_2(1, 0)
837  + t3dg_2(2, 2, 1) * t2s_2(2, 0)),
838  "T3dg(i,j,k)*T2s(j,l)(2,1,0)");
839  test_for_zero((t3dg_2(i, j, k) * t2s_2(j, l))(2, 1, 1)
840  - (t3dg_2(2, 0, 1) * t2s_2(0, 1)
841  + t3dg_2(2, 1, 1) * t2s_2(1, 1)
842  + t3dg_2(2, 2, 1) * t2s_2(2, 1)),
843  "T3dg(i,j,k)*T2s(j,l)(2,1,1)");
844  test_for_zero((t3dg_2(i, j, k) * t2s_2(j, l))(2, 1, 2)
845  - (t3dg_2(2, 0, 1) * t2s_2(0, 2)
846  + t3dg_2(2, 1, 1) * t2s_2(1, 2)
847  + t3dg_2(2, 2, 1) * t2s_2(2, 2)),
848  "T3dg(i,j,k)*T2s(j,l)(2,1,2)");
849  test_for_zero((t3dg_2(i, j, k) * t2s_2(j, l))(2, 2, 0)
850  - (t3dg_2(2, 0, 2) * t2s_2(0, 0)
851  + t3dg_2(2, 1, 2) * t2s_2(1, 0)
852  + t3dg_2(2, 2, 2) * t2s_2(2, 0)),
853  "T3dg(i,j,k)*T2s(j,l)(2,2,0)");
854  test_for_zero((t3dg_2(i, j, k) * t2s_2(j, l))(2, 2, 1)
855  - (t3dg_2(2, 0, 2) * t2s_2(0, 1)
856  + t3dg_2(2, 1, 2) * t2s_2(1, 1)
857  + t3dg_2(2, 2, 2) * t2s_2(2, 1)),
858  "T3dg(i,j,k)*T2s(j,l)(2,2,1)");
859  test_for_zero((t3dg_2(i, j, k) * t2s_2(j, l))(2, 2, 2)
860  - (t3dg_2(2, 0, 2) * t2s_2(0, 2)
861  + t3dg_2(2, 1, 2) * t2s_2(1, 2)
862  + t3dg_2(2, 2, 2) * t2s_2(2, 2)),
863  "T3dg(i,j,k)*T2s(j,l)(2,2,2)");
864 
865  test_for_zero((t3dg_3(i, j, k) * t2s_2(l, j))(0, 0, 0)
866  - (t3dg_3(0, 0, 0) * t2s_2(0, 0)
867  + t3dg_3(0, 1, 0) * t2s_2(0, 1)
868  + t3dg_3(0, 2, 0) * t2s_2(0, 2)),
869  "T3dg(i,j,k)*T2s(l,j)(0,0,0)");
870  test_for_zero((t3dg_3(i, j, k) * t2s_2(l, j))(0, 0, 1)
871  - (t3dg_3(0, 0, 0) * t2s_2(1, 0)
872  + t3dg_3(0, 1, 0) * t2s_2(1, 1)
873  + t3dg_3(0, 2, 0) * t2s_2(1, 2)),
874  "T3dg(i,j,k)*T2s(l,j)(0,0,1)");
875  test_for_zero((t3dg_3(i, j, k) * t2s_2(l, j))(0, 0, 2)
876  - (t3dg_3(0, 0, 0) * t2s_2(2, 0)
877  + t3dg_3(0, 1, 0) * t2s_2(2, 1)
878  + t3dg_3(0, 2, 0) * t2s_2(2, 2)),
879  "T3dg(i,j,k)*T2s(l,j)(0,0,2)");
880  test_for_zero((t3dg_3(i, j, k) * t2s_2(l, j))(0, 1, 0)
881  - (t3dg_3(0, 0, 1) * t2s_2(0, 0)
882  + t3dg_3(0, 1, 1) * t2s_2(0, 1)
883  + t3dg_3(0, 2, 1) * t2s_2(0, 2)),
884  "T3dg(i,j,k)*T2s(l,j)(0,1,0)");
885  test_for_zero((t3dg_3(i, j, k) * t2s_2(l, j))(0, 1, 1)
886  - (t3dg_3(0, 0, 1) * t2s_2(1, 0)
887  + t3dg_3(0, 1, 1) * t2s_2(1, 1)
888  + t3dg_3(0, 2, 1) * t2s_2(1, 2)),
889  "T3dg(i,j,k)*T2s(l,j)(0,1,1)");
890  test_for_zero((t3dg_3(i, j, k) * t2s_2(l, j))(0, 1, 2)
891  - (t3dg_3(0, 0, 1) * t2s_2(2, 0)
892  + t3dg_3(0, 1, 1) * t2s_2(2, 1)
893  + t3dg_3(0, 2, 1) * t2s_2(2, 2)),
894  "T3dg(i,j,k)*T2s(l,j)(0,1,2)");
895  test_for_zero((t3dg_3(i, j, k) * t2s_2(l, j))(0, 2, 0)
896  - (t3dg_3(0, 0, 2) * t2s_2(0, 0)
897  + t3dg_3(0, 1, 2) * t2s_2(0, 1)
898  + t3dg_3(0, 2, 2) * t2s_2(0, 2)),
899  "T3dg(i,j,k)*T2s(l,j)(0,2,0)");
900  test_for_zero((t3dg_3(i, j, k) * t2s_2(l, j))(0, 2, 1)
901  - (t3dg_3(0, 0, 2) * t2s_2(1, 0)
902  + t3dg_3(0, 1, 2) * t2s_2(1, 1)
903  + t3dg_3(0, 2, 2) * t2s_2(1, 2)),
904  "T3dg(i,j,k)*T2s(l,j)(0,2,1)");
905  test_for_zero((t3dg_3(i, j, k) * t2s_2(l, j))(0, 2, 2)
906  - (t3dg_3(0, 0, 2) * t2s_2(2, 0)
907  + t3dg_3(0, 1, 2) * t2s_2(2, 1)
908  + t3dg_3(0, 2, 2) * t2s_2(2, 2)),
909  "T3dg(i,j,k)*T2s(l,j)(0,2,2)");
910  test_for_zero((t3dg_3(i, j, k) * t2s_2(l, j))(1, 0, 0)
911  - (t3dg_3(1, 0, 0) * t2s_2(0, 0)
912  + t3dg_3(1, 1, 0) * t2s_2(0, 1)
913  + t3dg_3(1, 2, 0) * t2s_2(0, 2)),
914  "T3dg(i,j,k)*T2s(l,j)(1,0,0)");
915  test_for_zero((t3dg_3(i, j, k) * t2s_2(l, j))(1, 0, 1)
916  - (t3dg_3(1, 0, 0) * t2s_2(1, 0)
917  + t3dg_3(1, 1, 0) * t2s_2(1, 1)
918  + t3dg_3(1, 2, 0) * t2s_2(1, 2)),
919  "T3dg(i,j,k)*T2s(l,j)(1,0,1)");
920  test_for_zero((t3dg_3(i, j, k) * t2s_2(l, j))(1, 0, 2)
921  - (t3dg_3(1, 0, 0) * t2s_2(2, 0)
922  + t3dg_3(1, 1, 0) * t2s_2(2, 1)
923  + t3dg_3(1, 2, 0) * t2s_2(2, 2)),
924  "T3dg(i,j,k)*T2s(l,j)(1,0,2)");
925  test_for_zero((t3dg_3(i, j, k) * t2s_2(l, j))(1, 1, 0)
926  - (t3dg_3(1, 0, 1) * t2s_2(0, 0)
927  + t3dg_3(1, 1, 1) * t2s_2(0, 1)
928  + t3dg_3(1, 2, 1) * t2s_2(0, 2)),
929  "T3dg(i,j,k)*T2s(l,j)(1,1,0)");
930  test_for_zero((t3dg_3(i, j, k) * t2s_2(l, j))(1, 1, 1)
931  - (t3dg_3(1, 0, 1) * t2s_2(1, 0)
932  + t3dg_3(1, 1, 1) * t2s_2(1, 1)
933  + t3dg_3(1, 2, 1) * t2s_2(1, 2)),
934  "T3dg(i,j,k)*T2s(l,j)(1,1,1)");
935  test_for_zero((t3dg_3(i, j, k) * t2s_2(l, j))(1, 1, 2)
936  - (t3dg_3(1, 0, 1) * t2s_2(2, 0)
937  + t3dg_3(1, 1, 1) * t2s_2(2, 1)
938  + t3dg_3(1, 2, 1) * t2s_2(2, 2)),
939  "T3dg(i,j,k)*T2s(l,j)(1,1,2)");
940  test_for_zero((t3dg_3(i, j, k) * t2s_2(l, j))(1, 2, 0)
941  - (t3dg_3(1, 0, 2) * t2s_2(0, 0)
942  + t3dg_3(1, 1, 2) * t2s_2(0, 1)
943  + t3dg_3(1, 2, 2) * t2s_2(0, 2)),
944  "T3dg(i,j,k)*T2s(l,j)(1,2,0)");
945  test_for_zero((t3dg_3(i, j, k) * t2s_2(l, j))(1, 2, 1)
946  - (t3dg_3(1, 0, 2) * t2s_2(1, 0)
947  + t3dg_3(1, 1, 2) * t2s_2(1, 1)
948  + t3dg_3(1, 2, 2) * t2s_2(1, 2)),
949  "T3dg(i,j,k)*T2s(l,j)(1,2,1)");
950  test_for_zero((t3dg_3(i, j, k) * t2s_2(l, j))(1, 2, 2)
951  - (t3dg_3(1, 0, 2) * t2s_2(2, 0)
952  + t3dg_3(1, 1, 2) * t2s_2(2, 1)
953  + t3dg_3(1, 2, 2) * t2s_2(2, 2)),
954  "T3dg(i,j,k)*T2s(l,j)(1,2,2)");
955  test_for_zero((t3dg_3(i, j, k) * t2s_2(l, j))(2, 0, 0)
956  - (t3dg_3(2, 0, 0) * t2s_2(0, 0)
957  + t3dg_3(2, 1, 0) * t2s_2(0, 1)
958  + t3dg_3(2, 2, 0) * t2s_2(0, 2)),
959  "T3dg(i,j,k)*T2s(l,j)(2,0,0)");
960  test_for_zero((t3dg_3(i, j, k) * t2s_2(l, j))(2, 0, 1)
961  - (t3dg_3(2, 0, 0) * t2s_2(1, 0)
962  + t3dg_3(2, 1, 0) * t2s_2(1, 1)
963  + t3dg_3(2, 2, 0) * t2s_2(1, 2)),
964  "T3dg(i,j,k)*T2s(l,j)(2,0,1)");
965  test_for_zero((t3dg_3(i, j, k) * t2s_2(l, j))(2, 0, 2)
966  - (t3dg_3(2, 0, 0) * t2s_2(2, 0)
967  + t3dg_3(2, 1, 0) * t2s_2(2, 1)
968  + t3dg_3(2, 2, 0) * t2s_2(2, 2)),
969  "T3dg(i,j,k)*T2s(l,j)(2,0,2)");
970  test_for_zero((t3dg_3(i, j, k) * t2s_2(l, j))(2, 1, 0)
971  - (t3dg_3(2, 0, 1) * t2s_2(0, 0)
972  + t3dg_3(2, 1, 1) * t2s_2(0, 1)
973  + t3dg_3(2, 2, 1) * t2s_2(0, 2)),
974  "T3dg(i,j,k)*T2s(l,j)(2,1,0)");
975  test_for_zero((t3dg_3(i, j, k) * t2s_2(l, j))(2, 1, 1)
976  - (t3dg_3(2, 0, 1) * t2s_2(1, 0)
977  + t3dg_3(2, 1, 1) * t2s_2(1, 1)
978  + t3dg_3(2, 2, 1) * t2s_2(1, 2)),
979  "T3dg(i,j,k)*T2s(l,j)(2,1,1)");
980  test_for_zero((t3dg_3(i, j, k) * t2s_2(l, j))(2, 1, 2)
981  - (t3dg_3(2, 0, 1) * t2s_2(2, 0)
982  + t3dg_3(2, 1, 1) * t2s_2(2, 1)
983  + t3dg_3(2, 2, 1) * t2s_2(2, 2)),
984  "T3dg(i,j,k)*T2s(l,j)(2,1,2)");
985  test_for_zero((t3dg_3(i, j, k) * t2s_2(l, j))(2, 2, 0)
986  - (t3dg_3(2, 0, 2) * t2s_2(0, 0)
987  + t3dg_3(2, 1, 2) * t2s_2(0, 1)
988  + t3dg_3(2, 2, 2) * t2s_2(0, 2)),
989  "T3dg(i,j,k)*T2s(l,j)(2,2,0)");
990  test_for_zero((t3dg_3(i, j, k) * t2s_2(l, j))(2, 2, 1)
991  - (t3dg_3(2, 0, 2) * t2s_2(1, 0)
992  + t3dg_3(2, 1, 2) * t2s_2(1, 1)
993  + t3dg_3(2, 2, 2) * t2s_2(1, 2)),
994  "T3dg(i,j,k)*T2s(l,j)(2,2,1)");
995  test_for_zero((t3dg_3(i, j, k) * t2s_2(l, j))(2, 2, 2)
996  - (t3dg_3(2, 0, 2) * t2s_2(2, 0)
997  + t3dg_3(2, 1, 2) * t2s_2(2, 1)
998  + t3dg_3(2, 2, 2) * t2s_2(2, 2)),
999  "T3dg(i,j,k)*T2s(l,j)(2,2,2)");
1000 
1001  test_for_zero((t2s_2(j, l) * t3dg_3(i, j, k))(0, 0, 0)
1002  - (t3dg_3(0, 0, 0) * t2s_2(0, 0)
1003  + t3dg_3(0, 1, 0) * t2s_2(1, 0)
1004  + t3dg_3(0, 2, 0) * t2s_2(2, 0)),
1005  "T2s(j,l)*T3dg(i,j,k)(0,0,0)");
1006  test_for_zero((t2s_2(j, l) * t3dg_3(i, j, k))(0, 0, 1)
1007  - (t3dg_3(0, 0, 0) * t2s_2(0, 1)
1008  + t3dg_3(0, 1, 0) * t2s_2(1, 1)
1009  + t3dg_3(0, 2, 0) * t2s_2(2, 1)),
1010  "T2s(j,l)*T3dg(i,j,k)(0,0,1)");
1011  test_for_zero((t2s_2(j, l) * t3dg_3(i, j, k))(0, 0, 2)
1012  - (t3dg_3(0, 0, 0) * t2s_2(0, 2)
1013  + t3dg_3(0, 1, 0) * t2s_2(1, 2)
1014  + t3dg_3(0, 2, 0) * t2s_2(2, 2)),
1015  "T2s(j,l)*T3dg(i,j,k)(0,0,2)");
1016  test_for_zero((t2s_2(j, l) * t3dg_3(i, j, k))(0, 1, 0)
1017  - (t3dg_3(0, 0, 1) * t2s_2(0, 0)
1018  + t3dg_3(0, 1, 1) * t2s_2(1, 0)
1019  + t3dg_3(0, 2, 1) * t2s_2(2, 0)),
1020  "T2s(j,l)*T3dg(i,j,k)(0,1,0)");
1021  test_for_zero((t2s_2(j, l) * t3dg_3(i, j, k))(0, 1, 1)
1022  - (t3dg_3(0, 0, 1) * t2s_2(0, 1)
1023  + t3dg_3(0, 1, 1) * t2s_2(1, 1)
1024  + t3dg_3(0, 2, 1) * t2s_2(2, 1)),
1025  "T2s(j,l)*T3dg(i,j,k)(0,1,1)");
1026  test_for_zero((t2s_2(j, l) * t3dg_3(i, j, k))(0, 1, 2)
1027  - (t3dg_3(0, 0, 1) * t2s_2(0, 2)
1028  + t3dg_3(0, 1, 1) * t2s_2(1, 2)
1029  + t3dg_3(0, 2, 1) * t2s_2(2, 2)),
1030  "T2s(j,l)*T3dg(i,j,k)(0,1,2)");
1031  test_for_zero((t2s_2(j, l) * t3dg_3(i, j, k))(0, 2, 0)
1032  - (t3dg_3(0, 0, 2) * t2s_2(0, 0)
1033  + t3dg_3(0, 1, 2) * t2s_2(1, 0)
1034  + t3dg_3(0, 2, 2) * t2s_2(2, 0)),
1035  "T2s(j,l)*T3dg(i,j,k)(0,2,0)");
1036  test_for_zero((t2s_2(j, l) * t3dg_3(i, j, k))(0, 2, 1)
1037  - (t3dg_3(0, 0, 2) * t2s_2(0, 1)
1038  + t3dg_3(0, 1, 2) * t2s_2(1, 1)
1039  + t3dg_3(0, 2, 2) * t2s_2(2, 1)),
1040  "T2s(j,l)*T3dg(i,j,k)(0,2,1)");
1041  test_for_zero((t2s_2(j, l) * t3dg_3(i, j, k))(0, 2, 2)
1042  - (t3dg_3(0, 0, 2) * t2s_2(0, 2)
1043  + t3dg_3(0, 1, 2) * t2s_2(1, 2)
1044  + t3dg_3(0, 2, 2) * t2s_2(2, 2)),
1045  "T2s(j,l)*T3dg(i,j,k)(0,2,2)");
1046  test_for_zero((t2s_2(j, l) * t3dg_3(i, j, k))(1, 0, 0)
1047  - (t3dg_3(1, 0, 0) * t2s_2(0, 0)
1048  + t3dg_3(1, 1, 0) * t2s_2(1, 0)
1049  + t3dg_3(1, 2, 0) * t2s_2(2, 0)),
1050  "T2s(j,l)*T3dg(i,j,k)(1,0,0)");
1051  test_for_zero((t2s_2(j, l) * t3dg_3(i, j, k))(1, 0, 1)
1052  - (t3dg_3(1, 0, 0) * t2s_2(0, 1)
1053  + t3dg_3(1, 1, 0) * t2s_2(1, 1)
1054  + t3dg_3(1, 2, 0) * t2s_2(2, 1)),
1055  "T2s(j,l)*T3dg(i,j,k)(1,0,1)");
1056  test_for_zero((t2s_2(j, l) * t3dg_3(i, j, k))(1, 0, 2)
1057  - (t3dg_3(1, 0, 0) * t2s_2(0, 2)
1058  + t3dg_3(1, 1, 0) * t2s_2(1, 2)
1059  + t3dg_3(1, 2, 0) * t2s_2(2, 2)),
1060  "T2s(j,l)*T3dg(i,j,k)(1,0,2)");
1061  test_for_zero((t2s_2(j, l) * t3dg_3(i, j, k))(1, 1, 0)
1062  - (t3dg_3(1, 0, 1) * t2s_2(0, 0)
1063  + t3dg_3(1, 1, 1) * t2s_2(1, 0)
1064  + t3dg_3(1, 2, 1) * t2s_2(2, 0)),
1065  "T2s(j,l)*T3dg(i,j,k)(1,1,0)");
1066  test_for_zero((t2s_2(j, l) * t3dg_3(i, j, k))(1, 1, 1)
1067  - (t3dg_3(1, 0, 1) * t2s_2(0, 1)
1068  + t3dg_3(1, 1, 1) * t2s_2(1, 1)
1069  + t3dg_3(1, 2, 1) * t2s_2(2, 1)),
1070  "T2s(j,l)*T3dg(i,j,k)(1,1,1)");
1071  test_for_zero((t2s_2(j, l) * t3dg_3(i, j, k))(1, 1, 2)
1072  - (t3dg_3(1, 0, 1) * t2s_2(0, 2)
1073  + t3dg_3(1, 1, 1) * t2s_2(1, 2)
1074  + t3dg_3(1, 2, 1) * t2s_2(2, 2)),
1075  "T2s(j,l)*T3dg(i,j,k)(1,1,2)");
1076  test_for_zero((t2s_2(j, l) * t3dg_3(i, j, k))(1, 2, 0)
1077  - (t3dg_3(1, 0, 2) * t2s_2(0, 0)
1078  + t3dg_3(1, 1, 2) * t2s_2(1, 0)
1079  + t3dg_3(1, 2, 2) * t2s_2(2, 0)),
1080  "T2s(j,l)*T3dg(i,j,k)(1,2,0)");
1081  test_for_zero((t2s_2(j, l) * t3dg_3(i, j, k))(1, 2, 1)
1082  - (t3dg_3(1, 0, 2) * t2s_2(0, 1)
1083  + t3dg_3(1, 1, 2) * t2s_2(1, 1)
1084  + t3dg_3(1, 2, 2) * t2s_2(2, 1)),
1085  "T2s(j,l)*T3dg(i,j,k)(1,2,1)");
1086  test_for_zero((t2s_2(j, l) * t3dg_3(i, j, k))(1, 2, 2)
1087  - (t3dg_3(1, 0, 2) * t2s_2(0, 2)
1088  + t3dg_3(1, 1, 2) * t2s_2(1, 2)
1089  + t3dg_3(1, 2, 2) * t2s_2(2, 2)),
1090  "T2s(j,l)*T3dg(i,j,k)(1,2,2)");
1091  test_for_zero((t2s_2(j, l) * t3dg_3(i, j, k))(2, 0, 0)
1092  - (t3dg_3(2, 0, 0) * t2s_2(0, 0)
1093  + t3dg_3(2, 1, 0) * t2s_2(1, 0)
1094  + t3dg_3(2, 2, 0) * t2s_2(2, 0)),
1095  "T2s(j,l)*T3dg(i,j,k)(2,0,0)");
1096  test_for_zero((t2s_2(j, l) * t3dg_3(i, j, k))(2, 0, 1)
1097  - (t3dg_3(2, 0, 0) * t2s_2(0, 1)
1098  + t3dg_3(2, 1, 0) * t2s_2(1, 1)
1099  + t3dg_3(2, 2, 0) * t2s_2(2, 1)),
1100  "T2s(j,l)*T3dg(i,j,k)(2,0,1)");
1101  test_for_zero((t2s_2(j, l) * t3dg_3(i, j, k))(2, 0, 2)
1102  - (t3dg_3(2, 0, 0) * t2s_2(0, 2)
1103  + t3dg_3(2, 1, 0) * t2s_2(1, 2)
1104  + t3dg_3(2, 2, 0) * t2s_2(2, 2)),
1105  "T2s(j,l)*T3dg(i,j,k)(2,0,2)");
1106  test_for_zero((t2s_2(j, l) * t3dg_3(i, j, k))(2, 1, 0)
1107  - (t3dg_3(2, 0, 1) * t2s_2(0, 0)
1108  + t3dg_3(2, 1, 1) * t2s_2(1, 0)
1109  + t3dg_3(2, 2, 1) * t2s_2(2, 0)),
1110  "T2s(j,l)*T3dg(i,j,k)(2,1,0)");
1111  test_for_zero((t2s_2(j, l) * t3dg_3(i, j, k))(2, 1, 1)
1112  - (t3dg_3(2, 0, 1) * t2s_2(0, 1)
1113  + t3dg_3(2, 1, 1) * t2s_2(1, 1)
1114  + t3dg_3(2, 2, 1) * t2s_2(2, 1)),
1115  "T2s(j,l)*T3dg(i,j,k)(2,1,1)");
1116  test_for_zero((t2s_2(j, l) * t3dg_3(i, j, k))(2, 1, 2)
1117  - (t3dg_3(2, 0, 1) * t2s_2(0, 2)
1118  + t3dg_3(2, 1, 1) * t2s_2(1, 2)
1119  + t3dg_3(2, 2, 1) * t2s_2(2, 2)),
1120  "T2s(j,l)*T3dg(i,j,k)(2,1,2)");
1121  test_for_zero((t2s_2(j, l) * t3dg_3(i, j, k))(2, 2, 0)
1122  - (t3dg_3(2, 0, 2) * t2s_2(0, 0)
1123  + t3dg_3(2, 1, 2) * t2s_2(1, 0)
1124  + t3dg_3(2, 2, 2) * t2s_2(2, 0)),
1125  "T2s(j,l)*T3dg(i,j,k)(2,2,0)");
1126  test_for_zero((t2s_2(j, l) * t3dg_3(i, j, k))(2, 2, 1)
1127  - (t3dg_3(2, 0, 2) * t2s_2(0, 1)
1128  + t3dg_3(2, 1, 2) * t2s_2(1, 1)
1129  + t3dg_3(2, 2, 2) * t2s_2(2, 1)),
1130  "T2s(j,l)*T3dg(i,j,k)(2,2,1)");
1131  test_for_zero((t2s_2(j, l) * t3dg_3(i, j, k))(2, 2, 2)
1132  - (t3dg_3(2, 0, 2) * t2s_2(0, 2)
1133  + t3dg_3(2, 1, 2) * t2s_2(1, 2)
1134  + t3dg_3(2, 2, 2) * t2s_2(2, 2)),
1135  "T2s(j,l)*T3dg(i,j,k)(2,2,2)");
1136 
1137  test_for_zero((t2s_2(l, j) * t3dg_3(i, j, k))(0, 0, 0)
1138  - (t3dg_3(0, 0, 0) * t2s_2(0, 0)
1139  + t3dg_3(0, 1, 0) * t2s_2(0, 1)
1140  + t3dg_3(0, 2, 0) * t2s_2(0, 2)),
1141  "T2s(l,j)*T3dg(i,j,k)(0,0,0)");
1142  test_for_zero((t2s_2(l, j) * t3dg_3(i, j, k))(0, 0, 1)
1143  - (t3dg_3(0, 0, 0) * t2s_2(1, 0)
1144  + t3dg_3(0, 1, 0) * t2s_2(1, 1)
1145  + t3dg_3(0, 2, 0) * t2s_2(1, 2)),
1146  "T2s(l,j)*T3dg(i,j,k)(0,0,1)");
1147  test_for_zero((t2s_2(l, j) * t3dg_3(i, j, k))(0, 0, 2)
1148  - (t3dg_3(0, 0, 0) * t2s_2(2, 0)
1149  + t3dg_3(0, 1, 0) * t2s_2(2, 1)
1150  + t3dg_3(0, 2, 0) * t2s_2(2, 2)),
1151  "T2s(l,j)*T3dg(i,j,k)(0,0,2)");
1152  test_for_zero((t2s_2(l, j) * t3dg_3(i, j, k))(0, 1, 0)
1153  - (t3dg_3(0, 0, 1) * t2s_2(0, 0)
1154  + t3dg_3(0, 1, 1) * t2s_2(0, 1)
1155  + t3dg_3(0, 2, 1) * t2s_2(0, 2)),
1156  "T2s(l,j)*T3dg(i,j,k)(0,1,0)");
1157  test_for_zero((t2s_2(l, j) * t3dg_3(i, j, k))(0, 1, 1)
1158  - (t3dg_3(0, 0, 1) * t2s_2(1, 0)
1159  + t3dg_3(0, 1, 1) * t2s_2(1, 1)
1160  + t3dg_3(0, 2, 1) * t2s_2(1, 2)),
1161  "T2s(l,j)*T3dg(i,j,k)(0,1,1)");
1162  test_for_zero((t2s_2(l, j) * t3dg_3(i, j, k))(0, 1, 2)
1163  - (t3dg_3(0, 0, 1) * t2s_2(2, 0)
1164  + t3dg_3(0, 1, 1) * t2s_2(2, 1)
1165  + t3dg_3(0, 2, 1) * t2s_2(2, 2)),
1166  "T2s(l,j)*T3dg(i,j,k)(0,1,2)");
1167  test_for_zero((t2s_2(l, j) * t3dg_3(i, j, k))(0, 2, 0)
1168  - (t3dg_3(0, 0, 2) * t2s_2(0, 0)
1169  + t3dg_3(0, 1, 2) * t2s_2(0, 1)
1170  + t3dg_3(0, 2, 2) * t2s_2(0, 2)),
1171  "T2s(l,j)*T3dg(i,j,k)(0,2,0)");
1172  test_for_zero((t2s_2(l, j) * t3dg_3(i, j, k))(0, 2, 1)
1173  - (t3dg_3(0, 0, 2) * t2s_2(1, 0)
1174  + t3dg_3(0, 1, 2) * t2s_2(1, 1)
1175  + t3dg_3(0, 2, 2) * t2s_2(1, 2)),
1176  "T2s(l,j)*T3dg(i,j,k)(0,2,1)");
1177  test_for_zero((t2s_2(l, j) * t3dg_3(i, j, k))(0, 2, 2)
1178  - (t3dg_3(0, 0, 2) * t2s_2(2, 0)
1179  + t3dg_3(0, 1, 2) * t2s_2(2, 1)
1180  + t3dg_3(0, 2, 2) * t2s_2(2, 2)),
1181  "T2s(l,j)*T3dg(i,j,k)(0,2,2)");
1182  test_for_zero((t2s_2(l, j) * t3dg_3(i, j, k))(1, 0, 0)
1183  - (t3dg_3(1, 0, 0) * t2s_2(0, 0)
1184  + t3dg_3(1, 1, 0) * t2s_2(0, 1)
1185  + t3dg_3(1, 2, 0) * t2s_2(0, 2)),
1186  "T2s(l,j)*T3dg(i,j,k)(1,0,0)");
1187  test_for_zero((t2s_2(l, j) * t3dg_3(i, j, k))(1, 0, 1)
1188  - (t3dg_3(1, 0, 0) * t2s_2(1, 0)
1189  + t3dg_3(1, 1, 0) * t2s_2(1, 1)
1190  + t3dg_3(1, 2, 0) * t2s_2(1, 2)),
1191  "T2s(l,j)*T3dg(i,j,k)(1,0,1)");
1192  test_for_zero((t2s_2(l, j) * t3dg_3(i, j, k))(1, 0, 2)
1193  - (t3dg_3(1, 0, 0) * t2s_2(2, 0)
1194  + t3dg_3(1, 1, 0) * t2s_2(2, 1)
1195  + t3dg_3(1, 2, 0) * t2s_2(2, 2)),
1196  "T2s(l,j)*T3dg(i,j,k)(1,0,2)");
1197  test_for_zero((t2s_2(l, j) * t3dg_3(i, j, k))(1, 1, 0)
1198  - (t3dg_3(1, 0, 1) * t2s_2(0, 0)
1199  + t3dg_3(1, 1, 1) * t2s_2(0, 1)
1200  + t3dg_3(1, 2, 1) * t2s_2(0, 2)),
1201  "T2s(l,j)*T3dg(i,j,k)(1,1,0)");
1202  test_for_zero((t2s_2(l, j) * t3dg_3(i, j, k))(1, 1, 1)
1203  - (t3dg_3(1, 0, 1) * t2s_2(1, 0)
1204  + t3dg_3(1, 1, 1) * t2s_2(1, 1)
1205  + t3dg_3(1, 2, 1) * t2s_2(1, 2)),
1206  "T2s(l,j)*T3dg(i,j,k)(1,1,1)");
1207  test_for_zero((t2s_2(l, j) * t3dg_3(i, j, k))(1, 1, 2)
1208  - (t3dg_3(1, 0, 1) * t2s_2(2, 0)
1209  + t3dg_3(1, 1, 1) * t2s_2(2, 1)
1210  + t3dg_3(1, 2, 1) * t2s_2(2, 2)),
1211  "T2s(l,j)*T3dg(i,j,k)(1,1,2)");
1212  test_for_zero((t2s_2(l, j) * t3dg_3(i, j, k))(1, 2, 0)
1213  - (t3dg_3(1, 0, 2) * t2s_2(0, 0)
1214  + t3dg_3(1, 1, 2) * t2s_2(0, 1)
1215  + t3dg_3(1, 2, 2) * t2s_2(0, 2)),
1216  "T2s(l,j)*T3dg(i,j,k)(1,2,0)");
1217  test_for_zero((t2s_2(l, j) * t3dg_3(i, j, k))(1, 2, 1)
1218  - (t3dg_3(1, 0, 2) * t2s_2(1, 0)
1219  + t3dg_3(1, 1, 2) * t2s_2(1, 1)
1220  + t3dg_3(1, 2, 2) * t2s_2(1, 2)),
1221  "T2s(l,j)*T3dg(i,j,k)(1,2,1)");
1222  test_for_zero((t2s_2(l, j) * t3dg_3(i, j, k))(1, 2, 2)
1223  - (t3dg_3(1, 0, 2) * t2s_2(2, 0)
1224  + t3dg_3(1, 1, 2) * t2s_2(2, 1)
1225  + t3dg_3(1, 2, 2) * t2s_2(2, 2)),
1226  "T2s(l,j)*T3dg(i,j,k)(1,2,2)");
1227  test_for_zero((t2s_2(l, j) * t3dg_3(i, j, k))(2, 0, 0)
1228  - (t3dg_3(2, 0, 0) * t2s_2(0, 0)
1229  + t3dg_3(2, 1, 0) * t2s_2(0, 1)
1230  + t3dg_3(2, 2, 0) * t2s_2(0, 2)),
1231  "T2s(l,j)*T3dg(i,j,k)(2,0,0)");
1232  test_for_zero((t2s_2(l, j) * t3dg_3(i, j, k))(2, 0, 1)
1233  - (t3dg_3(2, 0, 0) * t2s_2(1, 0)
1234  + t3dg_3(2, 1, 0) * t2s_2(1, 1)
1235  + t3dg_3(2, 2, 0) * t2s_2(1, 2)),
1236  "T2s(l,j)*T3dg(i,j,k)(2,0,1)");
1237  test_for_zero((t2s_2(l, j) * t3dg_3(i, j, k))(2, 0, 2)
1238  - (t3dg_3(2, 0, 0) * t2s_2(2, 0)
1239  + t3dg_3(2, 1, 0) * t2s_2(2, 1)
1240  + t3dg_3(2, 2, 0) * t2s_2(2, 2)),
1241  "T2s(l,j)*T3dg(i,j,k)(2,0,2)");
1242  test_for_zero((t2s_2(l, j) * t3dg_3(i, j, k))(2, 1, 0)
1243  - (t3dg_3(2, 0, 1) * t2s_2(0, 0)
1244  + t3dg_3(2, 1, 1) * t2s_2(0, 1)
1245  + t3dg_3(2, 2, 1) * t2s_2(0, 2)),
1246  "T2s(l,j)*T3dg(i,j,k)(2,1,0)");
1247  test_for_zero((t2s_2(l, j) * t3dg_3(i, j, k))(2, 1, 1)
1248  - (t3dg_3(2, 0, 1) * t2s_2(1, 0)
1249  + t3dg_3(2, 1, 1) * t2s_2(1, 1)
1250  + t3dg_3(2, 2, 1) * t2s_2(1, 2)),
1251  "T2s(l,j)*T3dg(i,j,k)(2,1,1)");
1252  test_for_zero((t2s_2(l, j) * t3dg_3(i, j, k))(2, 1, 2)
1253  - (t3dg_3(2, 0, 1) * t2s_2(2, 0)
1254  + t3dg_3(2, 1, 1) * t2s_2(2, 1)
1255  + t3dg_3(2, 2, 1) * t2s_2(2, 2)),
1256  "T2s(l,j)*T3dg(i,j,k)(2,1,2)");
1257  test_for_zero((t2s_2(l, j) * t3dg_3(i, j, k))(2, 2, 0)
1258  - (t3dg_3(2, 0, 2) * t2s_2(0, 0)
1259  + t3dg_3(2, 1, 2) * t2s_2(0, 1)
1260  + t3dg_3(2, 2, 2) * t2s_2(0, 2)),
1261  "T2s(l,j)*T3dg(i,j,k)(2,2,0)");
1262  test_for_zero((t2s_2(l, j) * t3dg_3(i, j, k))(2, 2, 1)
1263  - (t3dg_3(2, 0, 2) * t2s_2(1, 0)
1264  + t3dg_3(2, 1, 2) * t2s_2(1, 1)
1265  + t3dg_3(2, 2, 2) * t2s_2(1, 2)),
1266  "T2s(l,j)*T3dg(i,j,k)(2,2,1)");
1267  test_for_zero((t2s_2(l, j) * t3dg_3(i, j, k))(2, 2, 2)
1268  - (t3dg_3(2, 0, 2) * t2s_2(2, 0)
1269  + t3dg_3(2, 1, 2) * t2s_2(2, 1)
1270  + t3dg_3(2, 2, 2) * t2s_2(2, 2)),
1271  "T2s(l,j)*T3dg(i,j,k)(2,2,2)");
1272 }
FTensor
JSON compatible output.
Definition: Christof_constructor.hpp:6
test_T3dg_37
void test_T3dg_37(const Tensor2< double, 3, 3 > &t2_2, const Tensor2_symmetric< double, 3 > &t2s_2, const Dg< double, 3, 3 > &t3dg_2, const Dg< double, 3, 3 > &t3dg_3)
Definition: test_T3dg_37.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
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