test_that("Modified Macdonald polynomial mu (q, t) = Modified Macdonald polynomial mu' (t, q)", {
n <- 4
mu <- c(2, 1, 1)
mup <- c(3, 1)
macHpoly <- modifiedMacdonaldPol(n, mu)
macHpolyp <- changeParameters(
modifiedMacdonaldPol(n, mup),
list(qlone(2), qlone(1))
)
expect_true(macHpoly == macHpolyp)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.