R Under development (unstable) (2023-04-23 r84305 ucrt) -- "Unsuffered Consequences" Copyright (C) 2023 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 (64-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > > library(rcdd) If you want correct answers, use rational arithmetic. See the Warnings sections in help pages for functions that do computational geometry. > > # needed because now uses R RNG for matrix row permutation > set.seed(42) > > d <- 4 > qux <- makeH(- diag(d), rep(0, d), rep(1, d), 1) > print(qux) [,1] [,2] [,3] [,4] [,5] [,6] [1,] 1 1 -1 -1 -1 -1 [2,] 0 0 1 0 0 0 [3,] 0 0 0 1 0 0 [4,] 0 0 0 0 1 0 [5,] 0 0 0 0 0 1 attr(,"representation") [1] "H" > > out <- scdd(qux) > print(out) $output [,1] [,2] [,3] [,4] [,5] [,6] [1,] 0 1 0 0 0 1 [2,] 0 1 0 0 1 0 [3,] 0 1 0 1 0 0 [4,] 0 1 1 0 0 0 attr(,"representation") [1] "V" > > out <- scdd(out$output) > print(out) $output [,1] [,2] [,3] [,4] [,5] [,6] [1,] 0 1 -1 -1 -1 0 [2,] 0 0 1 0 0 0 [3,] 0 0 0 1 0 0 [4,] 0 0 0 0 1 0 [5,] 1 -1 1 1 1 1 attr(,"representation") [1] "H" > > out <- scdd(qux, roworder = "randomrow") > print(out) $output [,1] [,2] [,3] [,4] [,5] [,6] [1,] 0 1 0 0 1 0 [2,] 0 1 1 0 0 0 [3,] 0 1 0 1 0 0 [4,] 0 1 0 0 0 1 attr(,"representation") [1] "V" > > out <- scdd(qux, roworder = "maxcutoff") > print(out) $output [,1] [,2] [,3] [,4] [,5] [,6] [1,] 0 1 1 0 0 0 [2,] 0 1 0 1 0 0 [3,] 0 1 0 0 1 0 [4,] 0 1 0 0 0 1 attr(,"representation") [1] "V" > > > proc.time() user system elapsed 0.12 0.06 0.18