10{
   14 
   18 
   19  t3as_1(
i, 
j, 
k) = t3as_2(
i, 
j, 
k) + t3as_3(
i, 
k, 
j);
 
   20  test_for_zero(t3as_1(0, 0, 0) - (t3as_2(0, 0, 0) - t3as_3(0, 0, 0)),
 
   21                "T3as(i,j,k) + T3as(i,k,j)(0,0,0)");
   22  test_for_zero(t3as_1(0, 0, 1) - (t3as_2(0, 0, 1) - t3as_3(0, 0, 1)),
 
   23                "T3as(i,j,k) + T3as(i,k,j)(0,0,1)");
   24  test_for_zero(t3as_1(0, 0, 2) - (t3as_2(0, 0, 2) - t3as_3(0, 0, 2)),
 
   25                "T3as(i,j,k) + T3as(i,k,j)(0,0,2)");
   26  test_for_zero(t3as_1(0, 1, 0) - (t3as_2(0, 1, 0) - t3as_3(0, 1, 0)),
 
   27                "T3as(i,j,k) + T3as(i,k,j)(0,1,0)");
   28  test_for_zero(t3as_1(0, 1, 1) - (t3as_2(0, 1, 1) - t3as_3(0, 1, 1)),
 
   29                "T3as(i,j,k) + T3as(i,k,j)(0,1,1)");
   30  test_for_zero(t3as_1(0, 1, 2) - (t3as_2(0, 1, 2) - t3as_3(0, 1, 2)),
 
   31                "T3as(i,j,k) + T3as(i,k,j)(0,1,2)");
   32  test_for_zero(t3as_1(0, 2, 0) - (t3as_2(0, 2, 0) - t3as_3(0, 2, 0)),
 
   33                "T3as(i,j,k) + T3as(i,k,j)(0,2,0)");
   34  test_for_zero(t3as_1(0, 2, 1) - (t3as_2(0, 2, 1) - t3as_3(0, 2, 1)),
 
   35                "T3as(i,j,k) + T3as(i,k,j)(0,2,1)");
   36  test_for_zero(t3as_1(0, 2, 2) - (t3as_2(0, 2, 2) - t3as_3(0, 2, 2)),
 
   37                "T3as(i,j,k) + T3as(i,k,j)(0,2,2)");
   38  test_for_zero(t3as_1(1, 0, 0) - (t3as_2(1, 0, 0) - t3as_3(1, 0, 0)),
 
   39                "T3as(i,j,k) + T3as(i,k,j)(1,0,0)");
   40  test_for_zero(t3as_1(1, 0, 1) - (t3as_2(1, 0, 1) - t3as_3(1, 0, 1)),
 
   41                "T3as(i,j,k) + T3as(i,k,j)(1,0,1)");
   42  test_for_zero(t3as_1(1, 0, 2) - (t3as_2(1, 0, 2) - t3as_3(1, 0, 2)),
 
   43                "T3as(i,j,k) + T3as(i,k,j)(1,0,2)");
   44  test_for_zero(t3as_1(1, 1, 0) - (t3as_2(1, 1, 0) - t3as_3(1, 1, 0)),
 
   45                "T3as(i,j,k) + T3as(i,k,j)(1,1,0)");
   46  test_for_zero(t3as_1(1, 1, 1) - (t3as_2(1, 1, 1) - t3as_3(1, 1, 1)),
 
   47                "T3as(i,j,k) + T3as(i,k,j)(1,1,1)");
   48  test_for_zero(t3as_1(1, 1, 2) - (t3as_2(1, 1, 2) - t3as_3(1, 1, 2)),
 
   49                "T3as(i,j,k) + T3as(i,k,j)(1,1,2)");
   50  test_for_zero(t3as_1(1, 2, 0) - (t3as_2(1, 2, 0) - t3as_3(1, 2, 0)),
 
   51                "T3as(i,j,k) + T3as(i,k,j)(1,2,0)");
   52  test_for_zero(t3as_1(1, 2, 1) - (t3as_2(1, 2, 1) - t3as_3(1, 2, 1)),
 
   53                "T3as(i,j,k) + T3as(i,k,j)(1,2,1)");
   54  test_for_zero(t3as_1(1, 2, 2) - (t3as_2(1, 2, 2) - t3as_3(1, 2, 2)),
 
   55                "T3as(i,j,k) + T3as(i,k,j)(1,2,2)");
   56  test_for_zero(t3as_1(2, 0, 0) - (t3as_2(2, 0, 0) - t3as_3(2, 0, 0)),
 
   57                "T3as(i,j,k) + T3as(i,k,j)(2,0,0)");
   58  test_for_zero(t3as_1(2, 0, 1) - (t3as_2(2, 0, 1) - t3as_3(2, 0, 1)),
 
   59                "T3as(i,j,k) + T3as(i,k,j)(2,0,1)");
   60  test_for_zero(t3as_1(2, 0, 2) - (t3as_2(2, 0, 2) - t3as_3(2, 0, 2)),
 
   61                "T3as(i,j,k) + T3as(i,k,j)(2,0,2)");
   62  test_for_zero(t3as_1(2, 1, 0) - (t3as_2(2, 1, 0) - t3as_3(2, 1, 0)),
 
   63                "T3as(i,j,k) + T3as(i,k,j)(2,1,0)");
   64  test_for_zero(t3as_1(2, 1, 1) - (t3as_2(2, 1, 1) - t3as_3(2, 1, 1)),
 
   65                "T3as(i,j,k) + T3as(i,k,j)(2,1,1)");
   66  test_for_zero(t3as_1(2, 1, 2) - (t3as_2(2, 1, 2) - t3as_3(2, 1, 2)),
 
   67                "T3as(i,j,k) + T3as(i,k,j)(2,1,2)");
   68  test_for_zero(t3as_1(2, 2, 0) - (t3as_2(2, 2, 0) - t3as_3(2, 2, 0)),
 
   69                "T3as(i,j,k) + T3as(i,k,j)(2,2,0)");
   70  test_for_zero(t3as_1(2, 2, 1) - (t3as_2(2, 2, 1) - t3as_3(2, 2, 1)),
 
   71                "T3as(i,j,k) + T3as(i,k,j)(2,2,1)");
   72  test_for_zero(t3as_1(2, 2, 2) - (t3as_2(2, 2, 2) - t3as_3(2, 2, 2)),
 
   73                "T3as(i,j,k) + T3as(i,k,j)(2,2,2)");
   74}
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)