Nothing
context("Rounding test")
library(volesti)
testRound <- function(P, exactvol, tol, name_string, num_of_exps, algo, rotation,seed){
if (rotation) {
P = rotate_polytope(P)$P
listHpoly = round_polytope(P, settings = list("seed" = seed))
} else {
listHpoly = round_polytope(P, settings = list("seed" = seed))
}
vol = 0
for (j in 1:num_of_exps) {
if (algo == "CB") {
vol = vol + listHpoly$round_value * volume(listHpoly$P, settings=list("algorithm"="CB", "error"=0.1, "seed" = seed))
} else {
vol = vol + listHpoly$round_value * volume(listHpoly$P, settings=list("algorithm"="CG", "error"=0.1, "seed" = seed))
}
}
vol = vol / num_of_exps
error = abs(vol - exactvol) / exactvol
if (error >= tol){
res = 0
} else {
res = 1
}
return(res)
}
cran_only = TRUE
for (i in 1:2) {
num_of_exps = 10
test_that("Rounding H-skinny_cube10", {
seed=5
P = gen_skinny_cube(10)
res = testRound(P, 102400, 0.3, 'H-skinny_cube10', num_of_exps, 'CB', FALSE,seed)
expect_equal(res, 1)
})
}
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.