v0.15.0
Loading...
Searching...
No Matches
test_T3dg_29.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_T3dg_29(const Tensor2_symmetric<double, 3> &t2s_2,
8 Dg<double, 3, 3> &t3dg_1, const Dg<double, 3, 3> &t3dg_2,
9 const Dg<double, 3, 3> &t3dg_3)
10{
11 Index<'i', 3> i;
12 Index<'j', 3> j;
13 Index<'k', 3> k;
14 Index<'l', 3> l;
15
16 Number<0> N0;
17 Number<1> N1;
18 Number<2> N2;
19
20 /* Dg tests */
21
22 t3dg_1(i, j, k) = t3dg_2(i, j, l) * t2s_2(k, l);
23 test_for_zero(t3dg_1(0, 0, 0)
24 - (t3dg_2(0, 0, 0) * t2s_2(0, 0)
25 + t3dg_2(0, 0, 1) * t2s_2(0, 1)
26 + t3dg_2(0, 0, 2) * t2s_2(0, 2)),
27 "T3dg(i,j,l)*T2s(k,l)(0,0,0)");
28 test_for_zero(t3dg_1(0, 0, 1)
29 - (t3dg_2(0, 0, 0) * t2s_2(1, 0)
30 + t3dg_2(0, 0, 1) * t2s_2(1, 1)
31 + t3dg_2(0, 0, 2) * t2s_2(1, 2)),
32 "T3dg(i,j,l)*T2s(k,l)(0,0,1)");
33 test_for_zero(t3dg_1(0, 0, 2)
34 - (t3dg_2(0, 0, 0) * t2s_2(2, 0)
35 + t3dg_2(0, 0, 1) * t2s_2(2, 1)
36 + t3dg_2(0, 0, 2) * t2s_2(2, 2)),
37 "T3dg(i,j,l)*T2s(k,l)(0,0,2)");
38 test_for_zero(t3dg_1(0, 1, 0)
39 - (t3dg_2(0, 1, 0) * t2s_2(0, 0)
40 + t3dg_2(0, 1, 1) * t2s_2(0, 1)
41 + t3dg_2(0, 1, 2) * t2s_2(0, 2)),
42 "T3dg(i,j,l)*T2s(k,l)(0,1,0)");
43 test_for_zero(t3dg_1(0, 1, 1)
44 - (t3dg_2(0, 1, 0) * t2s_2(1, 0)
45 + t3dg_2(0, 1, 1) * t2s_2(1, 1)
46 + t3dg_2(0, 1, 2) * t2s_2(1, 2)),
47 "T3dg(i,j,l)*T2s(k,l)(0,1,1)");
48 test_for_zero(t3dg_1(0, 1, 2)
49 - (t3dg_2(0, 1, 0) * t2s_2(2, 0)
50 + t3dg_2(0, 1, 1) * t2s_2(2, 1)
51 + t3dg_2(0, 1, 2) * t2s_2(2, 2)),
52 "T3dg(i,j,l)*T2s(k,l)(0,1,2)");
53 test_for_zero(t3dg_1(0, 2, 0)
54 - (t3dg_2(0, 2, 0) * t2s_2(0, 0)
55 + t3dg_2(0, 2, 1) * t2s_2(0, 1)
56 + t3dg_2(0, 2, 2) * t2s_2(0, 2)),
57 "T3dg(i,j,l)*T2s(k,l)(0,2,0)");
58 test_for_zero(t3dg_1(0, 2, 1)
59 - (t3dg_2(0, 2, 0) * t2s_2(1, 0)
60 + t3dg_2(0, 2, 1) * t2s_2(1, 1)
61 + t3dg_2(0, 2, 2) * t2s_2(1, 2)),
62 "T3dg(i,j,l)*T2s(k,l)(0,2,1)");
63 test_for_zero(t3dg_1(0, 2, 2)
64 - (t3dg_2(0, 2, 0) * t2s_2(2, 0)
65 + t3dg_2(0, 2, 1) * t2s_2(2, 1)
66 + t3dg_2(0, 2, 2) * t2s_2(2, 2)),
67 "T3dg(i,j,l)*T2s(k,l)(0,2,2)");
68 test_for_zero(t3dg_1(1, 0, 0)
69 - (t3dg_2(1, 0, 0) * t2s_2(0, 0)
70 + t3dg_2(1, 0, 1) * t2s_2(0, 1)
71 + t3dg_2(1, 0, 2) * t2s_2(0, 2)),
72 "T3dg(i,j,l)*T2s(k,l)(1,0,0)");
73 test_for_zero(t3dg_1(1, 0, 1)
74 - (t3dg_2(1, 0, 0) * t2s_2(1, 0)
75 + t3dg_2(1, 0, 1) * t2s_2(1, 1)
76 + t3dg_2(1, 0, 2) * t2s_2(1, 2)),
77 "T3dg(i,j,l)*T2s(k,l)(1,0,1)");
78 test_for_zero(t3dg_1(1, 0, 2)
79 - (t3dg_2(1, 0, 0) * t2s_2(2, 0)
80 + t3dg_2(1, 0, 1) * t2s_2(2, 1)
81 + t3dg_2(1, 0, 2) * t2s_2(2, 2)),
82 "T3dg(i,j,l)*T2s(k,l)(1,0,2)");
83 test_for_zero(t3dg_1(1, 1, 0)
84 - (t3dg_2(1, 1, 0) * t2s_2(0, 0)
85 + t3dg_2(1, 1, 1) * t2s_2(0, 1)
86 + t3dg_2(1, 1, 2) * t2s_2(0, 2)),
87 "T3dg(i,j,l)*T2s(k,l)(1,1,0)");
88 test_for_zero(t3dg_1(1, 1, 1)
89 - (t3dg_2(1, 1, 0) * t2s_2(1, 0)
90 + t3dg_2(1, 1, 1) * t2s_2(1, 1)
91 + t3dg_2(1, 1, 2) * t2s_2(1, 2)),
92 "T3dg(i,j,l)*T2s(k,l)(1,1,1)");
93 test_for_zero(t3dg_1(1, 1, 2)
94 - (t3dg_2(1, 1, 0) * t2s_2(2, 0)
95 + t3dg_2(1, 1, 1) * t2s_2(2, 1)
96 + t3dg_2(1, 1, 2) * t2s_2(2, 2)),
97 "T3dg(i,j,l)*T2s(k,l)(1,1,2)");
98 test_for_zero(t3dg_1(1, 2, 0)
99 - (t3dg_2(1, 2, 0) * t2s_2(0, 0)
100 + t3dg_2(1, 2, 1) * t2s_2(0, 1)
101 + t3dg_2(1, 2, 2) * t2s_2(0, 2)),
102 "T3dg(i,j,l)*T2s(k,l)(1,2,0)");
103 test_for_zero(t3dg_1(1, 2, 1)
104 - (t3dg_2(1, 2, 0) * t2s_2(1, 0)
105 + t3dg_2(1, 2, 1) * t2s_2(1, 1)
106 + t3dg_2(1, 2, 2) * t2s_2(1, 2)),
107 "T3dg(i,j,l)*T2s(k,l)(1,2,1)");
108 test_for_zero(t3dg_1(1, 2, 2)
109 - (t3dg_2(1, 2, 0) * t2s_2(2, 0)
110 + t3dg_2(1, 2, 1) * t2s_2(2, 1)
111 + t3dg_2(1, 2, 2) * t2s_2(2, 2)),
112 "T3dg(i,j,l)*T2s(k,l)(1,2,2)");
113 test_for_zero(t3dg_1(2, 0, 0)
114 - (t3dg_2(2, 0, 0) * t2s_2(0, 0)
115 + t3dg_2(2, 0, 1) * t2s_2(0, 1)
116 + t3dg_2(2, 0, 2) * t2s_2(0, 2)),
117 "T3dg(i,j,l)*T2s(k,l)(2,0,0)");
118 test_for_zero(t3dg_1(2, 0, 1)
119 - (t3dg_2(2, 0, 0) * t2s_2(1, 0)
120 + t3dg_2(2, 0, 1) * t2s_2(1, 1)
121 + t3dg_2(2, 0, 2) * t2s_2(1, 2)),
122 "T3dg(i,j,l)*T2s(k,l)(2,0,1)");
123 test_for_zero(t3dg_1(2, 0, 2)
124 - (t3dg_2(2, 0, 0) * t2s_2(2, 0)
125 + t3dg_2(2, 0, 1) * t2s_2(2, 1)
126 + t3dg_2(2, 0, 2) * t2s_2(2, 2)),
127 "T3dg(i,j,l)*T2s(k,l)(2,0,2)");
128 test_for_zero(t3dg_1(2, 1, 0)
129 - (t3dg_2(2, 1, 0) * t2s_2(0, 0)
130 + t3dg_2(2, 1, 1) * t2s_2(0, 1)
131 + t3dg_2(2, 1, 2) * t2s_2(0, 2)),
132 "T3dg(i,j,l)*T2s(k,l)(2,1,0)");
133 test_for_zero(t3dg_1(2, 1, 1)
134 - (t3dg_2(2, 1, 0) * t2s_2(1, 0)
135 + t3dg_2(2, 1, 1) * t2s_2(1, 1)
136 + t3dg_2(2, 1, 2) * t2s_2(1, 2)),
137 "T3dg(i,j,l)*T2s(k,l)(2,1,1)");
138 test_for_zero(t3dg_1(2, 1, 2)
139 - (t3dg_2(2, 1, 0) * t2s_2(2, 0)
140 + t3dg_2(2, 1, 1) * t2s_2(2, 1)
141 + t3dg_2(2, 1, 2) * t2s_2(2, 2)),
142 "T3dg(i,j,l)*T2s(k,l)(2,1,2)");
143 test_for_zero(t3dg_1(2, 2, 0)
144 - (t3dg_2(2, 2, 0) * t2s_2(0, 0)
145 + t3dg_2(2, 2, 1) * t2s_2(0, 1)
146 + t3dg_2(2, 2, 2) * t2s_2(0, 2)),
147 "T3dg(i,j,l)*T2s(k,l)(2,2,0)");
148 test_for_zero(t3dg_1(2, 2, 1)
149 - (t3dg_2(2, 2, 0) * t2s_2(1, 0)
150 + t3dg_2(2, 2, 1) * t2s_2(1, 1)
151 + t3dg_2(2, 2, 2) * t2s_2(1, 2)),
152 "T3dg(i,j,l)*T2s(k,l)(2,2,1)");
153 test_for_zero(t3dg_1(2, 2, 2)
154 - (t3dg_2(2, 2, 0) * t2s_2(2, 0)
155 + t3dg_2(2, 2, 1) * t2s_2(2, 1)
156 + t3dg_2(2, 2, 2) * t2s_2(2, 2)),
157 "T3dg(i,j,l)*T2s(k,l)(2,2,2)");
158
159 t3dg_1(i, j, k) = t2s_2(k, l) * t3dg_3(i, j, l);
160 test_for_zero(t3dg_1(0, 0, 0)
161 - (t3dg_3(0, 0, 0) * t2s_2(0, 0)
162 + t3dg_3(0, 0, 1) * t2s_2(0, 1)
163 + t3dg_3(0, 0, 2) * t2s_2(0, 2)),
164 "T2s(k,l)*T3dg(i,j,l)(0,0,0)");
165 test_for_zero(t3dg_1(0, 0, 1)
166 - (t3dg_3(0, 0, 0) * t2s_2(1, 0)
167 + t3dg_3(0, 0, 1) * t2s_2(1, 1)
168 + t3dg_3(0, 0, 2) * t2s_2(1, 2)),
169 "T2s(k,l)*T3dg(i,j,l)(0,0,1)");
170 test_for_zero(t3dg_1(0, 0, 2)
171 - (t3dg_3(0, 0, 0) * t2s_2(2, 0)
172 + t3dg_3(0, 0, 1) * t2s_2(2, 1)
173 + t3dg_3(0, 0, 2) * t2s_2(2, 2)),
174 "T2s(k,l)*T3dg(i,j,l)(0,0,2)");
175 test_for_zero(t3dg_1(0, 1, 0)
176 - (t3dg_3(0, 1, 0) * t2s_2(0, 0)
177 + t3dg_3(0, 1, 1) * t2s_2(0, 1)
178 + t3dg_3(0, 1, 2) * t2s_2(0, 2)),
179 "T2s(k,l)*T3dg(i,j,l)(0,1,0)");
180 test_for_zero(t3dg_1(0, 1, 1)
181 - (t3dg_3(0, 1, 0) * t2s_2(1, 0)
182 + t3dg_3(0, 1, 1) * t2s_2(1, 1)
183 + t3dg_3(0, 1, 2) * t2s_2(1, 2)),
184 "T2s(k,l)*T3dg(i,j,l)(0,1,1)");
185 test_for_zero(t3dg_1(0, 1, 2)
186 - (t3dg_3(0, 1, 0) * t2s_2(2, 0)
187 + t3dg_3(0, 1, 1) * t2s_2(2, 1)
188 + t3dg_3(0, 1, 2) * t2s_2(2, 2)),
189 "T2s(k,l)*T3dg(i,j,l)(0,1,2)");
190 test_for_zero(t3dg_1(0, 2, 0)
191 - (t3dg_3(0, 2, 0) * t2s_2(0, 0)
192 + t3dg_3(0, 2, 1) * t2s_2(0, 1)
193 + t3dg_3(0, 2, 2) * t2s_2(0, 2)),
194 "T2s(k,l)*T3dg(i,j,l)(0,2,0)");
195 test_for_zero(t3dg_1(0, 2, 1)
196 - (t3dg_3(0, 2, 0) * t2s_2(1, 0)
197 + t3dg_3(0, 2, 1) * t2s_2(1, 1)
198 + t3dg_3(0, 2, 2) * t2s_2(1, 2)),
199 "T2s(k,l)*T3dg(i,j,l)(0,2,1)");
200 test_for_zero(t3dg_1(0, 2, 2)
201 - (t3dg_3(0, 2, 0) * t2s_2(2, 0)
202 + t3dg_3(0, 2, 1) * t2s_2(2, 1)
203 + t3dg_3(0, 2, 2) * t2s_2(2, 2)),
204 "T2s(k,l)*T3dg(i,j,l)(0,2,2)");
205 test_for_zero(t3dg_1(1, 0, 0)
206 - (t3dg_3(1, 0, 0) * t2s_2(0, 0)
207 + t3dg_3(1, 0, 1) * t2s_2(0, 1)
208 + t3dg_3(1, 0, 2) * t2s_2(0, 2)),
209 "T2s(k,l)*T3dg(i,j,l)(1,0,0)");
210 test_for_zero(t3dg_1(1, 0, 1)
211 - (t3dg_3(1, 0, 0) * t2s_2(1, 0)
212 + t3dg_3(1, 0, 1) * t2s_2(1, 1)
213 + t3dg_3(1, 0, 2) * t2s_2(1, 2)),
214 "T2s(k,l)*T3dg(i,j,l)(1,0,1)");
215 test_for_zero(t3dg_1(1, 0, 2)
216 - (t3dg_3(1, 0, 0) * t2s_2(2, 0)
217 + t3dg_3(1, 0, 1) * t2s_2(2, 1)
218 + t3dg_3(1, 0, 2) * t2s_2(2, 2)),
219 "T2s(k,l)*T3dg(i,j,l)(1,0,2)");
220 test_for_zero(t3dg_1(1, 1, 0)
221 - (t3dg_3(1, 1, 0) * t2s_2(0, 0)
222 + t3dg_3(1, 1, 1) * t2s_2(0, 1)
223 + t3dg_3(1, 1, 2) * t2s_2(0, 2)),
224 "T2s(k,l)*T3dg(i,j,l)(1,1,0)");
225 test_for_zero(t3dg_1(1, 1, 1)
226 - (t3dg_3(1, 1, 0) * t2s_2(1, 0)
227 + t3dg_3(1, 1, 1) * t2s_2(1, 1)
228 + t3dg_3(1, 1, 2) * t2s_2(1, 2)),
229 "T2s(k,l)*T3dg(i,j,l)(1,1,1)");
230 test_for_zero(t3dg_1(1, 1, 2)
231 - (t3dg_3(1, 1, 0) * t2s_2(2, 0)
232 + t3dg_3(1, 1, 1) * t2s_2(2, 1)
233 + t3dg_3(1, 1, 2) * t2s_2(2, 2)),
234 "T2s(k,l)*T3dg(i,j,l)(1,1,2)");
235 test_for_zero(t3dg_1(1, 2, 0)
236 - (t3dg_3(1, 2, 0) * t2s_2(0, 0)
237 + t3dg_3(1, 2, 1) * t2s_2(0, 1)
238 + t3dg_3(1, 2, 2) * t2s_2(0, 2)),
239 "T2s(k,l)*T3dg(i,j,l)(1,2,0)");
240 test_for_zero(t3dg_1(1, 2, 1)
241 - (t3dg_3(1, 2, 0) * t2s_2(1, 0)
242 + t3dg_3(1, 2, 1) * t2s_2(1, 1)
243 + t3dg_3(1, 2, 2) * t2s_2(1, 2)),
244 "T2s(k,l)*T3dg(i,j,l)(1,2,1)");
245 test_for_zero(t3dg_1(1, 2, 2)
246 - (t3dg_3(1, 2, 0) * t2s_2(2, 0)
247 + t3dg_3(1, 2, 1) * t2s_2(2, 1)
248 + t3dg_3(1, 2, 2) * t2s_2(2, 2)),
249 "T2s(k,l)*T3dg(i,j,l)(1,2,2)");
250 test_for_zero(t3dg_1(2, 0, 0)
251 - (t3dg_3(2, 0, 0) * t2s_2(0, 0)
252 + t3dg_3(2, 0, 1) * t2s_2(0, 1)
253 + t3dg_3(2, 0, 2) * t2s_2(0, 2)),
254 "T2s(k,l)*T3dg(i,j,l)(2,0,0)");
255 test_for_zero(t3dg_1(2, 0, 1)
256 - (t3dg_3(2, 0, 0) * t2s_2(1, 0)
257 + t3dg_3(2, 0, 1) * t2s_2(1, 1)
258 + t3dg_3(2, 0, 2) * t2s_2(1, 2)),
259 "T2s(k,l)*T3dg(i,j,l)(2,0,1)");
260 test_for_zero(t3dg_1(2, 0, 2)
261 - (t3dg_3(2, 0, 0) * t2s_2(2, 0)
262 + t3dg_3(2, 0, 1) * t2s_2(2, 1)
263 + t3dg_3(2, 0, 2) * t2s_2(2, 2)),
264 "T2s(k,l)*T3dg(i,j,l)(2,0,2)");
265 test_for_zero(t3dg_1(2, 1, 0)
266 - (t3dg_3(2, 1, 0) * t2s_2(0, 0)
267 + t3dg_3(2, 1, 1) * t2s_2(0, 1)
268 + t3dg_3(2, 1, 2) * t2s_2(0, 2)),
269 "T2s(k,l)*T3dg(i,j,l)(2,1,0)");
270 test_for_zero(t3dg_1(2, 1, 1)
271 - (t3dg_3(2, 1, 0) * t2s_2(1, 0)
272 + t3dg_3(2, 1, 1) * t2s_2(1, 1)
273 + t3dg_3(2, 1, 2) * t2s_2(1, 2)),
274 "T2s(k,l)*T3dg(i,j,l)(2,1,1)");
275 test_for_zero(t3dg_1(2, 1, 2)
276 - (t3dg_3(2, 1, 0) * t2s_2(2, 0)
277 + t3dg_3(2, 1, 1) * t2s_2(2, 1)
278 + t3dg_3(2, 1, 2) * t2s_2(2, 2)),
279 "T2s(k,l)*T3dg(i,j,l)(2,1,2)");
280 test_for_zero(t3dg_1(2, 2, 0)
281 - (t3dg_3(2, 2, 0) * t2s_2(0, 0)
282 + t3dg_3(2, 2, 1) * t2s_2(0, 1)
283 + t3dg_3(2, 2, 2) * t2s_2(0, 2)),
284 "T2s(k,l)*T3dg(i,j,l)(2,2,0)");
285 test_for_zero(t3dg_1(2, 2, 1)
286 - (t3dg_3(2, 2, 0) * t2s_2(1, 0)
287 + t3dg_3(2, 2, 1) * t2s_2(1, 1)
288 + t3dg_3(2, 2, 2) * t2s_2(1, 2)),
289 "T2s(k,l)*T3dg(i,j,l)(2,2,1)");
290 test_for_zero(t3dg_1(2, 2, 2)
291 - (t3dg_3(2, 2, 0) * t2s_2(2, 0)
292 + t3dg_3(2, 2, 1) * t2s_2(2, 1)
293 + t3dg_3(2, 2, 2) * t2s_2(2, 2)),
294 "T2s(k,l)*T3dg(i,j,l)(2,2,2)");
295}
FTensor::Tensor2_symmetric
Definition Tensor2_symmetric_value.hpp:14
i
FTensor::Index< 'i', SPACE_DIM > i
Definition hcurl_divergence_operator_2d.cpp:27
l
FTensor::Index< 'l', 3 > l
Definition matrix_function.cpp:21
j
FTensor::Index< 'j', 3 > j
Definition matrix_function.cpp:19
k
FTensor::Index< 'k', 3 > k
Definition matrix_function.cpp:20
FTensor
Tensors class implemented by Walter Landry.
Definition FTensor.hpp:51
test_T3dg_29
void test_T3dg_29(const Tensor2_symmetric< double, 3 > &t2s_2, Dg< double, 3, 3 > &t3dg_1, const Dg< double, 3, 3 > &t3dg_2, const Dg< double, 3, 3 > &t3dg_3)
Definition test_T3dg_29.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