13{
   18 
   22 
   23  
   24 
   26                  - t2_2(0, 0) * t1_2(0) * 10,
   27                "T3*T(0,0,0)");
   29                  - t2_2(0, 0) * t1_2(1) * 10,
   30                "T3*T(0,0,1)");
   32                  - t2_2(0, 0) * t1_2(2) * 10,
   33                "T3*T(0,0,2)");
   35                  - t2_2(0, 1) * t1_2(0) * 10,
   36                "T3*T(0,1,0)");
   38                  - t2_2(0, 1) * t1_2(1) * 10,
   39                "T3*T(0,1,1)");
   41                  - t2_2(0, 1) * t1_2(2) * 10,
   42                "T3*T(0,1,2)");
   44                  - t2_2(0, 2) * t1_2(0) * 10,
   45                "T3*T(0,2,0)");
   47                  - t2_2(0, 2) * t1_2(1) * 10,
   48                "T3*T(0,2,1)");
   50                  - t2_2(0, 2) * t1_2(2) * 10,
   51                "T3*T(0,2,2)");
   53                  - t2_2(1, 0) * t1_2(0) * 10,
   54                "T3*T(1,0,0)");
   56                  - t2_2(1, 0) * t1_2(1) * 10,
   57                "T3*T(1,0,1)");
   59                  - t2_2(1, 0) * t1_2(2) * 10,
   60                "T3*T(1,0,2)");
   62                  - t2_2(1, 1) * t1_2(0) * 10,
   63                "T3*T(1,1,0)");
   65                  - t2_2(1, 1) * t1_2(1) * 10,
   66                "T3*T(1,1,1)");
   68                  - t2_2(1, 1) * t1_2(2) * 10,
   69                "T3*T(1,1,2)");
   71                  - t2_2(1, 2) * t1_2(0) * 10,
   72                "T3*T(1,2,0)");
   74                  - t2_2(1, 2) * t1_2(1) * 10,
   75                "T3*T(1,2,1)");
   77                  - t2_2(1, 2) * t1_2(2) * 10,
   78                "T3*T(1,2,2)");
   80                  - t2_2(2, 0) * t1_2(0) * 10,
   81                "T3*T(2,0,0)");
   83                  - t2_2(2, 0) * t1_2(1) * 10,
   84                "T3*T(2,0,1)");
   86                  - t2_2(2, 0) * t1_2(2) * 10,
   87                "T3*T(2,0,2)");
   89                  - t2_2(2, 1) * t1_2(0) * 10,
   90                "T3*T(2,1,0)");
   92                  - t2_2(2, 1) * t1_2(1) * 10,
   93                "T3*T(2,1,1)");
   95                  - t2_2(2, 1) * t1_2(2) * 10,
   96                "T3*T(2,1,2)");
   98                  - t2_2(2, 2) * t1_2(0) * 10,
   99                "T3*T(2,2,0)");
  101                  - t2_2(2, 2) * t1_2(1) * 10,
  102                "T3*T(2,2,1)");
  104                  - t2_2(2, 2) * t1_2(2) * 10,
  105                "T3*T(2,2,2)");
  106 
  108                  - t2_2(0, 0) * t1_2(0) * 10,
  109                "T*T3(0,0,0)");
  111                  - t2_2(0, 0) * t1_2(1) * 10,
  112                "T*T3(0,0,1)");
  114                  - t2_2(0, 0) * t1_2(2) * 10,
  115                "T*T3(0,0,2)");
  117                  - t2_2(0, 1) * t1_2(0) * 10,
  118                "T*T3(0,1,0)");
  120                  - t2_2(0, 1) * t1_2(1) * 10,
  121                "T*T3(0,1,1)");
  123                  - t2_2(0, 1) * t1_2(2) * 10,
  124                "T*T3(0,1,2)");
  126                  - t2_2(0, 2) * t1_2(0) * 10,
  127                "T*T3(0,2,0)");
  129                  - t2_2(0, 2) * t1_2(1) * 10,
  130                "T*T3(0,2,1)");
  132                  - t2_2(0, 2) * t1_2(2) * 10,
  133                "T*T3(0,2,2)");
  135                  - t2_2(1, 0) * t1_2(0) * 10,
  136                "T*T3(1,0,0)");
  138                  - t2_2(1, 0) * t1_2(1) * 10,
  139                "T*T3(1,0,1)");
  141                  - t2_2(1, 0) * t1_2(2) * 10,
  142                "T*T3(1,0,2)");
  144                  - t2_2(1, 1) * t1_2(0) * 10,
  145                "T*T3(1,1,0)");
  147                  - t2_2(1, 1) * t1_2(1) * 10,
  148                "T*T3(1,1,1)");
  150                  - t2_2(1, 1) * t1_2(2) * 10,
  151                "T*T3(1,1,2)");
  153                  - t2_2(1, 2) * t1_2(0) * 10,
  154                "T*T3(1,2,0)");
  156                  - t2_2(1, 2) * t1_2(1) * 10,
  157                "T*T3(1,2,1)");
  159                  - t2_2(1, 2) * t1_2(2) * 10,
  160                "T*T3(1,2,2)");
  162                  - t2_2(2, 0) * t1_2(0) * 10,
  163                "T*T3(2,0,0)");
  165                  - t2_2(2, 0) * t1_2(1) * 10,
  166                "T*T3(2,0,1)");
  168                  - t2_2(2, 0) * t1_2(2) * 10,
  169                "T*T3(2,0,2)");
  171                  - t2_2(2, 1) * t1_2(0) * 10,
  172                "T*T3(2,1,0)");
  174                  - t2_2(2, 1) * t1_2(1) * 10,
  175                "T*T3(2,1,1)");
  177                  - t2_2(2, 1) * t1_2(2) * 10,
  178                "T*T3(2,1,2)");
  180                  - t2_2(2, 2) * t1_2(0) * 10,
  181                "T*T3(2,2,0)");
  183                  - t2_2(2, 2) * t1_2(1) * 10,
  184                "T*T3(2,2,1)");
  186                  - t2_2(2, 2) * t1_2(2) * 10,
  187                "T*T3(2,2,2)");
  188 
  190 
  191  t2(
i, 
k) = (t2_2(
i, 
j) * t1_2(
k)) * t1_1(
j);
 
  192 
  194                  - (t2_2(0, 0) * t1_2(0) * t1_1(0)
  195                     + t2_2(0, 1) * t1_2(0) * t1_1(1)
  196                     + t2_2(0, 2) * t1_2(0) * t1_1(2)),
  197                "T3(i,j,k)*T1(j)(0,0)");
  199                  - (t2_2(0, 0) * t1_2(1) * t1_1(0)
  200                     + t2_2(0, 1) * t1_2(1) * t1_1(1)
  201                     + t2_2(0, 2) * t1_2(1) * t1_1(2)),
  202                "T3(i,j,k)*T1(j)(0,1)");
  204                  - (t2_2(0, 0) * t1_2(2) * t1_1(0)
  205                     + t2_2(0, 1) * t1_2(2) * t1_1(1)
  206                     + t2_2(0, 2) * t1_2(2) * t1_1(2)),
  207                "T3(i,j,k)*T1(j)(0,2)");
  209                  - (t2_2(1, 0) * t1_2(0) * t1_1(0)
  210                     + t2_2(1, 1) * t1_2(0) * t1_1(1)
  211                     + t2_2(1, 2) * t1_2(0) * t1_1(2)),
  212                "T3(i,j,k)*T1(j)(1,0)");
  214                  - (t2_2(1, 0) * t1_2(1) * t1_1(0)
  215                     + t2_2(1, 1) * t1_2(1) * t1_1(1)
  216                     + t2_2(1, 2) * t1_2(1) * t1_1(2)),
  217                "T3(i,j,k)*T1(j)(1,1)");
  219                  - (t2_2(1, 0) * t1_2(2) * t1_1(0)
  220                     + t2_2(1, 1) * t1_2(2) * t1_1(1)
  221                     + t2_2(1, 2) * t1_2(2) * t1_1(2)),
  222                "T3(i,j,k)*T1(j)(1,2)");
  224                  - (t2_2(2, 0) * t1_2(0) * t1_1(0)
  225                     + t2_2(2, 1) * t1_2(0) * t1_1(1)
  226                     + t2_2(2, 2) * t1_2(0) * t1_1(2)),
  227                "T3(i,j,k)*T1(j)(2,0)");
  229                  - (t2_2(2, 0) * t1_2(1) * t1_1(0)
  230                     + t2_2(2, 1) * t1_2(1) * t1_1(1)
  231                     + t2_2(2, 2) * t1_2(1) * t1_1(2)),
  232                "T3(i,j,k)*T1(j)(2,1)");
  234                  - (t2_2(2, 0) * t1_2(2) * t1_1(0)
  235                     + t2_2(2, 1) * t1_2(2) * t1_1(1)
  236                     + t2_2(2, 2) * t1_2(2) * t1_1(2)),
  237                "T3(i,j,k)*T1(j)(2,2)");
  238 
  239  t2(
i, 
k) = t1_2(
j) * (t2_2(
i, 
j) * t1_1(
k));
 
  240 
  242                  - (t2_2(0, 0) * t1_1(0) * t1_2(0)
  243                     + t2_2(0, 1) * t1_1(0) * t1_2(1)
  244                     + t2_2(0, 2) * t1_1(0) * t1_2(2)),
  245                "T1(j)*T3(i,j,k)(0,0)");
  247                  - (t2_2(0, 0) * t1_1(1) * t1_2(0)
  248                     + t2_2(0, 1) * t1_1(1) * t1_2(1)
  249                     + t2_2(0, 2) * t1_1(1) * t1_2(2)),
  250                "T1(j)*T3(i,j,k)(0,1)");
  252                  - (t2_2(0, 0) * t1_1(2) * t1_2(0)
  253                     + t2_2(0, 1) * t1_1(2) * t1_2(1)
  254                     + t2_2(0, 2) * t1_1(2) * t1_2(2)),
  255                "T1(j)*T3(i,j,k)(0,2)");
  257                  - (t2_2(1, 0) * t1_1(0) * t1_2(0)
  258                     + t2_2(1, 1) * t1_1(0) * t1_2(1)
  259                     + t2_2(1, 2) * t1_1(0) * t1_2(2)),
  260                "T1(j)*T3(i,j,k)(1,0)");
  262                  - (t2_2(1, 0) * t1_1(1) * t1_2(0)
  263                     + t2_2(1, 1) * t1_1(1) * t1_2(1)
  264                     + t2_2(1, 2) * t1_1(1) * t1_2(2)),
  265                "T1(j)*T3(i,j,k)(1,1)");
  267                  - (t2_2(1, 0) * t1_1(2) * t1_2(0)
  268                     + t2_2(1, 1) * t1_1(2) * t1_2(1)
  269                     + t2_2(1, 2) * t1_1(2) * t1_2(2)),
  270                "T1(j)*T3(i,j,k)(1,2)");
  272                  - (t2_2(2, 0) * t1_1(0) * t1_2(0)
  273                     + t2_2(2, 1) * t1_1(0) * t1_2(1)
  274                     + t2_2(2, 2) * t1_1(0) * t1_2(2)),
  275                "T1(j)*T3(i,j,k)(2,0)");
  277                  - (t2_2(2, 0) * t1_1(1) * t1_2(0)
  278                     + t2_2(2, 1) * t1_1(1) * t1_2(1)
  279                     + t2_2(2, 2) * t1_1(1) * t1_2(2)),
  280                "T1(j)*T3(i,j,k)(2,1)");
  282                  - (t2_2(2, 0) * t1_1(2) * t1_2(0)
  283                     + t2_2(2, 1) * t1_1(2) * t1_2(1)
  284                     + t2_2(2, 2) * t1_1(2) * t1_2(2)),
  285                "T1(j)*T3(i,j,k)(2,2)");
  286 
  287  
  289  for(int ii = 0; ii < 3; ++ii)
  290    {
  292                      - (t2_2(ii, 0) * t1_2(0) * t2s_2(0, 0)
  293                         + t2_2(ii, 0) * t1_2(1) * t2s_2(0, 1)
  294                         + t2_2(ii, 0) * t1_2(2) * t2s_2(0, 2)
  295                         + t2_2(ii, 1) * t1_2(0) * t2s_2(1, 0)
  296                         + t2_2(ii, 1) * t1_2(1) * t2s_2(1, 1)
  297                         + t2_2(ii, 1) * t1_2(2) * t2s_2(1, 2)
  298                         + t2_2(ii, 2) * t1_2(0) * t2s_2(2, 0)
  299                         + t2_2(ii, 2) * t1_2(1) * t2s_2(2, 1)
  300                         + t2_2(ii, 2) * t1_2(2) * t2s_2(2, 2)),
  301                    "T3(i,j,k)*T2s(j,k)(" + std::to_string(ii) + ")");
  302 
  304                      - (t2_2(ii, 0) * t1_2(0) * t2s_3(0, 0)
  305                         + t2_2(ii, 0) * t1_2(1) * t2s_3(0, 1)
  306                         + t2_2(ii, 0) * t1_2(2) * t2s_3(0, 2)
  307                         + t2_2(ii, 1) * t1_2(0) * t2s_3(1, 0)
  308                         + t2_2(ii, 1) * t1_2(1) * t2s_3(1, 1)
  309                         + t2_2(ii, 1) * t1_2(2) * t2s_3(1, 2)
  310                         + t2_2(ii, 2) * t1_2(0) * t2s_3(2, 0)
  311                         + t2_2(ii, 2) * t1_2(1) * t2s_3(2, 1)
  312                         + t2_2(ii, 2) * t1_2(2) * t2s_3(2, 2)),
  313                    "T2s(j,k)*T3(i,j,k)(" + std::to_string(ii) + ")");
  314 
  316                      - (t2_2(0, ii) * t1_2(0) * t2s_2(0, 0)
  317                         + t2_2(0, ii) * t1_2(1) * t2s_2(0, 1)
  318                         + t2_2(0, ii) * t1_2(2) * t2s_2(0, 2)
  319                         + t2_2(1, ii) * t1_2(0) * t2s_2(1, 0)
  320                         + t2_2(1, ii) * t1_2(1) * t2s_2(1, 1)
  321                         + t2_2(1, ii) * t1_2(2) * t2s_2(1, 2)
  322                         + t2_2(2, ii) * t1_2(0) * t2s_2(2, 0)
  323                         + t2_2(2, ii) * t1_2(1) * t2s_2(2, 1)
  324                         + t2_2(2, ii) * t1_2(2) * t2s_2(2, 2)),
  325                    "T3(j,i,k)*T2s(j,k)(" + std::to_string(ii) + ")");
  326 
  328                      - (t2_2(0, ii) * t1_2(0) * t2s_3(0, 0)
  329                         + t2_2(0, ii) * t1_2(1) * t2s_3(0, 1)
  330                         + t2_2(0, ii) * t1_2(2) * t2s_3(0, 2)
  331                         + t2_2(1, ii) * t1_2(0) * t2s_3(1, 0)
  332                         + t2_2(1, ii) * t1_2(1) * t2s_3(1, 1)
  333                         + t2_2(1, ii) * t1_2(2) * t2s_3(1, 2)
  334                         + t2_2(2, ii) * t1_2(0) * t2s_3(2, 0)
  335                         + t2_2(2, ii) * t1_2(1) * t2s_3(2, 1)
  336                         + t2_2(2, ii) * t1_2(2) * t2s_3(2, 2)),
  337                    "T2s(j,k)*T3(j,i,k)(" + std::to_string(ii) + ")");
  338 
  340                      - (t2_2(0, 0) * t1_2(ii) * t2s_2(0, 0)
  341                         + t2_2(0, 1) * t1_2(ii) * t2s_2(0, 1)
  342                         + t2_2(0, 2) * t1_2(ii) * t2s_2(0, 2)
  343                         + t2_2(1, 0) * t1_2(ii) * t2s_2(1, 0)
  344                         + t2_2(1, 1) * t1_2(ii) * t2s_2(1, 1)
  345                         + t2_2(1, 2) * t1_2(ii) * t2s_2(1, 2)
  346                         + t2_2(2, 0) * t1_2(ii) * t2s_2(2, 0)
  347                         + t2_2(2, 1) * t1_2(ii) * t2s_2(2, 1)
  348                         + t2_2(2, 2) * t1_2(ii) * t2s_2(2, 2)),
  349                    "T3(j,k,i)*T2s(j,k)(" + std::to_string(ii) + ")");
  350 
  352                      - (t2_2(0, 0) * t1_2(ii) * t2s_3(0, 0)
  353                         + t2_2(0, 1) * t1_2(ii) * t2s_3(0, 1)
  354                         + t2_2(0, 2) * t1_2(ii) * t2s_3(0, 2)
  355                         + t2_2(1, 0) * t1_2(ii) * t2s_3(1, 0)
  356                         + t2_2(1, 1) * t1_2(ii) * t2s_3(1, 1)
  357                         + t2_2(1, 2) * t1_2(ii) * t2s_3(1, 2)
  358                         + t2_2(2, 0) * t1_2(ii) * t2s_3(2, 0)
  359                         + t2_2(2, 1) * t1_2(ii) * t2s_3(2, 1)
  360                         + t2_2(2, 2) * t1_2(ii) * t2s_3(2, 2)),
  361                    "T2s(j,k)*T3(j,k,i)(" + std::to_string(ii) + ")");
  362    }
  363 
  364  
  365  for(int ii = 0; ii < 3; ++ii)
  366    {
  368                      - (t2_2(ii, 0) * t1_2(0) * t2_2(0, 0)
  369                         + t2_2(ii, 0) * t1_2(1) * t2_2(0, 1)
  370                         + t2_2(ii, 0) * t1_2(2) * t2_2(0, 2)
  371                         + t2_2(ii, 1) * t1_2(0) * t2_2(1, 0)
  372                         + t2_2(ii, 1) * t1_2(1) * t2_2(1, 1)
  373                         + t2_2(ii, 1) * t1_2(2) * t2_2(1, 2)
  374                         + t2_2(ii, 2) * t1_2(0) * t2_2(2, 0)
  375                         + t2_2(ii, 2) * t1_2(1) * t2_2(2, 1)
  376                         + t2_2(ii, 2) * t1_2(2) * t2_2(2, 2)),
  377                    "T3(i,j,k)*T2(j,k)(" + std::to_string(ii) + ")");
  378 
  380                      - (t2_2(ii, 0) * t1_2(0) * t2_2(0, 0)
  381                         + t2_2(ii, 0) * t1_2(1) * t2_2(0, 1)
  382                         + t2_2(ii, 0) * t1_2(2) * t2_2(0, 2)
  383                         + t2_2(ii, 1) * t1_2(0) * t2_2(1, 0)
  384                         + t2_2(ii, 1) * t1_2(1) * t2_2(1, 1)
  385                         + t2_2(ii, 1) * t1_2(2) * t2_2(1, 2)
  386                         + t2_2(ii, 2) * t1_2(0) * t2_2(2, 0)
  387                         + t2_2(ii, 2) * t1_2(1) * t2_2(2, 1)
  388                         + t2_2(ii, 2) * t1_2(2) * t2_2(2, 2)),
  389                    "T2(j,k)*T3(i,j,k)(" + std::to_string(ii) + ")");
  390 
  392                      - (t2_2(ii, 0) * t1_2(0) * t2_2(0, 0)
  393                         + t2_2(ii, 0) * t1_2(1) * t2_2(1, 0)
  394                         + t2_2(ii, 0) * t1_2(2) * t2_2(2, 0)
  395                         + t2_2(ii, 1) * t1_2(0) * t2_2(0, 1)
  396                         + t2_2(ii, 1) * t1_2(1) * t2_2(1, 1)
  397                         + t2_2(ii, 1) * t1_2(2) * t2_2(2, 1)
  398                         + t2_2(ii, 2) * t1_2(0) * t2_2(0, 2)
  399                         + t2_2(ii, 2) * t1_2(1) * t2_2(1, 2)
  400                         + t2_2(ii, 2) * t1_2(2) * t2_2(2, 2)),
  401                    "T3(i,j,k)*T2(k,j)(" + std::to_string(ii) + ")");
  402 
  404                      - (t2_2(ii, 0) * t1_2(0) * t2_2(0, 0)
  405                         + t2_2(ii, 0) * t1_2(1) * t2_2(1, 0)
  406                         + t2_2(ii, 0) * t1_2(2) * t2_2(2, 0)
  407                         + t2_2(ii, 1) * t1_2(0) * t2_2(0, 1)
  408                         + t2_2(ii, 1) * t1_2(1) * t2_2(1, 1)
  409                         + t2_2(ii, 1) * t1_2(2) * t2_2(2, 1)
  410                         + t2_2(ii, 2) * t1_2(0) * t2_2(0, 2)
  411                         + t2_2(ii, 2) * t1_2(1) * t2_2(1, 2)
  412                         + t2_2(ii, 2) * t1_2(2) * t2_2(2, 2)),
  413                    "T2(k,j)*T3(i,j,k)(" + std::to_string(ii) + ")");
  414 
  416                      - (t2_2(0, ii) * t1_2(0) * t2_2(0, 0)
  417                         + t2_2(0, ii) * t1_2(1) * t2_2(0, 1)
  418                         + t2_2(0, ii) * t1_2(2) * t2_2(0, 2)
  419                         + t2_2(1, ii) * t1_2(0) * t2_2(1, 0)
  420                         + t2_2(1, ii) * t1_2(1) * t2_2(1, 1)
  421                         + t2_2(1, ii) * t1_2(2) * t2_2(1, 2)
  422                         + t2_2(2, ii) * t1_2(0) * t2_2(2, 0)
  423                         + t2_2(2, ii) * t1_2(1) * t2_2(2, 1)
  424                         + t2_2(2, ii) * t1_2(2) * t2_2(2, 2)),
  425                    "T3(j,i,k)*T2(j,k)(" + std::to_string(ii) + ")");
  426 
  428                      - (t2_2(0, ii) * t1_2(0) * t2_2(0, 0)
  429                         + t2_2(0, ii) * t1_2(1) * t2_2(0, 1)
  430                         + t2_2(0, ii) * t1_2(2) * t2_2(0, 2)
  431                         + t2_2(1, ii) * t1_2(0) * t2_2(1, 0)
  432                         + t2_2(1, ii) * t1_2(1) * t2_2(1, 1)
  433                         + t2_2(1, ii) * t1_2(2) * t2_2(1, 2)
  434                         + t2_2(2, ii) * t1_2(0) * t2_2(2, 0)
  435                         + t2_2(2, ii) * t1_2(1) * t2_2(2, 1)
  436                         + t2_2(2, ii) * t1_2(2) * t2_2(2, 2)),
  437                    "T2(j,k)*T3(j,i,k)(" + std::to_string(ii) + ")");
  438 
  440                      - (t2_2(0, ii) * t1_2(0) * t2_2(0, 0)
  441                         + t2_2(0, ii) * t1_2(1) * t2_2(1, 0)
  442                         + t2_2(0, ii) * t1_2(2) * t2_2(2, 0)
  443                         + t2_2(1, ii) * t1_2(0) * t2_2(0, 1)
  444                         + t2_2(1, ii) * t1_2(1) * t2_2(1, 1)
  445                         + t2_2(1, ii) * t1_2(2) * t2_2(2, 1)
  446                         + t2_2(2, ii) * t1_2(0) * t2_2(0, 2)
  447                         + t2_2(2, ii) * t1_2(1) * t2_2(1, 2)
  448                         + t2_2(2, ii) * t1_2(2) * t2_2(2, 2)),
  449                    "T3(j,i,k)*T2(k,j)(" + std::to_string(ii) + ")");
  450 
  452                      - (t2_2(0, ii) * t1_2(0) * t2_2(0, 0)
  453                         + t2_2(0, ii) * t1_2(1) * t2_2(1, 0)
  454                         + t2_2(0, ii) * t1_2(2) * t2_2(2, 0)
  455                         + t2_2(1, ii) * t1_2(0) * t2_2(0, 1)
  456                         + t2_2(1, ii) * t1_2(1) * t2_2(1, 1)
  457                         + t2_2(1, ii) * t1_2(2) * t2_2(2, 1)
  458                         + t2_2(2, ii) * t1_2(0) * t2_2(0, 2)
  459                         + t2_2(2, ii) * t1_2(1) * t2_2(1, 2)
  460                         + t2_2(2, ii) * t1_2(2) * t2_2(2, 2)),
  461                    "T2(k,j)*T3(j,i,k)(" + std::to_string(ii) + ")");
  462 
  464                      - (t2_2(0, 0) * t1_2(ii) * t2_2(0, 0)
  465                         + t2_2(0, 1) * t1_2(ii) * t2_2(0, 1)
  466                         + t2_2(0, 2) * t1_2(ii) * t2_2(0, 2)
  467                         + t2_2(1, 0) * t1_2(ii) * t2_2(1, 0)
  468                         + t2_2(1, 1) * t1_2(ii) * t2_2(1, 1)
  469                         + t2_2(1, 2) * t1_2(ii) * t2_2(1, 2)
  470                         + t2_2(2, 0) * t1_2(ii) * t2_2(2, 0)
  471                         + t2_2(2, 1) * t1_2(ii) * t2_2(2, 1)
  472                         + t2_2(2, 2) * t1_2(ii) * t2_2(2, 2)),
  473                    "T3(j,k,i)*T2(j,k)(" + std::to_string(ii) + ")");
  474 
  476                      - (t2_2(0, 0) * t1_2(ii) * t2_2(0, 0)
  477                         + t2_2(0, 1) * t1_2(ii) * t2_2(0, 1)
  478                         + t2_2(0, 2) * t1_2(ii) * t2_2(0, 2)
  479                         + t2_2(1, 0) * t1_2(ii) * t2_2(1, 0)
  480                         + t2_2(1, 1) * t1_2(ii) * t2_2(1, 1)
  481                         + t2_2(1, 2) * t1_2(ii) * t2_2(1, 2)
  482                         + t2_2(2, 0) * t1_2(ii) * t2_2(2, 0)
  483                         + t2_2(2, 1) * t1_2(ii) * t2_2(2, 1)
  484                         + t2_2(2, 2) * t1_2(ii) * t2_2(2, 2)),
  485                    "T2(j,k)*T3(j,k,i)(" + std::to_string(ii) + ")");
  486 
  488                      - (t2_2(0, 0) * t1_2(ii) * t2_2(0, 0)
  489                         + t2_2(0, 1) * t1_2(ii) * t2_2(1, 0)
  490                         + t2_2(0, 2) * t1_2(ii) * t2_2(2, 0)
  491                         + t2_2(1, 0) * t1_2(ii) * t2_2(0, 1)
  492                         + t2_2(1, 1) * t1_2(ii) * t2_2(1, 1)
  493                         + t2_2(1, 2) * t1_2(ii) * t2_2(2, 1)
  494                         + t2_2(2, 0) * t1_2(ii) * t2_2(0, 2)
  495                         + t2_2(2, 1) * t1_2(ii) * t2_2(1, 2)
  496                         + t2_2(2, 2) * t1_2(ii) * t2_2(2, 2)),
  497                    "T3(j,k,i)*T2(k,j)(" + std::to_string(ii) + ")");
  498 
  500                      - (t2_2(0, 0) * t1_2(ii) * t2_2(0, 0)
  501                         + t2_2(0, 1) * t1_2(ii) * t2_2(1, 0)
  502                         + t2_2(0, 2) * t1_2(ii) * t2_2(2, 0)
  503                         + t2_2(1, 0) * t1_2(ii) * t2_2(0, 1)
  504                         + t2_2(1, 1) * t1_2(ii) * t2_2(1, 1)
  505                         + t2_2(1, 2) * t1_2(ii) * t2_2(2, 1)
  506                         + t2_2(2, 0) * t1_2(ii) * t2_2(0, 2)
  507                         + t2_2(2, 1) * t1_2(ii) * t2_2(1, 2)
  508                         + t2_2(2, 2) * t1_2(ii) * t2_2(2, 2)),
  509                    "T2(k,j)*T3(j,k,i)(" + std::to_string(ii) + ")");
  510    }
  511 
  513                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(0, 0, 0),
 
  514                "T3(i,j,k)*T2(k,l)(0,0,0)");
  516                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(0, 0, 1),
 
  517                "T3(i,j,k)*T2(k,l)(0,0,1)");
  519                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(0, 0, 2),
 
  520                "T3(i,j,k)*T2(k,l)(0,0,2)");
  522                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(0, 1, 0),
 
  523                "T3(i,j,k)*T2(k,l)(0,1,0)");
  525                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(0, 1, 1),
 
  526                "T3(i,j,k)*T2(k,l)(0,1,1)");
  528                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(0, 1, 2),
 
  529                "T3(i,j,k)*T2(k,l)(0,1,2)");
  531                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(0, 2, 0),
 
  532                "T3(i,j,k)*T2(k,l)(0,2,0)");
  534                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(0, 2, 1),
 
  535                "T3(i,j,k)*T2(k,l)(0,2,1)");
  537                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(0, 2, 2),
 
  538                "T3(i,j,k)*T2(k,l)(0,2,2)");
  540                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(1, 0, 0),
 
  541                "T3(i,j,k)*T2(k,l)(1,0,0)");
  543                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(1, 0, 1),
 
  544                "T3(i,j,k)*T2(k,l)(1,0,1)");
  546                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(1, 0, 2),
 
  547                "T3(i,j,k)*T2(k,l)(1,0,2)");
  549                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(1, 1, 0),
 
  550                "T3(i,j,k)*T2(k,l)(1,1,0)");
  552                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(1, 1, 1),
 
  553                "T3(i,j,k)*T2(k,l)(1,1,1)");
  555                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(1, 1, 2),
 
  556                "T3(i,j,k)*T2(k,l)(1,1,2)");
  558                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(1, 2, 0),
 
  559                "T3(i,j,k)*T2(k,l)(1,2,0)");
  561                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(1, 2, 1),
 
  562                "T3(i,j,k)*T2(k,l)(1,2,1)");
  564                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(1, 2, 2),
 
  565                "T3(i,j,k)*T2(k,l)(1,2,2)");
  567                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(2, 0, 0),
 
  568                "T3(i,j,k)*T2(k,l)(2,0,0)");
  570                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(2, 0, 1),
 
  571                "T3(i,j,k)*T2(k,l)(2,0,1)");
  573                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(2, 0, 2),
 
  574                "T3(i,j,k)*T2(k,l)(2,0,2)");
  576                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(2, 1, 0),
 
  577                "T3(i,j,k)*T2(k,l)(2,1,0)");
  579                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(2, 1, 1),
 
  580                "T3(i,j,k)*T2(k,l)(2,1,1)");
  582                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(2, 1, 2),
 
  583                "T3(i,j,k)*T2(k,l)(2,1,2)");
  585                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(2, 2, 0),
 
  586                "T3(i,j,k)*T2(k,l)(2,2,0)");
  588                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(2, 2, 1),
 
  589                "T3(i,j,k)*T2(k,l)(2,2,1)");
  591                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(2, 2, 2),
 
  592                "T3(i,j,k)*T2(k,l)(2,2,2)");
  593 
  595                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(0, 0, 0),
 
  596                "T3(i,j,k)*T2(l,k)(0,0,0)");
  598                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(0, 0, 1),
 
  599                "T3(i,j,k)*T2(l,k)(0,0,1)");
  601                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(0, 0, 2),
 
  602                "T3(i,j,k)*T2(l,k)(0,0,2)");
  604                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(0, 1, 0),
 
  605                "T3(i,j,k)*T2(l,k)(0,1,0)");
  607                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(0, 1, 1),
 
  608                "T3(i,j,k)*T2(l,k)(0,1,1)");
  610                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(0, 1, 2),
 
  611                "T3(i,j,k)*T2(l,k)(0,1,2)");
  613                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(0, 2, 0),
 
  614                "T3(i,j,k)*T2(l,k)(0,2,0)");
  616                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(0, 2, 1),
 
  617                "T3(i,j,k)*T2(l,k)(0,2,1)");
  619                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(0, 2, 2),
 
  620                "T3(i,j,k)*T2(l,k)(0,2,2)");
  622                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(1, 0, 0),
 
  623                "T3(i,j,k)*T2(l,k)(1,0,0)");
  625                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(1, 0, 1),
 
  626                "T3(i,j,k)*T2(l,k)(1,0,1)");
  628                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(1, 0, 2),
 
  629                "T3(i,j,k)*T2(l,k)(1,0,2)");
  631                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(1, 1, 0),
 
  632                "T3(i,j,k)*T2(l,k)(1,1,0)");
  634                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(1, 1, 1),
 
  635                "T3(i,j,k)*T2(l,k)(1,1,1)");
  637                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(1, 1, 2),
 
  638                "T3(i,j,k)*T2(l,k)(1,1,2)");
  640                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(1, 2, 0),
 
  641                "T3(i,j,k)*T2(l,k)(1,2,0)");
  643                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(1, 2, 1),
 
  644                "T3(i,j,k)*T2(l,k)(1,2,1)");
  646                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(1, 2, 2),
 
  647                "T3(i,j,k)*T2(l,k)(1,2,2)");
  649                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(2, 0, 0),
 
  650                "T3(i,j,k)*T2(l,k)(2,0,0)");
  652                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(2, 0, 1),
 
  653                "T3(i,j,k)*T2(l,k)(2,0,1)");
  655                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(2, 0, 2),
 
  656                "T3(i,j,k)*T2(l,k)(2,0,2)");
  658                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(2, 1, 0),
 
  659                "T3(i,j,k)*T2(l,k)(2,1,0)");
  661                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(2, 1, 1),
 
  662                "T3(i,j,k)*T2(l,k)(2,1,1)");
  664                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(2, 1, 2),
 
  665                "T3(i,j,k)*T2(l,k)(2,1,2)");
  667                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(2, 2, 0),
 
  668                "T3(i,j,k)*T2(l,k)(2,2,0)");
  670                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(2, 2, 1),
 
  671                "T3(i,j,k)*T2(l,k)(2,2,1)");
  673                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(2, 2, 2),
 
  674                "T3(i,j,k)*T2(l,k)(2,2,2)");
  675 
  677                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(0, 0, 0),
 
  678                "T2(k,l)*T3(i,j,k)(0,0,0)");
  680                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(0, 0, 1),
 
  681                "T2(k,l)*T3(i,j,k)(0,0,1)");
  683                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(0, 0, 2),
 
  684                "T2(k,l)*T3(i,j,k)(0,0,2)");
  686                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(0, 1, 0),
 
  687                "T2(k,l)*T3(i,j,k)(0,1,0)");
  689                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(0, 1, 1),
 
  690                "T2(k,l)*T3(i,j,k)(0,1,1)");
  692                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(0, 1, 2),
 
  693                "T2(k,l)*T3(i,j,k)(0,1,2)");
  695                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(0, 2, 0),
 
  696                "T2(k,l)*T3(i,j,k)(0,2,0)");
  698                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(0, 2, 1),
 
  699                "T2(k,l)*T3(i,j,k)(0,2,1)");
  701                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(0, 2, 2),
 
  702                "T2(k,l)*T3(i,j,k)(0,2,2)");
  704                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(1, 0, 0),
 
  705                "T2(k,l)*T3(i,j,k)(1,0,0)");
  707                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(1, 0, 1),
 
  708                "T2(k,l)*T3(i,j,k)(1,0,1)");
  710                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(1, 0, 2),
 
  711                "T2(k,l)*T3(i,j,k)(1,0,2)");
  713                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(1, 1, 0),
 
  714                "T2(k,l)*T3(i,j,k)(1,1,0)");
  716                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(1, 1, 1),
 
  717                "T2(k,l)*T3(i,j,k)(1,1,1)");
  719                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(1, 1, 2),
 
  720                "T2(k,l)*T3(i,j,k)(1,1,2)");
  722                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(1, 2, 0),
 
  723                "T2(k,l)*T3(i,j,k)(1,2,0)");
  725                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(1, 2, 1),
 
  726                "T2(k,l)*T3(i,j,k)(1,2,1)");
  728                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(1, 2, 2),
 
  729                "T2(k,l)*T3(i,j,k)(1,2,2)");
  731                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(2, 0, 0),
 
  732                "T2(k,l)*T3(i,j,k)(2,0,0)");
  734                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(2, 0, 1),
 
  735                "T2(k,l)*T3(i,j,k)(2,0,1)");
  737                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(2, 0, 2),
 
  738                "T2(k,l)*T3(i,j,k)(2,0,2)");
  740                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(2, 1, 0),
 
  741                "T2(k,l)*T3(i,j,k)(2,1,0)");
  743                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(2, 1, 1),
 
  744                "T2(k,l)*T3(i,j,k)(2,1,1)");
  746                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(2, 1, 2),
 
  747                "T2(k,l)*T3(i,j,k)(2,1,2)");
  749                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(2, 2, 0),
 
  750                "T2(k,l)*T3(i,j,k)(2,2,0)");
  752                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(2, 2, 1),
 
  753                "T2(k,l)*T3(i,j,k)(2,2,1)");
  755                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
