10{
   13 
   17 
   18  
   19 
   20  
   21 
   22  
   23 
   24  t2s_1(
j, 
i) += t2s_2(
i, 
j);
 
   25  test_for_zero(t2s_1(0, 0) - (2 * t2s_2(0, 0) + t1_2(0) * t1_2(0)),
 
   26                "T2s(j,i)+=T2s(i,j)(0,0)");
   27  test_for_zero(t2s_1(1, 0) - (2 * t2s_2(0, 1) + t1_2(0) * t1_2(1)),
 
   28                "T2s(j,i)+=T2s(i,j)(0,1)");
   29  test_for_zero(t2s_1(2, 0) - (2 * t2s_2(0, 2) + t1_2(0) * t1_2(2)),
 
   30                "T2s(j,i)+=T2s(i,j)(0,2)");
   31  test_for_zero(t2s_1(0, 1) - (2 * t2s_2(1, 0) + t1_2(1) * t1_2(0)),
 
   32                "T2s(j,i)+=T2s(i,j)(1,0)");
   33  test_for_zero(t2s_1(1, 1) - (2 * t2s_2(1, 1) + t1_2(1) * t1_2(1)),
 
   34                "T2s(j,i)+=T2s(i,j)(1,1)");
   35  test_for_zero(t2s_1(2, 1) - (2 * t2s_2(1, 2) + t1_2(1) * t1_2(2)),
 
   36                "T2s(j,i)+=T2s(i,j)(1,2)");
   37  test_for_zero(t2s_1(0, 2) - (2 * t2s_2(2, 0) + t1_2(2) * t1_2(0)),
 
   38                "T2s(j,i)+=T2s(i,j)(2,0)");
   39  test_for_zero(t2s_1(1, 2) - (2 * t2s_2(2, 1) + t1_2(2) * t1_2(1)),
 
   40                "T2s(j,i)+=T2s(i,j)(2,1)");
   41  test_for_zero(t2s_1(2, 2) - (2 * t2s_2(2, 2) + t1_2(2) * t1_2(2)),
 
   42                "T2s(j,i)+=T2s(i,j)(2,2)");
   43}
FTensor::Index< 'i', SPACE_DIM > i
FTensor::Index< 'j', 3 > j
void test_for_zero(const T &t, const std::string &s)