Line data    Source code

       1             : #include "Set/Set.H"
       2             : namespace Test
       3             : {
       4             : namespace Set
       5             : {
       6             : template <int dim, ::Set::Sym sym>
       7             : class Matrix4
       8             : {
       9             : public:
      10           6 :     int SymmetryTest(int /*verbose*/)
      11             :     {
      12           6 :         ::Set::Matrix4<dim,sym> matrix = ::Set::Matrix4<dim,sym>::Randomize();
      13             : 
      14             :         //Util::Message(INFO,matrix);
      15             :         //std::cout << matrix << std::endl;
      16             :         if (sym == ::Set::Sym::None)
      17             :             return 0;
      18             :         if (sym == ::Set::Sym::MajorMinor)
      19             :         {
      20           8 :             for (int i = 0; i < dim; i++)
      21          24 :             for (int j = 0; j < dim; j++)
      22          72 :             for (int k = 0; k < dim; k++)
      23         216 :             for (int l = 0; l < dim; l++)
      24             :             {
      25         324 :                 if (matrix(i,j,k,l) != matrix(k,l,i,j)) return 1; // breaks major symmetry
      26         324 :                 if (matrix(i,j,k,l) != matrix(i,j,l,k)) return 1; // breaks minor symmetry
      27         324 :                 if (matrix(i,j,k,l) != matrix(j,i,k,l)) return 1; // breaks minor symmetry
      28             :             }
      29           2 :             return 0;
      30             :         }
      31             :         if (sym == ::Set::Sym::Full)
      32             :         {
      33          14 :             for (int i = 0; i < dim; i++)
      34          36 :             for (int j = 0; j < dim; j++)
      35          96 :             for (int k = 0; k < dim; k++)
      36         264 :             for (int l = 0; l < dim; l++)
      37             :             {
      38         388 :                 if (matrix(i,j,k,l) != matrix(i,j,l,k)) return 1;
      39         388 :                 if (matrix(i,j,k,l) != matrix(i,k,j,l)) return 1;
      40         388 :                 if (matrix(i,j,k,l) != matrix(i,k,l,j)) return 1;
      41         388 :                 if (matrix(i,j,k,l) != matrix(i,l,j,k)) return 1;
      42         388 :                 if (matrix(i,j,k,l) != matrix(i,l,k,j)) return 1;
      43         388 :                 if (matrix(i,j,k,l) != matrix(j,i,k,l)) return 1;
      44         388 :                 if (matrix(i,j,k,l) != matrix(j,i,l,k)) return 1;
      45         388 :                 if (matrix(i,j,k,l) != matrix(j,k,i,l)) return 1;
      46         388 :                 if (matrix(i,j,k,l) != matrix(j,k,l,i)) return 1;
      47         388 :                 if (matrix(i,j,k,l) != matrix(j,l,i,k)) return 1;
      48         388 :                 if (matrix(i,j,k,l) != matrix(j,l,k,i)) return 1;
      49         388 :                 if (matrix(i,j,k,l) != matrix(k,i,j,l)) return 1;
      50         388 :                 if (matrix(i,j,k,l) != matrix(k,i,l,j)) return 1;
      51         388 :                 if (matrix(i,j,k,l) != matrix(k,j,i,l)) return 1;
      52         388 :                 if (matrix(i,j,k,l) != matrix(k,j,l,i)) return 1;
      53         388 :                 if (matrix(i,j,k,l) != matrix(k,l,i,j)) return 1;
      54         388 :                 if (matrix(i,j,k,l) != matrix(k,l,j,i)) return 1;
      55         388 :                 if (matrix(i,j,k,l) != matrix(l,i,j,k)) return 1;
      56         388 :                 if (matrix(i,j,k,l) != matrix(l,i,k,j)) return 1;
      57         388 :                 if (matrix(i,j,k,l) != matrix(l,j,i,k)) return 1;
      58         388 :                 if (matrix(i,j,k,l) != matrix(l,j,k,i)) return 1;
      59         388 :                 if (matrix(i,j,k,l) != matrix(l,k,i,j)) return 1;
      60         388 :                 if (matrix(i,j,k,l) != matrix(l,k,j,i)) return 1;
      61             :             }
      62           4 :             return 0;
      63             :         }
      64             : 
      65             : 
      66             :         return 1;
      67             :     }
      68             : };
      69             : }
      70             : }