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