9{
13
15 t2_a(
j,
i) = t2_4(
i,
j);
16
18 t2(
j,
k) = t2_a(
j,
i) * t2_5(
k,
i);
19
20 for(int jj = 0; jj < 3; ++jj)
21 for(int kk = 0; kk < 3; ++kk)
22 {
24 t2(jj, kk)
25 - (t2_a(jj, 0) * t2_5(kk, 0) + t2_a(jj, 1) * t2_5(kk, 1)
26 + t2_a(jj, 2) * t2_5(kk, 2) + t2_a(jj, 3) * t2_5(kk, 3)),
27 "T2(i,j)*T2(k,j)(" + std::to_string(jj) + "," + std::to_string(kk)
28 + ")");
29 }
30}
FTensor::Index< 'i', SPACE_DIM > i
FTensor::Index< 'j', 3 > j
FTensor::Index< 'k', 3 > k
void test_for_zero(const T &t, const std::string &s)