BOOST_AUTO_TEST_CASE_TEMPLATE( pow, Permutation, ptypes ) { // const auto& p = fixture<Permutation>().p; // const auto& I = fixture<Permutation>().I; // const auto& p102 = fixture<Permutation>().p102; // const auto& p120 = fixture<Permutation>().p120; // const auto& p021 = fixture<Permutation>().p021; // const auto& p201 = fixture<Permutation>().p201; // const auto& p210 = fixture<Permutation>().p210; const Permutation p0231({0,2,3,1}); // Check that powers of permutations are computed correctly BOOST_CHECK_EQUAL(p0231.pow(0), p0231.identity()); BOOST_CHECK_EQUAL(p0231.pow(1), p0231); BOOST_CHECK_EQUAL(p0231.pow(2), p0231 * p0231); BOOST_CHECK_EQUAL(p0231.pow(3), p0231 * p0231 * p0231); BOOST_CHECK_EQUAL(p0231.pow(4), p0231 * p0231 * p0231 * p0231); BOOST_CHECK_EQUAL(p0231.pow(5), p0231 * p0231 * p0231 * p0231 * p0231); BOOST_CHECK_EQUAL(p0231.pow(6), p0231 * p0231 * p0231 * p0231 * p0231 * p0231); BOOST_CHECK_EQUAL(p0231.pow(7), p0231 * p0231 * p0231 * p0231 * p0231 * p0231 * p0231); BOOST_CHECK_EQUAL(p0231.pow(8), p0231 * p0231 * p0231 * p0231 * p0231 * p0231 * p0231 * p0231); // Check that inverse powers of permutations are computed correctly Permutation I0231 = p0231.inv(); BOOST_CHECK_EQUAL(p0231.pow(-1), I0231); BOOST_CHECK_EQUAL(p0231.pow(-2), I0231 * I0231); BOOST_CHECK_EQUAL(p0231.pow(-3), I0231 * I0231 * I0231); BOOST_CHECK_EQUAL(p0231.pow(-4), I0231 * I0231 * I0231 * I0231); BOOST_CHECK_EQUAL(p0231.pow(-5), I0231 * I0231 * I0231 * I0231 * I0231); BOOST_CHECK_EQUAL(p0231.pow(-6), I0231 * I0231 * I0231 * I0231 * I0231 * I0231); BOOST_CHECK_EQUAL(p0231.pow(-7), I0231 * I0231 * I0231 * I0231 * I0231 * I0231 * I0231); BOOST_CHECK_EQUAL(p0231.pow(-8), I0231 * I0231 * I0231 * I0231 * I0231 * I0231 * I0231 * I0231); }