v0.14.0
Loading...
Searching...
No Matches
test_T2s_43.cpp
Go to the documentation of this file.
3#include <iostream>
4using namespace FTensor;
5using namespace std;
6
10{
11 Index<'i', 3> i;
12 Index<'j', 3> j;
13 Index<'k', 3> k;
14
18
19 /* Tensor2_symmetric tests */
20
21 t2_1(0, 0) = t2s_2(0, 0);
22 t2_1(0, 1) = t2s_2(0, 1);
23 t2_1(0, 2) = t2s_2(0, 2);
24 t2_1(1, 0) = t2s_2(1, 0);
25 t2_1(1, 1) = t2s_2(1, 1);
26 t2_1(1, 2) = t2s_2(1, 2);
27 t2_1(2, 0) = t2s_2(2, 0);
28 t2_1(2, 1) = t2s_2(2, 1);
29 t2_1(2, 2) = t2s_2(2, 2);
30
31 t2s_1(i, j) = t2s_2(k, i) ^ t2_1(k, j);
32 test_for_zero(t2s_1(0, 0)
33 - (t2s_2(0, 0) * t2_1(0, 0) + t2s_2(1, 0) * t2_1(1, 0)
34 + t2s_2(2, 0) * t2_1(2, 0)),
35 "T2s(k,i)^T2(k,j)(0,0)");
36 test_for_zero(t2s_1(0, 1)
37 - (t2s_2(0, 0) * t2_1(0, 1) + t2s_2(1, 0) * t2_1(1, 1)
38 + t2s_2(2, 0) * t2_1(2, 1)),
39 "T2s(k,i)^T2(k,j)(0,1)");
40 test_for_zero(t2s_1(0, 2)
41 - (t2s_2(0, 0) * t2_1(0, 2) + t2s_2(1, 0) * t2_1(1, 2)
42 + t2s_2(2, 0) * t2_1(2, 2)),
43 "T2s(k,i)^T2(k,j)(0,2)");
44 test_for_zero(t2s_1(1, 0)
45 - (t2s_2(0, 1) * t2_1(0, 0) + t2s_2(1, 1) * t2_1(1, 0)
46 + t2s_2(2, 1) * t2_1(2, 0)),
47 "T2s(k,i)^T2(k,j)(1,0)");
48 test_for_zero(t2s_1(1, 1)
49 - (t2s_2(0, 1) * t2_1(0, 1) + t2s_2(1, 1) * t2_1(1, 1)
50 + t2s_2(2, 1) * t2_1(2, 1)),
51 "T2s(k,i)^T2(k,j)(1,1)");
52 test_for_zero(t2s_1(1, 2)
53 - (t2s_2(0, 1) * t2_1(0, 2) + t2s_2(1, 1) * t2_1(1, 2)
54 + t2s_2(2, 1) * t2_1(2, 2)),
55 "T2s(k,i)^T2(k,j)(1,2)");
56 test_for_zero(t2s_1(2, 0)
57 - (t2s_2(0, 2) * t2_1(0, 0) + t2s_2(1, 2) * t2_1(1, 0)
58 + t2s_2(2, 2) * t2_1(2, 0)),
59 "T2s(k,i)^T2(k,j)(2,0)");
60 test_for_zero(t2s_1(2, 1)
61 - (t2s_2(0, 2) * t2_1(0, 1) + t2s_2(1, 2) * t2_1(1, 1)
62 + t2s_2(2, 2) * t2_1(2, 1)),
63 "T2s(k,i)^T2(k,j)(2,1)");
64 test_for_zero(t2s_1(2, 2)
65 - (t2s_2(0, 2) * t2_1(0, 2) + t2s_2(1, 2) * t2_1(1, 2)
66 + t2s_2(2, 2) * t2_1(2, 2)),
67 "T2s(k,i)^T2(k,j)(2,2)");
68 t2s_1(i, j) = t2_1(k, j) ^ t2s_2(k, i);
69 test_for_zero(t2s_1(0, 0)
70 - (t2s_2(0, 0) * t2_1(0, 0) + t2s_2(1, 0) * t2_1(1, 0)
71 + t2s_2(2, 0) * t2_1(2, 0)),
72 "T2(k,j)^T2s(k,i)(0,0)");
73 test_for_zero(t2s_1(0, 1)
74 - (t2s_2(0, 0) * t2_1(0, 1) + t2s_2(1, 0) * t2_1(1, 1)
75 + t2s_2(2, 0) * t2_1(2, 1)),
76 "T2(k,j)^T2s(k,i)(0,1)");
77 test_for_zero(t2s_1(0, 2)
78 - (t2s_2(0, 0) * t2_1(0, 2) + t2s_2(1, 0) * t2_1(1, 2)
79 + t2s_2(2, 0) * t2_1(2, 2)),
80 "T2(k,j)^T2s(k,i)(0,2)");
81 test_for_zero(t2s_1(1, 0)
82 - (t2s_2(0, 1) * t2_1(0, 0) + t2s_2(1, 1) * t2_1(1, 0)
83 + t2s_2(2, 1) * t2_1(2, 0)),
84 "T2(k,j)^T2s(k,i)(1,0)");
85 test_for_zero(t2s_1(1, 1)
86 - (t2s_2(0, 1) * t2_1(0, 1) + t2s_2(1, 1) * t2_1(1, 1)
87 + t2s_2(2, 1) * t2_1(2, 1)),
88 "T2(k,j)^T2s(k,i)(1,1)");
89 test_for_zero(t2s_1(1, 2)
90 - (t2s_2(0, 1) * t2_1(0, 2) + t2s_2(1, 1) * t2_1(1, 2)
91 + t2s_2(2, 1) * t2_1(2, 2)),
92 "T2(k,j)^T2s(k,i)(1,2)");
93 test_for_zero(t2s_1(2, 0)
94 - (t2s_2(0, 2) * t2_1(0, 0) + t2s_2(1, 2) * t2_1(1, 0)
95 + t2s_2(2, 2) * t2_1(2, 0)),
96 "T2(k,j)^T2s(k,i)(2,0)");
97 test_for_zero(t2s_1(2, 1)
98 - (t2s_2(0, 2) * t2_1(0, 1) + t2s_2(1, 2) * t2_1(1, 1)
99 + t2s_2(2, 2) * t2_1(2, 1)),
100 "T2(k,j)^T2s(k,i)(2,1)");
101 test_for_zero(t2s_1(2, 2)
102 - (t2s_2(0, 2) * t2_1(0, 2) + t2s_2(1, 2) * t2_1(1, 2)
103 + t2s_2(2, 2) * t2_1(2, 2)),
104 "T2(k,j)^T2s(k,i)(2,2)");
105}
static Number< 2 > N2
static Number< 1 > N1
static Number< 0 > N0
Tensors class implemented by Walter Landry.
FTensor::Index< 'i', SPACE_DIM > i
FTensor::Index< 'j', 3 > j
FTensor::Index< 'k', 3 > k
Tensors class implemented by Walter Landry.
Definition FTensor.hpp:51
void test_T2s_43(Tensor2< double, 3, 3 > &t2_1, Tensor2_symmetric< double, 3 > &t2s_1, const Tensor2_symmetric< double, 3 > &t2s_2)
void test_for_zero(const T &t, const std::string &s)