9{
12
16
17
19
20
21
22 t2(
i,
j) = t2_2(
i,
j) - t1_1(
i) * t1_2(
j);
24 "T2(i,j)-T2(i,j)(0,0)");
26 "T2(i,j)-T2(i,j)(0,1)");
28 "T2(i,j)-T2(i,j)(0,2)");
30 "T2(i,j)-T2(i,j)(1,0)");
32 "T2(i,j)-T2(i,j)(1,1)");
34 "T2(i,j)-T2(i,j)(1,2)");
36 "T2(i,j)-T2(i,j)(2,0)");
38 "T2(i,j)-T2(i,j)(2,1)");
40 "T2(i,j)-T2(i,j)(2,2)");
41
42 t2(
i,
j) = t2_2(
i,
j) - t1_1(
j) * t1_2(
i);
44 "T2(i,j)-T2(j,i)(0,0)");
46 "T2(i,j)-T2(j,i)(0,1)");
48 "T2(i,j)-T2(j,i)(0,2)");
50 "T2(i,j)-T2(j,i)(1,0)");
52 "T2(i,j)-T2(j,i)(1,1)");
54 "T2(i,j)-T2(j,i)(1,2)");
56 "T2(i,j)-T2(j,i)(2,0)");
58 "T2(i,j)-T2(j,i)(2,1)");
60 "T2(i,j)-T2(j,i)(2,2)");
61}
FTensor::Index< 'i', SPACE_DIM > i
FTensor::Index< 'j', 3 > j
void test_for_zero(const T &t, const std::string &s)