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