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 <http://www.gnu.org/licenses/>.
> 
> 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 array in MAT5 format:
> ## 1) without compression
> ## 2) with compression
> ##
> 
> ##
> ## logical: case-1
> ##
> a1_exp <- array(c(TRUE,  TRUE,  TRUE,  TRUE,  TRUE,
+                   FALSE, TRUE,  TRUE,  TRUE,  TRUE,
+                   FALSE, FALSE, TRUE,  TRUE,  TRUE,
+                   FALSE, FALSE, FALSE, TRUE,  TRUE,
+                   FALSE, FALSE, FALSE, FALSE, TRUE),
+                 c(5L, 5L))
> filename <- tempfile(fileext = ".mat")
> write.mat(list(a = a1_exp), filename = filename, compression = FALSE,
+           version = "MAT5")
> a1_obs <- read.mat(filename)[["a"]]
> unlink(filename)
> str(a1_obs)
 logi [1:5, 1:5] TRUE TRUE TRUE TRUE TRUE FALSE ...
> 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)
> str(a1_zlib_obs)
 logi [1:5, 1:5] TRUE TRUE TRUE TRUE TRUE FALSE ...
> stopifnot(identical(a1_zlib_obs, a1_exp))
> 
> ##
> ## logical: case-2
> ##
> a2_exp <- new("lgCMatrix",
+               i = c(6L, 0L, 4L, 5L, 0L, 2L, 4L, 1L, 4L, 6L, 7L, 3L,
+                     4L),
+               p = c(0L, 1L, 4L, 7L, 11L, 13L),
+               Dim = c(10L, 5L),
+               Dimnames = list(NULL, NULL),
+               x = c(TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE,
+                     TRUE, TRUE, TRUE, TRUE, TRUE),
+               factors = list())
> 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)
Formal class 'lgCMatrix' [package "Matrix"] with 6 slots
  ..@ i       : int [1:13] 6 0 4 5 0 2 4 1 4 6 ...
  ..@ p       : int [1:6] 0 1 4 7 11 13
  ..@ Dim     : int [1:2] 10 5
  ..@ Dimnames:List of 2
  .. ..$ : NULL
  .. ..$ : NULL
  ..@ x       : logi [1:13] TRUE TRUE TRUE TRUE TRUE TRUE ...
  ..@ factors : list()
> 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)
Formal class 'lgCMatrix' [package "Matrix"] with 6 slots
  ..@ i       : int [1:13] 6 0 4 5 0 2 4 1 4 6 ...
  ..@ p       : int [1:6] 0 1 4 7 11 13
  ..@ Dim     : int [1:2] 10 5
  ..@ Dimnames:List of 2
  .. ..$ : NULL
  .. ..$ : NULL
  ..@ x       : logi [1:13] TRUE TRUE TRUE TRUE TRUE TRUE ...
  ..@ factors : list()
> stopifnot(identical(a2_zlib_obs, a2_exp))
> 
> proc.time()
   user  system elapsed 
   0.78    0.18    0.95