R Under development (unstable) (2023-12-12 r85669 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. > ## rmatio, a R interface to the C library matio, MAT File I/O Library. > ## Copyright (C) 2013-2023 Stefan Widgren > ## > ## This program is free software: you can redistribute it and/or modify > ## it under the terms of the GNU General Public License as published by > ## the Free Software Foundation, either version 3 of the License, or > ## (at your option) any later version. > ## > ## rmatio is distributed in the hope that it will be useful, > ## but WITHOUT ANY WARRANTY; without even the implied warranty of > ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the > ## GNU General Public License for more details. > ## > ## You should have received a copy of the GNU General Public License > ## along with this program. If not, see . > > library(rmatio) > > ## For debugging > sessionInfo() R Under development (unstable) (2023-12-12 r85669 ucrt) Platform: x86_64-w64-mingw32/x64 Running under: Windows Server 2022 x64 (build 20348) Matrix products: default locale: [1] LC_COLLATE=C LC_CTYPE=German_Germany.utf8 [3] LC_MONETARY=C LC_NUMERIC=C [5] LC_TIME=C time zone: Europe/Berlin tzcode source: internal attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] rmatio_0.19.0 loaded via a namespace (and not attached): [1] compiler_4.4.0 Matrix_1.6-4 grid_4.4.0 lattice_0.22-5 > > ## > ## Check write and read of vector in MAT5 format: > ## 1) without compression > ## 2) with compression > ## > > ## > ## vector: case-1 > ## > a1_exp <- 1:5 > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a1_exp), + filename = filename, + compression = FALSE, + version = "MAT5") > a1_obs <- read.mat(filename)[["a"]] > unlink(filename) > storage.mode(a1_obs) <- "integer" > str(a1_obs) int [1:5] 1 2 3 4 5 > stopifnot(identical(a1_obs, a1_exp)) > > # Run the same test with compression > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a1_exp), + filename = filename, + compression = TRUE, + version = "MAT5") > a1_zlib_obs <- read.mat(filename)[["a"]] > unlink(filename) > storage.mode(a1_zlib_obs) <- "integer" > str(a1_zlib_obs) int [1:5] 1 2 3 4 5 > stopifnot(identical(a1_zlib_obs, a1_exp)) > > ## > ## vector: case-2 > ## > a2_exp <- c(1, 2, 3, 4, 5) > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a2_exp), + filename = filename, + compression = FALSE, + version = "MAT5") > a2_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a2_obs) num [1:5] 1 2 3 4 5 > stopifnot(identical(a2_obs, a2_exp)) > > # Run the same test with compression > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a2_exp), + filename = filename, + compression = TRUE, + version = "MAT5") > a2_zlib_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a2_zlib_obs) num [1:5] 1 2 3 4 5 > stopifnot(identical(a2_zlib_obs, a2_exp)) > > ## > ## vector: case-3 > ## > a3_exp <- 1 > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a3_exp), + filename = filename, + compression = FALSE, + version = "MAT5") > a3_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a3_obs) num 1 > stopifnot(identical(a3_obs, a3_exp)) > > # Run the same test with compression > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a3_exp), + filename = filename, + compression = TRUE, + version = "MAT5") > a3_zlib_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a3_zlib_obs) num 1 > stopifnot(identical(a3_zlib_obs, a3_exp)) > > proc.time() user system elapsed 0.79 0.14 0.93