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;
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;
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) +
")");
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.
void test_for_zero(const T &t, const std::string &s)