v0.14.0
Loading...
Searching...
No Matches
test_T4_008.cpp
Go to the documentation of this file.
3#include <iostream>
4using namespace FTensor;
5using namespace std;
6
7void test_T4_008() {
8
9 Index<'i', 3> i;
10 Index<'j', 3> j;
11 Index<'k', 3> k;
12 Index<'l', 3> l;
13
15
16 for (int ii = 0; ii != 3;++ii)
17 for (int jj = 0; jj != 3;++jj)
18 for (int kk = 0; kk != 3;++kk)
19 for (int ll = 0; ll != 3;++ll)
20 t_4(ii, jj, kk, ll) = 1 + ii + 10 * jj + 100 * kk + 1000 * ll;
21
24 for (int ii = 0; ii != 3;++ii)
25 for (int jj = ii; jj != 3;++jj) {
26 ts_2(ii, jj) = 1 + ii + 10 * jj;
27 t_2(ii, jj) = 1 + ii + 10 * jj;
28 t_2(jj, ii) = 1 + ii + 10 * jj;
29 }
30
32 t_2_2(k, l) = t_4(i, j, k, l) * ts_2(i, j);
33 t_2_2(k, l) -= t_4(i, j, k, l) * t_2(i, j);
34 for (int ii = 0; ii != 3; ++ii)
35 for (int jj = 0; jj != 3; ++jj) {
36 test_for_zero(t_2_2(ii, jj), "T4(i,j,k,l)*Ts2(i,j)(" + to_string(ii) +
37 "," + to_string(jj) + ")");
38 }
39
40
41}
Tensors class implemented by Walter Landry.
FTensor::Index< 'i', SPACE_DIM > i
FTensor::Index< 'l', 3 > l
FTensor::Index< 'j', 3 > j
FTensor::Index< 'k', 3 > k
Tensors class implemented by Walter Landry.
Definition FTensor.hpp:51
void test_T4_008()
void test_for_zero(const T &t, const std::string &s)