k, 
l)))(2, 2, 2),
 
  756                "T2(k,l)*T3(i,j,k)(2,2,2)");
  757 
  759                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(0, 0, 0),
 
  760                "T2(l,k)*T3(i,j,k)(0,0,0)");
  762                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(0, 0, 1),
 
  763                "T2(l,k)*T3(i,j,k)(0,0,1)");
  765                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(0, 0, 2),
 
  766                "T2(l,k)*T3(i,j,k)(0,0,2)");
  768                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(0, 1, 0),
 
  769                "T2(l,k)*T3(i,j,k)(0,1,0)");
  771                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(0, 1, 1),
 
  772                "T2(l,k)*T3(i,j,k)(0,1,1)");
  774                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(0, 1, 2),
 
  775                "T2(l,k)*T3(i,j,k)(0,1,2)");
  777                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(0, 2, 0),
 
  778                "T2(l,k)*T3(i,j,k)(0,2,0)");
  780                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(0, 2, 1),
 
  781                "T2(l,k)*T3(i,j,k)(0,2,1)");
  783                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(0, 2, 2),
 
  784                "T2(l,k)*T3(i,j,k)(0,2,2)");
  786                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(1, 0, 0),
 
  787                "T2(l,k)*T3(i,j,k)(1,0,0)");
  789                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(1, 0, 1),
 
  790                "T2(l,k)*T3(i,j,k)(1,0,1)");
  792                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(1, 0, 2),
 
  793                "T2(l,k)*T3(i,j,k)(1,0,2)");
  795                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(1, 1, 0),
 
  796                "T2(l,k)*T3(i,j,k)(1,1,0)");
  798                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(1, 1, 1),
 
  799                "T2(l,k)*T3(i,j,k)(1,1,1)");
  801                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(1, 1, 2),
 
  802                "T2(l,k)*T3(i,j,k)(1,1,2)");
  804                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(1, 2, 0),
 
  805                "T2(l,k)*T3(i,j,k)(1,2,0)");
  807                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(1, 2, 1),
 
  808                "T2(l,k)*T3(i,j,k)(1,2,1)");
  810                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(1, 2, 2),
 
  811                "T2(l,k)*T3(i,j,k)(1,2,2)");
  813                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(2, 0, 0),
 
  814                "T2(l,k)*T3(i,j,k)(2,0,0)");
  816                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(2, 0, 1),
 
  817                "T2(l,k)*T3(i,j,k)(2,0,1)");
  819                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(2, 0, 2),
 
  820                "T2(l,k)*T3(i,j,k)(2,0,2)");
  822                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(2, 1, 0),
 
  823                "T2(l,k)*T3(i,j,k)(2,1,0)");
  825                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(2, 1, 1),
 
  826                "T2(l,k)*T3(i,j,k)(2,1,1)");
  828                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(2, 1, 2),
 
  829                "T2(l,k)*T3(i,j,k)(2,1,2)");
  831                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(2, 2, 0),
 
  832                "T2(l,k)*T3(i,j,k)(2,2,0)");
  834                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(2, 2, 1),
 
  835                "T2(l,k)*T3(i,j,k)(2,2,1)");
  837                  - (t2_2(
i, 
j) * (t1_2(
k) * t2_3(
l, 
k)))(2, 2, 2),
 
  838                "T2(l,k)*T3(i,j,k)(2,2,2)");
  839 
  841                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(0, 0, 0),
 
  842                "T3(i,j,l)*T2(j,k)(0,0,0)");
  844                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(0, 0, 1),
 
  845                "T3(i,j,l)*T2(j,k)(0,0,1)");
  847                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(0, 0, 2),
 
  848                "T3(i,j,l)*T2(j,k)(0,0,2)");
  850                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(0, 1, 0),
 
  851                "T3(i,j,l)*T2(j,k)(0,1,0)");
  853                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(0, 1, 1),
 
  854                "T3(i,j,l)*T2(j,k)(0,1,1)");
  856                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(0, 1, 2),
 
  857                "T3(i,j,l)*T2(j,k)(0,1,2)");
  859                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(0, 2, 0),
 
  860                "T3(i,j,l)*T2(j,k)(0,2,0)");
  862                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(0, 2, 1),
 
  863                "T3(i,j,l)*T2(j,k)(0,2,1)");
  865                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(0, 2, 2),
 
  866                "T3(i,j,l)*T2(j,k)(0,2,2)");
  868                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(1, 0, 0),
 
  869                "T3(i,j,l)*T2(j,k)(1,0,0)");
  871                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(1, 0, 1),
 
  872                "T3(i,j,l)*T2(j,k)(1,0,1)");
  874                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(1, 0, 2),
 
  875                "T3(i,j,l)*T2(j,k)(1,0,2)");
  877                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(1, 1, 0),
 
  878                "T3(i,j,l)*T2(j,k)(1,1,0)");
  880                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(1, 1, 1),
 
  881                "T3(i,j,l)*T2(j,k)(1,1,1)");
  883                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(1, 1, 2),
 
  884                "T3(i,j,l)*T2(j,k)(1,1,2)");
  886                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(1, 2, 0),
 
  887                "T3(i,j,l)*T2(j,k)(1,2,0)");
  889                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(1, 2, 1),
 
  890                "T3(i,j,l)*T2(j,k)(1,2,1)");
  892                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(1, 2, 2),
 
  893                "T3(i,j,l)*T2(j,k)(1,2,2)");
  895                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(2, 0, 0),
 
  896                "T3(i,j,l)*T2(j,k)(2,0,0)");
  898                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(2, 0, 1),
 
  899                "T3(i,j,l)*T2(j,k)(2,0,1)");
  901                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(2, 0, 2),
 
  902                "T3(i,j,l)*T2(j,k)(2,0,2)");
  904                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(2, 1, 0),
 
  905                "T3(i,j,l)*T2(j,k)(2,1,0)");
  907                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(2, 1, 1),
 
  908                "T3(i,j,l)*T2(j,k)(2,1,1)");
  910                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(2, 1, 2),
 
  911                "T3(i,j,l)*T2(j,k)(2,1,2)");
  913                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(2, 2, 0),
 
  914                "T3(i,j,l)*T2(j,k)(2,2,0)");
  916                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(2, 2, 1),
 
  917                "T3(i,j,l)*T2(j,k)(2,2,1)");
  919                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(2, 2, 2),
 
  920                "T3(i,j,l)*T2(j,k)(2,2,2)");
  921 
  923                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(0, 0, 0),
 
  924                "T3(i,j,l)*T2(k,j)(0,0,0)");
  926                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(0, 0, 1),
 
  927                "T3(i,j,l)*T2(k,j)(0,0,1)");
  929                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(0, 0, 2),
 
  930                "T3(i,j,l)*T2(k,j)(0,0,2)");
  932                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(0, 1, 0),
 
  933                "T3(i,j,l)*T2(k,j)(0,1,0)");
  935                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(0, 1, 1),
 
  936                "T3(i,j,l)*T2(k,j)(0,1,1)");
  938                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(0, 1, 2),
 
  939                "T3(i,j,l)*T2(k,j)(0,1,2)");
  941                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(0, 2, 0),
 
  942                "T3(i,j,l)*T2(k,j)(0,2,0)");
  944                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(0, 2, 1),
 
  945                "T3(i,j,l)*T2(k,j)(0,2,1)");
  947                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(0, 2, 2),
 
  948                "T3(i,j,l)*T2(k,j)(0,2,2)");
  950                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(1, 0, 0),
 
  951                "T3(i,j,l)*T2(k,j)(1,0,0)");
  953                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(1, 0, 1),
 
  954                "T3(i,j,l)*T2(k,j)(1,0,1)");
  956                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(1, 0, 2),
 
  957                "T3(i,j,l)*T2(k,j)(1,0,2)");
  959                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(1, 1, 0),
 
  960                "T3(i,j,l)*T2(k,j)(1,1,0)");
  962                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(1, 1, 1),
 
  963                "T3(i,j,l)*T2(k,j)(1,1,1)");
  965                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(1, 1, 2),
 
  966                "T3(i,j,l)*T2(k,j)(1,1,2)");
  968                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(1, 2, 0),
 
  969                "T3(i,j,l)*T2(k,j)(1,2,0)");
  971                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(1, 2, 1),
 
  972                "T3(i,j,l)*T2(k,j)(1,2,1)");
  974                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(1, 2, 2),
 
  975                "T3(i,j,l)*T2(k,j)(1,2,2)");
  977                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(2, 0, 0),
 
  978                "T3(i,j,l)*T2(k,j)(2,0,0)");
  980                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(2, 0, 1),
 
  981                "T3(i,j,l)*T2(k,j)(2,0,1)");
  983                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(2, 0, 2),
 
  984                "T3(i,j,l)*T2(k,j)(2,0,2)");
  986                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(2, 1, 0),
 
  987                "T3(i,j,l)*T2(k,j)(2,1,0)");
  989                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(2, 1, 1),
 
  990                "T3(i,j,l)*T2(k,j)(2,1,1)");
  992                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(2, 1, 2),
 
  993                "T3(i,j,l)*T2(k,j)(2,1,2)");
  995                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(2, 2, 0),
 
  996                "T3(i,j,l)*T2(k,j)(2,2,0)");
  998                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(2, 2, 1),
 
  999                "T3(i,j,l)*T2(k,j)(2,2,1)");
 1001                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(2, 2, 2),
 
 1002                "T3(i,j,l)*T2(k,j)(2,2,2)");
 1003 
 1005                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(0, 0, 0),
 
 1006                "T2(j,k)*T3(i,j,l)(0,0,0)");
 1008                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(0, 0, 1),
 
 1009                "T2(j,k)*T3(i,j,l)(0,0,1)");
 1011                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(0, 0, 2),
 
 1012                "T2(j,k)*T3(i,j,l)(0,0,2)");
 1014                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(0, 1, 0),
 
 1015                "T2(j,k)*T3(i,j,l)(0,1,0)");
 1017                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(0, 1, 1),
 
 1018                "T2(j,k)*T3(i,j,l)(0,1,1)");
 1020                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(0, 1, 2),
 
 1021                "T2(j,k)*T3(i,j,l)(0,1,2)");
 1023                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(0, 2, 0),
 
 1024                "T2(j,k)*T3(i,j,l)(0,2,0)");
 1026                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(0, 2, 1),
 
 1027                "T2(j,k)*T3(i,j,l)(0,2,1)");
 1029                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(0, 2, 2),
 
 1030                "T2(j,k)*T3(i,j,l)(0,2,2)");
 1032                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(1, 0, 0),
 
 1033                "T2(j,k)*T3(i,j,l)(1,0,0)");
 1035                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(1, 0, 1),
 
 1036                "T2(j,k)*T3(i,j,l)(1,0,1)");
 1038                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(1, 0, 2),
 
 1039                "T2(j,k)*T3(i,j,l)(1,0,2)");
 1041                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(1, 1, 0),
 
 1042                "T2(j,k)*T3(i,j,l)(1,1,0)");
 1044                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(1, 1, 1),
 
 1045                "T2(j,k)*T3(i,j,l)(1,1,1)");
 1047                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(1, 1, 2),
 
 1048                "T2(j,k)*T3(i,j,l)(1,1,2)");
 1050                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(1, 2, 0),
 
 1051                "T2(j,k)*T3(i,j,l)(1,2,0)");
 1053                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(1, 2, 1),
 
 1054                "T2(j,k)*T3(i,j,l)(1,2,1)");
 1056                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(1, 2, 2),
 
 1057                "T2(j,k)*T3(i,j,l)(1,2,2)");
 1059                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(2, 0, 0),
 
 1060                "T2(j,k)*T3(i,j,l)(2,0,0)");
 1062                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(2, 0, 1),
 
 1063                "T2(j,k)*T3(i,j,l)(2,0,1)");
 1065                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(2, 0, 2),
 
 1066                "T2(j,k)*T3(i,j,l)(2,0,2)");
 1068                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(2, 1, 0),
 
 1069                "T2(j,k)*T3(i,j,l)(2,1,0)");
 1071                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(2, 1, 1),
 
 1072                "T2(j,k)*T3(i,j,l)(2,1,1)");
 1074                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(2, 1, 2),
 
 1075                "T2(j,k)*T3(i,j,l)(2,1,2)");
 1077                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(2, 2, 0),
 
 1078                "T2(j,k)*T3(i,j,l)(2,2,0)");
 1080                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(2, 2, 1),
 
 1081                "T2(j,k)*T3(i,j,l)(2,2,1)");
 1083                  - ((t2_2(
i, 
j) * t2_3(
j, 
k)) * t1_2(
l))(2, 2, 2),
 
 1084                "T2(j,k)*T3(i,j,l)(2,2,2)");
 1085 
 1087                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(0, 0, 0),
 
 1088                "T2(k,j)*T3(i,j,l)(0,0,0)");
 1090                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(0, 0, 1),
 
 1091                "T2(k,j)*T3(i,j,l)(0,0,1)");
 1093                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(0, 0, 2),
 
 1094                "T2(k,j)*T3(i,j,l)(0,0,2)");
 1096                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(0, 1, 0),
 
 1097                "T2(k,j)*T3(i,j,l)(0,1,0)");
 1099                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(0, 1, 1),
 
 1100                "T2(k,j)*T3(i,j,l)(0,1,1)");
 1102                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(0, 1, 2),
 
 1103                "T2(k,j)*T3(i,j,l)(0,1,2)");
 1105                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(0, 2, 0),
 
 1106                "T2(k,j)*T3(i,j,l)(0,2,0)");
 1108                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(0, 2, 1),
 
 1109                "T2(k,j)*T3(i,j,l)(0,2,1)");
 1111                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(0, 2, 2),
 
 1112                "T2(k,j)*T3(i,j,l)(0,2,2)");
 1114                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(1, 0, 0),
 
 1115                "T2(k,j)*T3(i,j,l)(1,0,0)");
 1117                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(1, 0, 1),
 
 1118                "T2(k,j)*T3(i,j,l)(1,0,1)");
 1120                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(1, 0, 2),
 
 1121                "T2(k,j)*T3(i,j,l)(1,0,2)");
 1123                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(1, 1, 0),
 
 1124                "T2(k,j)*T3(i,j,l)(1,1,0)");
 1126                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(1, 1, 1),
 
 1127                "T2(k,j)*T3(i,j,l)(1,1,1)");
 1129                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(1, 1, 2),
 
 1130                "T2(k,j)*T3(i,j,l)(1,1,2)");
 1132                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(1, 2, 0),
 
 1133                "T2(k,j)*T3(i,j,l)(1,2,0)");
 1135                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(1, 2, 1),
 
 1136                "T2(k,j)*T3(i,j,l)(1,2,1)");
 1138                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(1, 2, 2),
 
 1139                "T2(k,j)*T3(i,j,l)(1,2,2)");
 1141                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(2, 0, 0),
 
 1142                "T2(k,j)*T3(i,j,l)(2,0,0)");
 1144                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(2, 0, 1),
 
 1145                "T2(k,j)*T3(i,j,l)(2,0,1)");
 1147                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(2, 0, 2),
 
 1148                "T2(k,j)*T3(i,j,l)(2,0,2)");
 1150                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(2, 1, 0),
 
 1151                "T2(k,j)*T3(i,j,l)(2,1,0)");
 1153                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(2, 1, 1),
 
 1154                "T2(k,j)*T3(i,j,l)(2,1,1)");
 1156                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(2, 1, 2),
 
 1157                "T2(k,j)*T3(i,j,l)(2,1,2)");
 1159                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(2, 2, 0),
 
 1160                "T2(k,j)*T3(i,j,l)(2,2,0)");
 1162                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(2, 2, 1),
 
 1163                "T2(k,j)*T3(i,j,l)(2,2,1)");
 1165                  - ((t2_2(
i, 
j) * t2_3(
k, 
j)) * t1_2(
l))(2, 2, 2),
 
 1166                "T2(k,j)*T3(i,j,l)(2,2,2)");
 1167 
 1169                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(0, 0, 0),
 
 1170                "T3(i,k,l)*T2(i,j)(0,0,0)");
 1172                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(0, 0, 1),
 
 1173                "T3(i,k,l)*T2(i,j)(0,0,1)");
 1175                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(0, 0, 2),
 
 1176                "T3(i,k,l)*T2(i,j)(0,0,2)");
 1178                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(0, 1, 0),
 
 1179                "T3(i,k,l)*T2(i,j)(0,1,0)");
 1181                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(0, 1, 1),
 
 1182                "T3(i,k,l)*T2(i,j)(0,1,1)");
 1184                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(0, 1, 2),
 
 1185                "T3(i,k,l)*T2(i,j)(0,1,2)");
 1187                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(0, 2, 0),
 
 1188                "T3(i,k,l)*T2(i,j)(0,2,0)");
 1190                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(0, 2, 1),
 
 1191                "T3(i,k,l)*T2(i,j)(0,2,1)");
 1193                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(0, 2, 2),
 
 1194                "T3(i,k,l)*T2(i,j)(0,2,2)");
 1196                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(1, 0, 0),
 
 1197                "T3(i,k,l)*T2(i,j)(1,0,0)");
 1199                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(1, 0, 1),
 
 1200                "T3(i,k,l)*T2(i,j)(1,0,1)");
 1202                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(1, 0, 2),
 
 1203                "T3(i,k,l)*T2(i,j)(1,0,2)");
 1205                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(1, 1, 0),
 
 1206                "T3(i,k,l)*T2(i,j)(1,1,0)");
 1208                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(1, 1, 1),
 
 1209                "T3(i,k,l)*T2(i,j)(1,1,1)");
 1211                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(1, 1, 2),
 
 1212                "T3(i,k,l)*T2(i,j)(1,1,2)");
 1214                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(1, 2, 0),
 
 1215                "T3(i,k,l)*T2(i,j)(1,2,0)");
 1217                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(1, 2, 1),
 
 1218                "T3(i,k,l)*T2(i,j)(1,2,1)");
 1220                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(1, 2, 2),
 
 1221                "T3(i,k,l)*T2(i,j)(1,2,2)");
 1223                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(2, 0, 0),
 
 1224                "T3(i,k,l)*T2(i,j)(2,0,0)");
 1226                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(2, 0, 1),
 
 1227                "T3(i,k,l)*T2(i,j)(2,0,1)");
 1229                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(2, 0, 2),
 
 1230                "T3(i,k,l)*T2(i,j)(2,0,2)");
 1232                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(2, 1, 0),
 
 1233                "T3(i,k,l)*T2(i,j)(2,1,0)");
 1235                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(2, 1, 1),
 
 1236                "T3(i,k,l)*T2(i,j)(2,1,1)");
 1238                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(2, 1, 2),
 
 1239                "T3(i,k,l)*T2(i,j)(2,1,2)");
 1241                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(2, 2, 0),
 
 1242                "T3(i,k,l)*T2(i,j)(2,2,0)");
 1244                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(2, 2, 1),
 
 1245                "T3(i,k,l)*T2(i,j)(2,2,1)");
 1247                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(2, 2, 2),
 
 1248                "T3(i,k,l)*T2(i,j)(2,2,2)");
 1249 
 1251                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(0, 0, 0),
 
 1252                "T3(i,k,l)*T2(j,i)(0,0,0)");
 1254                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(0, 0, 1),
 
 1255                "T3(i,k,l)*T2(j,i)(0,0,1)");
 1257                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(0, 0, 2),
 
 1258                "T3(i,k,l)*T2(j,i)(0,0,2)");
 1260                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(0, 1, 0),
 
 1261                "T3(i,k,l)*T2(j,i)(0,1,0)");
 1263                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(0, 1, 1),
 
 1264                "T3(i,k,l)*T2(j,i)(0,1,1)");
 1266                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(0, 1, 2),
 
 1267                "T3(i,k,l)*T2(j,i)(0,1,2)");
 1269                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(0, 2, 0),
 
 1270                "T3(i,k,l)*T2(j,i)(0,2,0)");
 1272                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(0, 2, 1),
 
 1273                "T3(i,k,l)*T2(j,i)(0,2,1)");
 1275                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(0, 2, 2),
 
 1276                "T3(i,k,l)*T2(j,i)(0,2,2)");
 1278                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(1, 0, 0),
 
 1279                "T3(i,k,l)*T2(j,i)(1,0,0)");
 1281                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(1, 0, 1),
 
 1282                "T3(i,k,l)*T2(j,i)(1,0,1)");
 1284                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(1, 0, 2),
 
 1285                "T3(i,k,l)*T2(j,i)(1,0,2)");
 1287                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(1, 1, 0),
 
 1288                "T3(i,k,l)*T2(j,i)(1,1,0)");
 1290                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(1, 1, 1),
 
 1291                "T3(i,k,l)*T2(j,i)(1,1,1)");
 1293                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(1, 1, 2),
 
 1294                "T3(i,k,l)*T2(j,i)(1,1,2)");
 1296                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(1, 2, 0),
 
 1297                "T3(i,k,l)*T2(j,i)(1,2,0)");
 1299                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(1, 2, 1),
 
 1300                "T3(i,k,l)*T2(j,i)(1,2,1)");
 1302                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(1, 2, 2),
 
 1303                "T3(i,k,l)*T2(j,i)(1,2,2)");
 1305                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(2, 0, 0),
 
 1306                "T3(i,k,l)*T2(j,i)(2,0,0)");
 1308                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(2, 0, 1),
 
 1309                "T3(i,k,l)*T2(j,i)(2,0,1)");
 1311                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(2, 0, 2),
 
 1312                "T3(i,k,l)*T2(j,i)(2,0,2)");
 1314                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(2, 1, 0),
 
 1315                "T3(i,k,l)*T2(j,i)(2,1,0)");
 1317                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(2, 1, 1),
 
 1318                "T3(i,k,l)*T2(j,i)(2,1,1)");
 1320                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(2, 1, 2),
 
 1321                "T3(i,k,l)*T2(j,i)(2,1,2)");
 1323                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(2, 2, 0),
 
 1324                "T3(i,k,l)*T2(j,i)(2,2,0)");
 1326                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(2, 2, 1),
 
 1327                "T3(i,k,l)*T2(j,i)(2,2,1)");
 1329                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(2, 2, 2),
 
 1330                "T3(i,k,l)*T2(j,i)(2,2,2)");
 1331 
 1333                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(0, 0, 0),
 
 1334                "T2(i,j)*T3(i,k,l)(0,0,0)");
 1336                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(0, 0, 1),
 
 1337                "T2(i,j)*T3(i,k,l)(0,0,1)");
 1339                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(0, 0, 2),
 
 1340                "T2(i,j)*T3(i,k,l)(0,0,2)");
 1342                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(0, 1, 0),
 
 1343                "T2(i,j)*T3(i,k,l)(0,1,0)");
 1345                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(0, 1, 1),
 
 1346                "T2(i,j)*T3(i,k,l)(0,1,1)");
 1348                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(0, 1, 2),
 
 1349                "T2(i,j)*T3(i,k,l)(0,1,2)");
 1351                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(0, 2, 0),
 
 1352                "T2(i,j)*T3(i,k,l)(0,2,0)");
 1354                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(0, 2, 1),
 
 1355                "T2(i,j)*T3(i,k,l)(0,2,1)");
 1357                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(0, 2, 2),
 
 1358                "T2(i,j)*T3(i,k,l)(0,2,2)");
 1360                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(1, 0, 0),
 
 1361                "T2(i,j)*T3(i,k,l)(1,0,0)");
 1363                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(1, 0, 1),
 
 1364                "T2(i,j)*T3(i,k,l)(1,0,1)");
 1366                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(1, 0, 2),
 
 1367                "T2(i,j)*T3(i,k,l)(1,0,2)");
 1369                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(1, 1, 0),
 
 1370                "T2(i,j)*T3(i,k,l)(1,1,0)");
 1372                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(1, 1, 1),
 
 1373                "T2(i,j)*T3(i,k,l)(1,1,1)");
 1375                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(1, 1, 2),
 
 1376                "T2(i,j)*T3(i,k,l)(1,1,2)");
 1378                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(1, 2, 0),
 
 1379                "T2(i,j)*T3(i,k,l)(1,2,0)");
 1381                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(1, 2, 1),
 
 1382                "T2(i,j)*T3(i,k,l)(1,2,1)");
 1384                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(1, 2, 2),
 
 1385                "T2(i,j)*T3(i,k,l)(1,2,2)");
 1387                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(2, 0, 0),
 
 1388                "T2(i,j)*T3(i,k,l)(2,0,0)");
 1390                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(2, 0, 1),
 
 1391                "T2(i,j)*T3(i,k,l)(2,0,1)");
 1393                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(2, 0, 2),
 
 1394                "T2(i,j)*T3(i,k,l)(2,0,2)");
 1396                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(2, 1, 0),
 
 1397                "T2(i,j)*T3(i,k,l)(2,1,0)");
 1399                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(2, 1, 1),
 
 1400                "T2(i,j)*T3(i,k,l)(2,1,1)");
 1402                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(2, 1, 2),
 
 1403                "T2(i,j)*T3(i,k,l)(2,1,2)");
 1405                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(2, 2, 0),
 
 1406                "T2(i,j)*T3(i,k,l)(2,2,0)");
 1408                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(2, 2, 1),
 
 1409                "T2(i,j)*T3(i,k,l)(2,2,1)");
 1411                  - ((t2_3(
i, 
j) * t2_2(
i, 
k)) * t1_2(
l))(2, 2, 2),
 
 1412                "T2(i,j)*T3(i,k,l)(2,2,2)");
 1413 
 1415                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(0, 0, 0),
 
 1416                "T2(j,i)*T3(i,k,l)(0,0,0)");
 1418                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(0, 0, 1),
 
 1419                "T2(j,i)*T3(i,k,l)(0,0,1)");
 1421                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(0, 0, 2),
 
 1422                "T2(j,i)*T3(i,k,l)(0,0,2)");
 1424                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(0, 1, 0),
 
 1425                "T2(j,i)*T3(i,k,l)(0,1,0)");
 1427                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(0, 1, 1),
 
 1428                "T2(j,i)*T3(i,k,l)(0,1,1)");
 1430                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(0, 1, 2),
 
 1431                "T2(j,i)*T3(i,k,l)(0,1,2)");
 1433                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(0, 2, 0),
 
 1434                "T2(j,i)*T3(i,k,l)(0,2,0)");
 1436                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(0, 2, 1),
 
 1437                "T2(j,i)*T3(i,k,l)(0,2,1)");
 1439                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(0, 2, 2),
 
 1440                "T2(j,i)*T3(i,k,l)(0,2,2)");
 1442                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(1, 0, 0),
 
 1443                "T2(j,i)*T3(i,k,l)(1,0,0)");
 1445                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(1, 0, 1),
 
 1446                "T2(j,i)*T3(i,k,l)(1,0,1)");
 1448                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(1, 0, 2),
 
 1449                "T2(j,i)*T3(i,k,l)(1,0,2)");
 1451                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(1, 1, 0),
 
 1452                "T2(j,i)*T3(i,k,l)(1,1,0)");
 1454                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(1, 1, 1),
 
 1455                "T2(j,i)*T3(i,k,l)(1,1,1)");
 1457                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(1, 1, 2),
 
 1458                "T2(j,i)*T3(i,k,l)(1,1,2)");
 1460                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(1, 2, 0),
 
 1461                "T2(j,i)*T3(i,k,l)(1,2,0)");
 1463                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(1, 2, 1),
 
 1464                "T2(j,i)*T3(i,k,l)(1,2,1)");
 1466                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(1, 2, 2),
 
 1467                "T2(j,i)*T3(i,k,l)(1,2,2)");
 1469                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(2, 0, 0),
 
 1470                "T2(j,i)*T3(i,k,l)(2,0,0)");
 1472                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(2, 0, 1),
 
 1473                "T2(j,i)*T3(i,k,l)(2,0,1)");
 1475                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(2, 0, 2),
 
 1476                "T2(j,i)*T3(i,k,l)(2,0,2)");
 1478                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(2, 1, 0),
 
 1479                "T2(j,i)*T3(i,k,l)(2,1,0)");
 1481                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(2, 1, 1),
 
 1482                "T2(j,i)*T3(i,k,l)(2,1,1)");
 1484                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(2, 1, 2),
 
 1485                "T2(j,i)*T3(i,k,l)(2,1,2)");
 1487                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(2, 2, 0),
 
 1488                "T2(j,i)*T3(i,k,l)(2,2,0)");
 1490                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(2, 2, 1),
 
 1491                "T2(j,i)*T3(i,k,l)(2,2,1)");
 1493                  - ((t2_3(
j, 
i) * t2_2(
i, 
k)) * t1_2(
l))(2, 2, 2),
 
 1494                "T2(j,i)*T3(i,k,l)(2,2,2)");
 1495 
 1497                  - (t2_2(0, 0) * t1_2(0) * t2_3(0, 0) * t1_1(0)
 1498                     + t2_2(1, 0) * t1_2(0) * t2_3(1, 0) * t1_1(0)
 1499                     + t2_2(2, 0) * t1_2(0) * t2_3(2, 0) * t1_1(0)
 1500                     + t2_2(0, 1) * t1_2(0) * t2_3(0, 1) * t1_1(0)
 1501                     + t2_2(1, 1) * t1_2(0) * t2_3(1, 1) * t1_1(0)
 1502                     + t2_2(2, 1) * t1_2(0) * t2_3(2, 1) * t1_1(0)
 1503                     + t2_2(0, 2) * t1_2(0) * t2_3(0, 2) * t1_1(0)
 1504                     + t2_2(1, 2) * t1_2(0) * t2_3(1, 2) * t1_1(0)
 1505                     + t2_2(2, 2) * t1_2(0) * t2_3(2, 2) * t1_1(0)
 1506                     + t2_2(0, 0) * t1_2(1) * t2_3(0, 0) * t1_1(1)
 1507                     + t2_2(1, 0) * t1_2(1) * t2_3(1, 0) * t1_1(1)
 1508                     + t2_2(2, 0) * t1_2(1) * t2_3(2, 0) * t1_1(1)
 1509                     + t2_2(0, 1) * t1_2(1) * t2_3(0, 1) * t1_1(1)
 1510                     + t2_2(1, 1) * t1_2(1) * t2_3(1, 1) * t1_1(1)
 1511                     + t2_2(2, 1) * t1_2(1) * t2_3(2, 1) * t1_1(1)
 1512                     + t2_2(0, 2) * t1_2(1) * t2_3(0, 2) * t1_1(1)
 1513                     + t2_2(1, 2) * t1_2(1) * t2_3(1, 2) * t1_1(1)
 1514                     + t2_2(2, 2) * t1_2(1) * t2_3(2, 2) * t1_1(1)
 1515                     + t2_2(0, 0) * t1_2(2) * t2_3(0, 0) * t1_1(2)
 1516                     + t2_2(1, 0) * t1_2(2) * t2_3(1, 0) * t1_1(2)
 1517                     + t2_2(2, 0) * t1_2(2) * t2_3(2, 0) * t1_1(2)
 1518                     + t2_2(0, 1) * t1_2(2) * t2_3(0, 1) * t1_1(2)
 1519                     + t2_2(1, 1) * t1_2(2) * t2_3(1, 1) * t1_1(2)
 1520                     + t2_2(2, 1) * t1_2(2) * t2_3(2, 1) * t1_1(2)
 1521                     + t2_2(0, 2) * t1_2(2) * t2_3(0, 2) * t1_1(2)
 1522                     + t2_2(1, 2) * t1_2(2) * t2_3(1, 2) * t1_1(2)
 1523                     + t2_2(2, 2) * t1_2(2) * t2_3(2, 2) * t1_1(2)),
 1524                "T3(i,j,k)*T3(i,j,k)");
 1525 
 1527                  - (t2_2(0, 0) * t1_2(0) * t2_3(0, 0) * t1_1(0)
 1528                     + t2_2(1, 0) * t1_2(0) * t2_3(0, 1) * t1_1(0)
 1529                     + t2_2(2, 0) * t1_2(0) * t2_3(0, 2) * t1_1(0)
 1530                     + t2_2(0, 1) * t1_2(0) * t2_3(1, 0) * t1_1(0)
 1531                     + t2_2(1, 1) * t1_2(0) * t2_3(1, 1) * t1_1(0)
 1532                     + t2_2(2, 1) * t1_2(0) * t2_3(1, 2) * t1_1(0)
 1533                     + t2_2(0, 2) * t1_2(0) * t2_3(2, 0) * t1_1(0)
 1534                     + t2_2(1, 2) * t1_2(0) * t2_3(2, 1) * t1_1(0)
 1535                     + t2_2(2, 2) * t1_2(0) * t2_3(2, 2) * t1_1(0)
 1536                     + t2_2(0, 0) * t1_2(1) * t2_3(0, 0) * t1_1(1)
 1537                     + t2_2(1, 0) * t1_2(1) * t2_3(0, 1) * t1_1(1)
 1538                     + t2_2(2, 0) * t1_2(1) * t2_3(0, 2) * t1_1(1)
 1539                     + t2_2(0, 1) * t1_2(1) * t2_3(1, 0) * t1_1(1)
 1540                     + t2_2(1, 1) * t1_2(1) * t2_3(1, 1) * t1_1(1)
 1541                     + t2_2(2, 1) * t1_2(1) * t2_3(1, 2) * t1_1(1)
 1542                     + t2_2(0, 2) * t1_2(1) * t2_3(2, 0) * t1_1(1)
 1543                     + t2_2(1, 2) * t1_2(1) * t2_3(2, 1) * t1_1(1)
 1544                     + t2_2(2, 2) * t1_2(1) * t2_3(2, 2) * t1_1(1)
 1545                     + t2_2(0, 0) * t1_2(2) * t2_3(0, 0) * t1_1(2)
 1546                     + t2_2(1, 0) * t1_2(2) * t2_3(0, 1) * t1_1(2)
 1547                     + t2_2(2, 0) * t1_2(2) * t2_3(0, 2) * t1_1(2)
 1548                     + t2_2(0, 1) * t1_2(2) * t2_3(1, 0) * t1_1(2)
 1549                     + t2_2(1, 1) * t1_2(2) * t2_3(1, 1) * t1_1(2)
 1550                     + t2_2(2, 1) * t1_2(2) * t2_3(1, 2) * t1_1(2)
 1551                     + t2_2(0, 2) * t1_2(2) * t2_3(2, 0) * t1_1(2)
 1552                     + t2_2(1, 2) * t1_2(2) * t2_3(2, 1) * t1_1(2)
 1553                     + t2_2(2, 2) * t1_2(2) * t2_3(2, 2) * t1_1(2)),
 1554                "T3(i,j,k)*T3(j,i,k)");
 1555 
 1557                  - (t2_2(0, 0) * t1_2(0) * t2_3(0, 0) * t1_1(0)
 1558                     + t2_2(1, 0) * t1_2(0) * t2_3(1, 0) * t1_1(0)
 1559                     + t2_2(2, 0) * t1_2(0) * t2_3(2, 0) * t1_1(0)
 1560                     + t2_2(0, 1) * t1_2(0) * t2_3(0, 0) * t1_1(1)
 1561                     + t2_2(1, 1) * t1_2(0) * t2_3(1, 0) * t1_1(1)
 1562                     + t2_2(2, 1) * t1_2(0) * t2_3(2, 0) * t1_1(1)
 1563                     + t2_2(0, 2) * t1_2(0) * t2_3(0, 0) * t1_1(2)
 1564                     + t2_2(1, 2) * t1_2(0) * t2_3(1, 0) * t1_1(2)
 1565                     + t2_2(2, 2) * t1_2(0) * t2_3(2, 0) * t1_1(2)
 1566                     + t2_2(0, 0) * t1_2(1) * t2_3(0, 1) * t1_1(0)
 1567                     + t2_2(1, 0) * t1_2(1) * t2_3(1, 1) * t1_1(0)
 1568                     + t2_2(2, 0) * t1_2(1) * t2_3(2, 1) * t1_1(0)
 1569                     + t2_2(0, 1) * t1_2(1) * t2_3(0, 1) * t1_1(1)
 1570                     + t2_2(1, 1) * t1_2(1) * t2_3(1, 1) * t1_1(1)
 1571                     + t2_2(2, 1) * t1_2(1) * t2_3(2, 1) * t1_1(1)
 1572                     + t2_2(0, 2) * t1_2(1) * t2_3(0, 1) * t1_1(2)
 1573                     + t2_2(1, 2) * t1_2(1) * t2_3(1, 1) * t1_1(2)
 1574                     + t2_2(2, 2) * t1_2(1) * t2_3(2, 1) * t1_1(2)
 1575                     + t2_2(0, 0) * t1_2(2) * t2_3(0, 2) * t1_1(0)
 1576                     + t2_2(1, 0) * t1_2(2) * t2_3(1, 2) * t1_1(0)
 1577                     + t2_2(2, 0) * t1_2(2) * t2_3(2, 2) * t1_1(0)
 1578                     + t2_2(0, 1) * t1_2(2) * t2_3(0, 2) * t1_1(1)
 1579                     + t2_2(1, 1) * t1_2(2) * t2_3(1, 2) * t1_1(1)
 1580                     + t2_2(2, 1) * t1_2(2) * t2_3(2, 2) * t1_1(1)
 1581                     + t2_2(0, 2) * t1_2(2) * t2_3(0, 2) * t1_1(2)
 1582                     + t2_2(1, 2) * t1_2(2) * t2_3(1, 2) * t1_1(2)
 1583                     + t2_2(2, 2) * t1_2(2) * t2_3(2, 2) * t1_1(2)),
 1584                "T3(i,j,k)*T3(i,k,j)");
 1585 
 1587                  - (t2_2(0, 0) * t1_2(0) * t2_3(0, 0) * t1_1(0)
 1588                     + t2_2(1, 0) * t1_2(0) * t2_3(0, 0) * t1_1(1)
 1589                     + t2_2(2, 0) * t1_2(0) * t2_3(0, 0) * t1_1(2)
 1590                     + t2_2(0, 1) * t1_2(0) * t2_3(1, 0) * t1_1(0)
 1591                     + t2_2(1, 1) * t1_2(0) * t2_3(1, 0) * t1_1(1)
 1592                     + t2_2(2, 1) * t1_2(0) * t2_3(1, 0) * t1_1(2)
 1593                     + t2_2(0, 2) * t1_2(0) * t2_3(2, 0) * t1_1(0)
 1594                     + t2_2(1, 2) * t1_2(0) * t2_3(2, 0) * t1_1(1)
 1595                     + t2_2(2, 2) * t1_2(0) * t2_3(2, 0) * t1_1(2)
 1596                     + t2_2(0, 0) * t1_2(1) * t2_3(0, 1) * t1_1(0)
 1597                     + t2_2(1, 0) * t1_2(1) * t2_3(0, 1) * t1_1(1)
 1598                     + t2_2(2, 0) * t1_2(1) * t2_3(0, 1) * t1_1(2)
 1599                     + t2_2(0, 1) * t1_2(1) * t2_3(1, 1) * t1_1(0)
 1600                     + t2_2(1, 1) * t1_2(1) * t2_3(1, 1) * t1_1(1)
 1601                     + t2_2(2, 1) * t1_2(1) * t2_3(1, 1) * t1_1(2)
 1602                     + t2_2(0, 2) * t1_2(1) * t2_3(2, 1) * t1_1(0)
 1603                     + t2_2(1, 2) * t1_2(1) * t2_3(2, 1) * t1_1(1)
 1604                     + t2_2(2, 2) * t1_2(1) * t2_3(2, 1) * t1_1(2)
 1605                     + t2_2(0, 0) * t1_2(2) * t2_3(0, 2) * t1_1(0)
 1606                     + t2_2(1, 0) * t1_2(2) * t2_3(0, 2) * t1_1(1)
 1607                     + t2_2(2, 0) * t1_2(2) * t2_3(0, 2) * t1_1(2)
 1608                     + t2_2(0, 1) * t1_2(2) * t2_3(1, 2) * t1_1(0)
 1609                     + t2_2(1, 1) * t1_2(2) * t2_3(1, 2) * t1_1(1)
 1610                     + t2_2(2, 1) * t1_2(2) * t2_3(1, 2) * t1_1(2)
 1611                     + t2_2(0, 2) * t1_2(2) * t2_3(2, 2) * t1_1(0)
 1612                     + t2_2(1, 2) * t1_2(2) * t2_3(2, 2) * t1_1(1)
 1613                     + t2_2(2, 2) * t1_2(2) * t2_3(2, 2) * t1_1(2)),
 1614                "T3(i,j,k)*T3(j,k,i)");
 1615 
 1617                  - (t2_2(0, 0) * t1_2(0) * t2_3(0, 0) * t1_1(0)
 1618                     + t2_2(1, 0) * t1_2(0) * t2_3(0, 0) * t1_1(1)
 1619                     + t2_2(2, 0) * t1_2(0) * t2_3(0, 0) * t1_1(2)
 1620                     + t2_2(0, 1) * t1_2(0) * t2_3(0, 1) * t1_1(0)
 1621                     + t2_2(1, 1) * t1_2(0) * t2_3(0, 1) * t1_1(1)
 1622                     + t2_2(2, 1) * t1_2(0) * t2_3(0, 1) * t1_1(2)
 1623                     + t2_2(0, 2) * t1_2(0) * t2_3(0, 2) * t1_1(0)
 1624                     + t2_2(1, 2) * t1_2(0) * t2_3(0, 2) * t1_1(1)
 1625                     + t2_2(2, 2) * t1_2(0) * t2_3(0, 2) * t1_1(2)
 1626                     + t2_2(0, 0) * t1_2(1) * t2_3(1, 0) * t1_1(0)
 1627                     + t2_2(1, 0) * t1_2(1) * t2_3(1, 0) * t1_1(1)
 1628                     + t2_2(2, 0) * t1_2(1) * t2_3(1, 0) * t1_1(2)
 1629                     + t2_2(0, 1) * t1_2(1) * t2_3(1, 1) * t1_1(0)
 1630                     + t2_2(1, 1) * t1_2(1) * t2_3(1, 1) * t1_1(1)
 1631                     + t2_2(2, 1) * t1_2(1) * t2_3(1, 1) * t1_1(2)
 1632                     + t2_2(0, 2) * t1_2(1) * t2_3(1, 2) * t1_1(0)
 1633                     + t2_2(1, 2) * t1_2(1) * t2_3(1, 2) * t1_1(1)
 1634                     + t2_2(2, 2) * t1_2(1) * t2_3(1, 2) * t1_1(2)
 1635                     + t2_2(0, 0) * t1_2(2) * t2_3(2, 0) * t1_1(0)
 1636                     + t2_2(1, 0) * t1_2(2) * t2_3(2, 0) * t1_1(1)
 1637                     + t2_2(2, 0) * t1_2(2) * t2_3(2, 0) * t1_1(2)
 1638                     + t2_2(0, 1) * t1_2(2) * t2_3(2, 1) * t1_1(0)
 1639                     + t2_2(1, 1) * t1_2(2) * t2_3(2, 1) * t1_1(1)
 1640                     + t2_2(2, 1) * t1_2(2) * t2_3(2, 1) * t1_1(2)
 1641                     + t2_2(0, 2) * t1_2(2) * t2_3(2, 2) * t1_1(0)
 1642                     + t2_2(1, 2) * t1_2(2) * t2_3(2, 2) * t1_1(1)
 1643                     + t2_2(2, 2) * t1_2(2) * t2_3(2, 2) * t1_1(2)),
 1644                "T3(i,j,k)*T3(k,j,i)");
 1645 
 1647                  - (t2_2(0, 0) * t1_2(0) * t2_3(0, 0) * t1_1(0)
 1648                     + t2_2(1, 0) * t1_2(0) * t2_3(0, 1) * t1_1(0)
 1649                     + t2_2(2, 0) * t1_2(0) * t2_3(0, 2) * t1_1(0)
 1650                     + t2_2(0, 1) * t1_2(0) * t2_3(0, 0) * t1_1(1)
 1651                     + t2_2(1, 1) * t1_2(0) * t2_3(0, 1) * t1_1(1)
 1652                     + t2_2(2, 1) * t1_2(0) * t2_3(0, 2) * t1_1(1)
 1653                     + t2_2(0, 2) * t1_2(0) * t2_3(0, 0) * t1_1(2)
 1654                     + t2_2(1, 2) * t1_2(0) * t2_3(0, 1) * t1_1(2)
 1655                     + t2_2(2, 2) * t1_2(0) * t2_3(0, 2) * t1_1(2)
 1656                     + t2_2(0, 0) * t1_2(1) * t2_3(1, 0) * t1_1(0)
 1657                     + t2_2(1, 0) * t1_2(1) * t2_3(1, 1) * t1_1(0)
 1658                     + t2_2(2, 0) * t1_2(1) * t2_3(1, 2) * t1_1(0)
 1659                     + t2_2(0, 1) * t1_2(1) * t2_3(1, 0) * t1_1(1)
 1660                     + t2_2(1, 1) * t1_2(1) * t2_3(1, 1) * t1_1(1)
 1661                     + t2_2(2, 1) * t1_2(1) * t2_3(1, 2) * t1_1(1)
 1662                     + t2_2(0, 2) * t1_2(1) * t2_3(1, 0) * t1_1(2)
 1663                     + t2_2(1, 2) * t1_2(1) * t2_3(1, 1) * t1_1(2)
 1664                     + t2_2(2, 2) * t1_2(1) * t2_3(1, 2) * t1_1(2)
 1665                     + t2_2(0, 0) * t1_2(2) * t2_3(2, 0) * t1_1(0)
 1666                     + t2_2(1, 0) * t1_2(2) * t2_3(2, 1) * t1_1(0)
 1667                     + t2_2(2, 0) * t1_2(2) * t2_3(2, 2) * t1_1(0)
 1668                     + t2_2(0, 1) * t1_2(2) * t2_3(2, 0) * t1_1(1)
 1669                     + t2_2(1, 1) * t1_2(2) * t2_3(2, 1) * t1_1(1)
 1670                     + t2_2(2, 1) * t1_2(2) * t2_3(2, 2) * t1_1(1)
 1671                     + t2_2(0, 2) * t1_2(2) * t2_3(2, 0) * t1_1(2)
 1672                     + t2_2(1, 2) * t1_2(2) * t2_3(2, 1) * t1_1(2)
 1673                     + t2_2(2, 2) * t1_2(2) * t2_3(2, 2) * t1_1(2)),
 1674                "T3(i,j,k)*T3(k,i,j)");
 1675 
 1677                  - (t2_2(0, 0) * t1_2(0) * t3dg_2(0, 0, 0)
 1678                     + t2_2(1, 0) * t1_2(0) * t3dg_2(1, 0, 0)
 1679                     + t2_2(2, 0) * t1_2(0) * t3dg_2(2, 0, 0)
 1680                     + t2_2(0, 1) * t1_2(0) * t3dg_2(0, 1, 0)
 1681                     + t2_2(1, 1) * t1_2(0) * t3dg_2(1, 1, 0)
 1682                     + t2_2(2, 1) * t1_2(0) * t3dg_2(2, 1, 0)
 1683                     + t2_2(0, 2) * t1_2(0) * t3dg_2(0, 2, 0)
 1684                     + t2_2(1, 2) * t1_2(0) * t3dg_2(1, 2, 0)
 1685                     + t2_2(2, 2) * t1_2(0) * t3dg_2(2, 2, 0)
 1686                     + t2_2(0, 0) * t1_2(1) * t3dg_2(0, 0, 1)
 1687                     + t2_2(1, 0) * t1_2(1) * t3dg_2(1, 0, 1)
 1688                     + t2_2(2, 0) * t1_2(1) * t3dg_2(2, 0, 1)
 1689                     + t2_2(0, 1) * t1_2(1) * t3dg_2(0, 1, 1)
 1690                     + t2_2(1, 1) * t1_2(1) * t3dg_2(1, 1, 1)
 1691                     + t2_2(2, 1) * t1_2(1) * t3dg_2(2, 1, 1)
 1692                     + t2_2(0, 2) * t1_2(1) * t3dg_2(0, 2, 1)
 1693                     + t2_2(1, 2) * t1_2(1) * t3dg_2(1, 2, 1)
 1694                     + t2_2(2, 2) * t1_2(1) * t3dg_2(2, 2, 1)
 1695                     + t2_2(0, 0) * t1_2(2) * t3dg_2(0, 0, 2)
 1696                     + t2_2(1, 0) * t1_2(2) * t3dg_2(1, 0, 2)
 1697                     + t2_2(2, 0) * t1_2(2) * t3dg_2(2, 0, 2)
 1698                     + t2_2(0, 1) * t1_2(2) * t3dg_2(0, 1, 2)
 1699                     + t2_2(1, 1) * t1_2(2) * t3dg_2(1, 1, 2)
 1700                     + t2_2(2, 1) * t1_2(2) * t3dg_2(2, 1, 2)
 1701                     + t2_2(0, 2) * t1_2(2) * t3dg_2(0, 2, 2)
 1702                     + t2_2(1, 2) * t1_2(2) * t3dg_2(1, 2, 2)
 1703                     + t2_2(2, 2) * t1_2(2) * t3dg_2(2, 2, 2)),
 1704                "T3(i,j,k)*T3dg(i,j,k)");
 1705 
 1707                  - (t2_2(0, 0) * t1_2(0) * t3dg_2(0, 0, 0)
 1708                     + t2_2(1, 0) * t1_2(0) * t3dg_2(0, 1, 0)
 1709                     + t2_2(2, 0) * t1_2(0) * t3dg_2(0, 2, 0)
 1710                     + t2_2(0, 1) * t1_2(0) * t3dg_2(1, 0, 0)
 1711                     + t2_2(1, 1) * t1_2(0) * t3dg_2(1, 1, 0)
 1712                     + t2_2(2, 1) * t1_2(0) * t3dg_2(1, 2, 0)
 1713                     + t2_2(0, 2) * t1_2(0) * t3dg_2(2, 0, 0)
 1714                     + t2_2(1, 2) * t1_2(0) * t3dg_2(2, 1, 0)
 1715                     + t2_2(2, 2) * t1_2(0) * t3dg_2(2, 2, 0)
 1716                     + t2_2(0, 0) * t1_2(1) * t3dg_2(0, 0, 1)
 1717                     + t2_2(1, 0) * t1_2(1) * t3dg_2(0, 1, 1)
 1718                     + t2_2(2, 0) * t1_2(1) * t3dg_2(0, 2, 1)
 1719                     + t2_2(0, 1) * t1_2(1) * t3dg_2(1, 0, 1)
 1720                     + t2_2(1, 1) * t1_2(1) * t3dg_2(1, 1, 1)
 1721                     + t2_2(2, 1) * t1_2(1) * t3dg_2(1, 2, 1)
 1722                     + t2_2(0, 2) * t1_2(1) * t3dg_2(2, 0, 1)
 1723                     + t2_2(1, 2) * t1_2(1) * t3dg_2(2, 1, 1)
 1724                     + t2_2(2, 2) * t1_2(1) * t3dg_2(2, 2, 1)
 1725                     + t2_2(0, 0) * t1_2(2) * t3dg_2(0, 0, 2)
 1726                     + t2_2(1, 0) * t1_2(2) * t3dg_2(0, 1, 2)
 1727                     + t2_2(2, 0) * t1_2(2) * t3dg_2(0, 2, 2)
 1728                     + t2_2(0, 1) * t1_2(2) * t3dg_2(1, 0, 2)
 1729                     + t2_2(1, 1) * t1_2(2) * t3dg_2(1, 1, 2)
 1730                     + t2_2(2, 1) * t1_2(2) * t3dg_2(1, 2, 2)
 1731                     + t2_2(0, 2) * t1_2(2) * t3dg_2(2, 0, 2)
 1732                     + t2_2(1, 2) * t1_2(2) * t3dg_2(2, 1, 2)
 1733                     + t2_2(2, 2) * t1_2(2) * t3dg_2(2, 2, 2)),
 1734                "T3(i,j,k)*T3dg(j,i,k)");
 1735 
 1737                  - (t2_2(0, 0) * t1_2(0) * t3dg_2(0, 0, 0)
 1738                     + t2_2(1, 0) * t1_2(0) * t3dg_2(1, 0, 0)
 1739                     + t2_2(2, 0) * t1_2(0) * t3dg_2(2, 0, 0)
 1740                     + t2_2(0, 1) * t1_2(0) * t3dg_2(0, 0, 1)
 1741                     + t2_2(1, 1) * t1_2(0) * t3dg_2(1, 0, 1)
 1742                     + t2_2(2, 1) * t1_2(0) * t3dg_2(2, 0, 1)
 1743                     + t2_2(0, 2) * t1_2(0) * t3dg_2(0, 0, 2)
 1744                     + t2_2(1, 2) * t1_2(0) * t3dg_2(1, 0, 2)
 1745                     + t2_2(2, 2) * t1_2(0) * t3dg_2(2, 0, 2)
 1746                     + t2_2(0, 0) * t1_2(1) * t3dg_2(0, 1, 0)
 1747                     + t2_2(1, 0) * t1_2(1) * t3dg_2(1, 1, 0)
 1748                     + t2_2(2, 0) * t1_2(1) * t3dg_2(2, 1, 0)
 1749                     + t2_2(0, 1) * t1_2(1) * t3dg_2(0, 1, 1)
 1750                     + t2_2(1, 1) * t1_2(1) * t3dg_2(1, 1, 1)
 1751                     + t2_2(2, 1) * t1_2(1) * t3dg_2(2, 1, 1)
 1752                     + t2_2(0, 2) * t1_2(1) * t3dg_2(0, 1, 2)
 1753                     + t2_2(1, 2) * t1_2(1) * t3dg_2(1, 1, 2)
 1754                     + t2_2(2, 2) * t1_2(1) * t3dg_2(2, 1, 2)
 1755                     + t2_2(0, 0) * t1_2(2) * t3dg_2(0, 2, 0)
 1756                     + t2_2(1, 0) * t1_2(2) * t3dg_2(1, 2, 0)
 1757                     + t2_2(2, 0) * t1_2(2) * t3dg_2(2, 2, 0)
 1758                     + t2_2(0, 1) * t1_2(2) * t3dg_2(0, 2, 1)
 1759                     + t2_2(1, 1) * t1_2(2) * t3dg_2(1, 2, 1)
 1760                     + t2_2(2, 1) * t1_2(2) * t3dg_2(2, 2, 1)
 1761                     + t2_2(0, 2) * t1_2(2) * t3dg_2(0, 2, 2)
 1762                     + t2_2(1, 2) * t1_2(2) * t3dg_2(1, 2, 2)
 1763                     + t2_2(2, 2) * t1_2(2) * t3dg_2(2, 2, 2)),
 1764                "T3(i,j,k)*T3dg(i,k,j)");
 1765 
 1767                  - (t2_2(0, 0) * t1_2(0) * t3dg_2(0, 0, 0)
 1768                     + t2_2(1, 0) * t1_2(0) * t3dg_2(0, 0, 1)
 1769                     + t2_2(2, 0) * t1_2(0) * t3dg_2(0, 0, 2)
 1770                     + t2_2(0, 1) * t1_2(0) * t3dg_2(1, 0, 0)
 1771                     + t2_2(1, 1) * t1_2(0) * t3dg_2(1, 0, 1)
 1772                     + t2_2(2, 1) * t1_2(0) * t3dg_2(1, 0, 2)
 1773                     + t2_2(0, 2) * t1_2(0) * t3dg_2(2, 0, 0)
 1774                     + t2_2(1, 2) * t1_2(0) * t3dg_2(2, 0, 1)
 1775                     + t2_2(2, 2) * t1_2(0) * t3dg_2(2, 0, 2)
 1776                     + t2_2(0, 0) * t1_2(1) * t3dg_2(0, 1, 0)
 1777                     + t2_2(1, 0) * t1_2(1) * t3dg_2(0, 1, 1)
 1778                     + t2_2(2, 0) * t1_2(1) * t3dg_2(0, 1, 2)
 1779                     + t2_2(0, 1) * t1_2(1) * t3dg_2(1, 1, 0)
 1780                     + t2_2(1, 1) * t1_2(1) * t3dg_2(1, 1, 1)
 1781                     + t2_2(2, 1) * t1_2(1) * t3dg_2(1, 1, 2)
 1782                     + t2_2(0, 2) * t1_2(1) * t3dg_2(2, 1, 0)
 1783                     + t2_2(1, 2) * t1_2(1) * t3dg_2(2, 1, 1)
 1784                     + t2_2(2, 2) * t1_2(1) * t3dg_2(2, 1, 2)
 1785                     + t2_2(0, 0) * t1_2(2) * t3dg_2(0, 2, 0)
 1786                     + t2_2(1, 0) * t1_2(2) * t3dg_2(0, 2, 1)
 1787                     + t2_2(2, 0) * t1_2(2) * t3dg_2(0, 2, 2)
 1788                     + t2_2(0, 1) * t1_2(2) * t3dg_2(1, 2, 0)
 1789                     + t2_2(1, 1) * t1_2(2) * t3dg_2(1, 2, 1)
 1790                     + t2_2(2, 1) * t1_2(2) * t3dg_2(1, 2, 2)
 1791                     + t2_2(0, 2) * t1_2(2) * t3dg_2(2, 2, 0)
 1792                     + t2_2(1, 2) * t1_2(2) * t3dg_2(2, 2, 1)
 1793                     + t2_2(2, 2) * t1_2(2) * t3dg_2(2, 2, 2)),
 1794                "T3(i,j,k)*T3dg(j,k,i)");
 1795 
 1797                  - (t2_2(0, 0) * t1_2(0) * t3dg_2(0, 0, 0)
 1798                     + t2_2(1, 0) * t1_2(0) * t3dg_2(0, 0, 1)
 1799                     + t2_2(2, 0) * t1_2(0) * t3dg_2(0, 0, 2)
 1800                     + t2_2(0, 1) * t1_2(0) * t3dg_2(0, 1, 0)
 1801                     + t2_2(1, 1) * t1_2(0) * t3dg_2(0, 1, 1)
 1802                     + t2_2(2, 1) * t1_2(0) * t3dg_2(0, 1, 2)
 1803                     + t2_2(0, 2) * t1_2(0) * t3dg_2(0, 2, 0)
 1804                     + t2_2(1, 2) * t1_2(0) * t3dg_2(0, 2, 1)
 1805                     + t2_2(2, 2) * t1_2(0) * t3dg_2(0, 2, 2)
 1806                     + t2_2(0, 0) * t1_2(1) * t3dg_2(1, 0, 0)
 1807                     + t2_2(1, 0) * t1_2(1) * t3dg_2(1, 0, 1)
 1808                     + t2_2(2, 0) * t1_2(1) * t3dg_2(1, 0, 2)
 1809                     + t2_2(0, 1) * t1_2(1) * t3dg_2(1, 1, 0)
 1810                     + t2_2(1, 1) * t1_2(1) * t3dg_2(1, 1, 1)
 1811                     + t2_2(2, 1) * t1_2(1) * t3dg_2(1, 1, 2)
 1812                     + t2_2(0, 2) * t1_2(1) * t3dg_2(1, 2, 0)
 1813                     + t2_2(1, 2) * t1_2(1) * t3dg_2(1, 2, 1)
 1814                     + t2_2(2, 2) * t1_2(1) * t3dg_2(1, 2, 2)
 1815                     + t2_2(0, 0) * t1_2(2) * t3dg_2(2, 0, 0)
 1816                     + t2_2(1, 0) * t1_2(2) * t3dg_2(2, 0, 1)
 1817                     + t2_2(2, 0) * t1_2(2) * t3dg_2(2, 0, 2)
 1818                     + t2_2(0, 1) * t1_2(2) * t3dg_2(2, 1, 0)
 1819                     + t2_2(1, 1) * t1_2(2) * t3dg_2(2, 1, 1)
 1820                     + t2_2(2, 1) * t1_2(2) * t3dg_2(2, 1, 2)
 1821                     + t2_2(0, 2) * t1_2(2) * t3dg_2(2, 2, 0)
 1822                     + t2_2(1, 2) * t1_2(2) * t3dg_2(2, 2, 1)
 1823                     + t2_2(2, 2) * t1_2(2) * t3dg_2(2, 2, 2)),
 1824                "T3(i,j,k)*T3dg(k,j,i)");
 1825 
 1827                  - (t2_2(0, 0) * t1_2(0) * t3dg_2(0, 0, 0)
 1828                     + t2_2(1, 0) * t1_2(0) * t3dg_2(0, 1, 0)
 1829                     + t2_2(2, 0) * t1_2(0) * t3dg_2(0, 2, 0)
 1830                     + t2_2(0, 1) * t1_2(0) * t3dg_2(0, 0, 1)
 1831                     + t2_2(1, 1) * t1_2(0) * t3dg_2(0, 1, 1)
 1832                     + t2_2(2, 1) * t1_2(0) * t3dg_2(0, 2, 1)
 1833                     + t2_2(0, 2) * t1_2(0) * t3dg_2(0, 0, 2)
 1834                     + t2_2(1, 2) * t1_2(0) * t3dg_2(0, 1, 2)
 1835                     + t2_2(2, 2) * t1_2(0) * t3dg_2(0, 2, 2)
 1836                     + t2_2(0, 0) * t1_2(1) * t3dg_2(1, 0, 0)
 1837                     + t2_2(1, 0) * t1_2(1) * t3dg_2(1, 1, 0)
 1838                     + t2_2(2, 0) * t1_2(1) * t3dg_2(1, 2, 0)
 1839                     + t2_2(0, 1) * t1_2(1) * t3dg_2(1, 0, 1)
 1840                     + t2_2(1, 1) * t1_2(1) * t3dg_2(1, 1, 1)
 1841                     + t2_2(2, 1) * t1_2(1) * t3dg_2(1, 2, 1)
 1842                     + t2_2(0, 2) * t1_2(1) * t3dg_2(1, 0, 2)
 1843                     + t2_2(1, 2) * t1_2(1) * t3dg_2(1, 1, 2)
 1844                     + t2_2(2, 2) * t1_2(1) * t3dg_2(1, 2, 2)
 1845                     + t2_2(0, 0) * t1_2(2) * t3dg_2(2, 0, 0)
 1846                     + t2_2(1, 0) * t1_2(2) * t3dg_2(2, 1, 0)
 1847                     + t2_2(2, 0) * t1_2(2) * t3dg_2(2, 2, 0)
 1848                     + t2_2(0, 1) * t1_2(2) * t3dg_2(2, 0, 1)
 1849                     + t2_2(1, 1) * t1_2(2) * t3dg_2(2, 1, 1)
 1850                     + t2_2(2, 1) * t1_2(2) * t3dg_2(2, 2, 1)
 1851                     + t2_2(0, 2) * t1_2(2) * t3dg_2(2, 0, 2)
 1852                     + t2_2(1, 2) * t1_2(2) * t3dg_2(2, 1, 2)
 1853                     + t2_2(2, 2) * t1_2(2) * t3dg_2(2, 2, 2)),
 1854                "T3(i,j,k)*T3dg(k,i,j)");
 1855}
FTensor::Index< 'i', SPACE_DIM > i
FTensor::Index< 'l', 3 > l
FTensor::Index< 'j', 3 > j
FTensor::Index< 'k', 3 > k
void test_for_zero(const T &t, const std::string &s)