17    return a(1, 0, 0, 1) * b(1, 0, 0, 1) + 
a(2, 0, 0, 1) * b(2, 0, 0, 1)
 
   18           + 
a(0, 1, 0, 1) * b(0, 1, 0, 1) + 
a(2, 1, 0, 1) * b(2, 1, 0, 1)
 
   19           + 
a(0, 2, 0, 1) * b(0, 2, 0, 1) + 
a(1, 2, 0, 1) * b(1, 2, 0, 1)
 
   20           + 
a(1, 0, 0, 2) * b(1, 0, 0, 2) + 
a(2, 0, 0, 2) * b(2, 0, 0, 2)
 
   21           + 
a(0, 1, 0, 2) * b(0, 1, 0, 2) + 
a(2, 1, 0, 2) * b(2, 1, 0, 2)
 
   22           + 
a(0, 2, 0, 2) * b(0, 2, 0, 2) + 
a(1, 2, 0, 2) * b(1, 2, 0, 2)
 
   23           + 
a(1, 0, 1, 0) * b(1, 0, 1, 0) + 
a(2, 0, 1, 0) * b(2, 0, 1, 0)
 
   24           + 
a(0, 1, 1, 0) * b(0, 1, 1, 0) + 
a(2, 1, 1, 0) * b(2, 1, 1, 0)
 
   25           + 
a(0, 2, 1, 0) * b(0, 2, 1, 0) + 
a(1, 2, 1, 0) * b(1, 2, 1, 0)
 
   26           + 
a(1, 0, 1, 2) * b(1, 0, 1, 2) + 
a(2, 0, 1, 2) * b(2, 0, 1, 2)
 
   27           + 
a(0, 1, 1, 2) * b(0, 1, 1, 2) + 
a(2, 1, 1, 2) * b(2, 1, 1, 2)
 
   28           + 
a(0, 2, 1, 2) * b(0, 2, 1, 2) + 
a(1, 2, 1, 2) * b(1, 2, 1, 2)
 
   29           + 
a(1, 0, 2, 0) * b(1, 0, 2, 0) + 
a(2, 0, 2, 0) * b(2, 0, 2, 0)
 
   30           + 
a(0, 1, 2, 0) * b(0, 1, 2, 0) + 
a(2, 1, 2, 0) * b(2, 1, 2, 0)
 
   31           + 
a(0, 2, 2, 0) * b(0, 2, 2, 0) + 
a(1, 2, 2, 0) * b(1, 2, 2, 0)
 
   32           + 
a(1, 0, 2, 1) * b(1, 0, 2, 1) + 
a(2, 0, 2, 1) * b(2, 0, 2, 1)
 
   33           + 
a(0, 1, 2, 1) * b(0, 1, 2, 1) + 
a(2, 1, 2, 1) * b(2, 1, 2, 1)
 
   34           + 
a(0, 2, 2, 1) * b(0, 2, 2, 1) + 
a(1, 2, 2, 1) * b(1, 2, 2, 1);
 
 
   55    return a(1, 0, 0, 1) * b(1, 0, 0, 1) + 
a(2, 0, 0, 1) * b(2, 0, 0, 1)
 
   56           + 
a(0, 1, 0, 1) * b(0, 0, 1, 1) + 
a(2, 1, 0, 1) * b(2, 0, 1, 1)
 
   57           + 
a(0, 2, 0, 1) * b(0, 0, 2, 1) + 
a(1, 2, 0, 1) * b(1, 0, 2, 1)
 
   58           + 
a(1, 0, 0, 2) * b(1, 0, 0, 2) + 
a(2, 0, 0, 2) * b(2, 0, 0, 2)
 
   59           + 
a(0, 1, 0, 2) * b(0, 0, 1, 2) + 
a(2, 1, 0, 2) * b(2, 0, 1, 2)
 
   60           + 
a(0, 2, 0, 2) * b(0, 0, 2, 2) + 
a(1, 2, 0, 2) * b(1, 0, 2, 2)
 
   61           + 
a(1, 0, 1, 0) * b(1, 1, 0, 0) + 
a(2, 0, 1, 0) * b(2, 1, 0, 0)
 
   62           + 
a(0, 1, 1, 0) * b(0, 1, 1, 0) + 
a(2, 1, 1, 0) * b(2, 1, 1, 0)
 
   63           + 
a(0, 2, 1, 0) * b(0, 1, 2, 0) + 
a(1, 2, 1, 0) * b(1, 1, 2, 0)
 
   64           + 
a(1, 0, 1, 2) * b(1, 1, 0, 2) + 
a(2, 0, 1, 2) * b(2, 1, 0, 2)
 
   65           + 
a(0, 1, 1, 2) * b(0, 1, 1, 2) + 
a(2, 1, 1, 2) * b(2, 1, 1, 2)
 
   66           + 
a(0, 2, 1, 2) * b(0, 1, 2, 2) + 
a(1, 2, 1, 2) * b(1, 1, 2, 2)
 
   67           + 
a(1, 0, 2, 0) * b(1, 2, 0, 0) + 
a(2, 0, 2, 0) * b(2, 2, 0, 0)
 
   68           + 
a(0, 1, 2, 0) * b(0, 2, 1, 0) + 
a(2, 1, 2, 0) * b(2, 2, 1, 0)
 
   69           + 
a(0, 2, 2, 0) * b(0, 2, 2, 0) + 
a(1, 2, 2, 0) * b(1, 2, 2, 0)
 
   70           + 
a(1, 0, 2, 1) * b(1, 2, 0, 1) + 
a(2, 0, 2, 1) * b(2, 2, 0, 1)
 
   71           + 
a(0, 1, 2, 1) * b(0, 2, 1, 1) + 
a(2, 1, 2, 1) * b(2, 2, 1, 1)
 
   72           + 
a(0, 2, 2, 1) * b(0, 2, 2, 1) + 
a(1, 2, 2, 1) * b(1, 2, 2, 1);