8{
   15 
   19 
   21 
   22  t4ddg_2(
i, 
k, 
j, 
l) = t2s_2(
i, 
k) * t2s_2(
j, 
l);
 
   23  t4ddg_3(
i, 
l, 
j, 
k) = t2s_2(
i, 
l) * t2s_2(
j, 
k);
 
   24  t4ddg_1(
i, 
k, 
m, 
n) = (t4ddg_2(
i, 
j, 
k, 
l) ^ t4ddg_3(
j, 
l, 
m, 
n));
 
   26                  - (t4ddg_2(0, 0, 0, 0) * t4ddg_3(0, 0, 0, 0)
   27                     + t4ddg_2(0, 0, 0, 1) * t4ddg_3(0, 1, 0, 0)
   28                     + t4ddg_2(0, 0, 0, 2) * t4ddg_3(0, 2, 0, 0)
   29                     + t4ddg_2(0, 1, 0, 0) * t4ddg_3(1, 0, 0, 0)
   30                     + t4ddg_2(0, 1, 0, 1) * t4ddg_3(1, 1, 0, 0)
   31                     + t4ddg_2(0, 1, 0, 2) * t4ddg_3(1, 2, 0, 0)
   32                     + t4ddg_2(0, 2, 0, 0) * t4ddg_3(2, 0, 0, 0)
   33                     + t4ddg_2(0, 2, 0, 1) * t4ddg_3(2, 1, 0, 0)
   34                     + t4ddg_2(0, 2, 0, 2) * t4ddg_3(2, 2, 0, 0)),
   35                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,0,0,0)");
   37                  - (t4ddg_2(0, 0, 0, 0) * t4ddg_3(0, 0, 0, 1)
   38                     + t4ddg_2(0, 0, 0, 1) * t4ddg_3(0, 1, 0, 1)
   39                     + t4ddg_2(0, 0, 0, 2) * t4ddg_3(0, 2, 0, 1)
   40                     + t4ddg_2(0, 1, 0, 0) * t4ddg_3(1, 0, 0, 1)
   41                     + t4ddg_2(0, 1, 0, 1) * t4ddg_3(1, 1, 0, 1)
   42                     + t4ddg_2(0, 1, 0, 2) * t4ddg_3(1, 2, 0, 1)
   43                     + t4ddg_2(0, 2, 0, 0) * t4ddg_3(2, 0, 0, 1)
   44                     + t4ddg_2(0, 2, 0, 1) * t4ddg_3(2, 1, 0, 1)
   45                     + t4ddg_2(0, 2, 0, 2) * t4ddg_3(2, 2, 0, 1)),
   46                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,0,0,1)");
   48                  - (t4ddg_2(0, 0, 0, 0) * t4ddg_3(0, 0, 0, 2)
   49                     + t4ddg_2(0, 0, 0, 1) * t4ddg_3(0, 1, 0, 2)
   50                     + t4ddg_2(0, 0, 0, 2) * t4ddg_3(0, 2, 0, 2)
   51                     + t4ddg_2(0, 1, 0, 0) * t4ddg_3(1, 0, 0, 2)
   52                     + t4ddg_2(0, 1, 0, 1) * t4ddg_3(1, 1, 0, 2)
   53                     + t4ddg_2(0, 1, 0, 2) * t4ddg_3(1, 2, 0, 2)
   54                     + t4ddg_2(0, 2, 0, 0) * t4ddg_3(2, 0, 0, 2)
   55                     + t4ddg_2(0, 2, 0, 1) * t4ddg_3(2, 1, 0, 2)
   56                     + t4ddg_2(0, 2, 0, 2) * t4ddg_3(2, 2, 0, 2)),
   57                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,0,0,2)");
   59                  - (t4ddg_2(0, 0, 1, 0) * t4ddg_3(0, 0, 0, 0)
   60                     + t4ddg_2(0, 0, 1, 1) * t4ddg_3(0, 1, 0, 0)
   61                     + t4ddg_2(0, 0, 1, 2) * t4ddg_3(0, 2, 0, 0)
   62                     + t4ddg_2(0, 1, 1, 0) * t4ddg_3(1, 0, 0, 0)
   63                     + t4ddg_2(0, 1, 1, 1) * t4ddg_3(1, 1, 0, 0)
   64                     + t4ddg_2(0, 1, 1, 2) * t4ddg_3(1, 2, 0, 0)
   65                     + t4ddg_2(0, 2, 1, 0) * t4ddg_3(2, 0, 0, 0)
   66                     + t4ddg_2(0, 2, 1, 1) * t4ddg_3(2, 1, 0, 0)
   67                     + t4ddg_2(0, 2, 1, 2) * t4ddg_3(2, 2, 0, 0)),
   68                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,0,1,0)");
   70                  - (t4ddg_2(0, 0, 1, 0) * t4ddg_3(0, 0, 0, 1)
   71                     + t4ddg_2(0, 0, 1, 1) * t4ddg_3(0, 1, 0, 1)
   72                     + t4ddg_2(0, 0, 1, 2) * t4ddg_3(0, 2, 0, 1)
   73                     + t4ddg_2(0, 1, 1, 0) * t4ddg_3(1, 0, 0, 1)
   74                     + t4ddg_2(0, 1, 1, 1) * t4ddg_3(1, 1, 0, 1)
   75                     + t4ddg_2(0, 1, 1, 2) * t4ddg_3(1, 2, 0, 1)
   76                     + t4ddg_2(0, 2, 1, 0) * t4ddg_3(2, 0, 0, 1)
   77                     + t4ddg_2(0, 2, 1, 1) * t4ddg_3(2, 1, 0, 1)
   78                     + t4ddg_2(0, 2, 1, 2) * t4ddg_3(2, 2, 0, 1)),
   79                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,0,1,1)");
   81                  - (t4ddg_2(0, 0, 1, 0) * t4ddg_3(0, 0, 0, 2)
   82                     + t4ddg_2(0, 0, 1, 1) * t4ddg_3(0, 1, 0, 2)
   83                     + t4ddg_2(0, 0, 1, 2) * t4ddg_3(0, 2, 0, 2)
   84                     + t4ddg_2(0, 1, 1, 0) * t4ddg_3(1, 0, 0, 2)
   85                     + t4ddg_2(0, 1, 1, 1) * t4ddg_3(1, 1, 0, 2)
   86                     + t4ddg_2(0, 1, 1, 2) * t4ddg_3(1, 2, 0, 2)
   87                     + t4ddg_2(0, 2, 1, 0) * t4ddg_3(2, 0, 0, 2)
   88                     + t4ddg_2(0, 2, 1, 1) * t4ddg_3(2, 1, 0, 2)
   89                     + t4ddg_2(0, 2, 1, 2) * t4ddg_3(2, 2, 0, 2)),
   90                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,0,1,2)");
   92                  - (t4ddg_2(0, 0, 2, 0) * t4ddg_3(0, 0, 0, 0)
   93                     + t4ddg_2(0, 0, 2, 1) * t4ddg_3(0, 1, 0, 0)
   94                     + t4ddg_2(0, 0, 2, 2) * t4ddg_3(0, 2, 0, 0)
   95                     + t4ddg_2(0, 1, 2, 0) * t4ddg_3(1, 0, 0, 0)
   96                     + t4ddg_2(0, 1, 2, 1) * t4ddg_3(1, 1, 0, 0)
   97                     + t4ddg_2(0, 1, 2, 2) * t4ddg_3(1, 2, 0, 0)
   98                     + t4ddg_2(0, 2, 2, 0) * t4ddg_3(2, 0, 0, 0)
   99                     + t4ddg_2(0, 2, 2, 1) * t4ddg_3(2, 1, 0, 0)
  100                     + t4ddg_2(0, 2, 2, 2) * t4ddg_3(2, 2, 0, 0)),
  101                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,0,2,0)");
  103                  - (t4ddg_2(0, 0, 2, 0) * t4ddg_3(0, 0, 0, 1)
  104                     + t4ddg_2(0, 0, 2, 1) * t4ddg_3(0, 1, 0, 1)
  105                     + t4ddg_2(0, 0, 2, 2) * t4ddg_3(0, 2, 0, 1)
  106                     + t4ddg_2(0, 1, 2, 0) * t4ddg_3(1, 0, 0, 1)
  107                     + t4ddg_2(0, 1, 2, 1) * t4ddg_3(1, 1, 0, 1)
  108                     + t4ddg_2(0, 1, 2, 2) * t4ddg_3(1, 2, 0, 1)
  109                     + t4ddg_2(0, 2, 2, 0) * t4ddg_3(2, 0, 0, 1)
  110                     + t4ddg_2(0, 2, 2, 1) * t4ddg_3(2, 1, 0, 1)
  111                     + t4ddg_2(0, 2, 2, 2) * t4ddg_3(2, 2, 0, 1)),
  112                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,0,2,1)");
  114                  - (t4ddg_2(0, 0, 2, 0) * t4ddg_3(0, 0, 0, 2)
  115                     + t4ddg_2(0, 0, 2, 1) * t4ddg_3(0, 1, 0, 2)
  116                     + t4ddg_2(0, 0, 2, 2) * t4ddg_3(0, 2, 0, 2)
  117                     + t4ddg_2(0, 1, 2, 0) * t4ddg_3(1, 0, 0, 2)
  118                     + t4ddg_2(0, 1, 2, 1) * t4ddg_3(1, 1, 0, 2)
  119                     + t4ddg_2(0, 1, 2, 2) * t4ddg_3(1, 2, 0, 2)
  120                     + t4ddg_2(0, 2, 2, 0) * t4ddg_3(2, 0, 0, 2)
  121                     + t4ddg_2(0, 2, 2, 1) * t4ddg_3(2, 1, 0, 2)
  122                     + t4ddg_2(0, 2, 2, 2) * t4ddg_3(2, 2, 0, 2)),
  123                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,0,2,2)");
  125                  - (t4ddg_2(0, 0, 0, 0) * t4ddg_3(0, 0, 1, 0)
  126                     + t4ddg_2(0, 0, 0, 1) * t4ddg_3(0, 1, 1, 0)
  127                     + t4ddg_2(0, 0, 0, 2) * t4ddg_3(0, 2, 1, 0)
  128                     + t4ddg_2(0, 1, 0, 0) * t4ddg_3(1, 0, 1, 0)
  129                     + t4ddg_2(0, 1, 0, 1) * t4ddg_3(1, 1, 1, 0)
  130                     + t4ddg_2(0, 1, 0, 2) * t4ddg_3(1, 2, 1, 0)
  131                     + t4ddg_2(0, 2, 0, 0) * t4ddg_3(2, 0, 1, 0)
  132                     + t4ddg_2(0, 2, 0, 1) * t4ddg_3(2, 1, 1, 0)
  133                     + t4ddg_2(0, 2, 0, 2) * t4ddg_3(2, 2, 1, 0)),
  134                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,1,0,0)");
  136                  - (t4ddg_2(0, 0, 0, 0) * t4ddg_3(0, 0, 1, 1)
  137                     + t4ddg_2(0, 0, 0, 1) * t4ddg_3(0, 1, 1, 1)
  138                     + t4ddg_2(0, 0, 0, 2) * t4ddg_3(0, 2, 1, 1)
  139                     + t4ddg_2(0, 1, 0, 0) * t4ddg_3(1, 0, 1, 1)
  140                     + t4ddg_2(0, 1, 0, 1) * t4ddg_3(1, 1, 1, 1)
  141                     + t4ddg_2(0, 1, 0, 2) * t4ddg_3(1, 2, 1, 1)
  142                     + t4ddg_2(0, 2, 0, 0) * t4ddg_3(2, 0, 1, 1)
  143                     + t4ddg_2(0, 2, 0, 1) * t4ddg_3(2, 1, 1, 1)
  144                     + t4ddg_2(0, 2, 0, 2) * t4ddg_3(2, 2, 1, 1)),
  145                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,1,0,1)");
  147                  - (t4ddg_2(0, 0, 0, 0) * t4ddg_3(0, 0, 1, 2)
  148                     + t4ddg_2(0, 0, 0, 1) * t4ddg_3(0, 1, 1, 2)
  149                     + t4ddg_2(0, 0, 0, 2) * t4ddg_3(0, 2, 1, 2)
  150                     + t4ddg_2(0, 1, 0, 0) * t4ddg_3(1, 0, 1, 2)
  151                     + t4ddg_2(0, 1, 0, 1) * t4ddg_3(1, 1, 1, 2)
  152                     + t4ddg_2(0, 1, 0, 2) * t4ddg_3(1, 2, 1, 2)
  153                     + t4ddg_2(0, 2, 0, 0) * t4ddg_3(2, 0, 1, 2)
  154                     + t4ddg_2(0, 2, 0, 1) * t4ddg_3(2, 1, 1, 2)
  155                     + t4ddg_2(0, 2, 0, 2) * t4ddg_3(2, 2, 1, 2)),
  156                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,1,0,2)");
  158                  - (t4ddg_2(0, 0, 1, 0) * t4ddg_3(0, 0, 1, 0)
  159                     + t4ddg_2(0, 0, 1, 1) * t4ddg_3(0, 1, 1, 0)
  160                     + t4ddg_2(0, 0, 1, 2) * t4ddg_3(0, 2, 1, 0)
  161                     + t4ddg_2(0, 1, 1, 0) * t4ddg_3(1, 0, 1, 0)
  162                     + t4ddg_2(0, 1, 1, 1) * t4ddg_3(1, 1, 1, 0)
  163                     + t4ddg_2(0, 1, 1, 2) * t4ddg_3(1, 2, 1, 0)
  164                     + t4ddg_2(0, 2, 1, 0) * t4ddg_3(2, 0, 1, 0)
  165                     + t4ddg_2(0, 2, 1, 1) * t4ddg_3(2, 1, 1, 0)
  166                     + t4ddg_2(0, 2, 1, 2) * t4ddg_3(2, 2, 1, 0)),
  167                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,1,1,0)");
  169                  - (t4ddg_2(0, 0, 1, 0) * t4ddg_3(0, 0, 1, 1)
  170                     + t4ddg_2(0, 0, 1, 1) * t4ddg_3(0, 1, 1, 1)
  171                     + t4ddg_2(0, 0, 1, 2) * t4ddg_3(0, 2, 1, 1)
  172                     + t4ddg_2(0, 1, 1, 0) * t4ddg_3(1, 0, 1, 1)
  173                     + t4ddg_2(0, 1, 1, 1) * t4ddg_3(1, 1, 1, 1)
  174                     + t4ddg_2(0, 1, 1, 2) * t4ddg_3(1, 2, 1, 1)
  175                     + t4ddg_2(0, 2, 1, 0) * t4ddg_3(2, 0, 1, 1)
  176                     + t4ddg_2(0, 2, 1, 1) * t4ddg_3(2, 1, 1, 1)
  177                     + t4ddg_2(0, 2, 1, 2) * t4ddg_3(2, 2, 1, 1)),
  178                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,1,1,1)");
  180                  - (t4ddg_2(0, 0, 1, 0) * t4ddg_3(0, 0, 1, 2)
  181                     + t4ddg_2(0, 0, 1, 1) * t4ddg_3(0, 1, 1, 2)
  182                     + t4ddg_2(0, 0, 1, 2) * t4ddg_3(0, 2, 1, 2)
  183                     + t4ddg_2(0, 1, 1, 0) * t4ddg_3(1, 0, 1, 2)
  184                     + t4ddg_2(0, 1, 1, 1) * t4ddg_3(1, 1, 1, 2)
  185                     + t4ddg_2(0, 1, 1, 2) * t4ddg_3(1, 2, 1, 2)
  186                     + t4ddg_2(0, 2, 1, 0) * t4ddg_3(2, 0, 1, 2)
  187                     + t4ddg_2(0, 2, 1, 1) * t4ddg_3(2, 1, 1, 2)
  188                     + t4ddg_2(0, 2, 1, 2) * t4ddg_3(2, 2, 1, 2)),
  189                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,1,1,2)");
  191                  - (t4ddg_2(0, 0, 2, 0) * t4ddg_3(0, 0, 1, 0)
  192                     + t4ddg_2(0, 0, 2, 1) * t4ddg_3(0, 1, 1, 0)
  193                     + t4ddg_2(0, 0, 2, 2) * t4ddg_3(0, 2, 1, 0)
  194                     + t4ddg_2(0, 1, 2, 0) * t4ddg_3(1, 0, 1, 0)
  195                     + t4ddg_2(0, 1, 2, 1) * t4ddg_3(1, 1, 1, 0)
  196                     + t4ddg_2(0, 1, 2, 2) * t4ddg_3(1, 2, 1, 0)
  197                     + t4ddg_2(0, 2, 2, 0) * t4ddg_3(2, 0, 1, 0)
  198                     + t4ddg_2(0, 2, 2, 1) * t4ddg_3(2, 1, 1, 0)
  199                     + t4ddg_2(0, 2, 2, 2) * t4ddg_3(2, 2, 1, 0)),
  200                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,1,2,0)");
  202                  - (t4ddg_2(0, 0, 2, 0) * t4ddg_3(0, 0, 1, 1)
  203                     + t4ddg_2(0, 0, 2, 1) * t4ddg_3(0, 1, 1, 1)
  204                     + t4ddg_2(0, 0, 2, 2) * t4ddg_3(0, 2, 1, 1)
  205                     + t4ddg_2(0, 1, 2, 0) * t4ddg_3(1, 0, 1, 1)
  206                     + t4ddg_2(0, 1, 2, 1) * t4ddg_3(1, 1, 1, 1)
  207                     + t4ddg_2(0, 1, 2, 2) * t4ddg_3(1, 2, 1, 1)
  208                     + t4ddg_2(0, 2, 2, 0) * t4ddg_3(2, 0, 1, 1)
  209                     + t4ddg_2(0, 2, 2, 1) * t4ddg_3(2, 1, 1, 1)
  210                     + t4ddg_2(0, 2, 2, 2) * t4ddg_3(2, 2, 1, 1)),
  211                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,1,2,1)");
  213                  - (t4ddg_2(0, 0, 2, 0) * t4ddg_3(0, 0, 1, 2)
  214                     + t4ddg_2(0, 0, 2, 1) * t4ddg_3(0, 1, 1, 2)
  215                     + t4ddg_2(0, 0, 2, 2) * t4ddg_3(0, 2, 1, 2)
  216                     + t4ddg_2(0, 1, 2, 0) * t4ddg_3(1, 0, 1, 2)
  217                     + t4ddg_2(0, 1, 2, 1) * t4ddg_3(1, 1, 1, 2)
  218                     + t4ddg_2(0, 1, 2, 2) * t4ddg_3(1, 2, 1, 2)
  219                     + t4ddg_2(0, 2, 2, 0) * t4ddg_3(2, 0, 1, 2)
  220                     + t4ddg_2(0, 2, 2, 1) * t4ddg_3(2, 1, 1, 2)
  221                     + t4ddg_2(0, 2, 2, 2) * t4ddg_3(2, 2, 1, 2)),
  222                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,1,2,2)");
  224                  - (t4ddg_2(0, 0, 0, 0) * t4ddg_3(0, 0, 2, 0)
  225                     + t4ddg_2(0, 0, 0, 1) * t4ddg_3(0, 1, 2, 0)
  226                     + t4ddg_2(0, 0, 0, 2) * t4ddg_3(0, 2, 2, 0)
  227                     + t4ddg_2(0, 1, 0, 0) * t4ddg_3(1, 0, 2, 0)
  228                     + t4ddg_2(0, 1, 0, 1) * t4ddg_3(1, 1, 2, 0)
  229                     + t4ddg_2(0, 1, 0, 2) * t4ddg_3(1, 2, 2, 0)
  230                     + t4ddg_2(0, 2, 0, 0) * t4ddg_3(2, 0, 2, 0)
  231                     + t4ddg_2(0, 2, 0, 1) * t4ddg_3(2, 1, 2, 0)
  232                     + t4ddg_2(0, 2, 0, 2) * t4ddg_3(2, 2, 2, 0)),
  233                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,2,0,0)");
  235                  - (t4ddg_2(0, 0, 0, 0) * t4ddg_3(0, 0, 2, 1)
  236                     + t4ddg_2(0, 0, 0, 1) * t4ddg_3(0, 1, 2, 1)
  237                     + t4ddg_2(0, 0, 0, 2) * t4ddg_3(0, 2, 2, 1)
  238                     + t4ddg_2(0, 1, 0, 0) * t4ddg_3(1, 0, 2, 1)
  239                     + t4ddg_2(0, 1, 0, 1) * t4ddg_3(1, 1, 2, 1)
  240                     + t4ddg_2(0, 1, 0, 2) * t4ddg_3(1, 2, 2, 1)
  241                     + t4ddg_2(0, 2, 0, 0) * t4ddg_3(2, 0, 2, 1)
  242                     + t4ddg_2(0, 2, 0, 1) * t4ddg_3(2, 1, 2, 1)
  243                     + t4ddg_2(0, 2, 0, 2) * t4ddg_3(2, 2, 2, 1)),
  244                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,2,0,1)");
  246                  - (t4ddg_2(0, 0, 0, 0) * t4ddg_3(0, 0, 2, 2)
  247                     + t4ddg_2(0, 0, 0, 1) * t4ddg_3(0, 1, 2, 2)
  248                     + t4ddg_2(0, 0, 0, 2) * t4ddg_3(0, 2, 2, 2)
  249                     + t4ddg_2(0, 1, 0, 0) * t4ddg_3(1, 0, 2, 2)
  250                     + t4ddg_2(0, 1, 0, 1) * t4ddg_3(1, 1, 2, 2)
  251                     + t4ddg_2(0, 1, 0, 2) * t4ddg_3(1, 2, 2, 2)
  252                     + t4ddg_2(0, 2, 0, 0) * t4ddg_3(2, 0, 2, 2)
  253                     + t4ddg_2(0, 2, 0, 1) * t4ddg_3(2, 1, 2, 2)
  254                     + t4ddg_2(0, 2, 0, 2) * t4ddg_3(2, 2, 2, 2)),
  255                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,2,0,2)");
  257                  - (t4ddg_2(0, 0, 1, 0) * t4ddg_3(0, 0, 2, 0)
  258                     + t4ddg_2(0, 0, 1, 1) * t4ddg_3(0, 1, 2, 0)
  259                     + t4ddg_2(0, 0, 1, 2) * t4ddg_3(0, 2, 2, 0)
  260                     + t4ddg_2(0, 1, 1, 0) * t4ddg_3(1, 0, 2, 0)
  261                     + t4ddg_2(0, 1, 1, 1) * t4ddg_3(1, 1, 2, 0)
  262                     + t4ddg_2(0, 1, 1, 2) * t4ddg_3(1, 2, 2, 0)
  263                     + t4ddg_2(0, 2, 1, 0) * t4ddg_3(2, 0, 2, 0)
  264                     + t4ddg_2(0, 2, 1, 1) * t4ddg_3(2, 1, 2, 0)
  265                     + t4ddg_2(0, 2, 1, 2) * t4ddg_3(2, 2, 2, 0)),
  266                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,2,1,0)");
  268                  - (t4ddg_2(0, 0, 1, 0) * t4ddg_3(0, 0, 2, 1)
  269                     + t4ddg_2(0, 0, 1, 1) * t4ddg_3(0, 1, 2, 1)
  270                     + t4ddg_2(0, 0, 1, 2) * t4ddg_3(0, 2, 2, 1)
  271                     + t4ddg_2(0, 1, 1, 0) * t4ddg_3(1, 0, 2, 1)
  272                     + t4ddg_2(0, 1, 1, 1) * t4ddg_3(1, 1, 2, 1)
  273                     + t4ddg_2(0, 1, 1, 2) * t4ddg_3(1, 2, 2, 1)
  274                     + t4ddg_2(0, 2, 1, 0) * t4ddg_3(2, 0, 2, 1)
  275                     + t4ddg_2(0, 2, 1, 1) * t4ddg_3(2, 1, 2, 1)
  276                     + t4ddg_2(0, 2, 1, 2) * t4ddg_3(2, 2, 2, 1)),
  277                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,2,1,1)");
  279                  - (t4ddg_2(0, 0, 1, 0) * t4ddg_3(0, 0, 2, 2)
  280                     + t4ddg_2(0, 0, 1, 1) * t4ddg_3(0, 1, 2, 2)
  281                     + t4ddg_2(0, 0, 1, 2) * t4ddg_3(0, 2, 2, 2)
  282                     + t4ddg_2(0, 1, 1, 0) * t4ddg_3(1, 0, 2, 2)
  283                     + t4ddg_2(0, 1, 1, 1) * t4ddg_3(1, 1, 2, 2)
  284                     + t4ddg_2(0, 1, 1, 2) * t4ddg_3(1, 2, 2, 2)
  285                     + t4ddg_2(0, 2, 1, 0) * t4ddg_3(2, 0, 2, 2)
  286                     + t4ddg_2(0, 2, 1, 1) * t4ddg_3(2, 1, 2, 2)
  287                     + t4ddg_2(0, 2, 1, 2) * t4ddg_3(2, 2, 2, 2)),
  288                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,2,1,2)");
  290                  - (t4ddg_2(0, 0, 2, 0) * t4ddg_3(0, 0, 2, 0)
  291                     + t4ddg_2(0, 0, 2, 1) * t4ddg_3(0, 1, 2, 0)
  292                     + t4ddg_2(0, 0, 2, 2) * t4ddg_3(0, 2, 2, 0)
  293                     + t4ddg_2(0, 1, 2, 0) * t4ddg_3(1, 0, 2, 0)
  294                     + t4ddg_2(0, 1, 2, 1) * t4ddg_3(1, 1, 2, 0)
  295                     + t4ddg_2(0, 1, 2, 2) * t4ddg_3(1, 2, 2, 0)
  296                     + t4ddg_2(0, 2, 2, 0) * t4ddg_3(2, 0, 2, 0)
  297                     + t4ddg_2(0, 2, 2, 1) * t4ddg_3(2, 1, 2, 0)
  298                     + t4ddg_2(0, 2, 2, 2) * t4ddg_3(2, 2, 2, 0)),
  299                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,2,2,0)");
  301                  - (t4ddg_2(0, 0, 2, 0) * t4ddg_3(0, 0, 2, 1)
  302                     + t4ddg_2(0, 0, 2, 1) * t4ddg_3(0, 1, 2, 1)
  303                     + t4ddg_2(0, 0, 2, 2) * t4ddg_3(0, 2, 2, 1)
  304                     + t4ddg_2(0, 1, 2, 0) * t4ddg_3(1, 0, 2, 1)
  305                     + t4ddg_2(0, 1, 2, 1) * t4ddg_3(1, 1, 2, 1)
  306                     + t4ddg_2(0, 1, 2, 2) * t4ddg_3(1, 2, 2, 1)
  307                     + t4ddg_2(0, 2, 2, 0) * t4ddg_3(2, 0, 2, 1)
  308                     + t4ddg_2(0, 2, 2, 1) * t4ddg_3(2, 1, 2, 1)
  309                     + t4ddg_2(0, 2, 2, 2) * t4ddg_3(2, 2, 2, 1)),
  310                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,2,2,1)");
  312                  - (t4ddg_2(0, 0, 2, 0) * t4ddg_3(0, 0, 2, 2)
  313                     + t4ddg_2(0, 0, 2, 1) * t4ddg_3(0, 1, 2, 2)
  314                     + t4ddg_2(0, 0, 2, 2) * t4ddg_3(0, 2, 2, 2)
  315                     + t4ddg_2(0, 1, 2, 0) * t4ddg_3(1, 0, 2, 2)
  316                     + t4ddg_2(0, 1, 2, 1) * t4ddg_3(1, 1, 2, 2)
  317                     + t4ddg_2(0, 1, 2, 2) * t4ddg_3(1, 2, 2, 2)
  318                     + t4ddg_2(0, 2, 2, 0) * t4ddg_3(2, 0, 2, 2)
  319                     + t4ddg_2(0, 2, 2, 1) * t4ddg_3(2, 1, 2, 2)
  320                     + t4ddg_2(0, 2, 2, 2) * t4ddg_3(2, 2, 2, 2)),
  321                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(0,2,2,2)");
  323                  - (t4ddg_2(1, 0, 0, 0) * t4ddg_3(0, 0, 0, 0)
  324                     + t4ddg_2(1, 0, 0, 1) * t4ddg_3(0, 1, 0, 0)
  325                     + t4ddg_2(1, 0, 0, 2) * t4ddg_3(0, 2, 0, 0)
  326                     + t4ddg_2(1, 1, 0, 0) * t4ddg_3(1, 0, 0, 0)
  327                     + t4ddg_2(1, 1, 0, 1) * t4ddg_3(1, 1, 0, 0)
  328                     + t4ddg_2(1, 1, 0, 2) * t4ddg_3(1, 2, 0, 0)
  329                     + t4ddg_2(1, 2, 0, 0) * t4ddg_3(2, 0, 0, 0)
  330                     + t4ddg_2(1, 2, 0, 1) * t4ddg_3(2, 1, 0, 0)
  331                     + t4ddg_2(1, 2, 0, 2) * t4ddg_3(2, 2, 0, 0)),
  332                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,0,0,0)");
  334                  - (t4ddg_2(1, 0, 0, 0) * t4ddg_3(0, 0, 0, 1)
  335                     + t4ddg_2(1, 0, 0, 1) * t4ddg_3(0, 1, 0, 1)
  336                     + t4ddg_2(1, 0, 0, 2) * t4ddg_3(0, 2, 0, 1)
  337                     + t4ddg_2(1, 1, 0, 0) * t4ddg_3(1, 0, 0, 1)
  338                     + t4ddg_2(1, 1, 0, 1) * t4ddg_3(1, 1, 0, 1)
  339                     + t4ddg_2(1, 1, 0, 2) * t4ddg_3(1, 2, 0, 1)
  340                     + t4ddg_2(1, 2, 0, 0) * t4ddg_3(2, 0, 0, 1)
  341                     + t4ddg_2(1, 2, 0, 1) * t4ddg_3(2, 1, 0, 1)
  342                     + t4ddg_2(1, 2, 0, 2) * t4ddg_3(2, 2, 0, 1)),
  343                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,0,0,1)");
  345                  - (t4ddg_2(1, 0, 0, 0) * t4ddg_3(0, 0, 0, 2)
  346                     + t4ddg_2(1, 0, 0, 1) * t4ddg_3(0, 1, 0, 2)
  347                     + t4ddg_2(1, 0, 0, 2) * t4ddg_3(0, 2, 0, 2)
  348                     + t4ddg_2(1, 1, 0, 0) * t4ddg_3(1, 0, 0, 2)
  349                     + t4ddg_2(1, 1, 0, 1) * t4ddg_3(1, 1, 0, 2)
  350                     + t4ddg_2(1, 1, 0, 2) * t4ddg_3(1, 2, 0, 2)
  351                     + t4ddg_2(1, 2, 0, 0) * t4ddg_3(2, 0, 0, 2)
  352                     + t4ddg_2(1, 2, 0, 1) * t4ddg_3(2, 1, 0, 2)
  353                     + t4ddg_2(1, 2, 0, 2) * t4ddg_3(2, 2, 0, 2)),
  354                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,0,0,2)");
  356                  - (t4ddg_2(1, 0, 1, 0) * t4ddg_3(0, 0, 0, 0)
  357                     + t4ddg_2(1, 0, 1, 1) * t4ddg_3(0, 1, 0, 0)
  358                     + t4ddg_2(1, 0, 1, 2) * t4ddg_3(0, 2, 0, 0)
  359                     + t4ddg_2(1, 1, 1, 0) * t4ddg_3(1, 0, 0, 0)
  360                     + t4ddg_2(1, 1, 1, 1) * t4ddg_3(1, 1, 0, 0)
  361                     + t4ddg_2(1, 1, 1, 2) * t4ddg_3(1, 2, 0, 0)
  362                     + t4ddg_2(1, 2, 1, 0) * t4ddg_3(2, 0, 0, 0)
  363                     + t4ddg_2(1, 2, 1, 1) * t4ddg_3(2, 1, 0, 0)
  364                     + t4ddg_2(1, 2, 1, 2) * t4ddg_3(2, 2, 0, 0)),
  365                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,0,1,0)");
  367                  - (t4ddg_2(1, 0, 1, 0) * t4ddg_3(0, 0, 0, 1)
  368                     + t4ddg_2(1, 0, 1, 1) * t4ddg_3(0, 1, 0, 1)
  369                     + t4ddg_2(1, 0, 1, 2) * t4ddg_3(0, 2, 0, 1)
  370                     + t4ddg_2(1, 1, 1, 0) * t4ddg_3(1, 0, 0, 1)
  371                     + t4ddg_2(1, 1, 1, 1) * t4ddg_3(1, 1, 0, 1)
  372                     + t4ddg_2(1, 1, 1, 2) * t4ddg_3(1, 2, 0, 1)
  373                     + t4ddg_2(1, 2, 1, 0) * t4ddg_3(2, 0, 0, 1)
  374                     + t4ddg_2(1, 2, 1, 1) * t4ddg_3(2, 1, 0, 1)
  375                     + t4ddg_2(1, 2, 1, 2) * t4ddg_3(2, 2, 0, 1)),
  376                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,0,1,1)");
  378                  - (t4ddg_2(1, 0, 1, 0) * t4ddg_3(0, 0, 0, 2)
  379                     + t4ddg_2(1, 0, 1, 1) * t4ddg_3(0, 1, 0, 2)
  380                     + t4ddg_2(1, 0, 1, 2) * t4ddg_3(0, 2, 0, 2)
  381                     + t4ddg_2(1, 1, 1, 0) * t4ddg_3(1, 0, 0, 2)
  382                     + t4ddg_2(1, 1, 1, 1) * t4ddg_3(1, 1, 0, 2)
  383                     + t4ddg_2(1, 1, 1, 2) * t4ddg_3(1, 2, 0, 2)
  384                     + t4ddg_2(1, 2, 1, 0) * t4ddg_3(2, 0, 0, 2)
  385                     + t4ddg_2(1, 2, 1, 1) * t4ddg_3(2, 1, 0, 2)
  386                     + t4ddg_2(1, 2, 1, 2) * t4ddg_3(2, 2, 0, 2)),
  387                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,0,1,2)");
  389                  - (t4ddg_2(1, 0, 2, 0) * t4ddg_3(0, 0, 0, 0)
  390                     + t4ddg_2(1, 0, 2, 1) * t4ddg_3(0, 1, 0, 0)
  391                     + t4ddg_2(1, 0, 2, 2) * t4ddg_3(0, 2, 0, 0)
  392                     + t4ddg_2(1, 1, 2, 0) * t4ddg_3(1, 0, 0, 0)
  393                     + t4ddg_2(1, 1, 2, 1) * t4ddg_3(1, 1, 0, 0)
  394                     + t4ddg_2(1, 1, 2, 2) * t4ddg_3(1, 2, 0, 0)
  395                     + t4ddg_2(1, 2, 2, 0) * t4ddg_3(2, 0, 0, 0)
  396                     + t4ddg_2(1, 2, 2, 1) * t4ddg_3(2, 1, 0, 0)
  397                     + t4ddg_2(1, 2, 2, 2) * t4ddg_3(2, 2, 0, 0)),
  398                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,0,2,0)");
  400                  - (t4ddg_2(1, 0, 2, 0) * t4ddg_3(0, 0, 0, 1)
  401                     + t4ddg_2(1, 0, 2, 1) * t4ddg_3(0, 1, 0, 1)
  402                     + t4ddg_2(1, 0, 2, 2) * t4ddg_3(0, 2, 0, 1)
  403                     + t4ddg_2(1, 1, 2, 0) * t4ddg_3(1, 0, 0, 1)
  404                     + t4ddg_2(1, 1, 2, 1) * t4ddg_3(1, 1, 0, 1)
  405                     + t4ddg_2(1, 1, 2, 2) * t4ddg_3(1, 2, 0, 1)
  406                     + t4ddg_2(1, 2, 2, 0) * t4ddg_3(2, 0, 0, 1)
  407                     + t4ddg_2(1, 2, 2, 1) * t4ddg_3(2, 1, 0, 1)
  408                     + t4ddg_2(1, 2, 2, 2) * t4ddg_3(2, 2, 0, 1)),
  409                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,0,2,1)");
  411                  - (t4ddg_2(1, 0, 2, 0) * t4ddg_3(0, 0, 0, 2)
  412                     + t4ddg_2(1, 0, 2, 1) * t4ddg_3(0, 1, 0, 2)
  413                     + t4ddg_2(1, 0, 2, 2) * t4ddg_3(0, 2, 0, 2)
  414                     + t4ddg_2(1, 1, 2, 0) * t4ddg_3(1, 0, 0, 2)
  415                     + t4ddg_2(1, 1, 2, 1) * t4ddg_3(1, 1, 0, 2)
  416                     + t4ddg_2(1, 1, 2, 2) * t4ddg_3(1, 2, 0, 2)
  417                     + t4ddg_2(1, 2, 2, 0) * t4ddg_3(2, 0, 0, 2)
  418                     + t4ddg_2(1, 2, 2, 1) * t4ddg_3(2, 1, 0, 2)
  419                     + t4ddg_2(1, 2, 2, 2) * t4ddg_3(2, 2, 0, 2)),
  420                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,0,2,2)");
  422                  - (t4ddg_2(1, 0, 0, 0) * t4ddg_3(0, 0, 1, 0)
  423                     + t4ddg_2(1, 0, 0, 1) * t4ddg_3(0, 1, 1, 0)
  424                     + t4ddg_2(1, 0, 0, 2) * t4ddg_3(0, 2, 1, 0)
  425                     + t4ddg_2(1, 1, 0, 0) * t4ddg_3(1, 0, 1, 0)
  426                     + t4ddg_2(1, 1, 0, 1) * t4ddg_3(1, 1, 1, 0)
  427                     + t4ddg_2(1, 1, 0, 2) * t4ddg_3(1, 2, 1, 0)
  428                     + t4ddg_2(1, 2, 0, 0) * t4ddg_3(2, 0, 1, 0)
  429                     + t4ddg_2(1, 2, 0, 1) * t4ddg_3(2, 1, 1, 0)
  430                     + t4ddg_2(1, 2, 0, 2) * t4ddg_3(2, 2, 1, 0)),
  431                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,1,0,0)");
  433                  - (t4ddg_2(1, 0, 0, 0) * t4ddg_3(0, 0, 1, 1)
  434                     + t4ddg_2(1, 0, 0, 1) * t4ddg_3(0, 1, 1, 1)
  435                     + t4ddg_2(1, 0, 0, 2) * t4ddg_3(0, 2, 1, 1)
  436                     + t4ddg_2(1, 1, 0, 0) * t4ddg_3(1, 0, 1, 1)
  437                     + t4ddg_2(1, 1, 0, 1) * t4ddg_3(1, 1, 1, 1)
  438                     + t4ddg_2(1, 1, 0, 2) * t4ddg_3(1, 2, 1, 1)
  439                     + t4ddg_2(1, 2, 0, 0) * t4ddg_3(2, 0, 1, 1)
  440                     + t4ddg_2(1, 2, 0, 1) * t4ddg_3(2, 1, 1, 1)
  441                     + t4ddg_2(1, 2, 0, 2) * t4ddg_3(2, 2, 1, 1)),
  442                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,1,0,1)");
  444                  - (t4ddg_2(1, 0, 0, 0) * t4ddg_3(0, 0, 1, 2)
  445                     + t4ddg_2(1, 0, 0, 1) * t4ddg_3(0, 1, 1, 2)
  446                     + t4ddg_2(1, 0, 0, 2) * t4ddg_3(0, 2, 1, 2)
  447                     + t4ddg_2(1, 1, 0, 0) * t4ddg_3(1, 0, 1, 2)
  448                     + t4ddg_2(1, 1, 0, 1) * t4ddg_3(1, 1, 1, 2)
  449                     + t4ddg_2(1, 1, 0, 2) * t4ddg_3(1, 2, 1, 2)
  450                     + t4ddg_2(1, 2, 0, 0) * t4ddg_3(2, 0, 1, 2)
  451                     + t4ddg_2(1, 2, 0, 1) * t4ddg_3(2, 1, 1, 2)
  452                     + t4ddg_2(1, 2, 0, 2) * t4ddg_3(2, 2, 1, 2)),
  453                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,1,0,2)");
  455                  - (t4ddg_2(1, 0, 1, 0) * t4ddg_3(0, 0, 1, 0)
  456                     + t4ddg_2(1, 0, 1, 1) * t4ddg_3(0, 1, 1, 0)
  457                     + t4ddg_2(1, 0, 1, 2) * t4ddg_3(0, 2, 1, 0)
  458                     + t4ddg_2(1, 1, 1, 0) * t4ddg_3(1, 0, 1, 0)
  459                     + t4ddg_2(1, 1, 1, 1) * t4ddg_3(1, 1, 1, 0)
  460                     + t4ddg_2(1, 1, 1, 2) * t4ddg_3(1, 2, 1, 0)
  461                     + t4ddg_2(1, 2, 1, 0) * t4ddg_3(2, 0, 1, 0)
  462                     + t4ddg_2(1, 2, 1, 1) * t4ddg_3(2, 1, 1, 0)
  463                     + t4ddg_2(1, 2, 1, 2) * t4ddg_3(2, 2, 1, 0)),
  464                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,1,1,0)");
  466                  - (t4ddg_2(1, 0, 1, 0) * t4ddg_3(0, 0, 1, 1)
  467                     + t4ddg_2(1, 0, 1, 1) * t4ddg_3(0, 1, 1, 1)
  468                     + t4ddg_2(1, 0, 1, 2) * t4ddg_3(0, 2, 1, 1)
  469                     + t4ddg_2(1, 1, 1, 0) * t4ddg_3(1, 0, 1, 1)
  470                     + t4ddg_2(1, 1, 1, 1) * t4ddg_3(1, 1, 1, 1)
  471                     + t4ddg_2(1, 1, 1, 2) * t4ddg_3(1, 2, 1, 1)
  472                     + t4ddg_2(1, 2, 1, 0) * t4ddg_3(2, 0, 1, 1)
  473                     + t4ddg_2(1, 2, 1, 1) * t4ddg_3(2, 1, 1, 1)
  474                     + t4ddg_2(1, 2, 1, 2) * t4ddg_3(2, 2, 1, 1)),
  475                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,1,1,1)");
  477                  - (t4ddg_2(1, 0, 1, 0) * t4ddg_3(0, 0, 1, 2)
  478                     + t4ddg_2(1, 0, 1, 1) * t4ddg_3(0, 1, 1, 2)
  479                     + t4ddg_2(1, 0, 1, 2) * t4ddg_3(0, 2, 1, 2)
  480                     + t4ddg_2(1, 1, 1, 0) * t4ddg_3(1, 0, 1, 2)
  481                     + t4ddg_2(1, 1, 1, 1) * t4ddg_3(1, 1, 1, 2)
  482                     + t4ddg_2(1, 1, 1, 2) * t4ddg_3(1, 2, 1, 2)
  483                     + t4ddg_2(1, 2, 1, 0) * t4ddg_3(2, 0, 1, 2)
  484                     + t4ddg_2(1, 2, 1, 1) * t4ddg_3(2, 1, 1, 2)
  485                     + t4ddg_2(1, 2, 1, 2) * t4ddg_3(2, 2, 1, 2)),
  486                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,1,1,2)");
  488                  - (t4ddg_2(1, 0, 2, 0) * t4ddg_3(0, 0, 1, 0)
  489                     + t4ddg_2(1, 0, 2, 1) * t4ddg_3(0, 1, 1, 0)
  490                     + t4ddg_2(1, 0, 2, 2) * t4ddg_3(0, 2, 1, 0)
  491                     + t4ddg_2(1, 1, 2, 0) * t4ddg_3(1, 0, 1, 0)
  492                     + t4ddg_2(1, 1, 2, 1) * t4ddg_3(1, 1, 1, 0)
  493                     + t4ddg_2(1, 1, 2, 2) * t4ddg_3(1, 2, 1, 0)
  494                     + t4ddg_2(1, 2, 2, 0) * t4ddg_3(2, 0, 1, 0)
  495                     + t4ddg_2(1, 2, 2, 1) * t4ddg_3(2, 1, 1, 0)
  496                     + t4ddg_2(1, 2, 2, 2) * t4ddg_3(2, 2, 1, 0)),
  497                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,1,2,0)");
  499                  - (t4ddg_2(1, 0, 2, 0) * t4ddg_3(0, 0, 1, 1)
  500                     + t4ddg_2(1, 0, 2, 1) * t4ddg_3(0, 1, 1, 1)
  501                     + t4ddg_2(1, 0, 2, 2) * t4ddg_3(0, 2, 1, 1)
  502                     + t4ddg_2(1, 1, 2, 0) * t4ddg_3(1, 0, 1, 1)
  503                     + t4ddg_2(1, 1, 2, 1) * t4ddg_3(1, 1, 1, 1)
  504                     + t4ddg_2(1, 1, 2, 2) * t4ddg_3(1, 2, 1, 1)
  505                     + t4ddg_2(1, 2, 2, 0) * t4ddg_3(2, 0, 1, 1)
  506                     + t4ddg_2(1, 2, 2, 1) * t4ddg_3(2, 1, 1, 1)
  507                     + t4ddg_2(1, 2, 2, 2) * t4ddg_3(2, 2, 1, 1)),
  508                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,1,2,1)");
  510                  - (t4ddg_2(1, 0, 2, 0) * t4ddg_3(0, 0, 1, 2)
  511                     + t4ddg_2(1, 0, 2, 1) * t4ddg_3(0, 1, 1, 2)
  512                     + t4ddg_2(1, 0, 2, 2) * t4ddg_3(0, 2, 1, 2)
  513                     + t4ddg_2(1, 1, 2, 0) * t4ddg_3(1, 0, 1, 2)
  514                     + t4ddg_2(1, 1, 2, 1) * t4ddg_3(1, 1, 1, 2)
  515                     + t4ddg_2(1, 1, 2, 2) * t4ddg_3(1, 2, 1, 2)
  516                     + t4ddg_2(1, 2, 2, 0) * t4ddg_3(2, 0, 1, 2)
  517                     + t4ddg_2(1, 2, 2, 1) * t4ddg_3(2, 1, 1, 2)
  518                     + t4ddg_2(1, 2, 2, 2) * t4ddg_3(2, 2, 1, 2)),
  519                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,1,2,2)");
  521                  - (t4ddg_2(1, 0, 0, 0) * t4ddg_3(0, 0, 2, 0)
  522                     + t4ddg_2(1, 0, 0, 1) * t4ddg_3(0, 1, 2, 0)
  523                     + t4ddg_2(1, 0, 0, 2) * t4ddg_3(0, 2, 2, 0)
  524                     + t4ddg_2(1, 1, 0, 0) * t4ddg_3(1, 0, 2, 0)
  525                     + t4ddg_2(1, 1, 0, 1) * t4ddg_3(1, 1, 2, 0)
  526                     + t4ddg_2(1, 1, 0, 2) * t4ddg_3(1, 2, 2, 0)
  527                     + t4ddg_2(1, 2, 0, 0) * t4ddg_3(2, 0, 2, 0)
  528                     + t4ddg_2(1, 2, 0, 1) * t4ddg_3(2, 1, 2, 0)
  529                     + t4ddg_2(1, 2, 0, 2) * t4ddg_3(2, 2, 2, 0)),
  530                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,2,0,0)");
  532                  - (t4ddg_2(1, 0, 0, 0) * t4ddg_3(0, 0, 2, 1)
  533                     + t4ddg_2(1, 0, 0, 1) * t4ddg_3(0, 1, 2, 1)
  534                     + t4ddg_2(1, 0, 0, 2) * t4ddg_3(0, 2, 2, 1)
  535                     + t4ddg_2(1, 1, 0, 0) * t4ddg_3(1, 0, 2, 1)
  536                     + t4ddg_2(1, 1, 0, 1) * t4ddg_3(1, 1, 2, 1)
  537                     + t4ddg_2(1, 1, 0, 2) * t4ddg_3(1, 2, 2, 1)
  538                     + t4ddg_2(1, 2, 0, 0) * t4ddg_3(2, 0, 2, 1)
  539                     + t4ddg_2(1, 2, 0, 1) * t4ddg_3(2, 1, 2, 1)
  540                     + t4ddg_2(1, 2, 0, 2) * t4ddg_3(2, 2, 2, 1)),
  541                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,2,0,1)");
  543                  - (t4ddg_2(1, 0, 0, 0) * t4ddg_3(0, 0, 2, 2)
  544                     + t4ddg_2(1, 0, 0, 1) * t4ddg_3(0, 1, 2, 2)
  545                     + t4ddg_2(1, 0, 0, 2) * t4ddg_3(0, 2, 2, 2)
  546                     + t4ddg_2(1, 1, 0, 0) * t4ddg_3(1, 0, 2, 2)
  547                     + t4ddg_2(1, 1, 0, 1) * t4ddg_3(1, 1, 2, 2)
  548                     + t4ddg_2(1, 1, 0, 2) * t4ddg_3(1, 2, 2, 2)
  549                     + t4ddg_2(1, 2, 0, 0) * t4ddg_3(2, 0, 2, 2)
  550                     + t4ddg_2(1, 2, 0, 1) * t4ddg_3(2, 1, 2, 2)
  551                     + t4ddg_2(1, 2, 0, 2) * t4ddg_3(2, 2, 2, 2)),
  552                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,2,0,2)");
  554                  - (t4ddg_2(1, 0, 1, 0) * t4ddg_3(0, 0, 2, 0)
  555                     + t4ddg_2(1, 0, 1, 1) * t4ddg_3(0, 1, 2, 0)
  556                     + t4ddg_2(1, 0, 1, 2) * t4ddg_3(0, 2, 2, 0)
  557                     + t4ddg_2(1, 1, 1, 0) * t4ddg_3(1, 0, 2, 0)
  558                     + t4ddg_2(1, 1, 1, 1) * t4ddg_3(1, 1, 2, 0)
  559                     + t4ddg_2(1, 1, 1, 2) * t4ddg_3(1, 2, 2, 0)
  560                     + t4ddg_2(1, 2, 1, 0) * t4ddg_3(2, 0, 2, 0)
  561                     + t4ddg_2(1, 2, 1, 1) * t4ddg_3(2, 1, 2, 0)
  562                     + t4ddg_2(1, 2, 1, 2) * t4ddg_3(2, 2, 2, 0)),
  563                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,2,1,0)");
  565                  - (t4ddg_2(1, 0, 1, 0) * t4ddg_3(0, 0, 2, 1)
  566                     + t4ddg_2(1, 0, 1, 1) * t4ddg_3(0, 1, 2, 1)
  567                     + t4ddg_2(1, 0, 1, 2) * t4ddg_3(0, 2, 2, 1)
  568                     + t4ddg_2(1, 1, 1, 0) * t4ddg_3(1, 0, 2, 1)
  569                     + t4ddg_2(1, 1, 1, 1) * t4ddg_3(1, 1, 2, 1)
  570                     + t4ddg_2(1, 1, 1, 2) * t4ddg_3(1, 2, 2, 1)
  571                     + t4ddg_2(1, 2, 1, 0) * t4ddg_3(2, 0, 2, 1)
  572                     + t4ddg_2(1, 2, 1, 1) * t4ddg_3(2, 1, 2, 1)
  573                     + t4ddg_2(1, 2, 1, 2) * t4ddg_3(2, 2, 2, 1)),
  574                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,2,1,1)");
  576                  - (t4ddg_2(1, 0, 1, 0) * t4ddg_3(0, 0, 2, 2)
  577                     + t4ddg_2(1, 0, 1, 1) * t4ddg_3(0, 1, 2, 2)
  578                     + t4ddg_2(1, 0, 1, 2) * t4ddg_3(0, 2, 2, 2)
  579                     + t4ddg_2(1, 1, 1, 0) * t4ddg_3(1, 0, 2, 2)
  580                     + t4ddg_2(1, 1, 1, 1) * t4ddg_3(1, 1, 2, 2)
  581                     + t4ddg_2(1, 1, 1, 2) * t4ddg_3(1, 2, 2, 2)
  582                     + t4ddg_2(1, 2, 1, 0) * t4ddg_3(2, 0, 2, 2)
  583                     + t4ddg_2(1, 2, 1, 1) * t4ddg_3(2, 1, 2, 2)
  584                     + t4ddg_2(1, 2, 1, 2) * t4ddg_3(2, 2, 2, 2)),
  585                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,2,1,2)");
  587                  - (t4ddg_2(1, 0, 2, 0) * t4ddg_3(0, 0, 2, 0)
  588                     + t4ddg_2(1, 0, 2, 1) * t4ddg_3(0, 1, 2, 0)
  589                     + t4ddg_2(1, 0, 2, 2) * t4ddg_3(0, 2, 2, 0)
  590                     + t4ddg_2(1, 1, 2, 0) * t4ddg_3(1, 0, 2, 0)
  591                     + t4ddg_2(1, 1, 2, 1) * t4ddg_3(1, 1, 2, 0)
  592                     + t4ddg_2(1, 1, 2, 2) * t4ddg_3(1, 2, 2, 0)
  593                     + t4ddg_2(1, 2, 2, 0) * t4ddg_3(2, 0, 2, 0)
  594                     + t4ddg_2(1, 2, 2, 1) * t4ddg_3(2, 1, 2, 0)
  595                     + t4ddg_2(1, 2, 2, 2) * t4ddg_3(2, 2, 2, 0)),
  596                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,2,2,0)");
  598                  - (t4ddg_2(1, 0, 2, 0) * t4ddg_3(0, 0, 2, 1)
  599                     + t4ddg_2(1, 0, 2, 1) * t4ddg_3(0, 1, 2, 1)
  600                     + t4ddg_2(1, 0, 2, 2) * t4ddg_3(0, 2, 2, 1)
  601                     + t4ddg_2(1, 1, 2, 0) * t4ddg_3(1, 0, 2, 1)
  602                     + t4ddg_2(1, 1, 2, 1) * t4ddg_3(1, 1, 2, 1)
  603                     + t4ddg_2(1, 1, 2, 2) * t4ddg_3(1, 2, 2, 1)
  604                     + t4ddg_2(1, 2, 2, 0) * t4ddg_3(2, 0, 2, 1)
  605                     + t4ddg_2(1, 2, 2, 1) * t4ddg_3(2, 1, 2, 1)
  606                     + t4ddg_2(1, 2, 2, 2) * t4ddg_3(2, 2, 2, 1)),
  607                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,2,2,1)");
  609                  - (t4ddg_2(1, 0, 2, 0) * t4ddg_3(0, 0, 2, 2)
  610                     + t4ddg_2(1, 0, 2, 1) * t4ddg_3(0, 1, 2, 2)
  611                     + t4ddg_2(1, 0, 2, 2) * t4ddg_3(0, 2, 2, 2)
  612                     + t4ddg_2(1, 1, 2, 0) * t4ddg_3(1, 0, 2, 2)
  613                     + t4ddg_2(1, 1, 2, 1) * t4ddg_3(1, 1, 2, 2)
  614                     + t4ddg_2(1, 1, 2, 2) * t4ddg_3(1, 2, 2, 2)
  615                     + t4ddg_2(1, 2, 2, 0) * t4ddg_3(2, 0, 2, 2)
  616                     + t4ddg_2(1, 2, 2, 1) * t4ddg_3(2, 1, 2, 2)
  617                     + t4ddg_2(1, 2, 2, 2) * t4ddg_3(2, 2, 2, 2)),
  618                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(1,2,2,2)");
  620                  - (t4ddg_2(2, 0, 0, 0) * t4ddg_3(0, 0, 0, 0)
  621                     + t4ddg_2(2, 0, 0, 1) * t4ddg_3(0, 1, 0, 0)
  622                     + t4ddg_2(2, 0, 0, 2) * t4ddg_3(0, 2, 0, 0)
  623                     + t4ddg_2(2, 1, 0, 0) * t4ddg_3(1, 0, 0, 0)
  624                     + t4ddg_2(2, 1, 0, 1) * t4ddg_3(1, 1, 0, 0)
  625                     + t4ddg_2(2, 1, 0, 2) * t4ddg_3(1, 2, 0, 0)
  626                     + t4ddg_2(2, 2, 0, 0) * t4ddg_3(2, 0, 0, 0)
  627                     + t4ddg_2(2, 2, 0, 1) * t4ddg_3(2, 1, 0, 0)
  628                     + t4ddg_2(2, 2, 0, 2) * t4ddg_3(2, 2, 0, 0)),
  629                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,0,0,0)");
  631                  - (t4ddg_2(2, 0, 0, 0) * t4ddg_3(0, 0, 0, 1)
  632                     + t4ddg_2(2, 0, 0, 1) * t4ddg_3(0, 1, 0, 1)
  633                     + t4ddg_2(2, 0, 0, 2) * t4ddg_3(0, 2, 0, 1)
  634                     + t4ddg_2(2, 1, 0, 0) * t4ddg_3(1, 0, 0, 1)
  635                     + t4ddg_2(2, 1, 0, 1) * t4ddg_3(1, 1, 0, 1)
  636                     + t4ddg_2(2, 1, 0, 2) * t4ddg_3(1, 2, 0, 1)
  637                     + t4ddg_2(2, 2, 0, 0) * t4ddg_3(2, 0, 0, 1)
  638                     + t4ddg_2(2, 2, 0, 1) * t4ddg_3(2, 1, 0, 1)
  639                     + t4ddg_2(2, 2, 0, 2) * t4ddg_3(2, 2, 0, 1)),
  640                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,0,0,1)");
  642                  - (t4ddg_2(2, 0, 0, 0) * t4ddg_3(0, 0, 0, 2)
  643                     + t4ddg_2(2, 0, 0, 1) * t4ddg_3(0, 1, 0, 2)
  644                     + t4ddg_2(2, 0, 0, 2) * t4ddg_3(0, 2, 0, 2)
  645                     + t4ddg_2(2, 1, 0, 0) * t4ddg_3(1, 0, 0, 2)
  646                     + t4ddg_2(2, 1, 0, 1) * t4ddg_3(1, 1, 0, 2)
  647                     + t4ddg_2(2, 1, 0, 2) * t4ddg_3(1, 2, 0, 2)
  648                     + t4ddg_2(2, 2, 0, 0) * t4ddg_3(2, 0, 0, 2)
  649                     + t4ddg_2(2, 2, 0, 1) * t4ddg_3(2, 1, 0, 2)
  650                     + t4ddg_2(2, 2, 0, 2) * t4ddg_3(2, 2, 0, 2)),
  651                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,0,0,2)");
  653                  - (t4ddg_2(2, 0, 1, 0) * t4ddg_3(0, 0, 0, 0)
  654                     + t4ddg_2(2, 0, 1, 1) * t4ddg_3(0, 1, 0, 0)
  655                     + t4ddg_2(2, 0, 1, 2) * t4ddg_3(0, 2, 0, 0)
  656                     + t4ddg_2(2, 1, 1, 0) * t4ddg_3(1, 0, 0, 0)
  657                     + t4ddg_2(2, 1, 1, 1) * t4ddg_3(1, 1, 0, 0)
  658                     + t4ddg_2(2, 1, 1, 2) * t4ddg_3(1, 2, 0, 0)
  659                     + t4ddg_2(2, 2, 1, 0) * t4ddg_3(2, 0, 0, 0)
  660                     + t4ddg_2(2, 2, 1, 1) * t4ddg_3(2, 1, 0, 0)
  661                     + t4ddg_2(2, 2, 1, 2) * t4ddg_3(2, 2, 0, 0)),
  662                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,0,1,0)");
  664                  - (t4ddg_2(2, 0, 1, 0) * t4ddg_3(0, 0, 0, 1)
  665                     + t4ddg_2(2, 0, 1, 1) * t4ddg_3(0, 1, 0, 1)
  666                     + t4ddg_2(2, 0, 1, 2) * t4ddg_3(0, 2, 0, 1)
  667                     + t4ddg_2(2, 1, 1, 0) * t4ddg_3(1, 0, 0, 1)
  668                     + t4ddg_2(2, 1, 1, 1) * t4ddg_3(1, 1, 0, 1)
  669                     + t4ddg_2(2, 1, 1, 2) * t4ddg_3(1, 2, 0, 1)
  670                     + t4ddg_2(2, 2, 1, 0) * t4ddg_3(2, 0, 0, 1)
  671                     + t4ddg_2(2, 2, 1, 1) * t4ddg_3(2, 1, 0, 1)
  672                     + t4ddg_2(2, 2, 1, 2) * t4ddg_3(2, 2, 0, 1)),
  673                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,0,1,1)");
  675                  - (t4ddg_2(2, 0, 1, 0) * t4ddg_3(0, 0, 0, 2)
  676                     + t4ddg_2(2, 0, 1, 1) * t4ddg_3(0, 1, 0, 2)
  677                     + t4ddg_2(2, 0, 1, 2) * t4ddg_3(0, 2, 0, 2)
  678                     + t4ddg_2(2, 1, 1, 0) * t4ddg_3(1, 0, 0, 2)
  679                     + t4ddg_2(2, 1, 1, 1) * t4ddg_3(1, 1, 0, 2)
  680                     + t4ddg_2(2, 1, 1, 2) * t4ddg_3(1, 2, 0, 2)
  681                     + t4ddg_2(2, 2, 1, 0) * t4ddg_3(2, 0, 0, 2)
  682                     + t4ddg_2(2, 2, 1, 1) * t4ddg_3(2, 1, 0, 2)
  683                     + t4ddg_2(2, 2, 1, 2) * t4ddg_3(2, 2, 0, 2)),
  684                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,0,1,2)");
  686                  - (t4ddg_2(2, 0, 2, 0) * t4ddg_3(0, 0, 0, 0)
  687                     + t4ddg_2(2, 0, 2, 1) * t4ddg_3(0, 1, 0, 0)
  688                     + t4ddg_2(2, 0, 2, 2) * t4ddg_3(0, 2, 0, 0)
  689                     + t4ddg_2(2, 1, 2, 0) * t4ddg_3(1, 0, 0, 0)
  690                     + t4ddg_2(2, 1, 2, 1) * t4ddg_3(1, 1, 0, 0)
  691                     + t4ddg_2(2, 1, 2, 2) * t4ddg_3(1, 2, 0, 0)
  692                     + t4ddg_2(2, 2, 2, 0) * t4ddg_3(2, 0, 0, 0)
  693                     + t4ddg_2(2, 2, 2, 1) * t4ddg_3(2, 1, 0, 0)
  694                     + t4ddg_2(2, 2, 2, 2) * t4ddg_3(2, 2, 0, 0)),
  695                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,0,2,0)");
  697                  - (t4ddg_2(2, 0, 2, 0) * t4ddg_3(0, 0, 0, 1)
  698                     + t4ddg_2(2, 0, 2, 1) * t4ddg_3(0, 1, 0, 1)
  699                     + t4ddg_2(2, 0, 2, 2) * t4ddg_3(0, 2, 0, 1)
  700                     + t4ddg_2(2, 1, 2, 0) * t4ddg_3(1, 0, 0, 1)
  701                     + t4ddg_2(2, 1, 2, 1) * t4ddg_3(1, 1, 0, 1)
  702                     + t4ddg_2(2, 1, 2, 2) * t4ddg_3(1, 2, 0, 1)
  703                     + t4ddg_2(2, 2, 2, 0) * t4ddg_3(2, 0, 0, 1)
  704                     + t4ddg_2(2, 2, 2, 1) * t4ddg_3(2, 1, 0, 1)
  705                     + t4ddg_2(2, 2, 2, 2) * t4ddg_3(2, 2, 0, 1)),
  706                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,0,2,1)");
  708                  - (t4ddg_2(2, 0, 2, 0) * t4ddg_3(0, 0, 0, 2)
  709                     + t4ddg_2(2, 0, 2, 1) * t4ddg_3(0, 1, 0, 2)
  710                     + t4ddg_2(2, 0, 2, 2) * t4ddg_3(0, 2, 0, 2)
  711                     + t4ddg_2(2, 1, 2, 0) * t4ddg_3(1, 0, 0, 2)
  712                     + t4ddg_2(2, 1, 2, 1) * t4ddg_3(1, 1, 0, 2)
  713                     + t4ddg_2(2, 1, 2, 2) * t4ddg_3(1, 2, 0, 2)
  714                     + t4ddg_2(2, 2, 2, 0) * t4ddg_3(2, 0, 0, 2)
  715                     + t4ddg_2(2, 2, 2, 1) * t4ddg_3(2, 1, 0, 2)
  716                     + t4ddg_2(2, 2, 2, 2) * t4ddg_3(2, 2, 0, 2)),
  717                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,0,2,2)");
  719                  - (t4ddg_2(2, 0, 0, 0) * t4ddg_3(0, 0, 1, 0)
  720                     + t4ddg_2(2, 0, 0, 1) * t4ddg_3(0, 1, 1, 0)
  721                     + t4ddg_2(2, 0, 0, 2) * t4ddg_3(0, 2, 1, 0)
  722                     + t4ddg_2(2, 1, 0, 0) * t4ddg_3(1, 0, 1, 0)
  723                     + t4ddg_2(2, 1, 0, 1) * t4ddg_3(1, 1, 1, 0)
  724                     + t4ddg_2(2, 1, 0, 2) * t4ddg_3(1, 2, 1, 0)
  725                     + t4ddg_2(2, 2, 0, 0) * t4ddg_3(2, 0, 1, 0)
  726                     + t4ddg_2(2, 2, 0, 1) * t4ddg_3(2, 1, 1, 0)
  727                     + t4ddg_2(2, 2, 0, 2) * t4ddg_3(2, 2, 1, 0)),
  728                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,1,0,0)");
  730                  - (t4ddg_2(2, 0, 0, 0) * t4ddg_3(0, 0, 1, 1)
  731                     + t4ddg_2(2, 0, 0, 1) * t4ddg_3(0, 1, 1, 1)
  732                     + t4ddg_2(2, 0, 0, 2) * t4ddg_3(0, 2, 1, 1)
  733                     + t4ddg_2(2, 1, 0, 0) * t4ddg_3(1, 0, 1, 1)
  734                     + t4ddg_2(2, 1, 0, 1) * t4ddg_3(1, 1, 1, 1)
  735                     + t4ddg_2(2, 1, 0, 2) * t4ddg_3(1, 2, 1, 1)
  736                     + t4ddg_2(2, 2, 0, 0) * t4ddg_3(2, 0, 1, 1)
  737                     + t4ddg_2(2, 2, 0, 1) * t4ddg_3(2, 1, 1, 1)
  738                     + t4ddg_2(2, 2, 0, 2) * t4ddg_3(2, 2, 1, 1)),
  739                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,1,0,1)");
  741                  - (t4ddg_2(2, 0, 0, 0) * t4ddg_3(0, 0, 1, 2)
  742                     + t4ddg_2(2, 0, 0, 1) * t4ddg_3(0, 1, 1, 2)
  743                     + t4ddg_2(2, 0, 0, 2) * t4ddg_3(0, 2, 1, 2)
  744                     + t4ddg_2(2, 1, 0, 0) * t4ddg_3(1, 0, 1, 2)
  745                     + t4ddg_2(2, 1, 0, 1) * t4ddg_3(1, 1, 1, 2)
  746                     + t4ddg_2(2, 1, 0, 2) * t4ddg_3(1, 2, 1, 2)
  747                     + t4ddg_2(2, 2, 0, 0) * t4ddg_3(2, 0, 1, 2)
  748                     + t4ddg_2(2, 2, 0, 1) * t4ddg_3(2, 1, 1, 2)
  749                     + t4ddg_2(2, 2, 0, 2) * t4ddg_3(2, 2, 1, 2)),
  750                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,1,0,2)");
  752                  - (t4ddg_2(2, 0, 1, 0) * t4ddg_3(0, 0, 1, 0)
  753                     + t4ddg_2(2, 0, 1, 1) * t4ddg_3(0, 1, 1, 0)
  754                     + t4ddg_2(2, 0, 1, 2) * t4ddg_3(0, 2, 1, 0)
  755                     + t4ddg_2(2, 1, 1, 0) * t4ddg_3(1, 0, 1, 0)
  756                     + t4ddg_2(2, 1, 1, 1) * t4ddg_3(1, 1, 1, 0)
  757                     + t4ddg_2(2, 1, 1, 2) * t4ddg_3(1, 2, 1, 0)
  758                     + t4ddg_2(2, 2, 1, 0) * t4ddg_3(2, 0, 1, 0)
  759                     + t4ddg_2(2, 2, 1, 1) * t4ddg_3(2, 1, 1, 0)
  760                     + t4ddg_2(2, 2, 1, 2) * t4ddg_3(2, 2, 1, 0)),
  761                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,1,1,0)");
  763                  - (t4ddg_2(2, 0, 1, 0) * t4ddg_3(0, 0, 1, 1)
  764                     + t4ddg_2(2, 0, 1, 1) * t4ddg_3(0, 1, 1, 1)
  765                     + t4ddg_2(2, 0, 1, 2) * t4ddg_3(0, 2, 1, 1)
  766                     + t4ddg_2(2, 1, 1, 0) * t4ddg_3(1, 0, 1, 1)
  767                     + t4ddg_2(2, 1, 1, 1) * t4ddg_3(1, 1, 1, 1)
  768                     + t4ddg_2(2, 1, 1, 2) * t4ddg_3(1, 2, 1, 1)
  769                     + t4ddg_2(2, 2, 1, 0) * t4ddg_3(2, 0, 1, 1)
  770                     + t4ddg_2(2, 2, 1, 1) * t4ddg_3(2, 1, 1, 1)
  771                     + t4ddg_2(2, 2, 1, 2) * t4ddg_3(2, 2, 1, 1)),
  772                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,1,1,1)");
  774                  - (t4ddg_2(2, 0, 1, 0) * t4ddg_3(0, 0, 1, 2)
  775                     + t4ddg_2(2, 0, 1, 1) * t4ddg_3(0, 1, 1, 2)
  776                     + t4ddg_2(2, 0, 1, 2) * t4ddg_3(0, 2, 1, 2)
  777                     + t4ddg_2(2, 1, 1, 0) * t4ddg_3(1, 0, 1, 2)
  778                     + t4ddg_2(2, 1, 1, 1) * t4ddg_3(1, 1, 1, 2)
  779                     + t4ddg_2(2, 1, 1, 2) * t4ddg_3(1, 2, 1, 2)
  780                     + t4ddg_2(2, 2, 1, 0) * t4ddg_3(2, 0, 1, 2)
  781                     + t4ddg_2(2, 2, 1, 1) * t4ddg_3(2, 1, 1, 2)
  782                     + t4ddg_2(2, 2, 1, 2) * t4ddg_3(2, 2, 1, 2)),
  783                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,1,1,2)");
  785                  - (t4ddg_2(2, 0, 2, 0) * t4ddg_3(0, 0, 1, 0)
  786                     + t4ddg_2(2, 0, 2, 1) * t4ddg_3(0, 1, 1, 0)
  787                     + t4ddg_2(2, 0, 2, 2) * t4ddg_3(0, 2, 1, 0)
  788                     + t4ddg_2(2, 1, 2, 0) * t4ddg_3(1, 0, 1, 0)
  789                     + t4ddg_2(2, 1, 2, 1) * t4ddg_3(1, 1, 1, 0)
  790                     + t4ddg_2(2, 1, 2, 2) * t4ddg_3(1, 2, 1, 0)
  791                     + t4ddg_2(2, 2, 2, 0) * t4ddg_3(2, 0, 1, 0)
  792                     + t4ddg_2(2, 2, 2, 1) * t4ddg_3(2, 1, 1, 0)
  793                     + t4ddg_2(2, 2, 2, 2) * t4ddg_3(2, 2, 1, 0)),
  794                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,1,2,0)");
  796                  - (t4ddg_2(2, 0, 2, 0) * t4ddg_3(0, 0, 1, 1)
  797                     + t4ddg_2(2, 0, 2, 1) * t4ddg_3(0, 1, 1, 1)
  798                     + t4ddg_2(2, 0, 2, 2) * t4ddg_3(0, 2, 1, 1)
  799                     + t4ddg_2(2, 1, 2, 0) * t4ddg_3(1, 0, 1, 1)
  800                     + t4ddg_2(2, 1, 2, 1) * t4ddg_3(1, 1, 1, 1)
  801                     + t4ddg_2(2, 1, 2, 2) * t4ddg_3(1, 2, 1, 1)
  802                     + t4ddg_2(2, 2, 2, 0) * t4ddg_3(2, 0, 1, 1)
  803                     + t4ddg_2(2, 2, 2, 1) * t4ddg_3(2, 1, 1, 1)
  804                     + t4ddg_2(2, 2, 2, 2) * t4ddg_3(2, 2, 1, 1)),
  805                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,1,2,1)");
  807                  - (t4ddg_2(2, 0, 2, 0) * t4ddg_3(0, 0, 1, 2)
  808                     + t4ddg_2(2, 0, 2, 1) * t4ddg_3(0, 1, 1, 2)
  809                     + t4ddg_2(2, 0, 2, 2) * t4ddg_3(0, 2, 1, 2)
  810                     + t4ddg_2(2, 1, 2, 0) * t4ddg_3(1, 0, 1, 2)
  811                     + t4ddg_2(2, 1, 2, 1) * t4ddg_3(1, 1, 1, 2)
  812                     + t4ddg_2(2, 1, 2, 2) * t4ddg_3(1, 2, 1, 2)
  813                     + t4ddg_2(2, 2, 2, 0) * t4ddg_3(2, 0, 1, 2)
  814                     + t4ddg_2(2, 2, 2, 1) * t4ddg_3(2, 1, 1, 2)
  815                     + t4ddg_2(2, 2, 2, 2) * t4ddg_3(2, 2, 1, 2)),
  816                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,1,2,2)");
  818                  - (t4ddg_2(2, 0, 0, 0) * t4ddg_3(0, 0, 2, 0)
  819                     + t4ddg_2(2, 0, 0, 1) * t4ddg_3(0, 1, 2, 0)
  820                     + t4ddg_2(2, 0, 0, 2) * t4ddg_3(0, 2, 2, 0)
  821                     + t4ddg_2(2, 1, 0, 0) * t4ddg_3(1, 0, 2, 0)
  822                     + t4ddg_2(2, 1, 0, 1) * t4ddg_3(1, 1, 2, 0)
  823                     + t4ddg_2(2, 1, 0, 2) * t4ddg_3(1, 2, 2, 0)
  824                     + t4ddg_2(2, 2, 0, 0) * t4ddg_3(2, 0, 2, 0)
  825                     + t4ddg_2(2, 2, 0, 1) * t4ddg_3(2, 1, 2, 0)
  826                     + t4ddg_2(2, 2, 0, 2) * t4ddg_3(2, 2, 2, 0)),
  827                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,2,0,0)");
  829                  - (t4ddg_2(2, 0, 0, 0) * t4ddg_3(0, 0, 2, 1)
  830                     + t4ddg_2(2, 0, 0, 1) * t4ddg_3(0, 1, 2, 1)
  831                     + t4ddg_2(2, 0, 0, 2) * t4ddg_3(0, 2, 2, 1)
  832                     + t4ddg_2(2, 1, 0, 0) * t4ddg_3(1, 0, 2, 1)
  833                     + t4ddg_2(2, 1, 0, 1) * t4ddg_3(1, 1, 2, 1)
  834                     + t4ddg_2(2, 1, 0, 2) * t4ddg_3(1, 2, 2, 1)
  835                     + t4ddg_2(2, 2, 0, 0) * t4ddg_3(2, 0, 2, 1)
  836                     + t4ddg_2(2, 2, 0, 1) * t4ddg_3(2, 1, 2, 1)
  837                     + t4ddg_2(2, 2, 0, 2) * t4ddg_3(2, 2, 2, 1)),
  838                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,2,0,1)");
  840                  - (t4ddg_2(2, 0, 0, 0) * t4ddg_3(0, 0, 2, 2)
  841                     + t4ddg_2(2, 0, 0, 1) * t4ddg_3(0, 1, 2, 2)
  842                     + t4ddg_2(2, 0, 0, 2) * t4ddg_3(0, 2, 2, 2)
  843                     + t4ddg_2(2, 1, 0, 0) * t4ddg_3(1, 0, 2, 2)
  844                     + t4ddg_2(2, 1, 0, 1) * t4ddg_3(1, 1, 2, 2)
  845                     + t4ddg_2(2, 1, 0, 2) * t4ddg_3(1, 2, 2, 2)
  846                     + t4ddg_2(2, 2, 0, 0) * t4ddg_3(2, 0, 2, 2)
  847                     + t4ddg_2(2, 2, 0, 1) * t4ddg_3(2, 1, 2, 2)
  848                     + t4ddg_2(2, 2, 0, 2) * t4ddg_3(2, 2, 2, 2)),
  849                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,2,0,2)");
  851                  - (t4ddg_2(2, 0, 1, 0) * t4ddg_3(0, 0, 2, 0)
  852                     + t4ddg_2(2, 0, 1, 1) * t4ddg_3(0, 1, 2, 0)
  853                     + t4ddg_2(2, 0, 1, 2) * t4ddg_3(0, 2, 2, 0)
  854                     + t4ddg_2(2, 1, 1, 0) * t4ddg_3(1, 0, 2, 0)
  855                     + t4ddg_2(2, 1, 1, 1) * t4ddg_3(1, 1, 2, 0)
  856                     + t4ddg_2(2, 1, 1, 2) * t4ddg_3(1, 2, 2, 0)
  857                     + t4ddg_2(2, 2, 1, 0) * t4ddg_3(2, 0, 2, 0)
  858                     + t4ddg_2(2, 2, 1, 1) * t4ddg_3(2, 1, 2, 0)
  859                     + t4ddg_2(2, 2, 1, 2) * t4ddg_3(2, 2, 2, 0)),
  860                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,2,1,0)");
  862                  - (t4ddg_2(2, 0, 1, 0) * t4ddg_3(0, 0, 2, 1)
  863                     + t4ddg_2(2, 0, 1, 1) * t4ddg_3(0, 1, 2, 1)
  864                     + t4ddg_2(2, 0, 1, 2) * t4ddg_3(0, 2, 2, 1)
  865                     + t4ddg_2(2, 1, 1, 0) * t4ddg_3(1, 0, 2, 1)
  866                     + t4ddg_2(2, 1, 1, 1) * t4ddg_3(1, 1, 2, 1)
  867                     + t4ddg_2(2, 1, 1, 2) * t4ddg_3(1, 2, 2, 1)
  868                     + t4ddg_2(2, 2, 1, 0) * t4ddg_3(2, 0, 2, 1)
  869                     + t4ddg_2(2, 2, 1, 1) * t4ddg_3(2, 1, 2, 1)
  870                     + t4ddg_2(2, 2, 1, 2) * t4ddg_3(2, 2, 2, 1)),
  871                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,2,1,1)");
  873                  - (t4ddg_2(2, 0, 1, 0) * t4ddg_3(0, 0, 2, 2)
  874                     + t4ddg_2(2, 0, 1, 1) * t4ddg_3(0, 1, 2, 2)
  875                     + t4ddg_2(2, 0, 1, 2) * t4ddg_3(0, 2, 2, 2)
  876                     + t4ddg_2(2, 1, 1, 0) * t4ddg_3(1, 0, 2, 2)
  877                     + t4ddg_2(2, 1, 1, 1) * t4ddg_3(1, 1, 2, 2)
  878                     + t4ddg_2(2, 1, 1, 2) * t4ddg_3(1, 2, 2, 2)
  879                     + t4ddg_2(2, 2, 1, 0) * t4ddg_3(2, 0, 2, 2)
  880                     + t4ddg_2(2, 2, 1, 1) * t4ddg_3(2, 1, 2, 2)
  881                     + t4ddg_2(2, 2, 1, 2) * t4ddg_3(2, 2, 2, 2)),
  882                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,2,1,2)");
  884                  - (t4ddg_2(2, 0, 2, 0) * t4ddg_3(0, 0, 2, 0)
  885                     + t4ddg_2(2, 0, 2, 1) * t4ddg_3(0, 1, 2, 0)
  886                     + t4ddg_2(2, 0, 2, 2) * t4ddg_3(0, 2, 2, 0)
  887                     + t4ddg_2(2, 1, 2, 0) * t4ddg_3(1, 0, 2, 0)
  888                     + t4ddg_2(2, 1, 2, 1) * t4ddg_3(1, 1, 2, 0)
  889                     + t4ddg_2(2, 1, 2, 2) * t4ddg_3(1, 2, 2, 0)
  890                     + t4ddg_2(2, 2, 2, 0) * t4ddg_3(2, 0, 2, 0)
  891                     + t4ddg_2(2, 2, 2, 1) * t4ddg_3(2, 1, 2, 0)
  892                     + t4ddg_2(2, 2, 2, 2) * t4ddg_3(2, 2, 2, 0)),
  893                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,2,2,0)");
  895                  - (t4ddg_2(2, 0, 2, 0) * t4ddg_3(0, 0, 2, 1)
  896                     + t4ddg_2(2, 0, 2, 1) * t4ddg_3(0, 1, 2, 1)
  897                     + t4ddg_2(2, 0, 2, 2) * t4ddg_3(0, 2, 2, 1)
  898                     + t4ddg_2(2, 1, 2, 0) * t4ddg_3(1, 0, 2, 1)
  899                     + t4ddg_2(2, 1, 2, 1) * t4ddg_3(1, 1, 2, 1)
  900                     + t4ddg_2(2, 1, 2, 2) * t4ddg_3(1, 2, 2, 1)
  901                     + t4ddg_2(2, 2, 2, 0) * t4ddg_3(2, 0, 2, 1)
  902                     + t4ddg_2(2, 2, 2, 1) * t4ddg_3(2, 1, 2, 1)
  903                     + t4ddg_2(2, 2, 2, 2) * t4ddg_3(2, 2, 2, 1)),
  904                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,2,2,1)");
  906                  - (t4ddg_2(2, 0, 2, 0) * t4ddg_3(0, 0, 2, 2)
  907                     + t4ddg_2(2, 0, 2, 1) * t4ddg_3(0, 1, 2, 2)
  908                     + t4ddg_2(2, 0, 2, 2) * t4ddg_3(0, 2, 2, 2)
  909                     + t4ddg_2(2, 1, 2, 0) * t4ddg_3(1, 0, 2, 2)
  910                     + t4ddg_2(2, 1, 2, 1) * t4ddg_3(1, 1, 2, 2)
  911                     + t4ddg_2(2, 1, 2, 2) * t4ddg_3(1, 2, 2, 2)
  912                     + t4ddg_2(2, 2, 2, 0) * t4ddg_3(2, 0, 2, 2)
  913                     + t4ddg_2(2, 2, 2, 1) * t4ddg_3(2, 1, 2, 2)
  914                     + t4ddg_2(2, 2, 2, 2) * t4ddg_3(2, 2, 2, 2)),
  915                "T4ddg(i,j,k,l)^T4ddg(j,l,m,n)(2,2,2,2)");
  916}
FTensor::Index< 'i', SPACE_DIM > i
const double n
refractive index of diffusive medium
FTensor::Index< 'l', 3 > l
FTensor::Index< 'j', 3 > j
FTensor::Index< 'k', 3 > k
FTensor::Index< 'm', 3 > m
void test_for_zero(const T &t, const std::string &s)