R Under development (unstable) (2023-12-13 r85679 ucrt) -- "Unsuffered Consequences" Copyright (C) 2023 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 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. > > d <- 4 > qux <- makeH(- diag(d), rep(0, d), rep(1, d), 1) > qux <- d2q(qux) > 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" > > > proc.time() user system elapsed 0.14 0.03 0.15