TEST() [10/10]

TEST	(	ResMonoidSparse	,
		mult	)

Definition at line 176 of file ResTest.cpp.

{
  ResMonoidSparse M(4,
                    std::vector<int>{1, 1, 1, 1},
                    std::vector<int>{},
                    MonomialOrderingType::GRevLex);
 
  int exp1[]{0, 0, 0, 0};
  int exp2[]{0, 1, 0, 3};
  int exp3[]{0, 0, 0, 0};
  int mon1[]{0, 0, 0, 0, 0, 0};
  int mon2[]{0, 0, 0, 0, 0, 0, 0, 0, 0};
  int mon3[]{0, 0, 0, 0, 0, 0, 0, 0, 0};
  int mon[]{0, 0, 0, 0, 0, 0, 0, 0, 0};
  int comp;
  M.from_expvector(exp1, 0, mon1);
  M.from_expvector(exp2, 0, mon2);
  M.mult(mon1, mon2, mon);
  M.to_expvector(mon, exp3, comp);
  EXPECT_EQ(3, M.monomial_size(mon1));
  EXPECT_EQ(7, M.monomial_size(mon2));
  EXPECT_EQ(7, M.monomial_size(mon));
 
  for (int i = 0; i < M.n_vars(); i++)
    for (int j = i; j < M.n_vars(); j++)
      for (int k = 0; k < M.n_vars(); k++)
        {
          int exp1[]{0, 0, 0, 0};
          exp1[i]++;
          exp1[j]++;
          int exp2[]{0, 0, 0, 0};
          exp2[k]++;
          M.from_expvector(exp1, 0, mon1);
          M.from_expvector(exp2, 0, mon2);
          M.mult(mon1, mon2, mon);
          M.mult(mon1, mon2, mon3);
          EXPECT_TRUE(M.is_equal(mon, mon3));
          M.to_expvector(mon, exp3, comp);
          int exp3a[]{0, 0, 0, 0};
          exp3a[i]++;
          exp3a[j]++;
          exp3a[k]++;
          for (int ell = 0; ell < M.n_vars(); ell++)
            EXPECT_EQ(exp3[ell], exp3a[ell]);
          EXPECT_EQ(M.monomial_size(mon),
                    M.monomial_size(mon1) + M.monomial_size(mon2) - 3);
        }
}

References GRevLex.