9{
13
15 t2_a(
i,
j) = t2_5(
j,
i);
16
18 t2(
j,
k) = t2_4(
i,
j) * t2_a(
i,
k);
19
20 for(int jj = 0; jj < 3; ++jj)
21 for(int kk = 0; kk < 3; ++kk)
22 {
24 t2(jj, kk)
25 - (t2_4(0, jj) * t2_a(0, kk) + t2_4(1, jj) * t2_a(1, kk)
26 + t2_4(2, jj) * t2_a(2, kk) + t2_4(3, jj) * t2_a(3, kk)),
27 "T2(i,j)*T2(i,k)(" + 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)