Nothing
test_that(
"Zonal = 0 if l(lambda)>l(x)", {
# numeric
expect_equal(ZonalR(c(1,2), c(3,2,1)), 0)
expect_equal(ZonalR(c(1,2), c(3,2,1), algorithm = "naive"), 0)
# gmp
x <- as.bigq(c(1L,2L))
lambda <- c(3,2,1)
expect_identical(ZonalR(x, lambda), as.bigq(0L))
expect_identical(ZonalR(x, lambda, algorithm = "naive"), as.bigq(0L))
# polynomial
n <- 2
lambda <- c(3,2,1)
expect_true(ZonalPolR(n, lambda) == as.qspray(0))
expect_identical(ZonalPolR(n, lambda, algorithm = "naive"),
as.qspray(0))
expect_identical(ZonalPolR(n, lambda, exact = FALSE, algorithm = "naive"),
mvp::constant(0))
expect_identical(ZonalPolR(n, lambda, algorithm = "naive", basis = "MSF"),
as.qspray(0))
expect_identical(ZonalPolR(n, lambda, exact = FALSE, algorithm = "naive",
basis = "MSF"),
mvp::constant(0))
}
)
test_that(
"Zonal polynomials sum to the trace - gmp", {
x <- as.bigq(c(1L,2L,4L,7L), c(2L,3L,1L,2L))
expected <- sum(x)^3
obtained_DK <- ZonalR(x, 3) + ZonalR(x, c(2,1)) + ZonalR(x, c(1,1,1))
obtained_naive <- ZonalR(x, 3, algorithm = "naive") +
ZonalR(x, c(2,1), algorithm = "naive") +
ZonalR(x, c(1,1,1), algorithm = "naive")
expect_identical(obtained_DK, expected)
expect_identical(obtained_naive, expected)
}
)
test_that(
"Zonal polynomials sum to the trace - numeric", {
x <- c(1,2,4,7) / c(2,3,1,2)
expected <- sum(x)^3
obtained_DK <- ZonalR(x, 3) + ZonalR(x, c(2,1)) + ZonalR(x, c(1,1,1))
obtained_naive <- ZonalR(x, 3, algorithm = "naive") +
ZonalR(x, c(2,1), algorithm = "naive") +
ZonalR(x, c(1,1,1), algorithm = "naive")
expect_equal(obtained_DK, expected)
expect_equal(obtained_naive, expected)
}
)
test_that(
"Zonal polynomials sum to the trace - complex", {
x <- c(1i, 2+1i, 4-3i, 7-3i) / c(2,3,1,2)
expected <- sum(x)^3
obtained_DK <- ZonalR(x, 3) + ZonalR(x, c(2,1)) + ZonalR(x, c(1,1,1))
obtained_naive <- ZonalR(x, 3, algorithm = "naive") +
ZonalR(x, c(2,1), algorithm = "naive") +
ZonalR(x, c(1,1,1), algorithm = "naive")
expect_equal(obtained_DK, expected)
expect_equal(obtained_naive, expected)
}
)
test_that(
"ZonalPol is correct", {
lambda <- c(3,2)
pol <- ZonalPolR(4, lambda, algorithm = "naive")
x <- as.bigq(c(6L,-7L,8L,9L), c(1L,2L,3L,4L))
polEval <- qspray::evalQspray(pol, x)
expect_identical(polEval, ZonalR(as.bigq(x), lambda))
}
)
test_that(
"Zonal polynomials sum to the trace - polynomial", {
n <- 4
expected <- (qlone(1) + qlone(2) + qlone(3) + qlone(4))^3
obtained <- ZonalPolR(n, 3) + ZonalPolR(n, c(2,1)) + ZonalPolR(n, c(1,1,1))
expect_true(expected == obtained)
}
)
test_that(
"ZonalCPP is correct", {
x <- as.bigq(c(6L, -7L, 8L, 9L), c(1L, 2L, 3L, 4L))
lambda <- c(3, 2)
res <- Zonal(x, lambda)
expect_identical(res, ZonalR(x, lambda))
#
x <- c(6, -7, 8, 9) / c(1, 2, 3, 4)
lambda <- c(3, 2)
res <- Zonal(x, lambda)
expect_equal(res, ZonalR(x, lambda))
}
)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.