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) > library(Matrix) > > ## 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] Matrix_1.6-4 rmatio_0.19.0 loaded via a namespace (and not attached): [1] compiler_4.4.0 grid_4.4.0 lattice_0.22-5 > > ## > ## Check write and read of cell arrays in MAT5 format: > ## 1) without compression > ## 2) with compression > ## > ## In rmatio, cell arrays are mapped to an unnamed list > > ## > ## Note: > ## If the list contains elements of differents lengths i.e. > ## a14_in <- list(c("a", "bb"), c("c", "dd")) # nolint > ## then the expected result is not identical, since each element > ## is saved in a cell and the expected result of a14_in is therefore > ## a14_exp <- list(list("a", "bb"), list("c", "dd")) # nolint > ## > > ## > ## cell: case-1 > ## > a1_exp <- list() > 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) list() > 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) list() > stopifnot(identical(a1_zlib_obs, a1_exp)) > > ## > ## cell: case-2 > ## > a2_exp <- list(complex(0), logical(0), character(0), numeric(0), integer(0)) > 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) List of 5 $ : cplx(0) $ : logi(0) $ : chr(0) $ : num(0) $ : int(0) > 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) List of 5 $ : cplx(0) $ : logi(0) $ : chr(0) $ : num(0) $ : int(0) > stopifnot(identical(a2_zlib_obs, a2_exp)) > > ## > ## cell: case-3 > ## > a3_exp <- list(list(array(c(1, 3, 2, 4), c(2, 2)), + array(c(5, 8, 6, 9, 7, 10), c(2, 3)), + array(c(11, 15, 12, 16, 13, 17, 14, 18), c(2, 4))), + list(array(c(19, 21, 20, 22), c(2, 2)), + array(c(23, 25, 27, 24, 26, 28), c(3L, 2L)), + array(c(29, 31, 33, 35, 30, 32, 34, 36), c(4, 2)))) > 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) List of 2 $ :List of 3 ..$ : num [1:2, 1:2] 1 3 2 4 ..$ : num [1:2, 1:3] 5 8 6 9 7 10 ..$ : num [1:2, 1:4] 11 15 12 16 13 17 14 18 $ :List of 3 ..$ : num [1:2, 1:2] 19 21 20 22 ..$ : num [1:3, 1:2] 23 25 27 24 26 28 ..$ : num [1:4, 1:2] 29 31 33 35 30 32 34 36 > 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) List of 2 $ :List of 3 ..$ : num [1:2, 1:2] 1 3 2 4 ..$ : num [1:2, 1:3] 5 8 6 9 7 10 ..$ : num [1:2, 1:4] 11 15 12 16 13 17 14 18 $ :List of 3 ..$ : num [1:2, 1:2] 19 21 20 22 ..$ : num [1:3, 1:2] 23 25 27 24 26 28 ..$ : num [1:4, 1:2] 29 31 33 35 30 32 34 36 > stopifnot(identical(a3_zlib_obs, a3_exp)) > > ## > ## cell: case-4 > ## > a4_exp <- list(list(array(c(1L, 3L, 2L, 4L), c(2, 2)), + array(c(5L, 8L, 6L, 9L, 7L, 10L), c(2, 3)), + array(c(11L, 15L, 12L, 16L, 13L, 17L, 14L, 18L), c(2, 4))), + list(array(c(19L, 21L, 20L, 22L), c(2, 2)), + array(c(23L, 25L, 27L, 24L, 26L, 28L), c(3L, 2L)), + array(c(29L, 31L, 33L, 35L, 30L, 32L, 34L, 36L), c(4, 2)))) > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a4_exp), + filename = filename, + compression = FALSE, + version = "MAT5") > a4_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a4_obs) List of 2 $ :List of 3 ..$ : int [1:2, 1:2] 1 3 2 4 ..$ : int [1:2, 1:3] 5 8 6 9 7 10 ..$ : int [1:2, 1:4] 11 15 12 16 13 17 14 18 $ :List of 3 ..$ : int [1:2, 1:2] 19 21 20 22 ..$ : int [1:3, 1:2] 23 25 27 24 26 28 ..$ : int [1:4, 1:2] 29 31 33 35 30 32 34 36 > stopifnot(identical(a4_obs, a4_exp)) > > ## Run the same test with compression > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a4_exp), + filename = filename, + compression = TRUE, + version = "MAT5") > a4_zlib_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a4_zlib_obs) List of 2 $ :List of 3 ..$ : int [1:2, 1:2] 1 3 2 4 ..$ : int [1:2, 1:3] 5 8 6 9 7 10 ..$ : int [1:2, 1:4] 11 15 12 16 13 17 14 18 $ :List of 3 ..$ : int [1:2, 1:2] 19 21 20 22 ..$ : int [1:3, 1:2] 23 25 27 24 26 28 ..$ : int [1:4, 1:2] 29 31 33 35 30 32 34 36 > stopifnot(identical(a4_zlib_obs, a4_exp)) > > ## > ## cell: case-5 > ## > a5_exp <- list(list(triu(Matrix(1:20, nrow = 4, ncol = 5, sparse = TRUE)), + tril(Matrix(1:20, nrow = 5, ncol = 4, sparse = TRUE, + byrow = TRUE)))) > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a5_exp), + filename = filename, + compression = FALSE, + version = "MAT5") > a5_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a5_obs) List of 1 $ :List of 2 ..$ :Formal class 'dgCMatrix' [package "Matrix"] with 6 slots .. .. ..@ i : int [1:14] 0 0 1 0 1 2 0 1 2 3 ... .. .. ..@ p : int [1:6] 0 1 3 6 10 14 .. .. ..@ Dim : int [1:2] 4 5 .. .. ..@ Dimnames:List of 2 .. .. .. ..$ : NULL .. .. .. ..$ : NULL .. .. ..@ x : num [1:14] 1 5 6 9 10 11 13 14 15 16 ... .. .. ..@ factors : list() ..$ :Formal class 'dgCMatrix' [package "Matrix"] with 6 slots .. .. ..@ i : int [1:14] 0 1 2 3 4 1 2 3 4 2 ... .. .. ..@ p : int [1:5] 0 5 9 12 14 .. .. ..@ Dim : int [1:2] 5 4 .. .. ..@ Dimnames:List of 2 .. .. .. ..$ : NULL .. .. .. ..$ : NULL .. .. ..@ x : num [1:14] 1 5 9 13 17 6 10 14 18 11 ... .. .. ..@ factors : list() > stopifnot(identical(a5_obs, a5_exp)) > > ## Run the same test with compression > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a5_exp), + filename = filename, + compression = TRUE, + version = "MAT5") > a5_zlib_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a5_zlib_obs) List of 1 $ :List of 2 ..$ :Formal class 'dgCMatrix' [package "Matrix"] with 6 slots .. .. ..@ i : int [1:14] 0 0 1 0 1 2 0 1 2 3 ... .. .. ..@ p : int [1:6] 0 1 3 6 10 14 .. .. ..@ Dim : int [1:2] 4 5 .. .. ..@ Dimnames:List of 2 .. .. .. ..$ : NULL .. .. .. ..$ : NULL .. .. ..@ x : num [1:14] 1 5 6 9 10 11 13 14 15 16 ... .. .. ..@ factors : list() ..$ :Formal class 'dgCMatrix' [package "Matrix"] with 6 slots .. .. ..@ i : int [1:14] 0 1 2 3 4 1 2 3 4 2 ... .. .. ..@ p : int [1:5] 0 5 9 12 14 .. .. ..@ Dim : int [1:2] 5 4 .. .. ..@ Dimnames:List of 2 .. .. .. ..$ : NULL .. .. .. ..$ : NULL .. .. ..@ x : num [1:14] 1 5 9 13 17 6 10 14 18 11 ... .. .. ..@ factors : list() > stopifnot(identical(a5_zlib_obs, a5_exp)) > > ## > ## cell: case-6 > ## > a6_exp <- list(array(c(1 + 21i, 0 + 0i, 0 + 0i, 0 + 0i, 5 + 25i, + 6 + 26i, 0 + 0i, 0 + 0i, 9 + 29i, 10 + 30i, + 11 + 31i, 0 + 0i, 13 + 33i, 14 + 34i, 15 + 35i, + 16 + 36i, 17 + 37i, 18 + 38i, 19 + 39i, + 20 + 40i), + c(4, 5)), + array(c(1 - 21i, 5 - 25i, 9 - 29i, 13 - 33i, 17 - 37i, + 0 + 0i, 6 - 26i, 10 - 30i, 14 - 34i, 18 - 38i, + 0 + 0i, 0 + 0i, 11 - 31i, 15 - 35i, 19 - 39i, + 0 + 0i, 0 + 0i, 0 + 0i, 16 - 36i, 20 - 40i), + c(5, 4))) > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a6_exp), + filename = filename, + compression = FALSE, + version = "MAT5") > a6_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a6_obs) List of 2 $ : cplx [1:4, 1:5] 1+21i 0+0i 0+0i ... $ : cplx [1:5, 1:4] 1-21i 5-25i 9-29i ... > stopifnot(identical(a6_obs, a6_exp)) > > ## Run the same test with compression > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a6_exp), + filename = filename, + compression = TRUE, + version = "MAT5") > a6_zlib_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a6_zlib_obs) List of 2 $ : cplx [1:4, 1:5] 1+21i 0+0i 0+0i ... $ : cplx [1:5, 1:4] 1-21i 5-25i 9-29i ... > stopifnot(identical(a6_zlib_obs, a6_exp)) > > ## > ## cell: case-7 > ## > a7_exp <- list(list("abcdefghijklmnopqrstuvwxyz", + "1234567890!@#$%^&*()-_=+`~"), # + list("ABCDEFGHIJKLMNOPQRSTUVWXYZ", + "[{]}\\|;:'\",<.>/? ")) > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a7_exp), + filename = filename, + compression = FALSE, + version = "MAT5") > a7_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a7_obs) List of 2 $ :List of 2 ..$ : chr "abcdefghijklmnopqrstuvwxyz" ..$ : chr "1234567890!@#$%^&*()-_=+`~" $ :List of 2 ..$ : chr "ABCDEFGHIJKLMNOPQRSTUVWXYZ" ..$ : chr "[{]}\\|;:'\",<.>/? " > stopifnot(identical(a7_obs, a7_exp)) > > ## Run the same test with compression > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a7_exp), + filename = filename, + compression = TRUE, + version = "MAT5") > a7_zlib_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a7_zlib_obs) List of 2 $ :List of 2 ..$ : chr "abcdefghijklmnopqrstuvwxyz" ..$ : chr "1234567890!@#$%^&*()-_=+`~" $ :List of 2 ..$ : chr "ABCDEFGHIJKLMNOPQRSTUVWXYZ" ..$ : chr "[{]}\\|;:'\",<.>/? " > stopifnot(identical(a7_zlib_obs, a7_exp)) > > ## > ## cell: case-8 > ## > a8_exp <- list(structure(list(), .Names = character(0)), + list(), + structure(list(field1 = numeric(0), + field2 = character(0)), + .Names = c("field1", "field2"))) > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a8_exp), + filename = filename, + compression = FALSE, + version = "MAT5") > a8_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a8_obs) List of 3 $ : Named list() $ : list() $ :List of 2 ..$ field1: num(0) ..$ field2: chr(0) > stopifnot(identical(a8_obs, a8_exp)) > > ## Run the same test with compression > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a8_exp), + filename = filename, + compression = TRUE, + version = "MAT5") > a8_zlib_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a8_zlib_obs) List of 3 $ : Named list() $ : list() $ :List of 2 ..$ field1: num(0) ..$ field2: chr(0) > stopifnot(identical(a8_zlib_obs, a8_exp)) > > ## > ## cell: case-9 > ## > a9_exp <- list(list(structure(list( + field1 = list(1, 14), + field2 = list( + structure(c(2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13), + .Dim = 3:4), + structure(c(15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26), + .Dim = 3:4))), + .Names = c("field1", "field2")), + structure(list( + field1 = list(1, 14), + field2 = list(structure(c(2, 3, 4, 5, 6, 7, 8, + 9, 10, 11, 12, 13), + .Dim = 3:4), + structure(c(15, 16, 17, 18, 19, + 20, 21, 22, 23, 24, + 25, 26), + .Dim = 3:4))), + .Names = c("field1", "field2")), + structure(list( + field1 = list(1, 14), + field2 = list(structure(c(2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, + 13), + .Dim = 3:4), + structure(c(15, 16, 17, 18, 19, 20, 21, 22, 23, + 24, 25, 26), + .Dim = 3:4))), + .Names = c("field1", "field2")), + structure(list( + field1 = list(1, 14), + field2 = list(structure(c(2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, + 13), + .Dim = 3:4), + structure(c(15, 16, 17, 18, 19, 20, 21, 22, 23, + 24, 25, 26), + .Dim = 3:4))), + .Names = c("field1", "field2")), + structure(list( + field1 = list(1L, 14L), + field2 = list(structure(2:13, + .Dim = 3:4), + structure(15:26, + .Dim = 3:4))), + .Names = c("field1", "field2")), + structure(list( + field1 = list(1, 14), + field2 = list(structure(c(2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, + 13), + .Dim = 3:4), + structure(c(15, 16, 17, 18, 19, 20, 21, 22, 23, + 24, 25, 26), + .Dim = 3:4))), + .Names = c("field1", "field2")), + structure(list( + field1 = list(1L, 14L), + field2 = list(structure(2:13, + .Dim = 3:4), + structure(15:26, + .Dim = 3:4))), + .Names = c("field1", "field2")), + structure(list( + field1 = list(1L, 14L), + field2 = list(structure(2:13, + .Dim = 3:4), + structure(15:26, + .Dim = 3:4))), + .Names = c("field1", "field2")), + structure(list( + field1 = list(1L, 14L), + field2 = list(structure(2:13, + .Dim = 3:4), + structure(15:26, + .Dim = 3:4))), + .Names = c("field1", "field2")), + structure(list( + field1 = list(1L, 14L), + field2 = list(structure(2:13, + .Dim = 3:4), + structure(15:26, + .Dim = 3:4))), + .Names = c("field1", "field2"))), + list(structure(list( + field1 = list(1 + 51i, 14 + 64i), + field2 = list(structure(c(2 + 52i, 3 + 53i, 4 + 54i, 5 + 55i, + 6 + 56i, 7 + 57i, 8 + 58i, 9 + 59i, + 10 + 60i, 11 + 61i, 12 + 62i, + 13 + 63i), + .Dim = 3:4), + structure(c(15 + 65i, 16 + 66i, 17 + 67i, + 18 + 68i, 19 + 69i, 20 + 70i, + 21 + 71i, 22 + 72i, 23 + 73i, + 24 + 74i, 25 + 75i, 26 + 76i), + .Dim = 3:4))), + .Names = c("field1", "field2")), + structure(list( + field1 = list(1 + 51i, 14 + 64i), + field2 = list(structure(c(2 + 52i, 3 + 53i, 4 + 54i, + 5 + 55i, 6 + 56i, 7 + 57i, + 8 + 58i, 9 + 59i, 10 + 60i, + 11 + 61i, 12 + 62i, 13 + 63i), + .Dim = 3:4), + structure(c(15 + 65i, 16 + 66i, 17 + 67i, + 18 + 68i, 19 + 69i, 20 + 70i, + 21 + 71i, 22 + 72i, 23 + 73i, + 24 + 74i, 25 + 75i, 26 + 76i), + .Dim = 3:4))), + .Names = c("field1", "field2")), + structure(list( + field1 = list(1 + 51i, 14 + 64i), + field2 = list(structure(c(2 + 52i, 3 + 53i, 4 + 54i, + 5 + 55i, 6 + 56i, 7 + 57i, + 8 + 58i, 9 + 59i, 10 + 60i, + 11 + 61i, 12 + 62i, 13 + 63i), + .Dim = 3:4), + structure(c(15 + 65i, 16 + 66i, 17 + 67i, + 18 + 68i, 19 + 69i, 20 + 70i, + 21 + 71i, 22 + 72i, 23 + 73i, + 24 + 74i, 25 + 75i, 26 + 76i), + .Dim = 3:4))), + .Names = c("field1", "field2")), + structure(list( + field1 = list(1 + 51i, 14 + 64i), + field2 = list(structure(c(2 + 52i, 3 + 53i, 4 + 54i, + 5 + 55i, 6 + 56i, 7 + 57i, + 8 + 58i, 9 + 59i, 10 + 60i, + 11 + 61i, 12 + 62i, 13 + 63i), + .Dim = 3:4), + structure(c(15 + 65i, 16 + 66i, 17 + 67i, + 18 + 68i, 19 + 69i, 20 + 70i, + 21 + 71i, 22 + 72i, 23 + 73i, + 24 + 74i, 25 + 75i, 26 + 76i), + .Dim = 3:4))), + .Names = c("field1", "field2")), + structure(list( + field1 = list(1 + 51i, 14 + 64i), + field2 = list(structure(c(2 + 52i, 3 + 53i, 4 + 54i, + 5 + 55i, 6 + 56i, 7 + 57i, + 8 + 58i, 9 + 59i, 10 + 60i, + 11 + 61i, 12 + 62i, 13 + 63i), + .Dim = 3:4), + structure(c(15 + 65i, 16 + 66i, 17 + 67i, + 18 + 68i, 19 + 69i, 20 + 70i, + 21 + 71i, 22 + 72i, 23 + 73i, + 24 + 74i, 25 + 75i, 26 + 76i), + .Dim = 3:4))), + .Names = c("field1", "field2")), + structure(list( + field1 = list(1 + 51i, 14 + 64i), + field2 = list(structure(c(2 + 52i, 3 + 53i, 4 + 54i, + 5 + 55i, 6 + 56i, 7 + 57i, + 8 + 58i, 9 + 59i, 10 + 60i, + 11 + 61i, 12 + 62i, 13 + 63i), + .Dim = 3:4), + structure(c(15 + 65i, 16 + 66i, 17 + 67i, + 18 + 68i, 19 + 69i, 20 + 70i, + 21 + 71i, 22 + 72i, 23 + 73i, + 24 + 74i, 25 + 75i, 26 + 76i), + .Dim = 3:4))), + .Names = c("field1", "field2")), + structure(list( + field1 = list(1 + 51i, 14 + 64i), + field2 = list(structure(c(2 + 52i, 3 + 53i, 4 + 54i, + 5 + 55i, 6 + 56i, 7 + 57i, + 8 + 58i, 9 + 59i, 10 + 60i, + 11 + 61i, 12 + 62i, 13 + 63i), + .Dim = 3:4), + structure(c(15 + 65i, 16 + 66i, 17 + 67i, + 18 + 68i, 19 + 69i, 20 + 70i, + 21 + 71i, 22 + 72i, 23 + 73i, + 24 + 74i, 25 + 75i, 26 + 76i), + .Dim = 3:4))), + .Names = c("field1", "field2")), + structure(list( + field1 = list(1 + 51i, 14 + 64i), + field2 = list(structure(c(2 + 52i, 3 + 53i, 4 + 54i, + 5 + 55i, 6 + 56i, 7 + 57i, + 8 + 58i, 9 + 59i, 10 + 60i, + 11 + 61i, 12 + 62i, 13 + 63i), + .Dim = 3:4), + structure(c(15 + 65i, 16 + 66i, 17 + 67i, + 18 + 68i, 19 + 69i, 20 + 70i, + 21 + 71i, 22 + 72i, 23 + 73i, + 24 + 74i, 25 + 75i, 26 + 76i), + .Dim = 3:4))), + .Names = c("field1", "field2")), + structure(list( + field1 = list(1 + 51i, 14 + 64i), + field2 = list(structure(c(2 + 52i, 3 + 53i, 4 + 54i, + 5 + 55i, 6 + 56i, 7 + 57i, + 8 + 58i, 9 + 59i, 10 + 60i, + 11 + 61i, 12 + 62i, 13 + 63i), + .Dim = 3:4), + structure(c(15 + 65i, 16 + 66i, 17 + 67i, + 18 + 68i, 19 + 69i, 20 + 70i, + 21 + 71i, 22 + 72i, 23 + 73i, + 24 + 74i, 25 + 75i, 26 + 76i), + .Dim = 3:4))), + .Names = c("field1", "field2")), + structure(list( + field1 = list(1 + 51i, 14 + 64i), + field2 = list(structure(c(2 + 52i, 3 + 53i, 4 + 54i, + 5 + 55i, 6 + 56i, 7 + 57i, + 8 + 58i, 9 + 59i, 10 + 60i, + 11 + 61i, 12 + 62i, 13 + 63i), + .Dim = 3:4), + structure(c(15 + 65i, 16 + 66i, 17 + 67i, + 18 + 68i, 19 + 69i, 20 + 70i, + 21 + 71i, 22 + 72i, 23 + 73i, + 24 + 74i, 25 + 75i, 26 + 76i), + .Dim = 3:4))), + .Names = c("field1", "field2")))) > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a9_exp), + filename = filename, + compression = FALSE, + version = "MAT5") > a9_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a9_obs) List of 2 $ :List of 10 ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : num 1 .. .. ..$ : num 14 .. ..$ field2:List of 2 .. .. ..$ : num [1:3, 1:4] 2 3 4 5 6 7 8 9 10 11 ... .. .. ..$ : num [1:3, 1:4] 15 16 17 18 19 20 21 22 23 24 ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : num 1 .. .. ..$ : num 14 .. ..$ field2:List of 2 .. .. ..$ : num [1:3, 1:4] 2 3 4 5 6 7 8 9 10 11 ... .. .. ..$ : num [1:3, 1:4] 15 16 17 18 19 20 21 22 23 24 ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : num 1 .. .. ..$ : num 14 .. ..$ field2:List of 2 .. .. ..$ : num [1:3, 1:4] 2 3 4 5 6 7 8 9 10 11 ... .. .. ..$ : num [1:3, 1:4] 15 16 17 18 19 20 21 22 23 24 ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : num 1 .. .. ..$ : num 14 .. ..$ field2:List of 2 .. .. ..$ : num [1:3, 1:4] 2 3 4 5 6 7 8 9 10 11 ... .. .. ..$ : num [1:3, 1:4] 15 16 17 18 19 20 21 22 23 24 ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : int 1 .. .. ..$ : int 14 .. ..$ field2:List of 2 .. .. ..$ : int [1:3, 1:4] 2 3 4 5 6 7 8 9 10 11 ... .. .. ..$ : int [1:3, 1:4] 15 16 17 18 19 20 21 22 23 24 ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : num 1 .. .. ..$ : num 14 .. ..$ field2:List of 2 .. .. ..$ : num [1:3, 1:4] 2 3 4 5 6 7 8 9 10 11 ... .. .. ..$ : num [1:3, 1:4] 15 16 17 18 19 20 21 22 23 24 ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : int 1 .. .. ..$ : int 14 .. ..$ field2:List of 2 .. .. ..$ : int [1:3, 1:4] 2 3 4 5 6 7 8 9 10 11 ... .. .. ..$ : int [1:3, 1:4] 15 16 17 18 19 20 21 22 23 24 ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : int 1 .. .. ..$ : int 14 .. ..$ field2:List of 2 .. .. ..$ : int [1:3, 1:4] 2 3 4 5 6 7 8 9 10 11 ... .. .. ..$ : int [1:3, 1:4] 15 16 17 18 19 20 21 22 23 24 ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : int 1 .. .. ..$ : int 14 .. ..$ field2:List of 2 .. .. ..$ : int [1:3, 1:4] 2 3 4 5 6 7 8 9 10 11 ... .. .. ..$ : int [1:3, 1:4] 15 16 17 18 19 20 21 22 23 24 ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : int 1 .. .. ..$ : int 14 .. ..$ field2:List of 2 .. .. ..$ : int [1:3, 1:4] 2 3 4 5 6 7 8 9 10 11 ... .. .. ..$ : int [1:3, 1:4] 15 16 17 18 19 20 21 22 23 24 ... $ :List of 10 ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : cplx 1+51i .. .. ..$ : cplx 14+64i .. ..$ field2:List of 2 .. .. ..$ : cplx [1:3, 1:4] 2+52i 3+53i 4+54i ... .. .. ..$ : cplx [1:3, 1:4] 15+65i 16+66i 17+67i ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : cplx 1+51i .. .. ..$ : cplx 14+64i .. ..$ field2:List of 2 .. .. ..$ : cplx [1:3, 1:4] 2+52i 3+53i 4+54i ... .. .. ..$ : cplx [1:3, 1:4] 15+65i 16+66i 17+67i ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : cplx 1+51i .. .. ..$ : cplx 14+64i .. ..$ field2:List of 2 .. .. ..$ : cplx [1:3, 1:4] 2+52i 3+53i 4+54i ... .. .. ..$ : cplx [1:3, 1:4] 15+65i 16+66i 17+67i ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : cplx 1+51i .. .. ..$ : cplx 14+64i .. ..$ field2:List of 2 .. .. ..$ : cplx [1:3, 1:4] 2+52i 3+53i 4+54i ... .. .. ..$ : cplx [1:3, 1:4] 15+65i 16+66i 17+67i ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : cplx 1+51i .. .. ..$ : cplx 14+64i .. ..$ field2:List of 2 .. .. ..$ : cplx [1:3, 1:4] 2+52i 3+53i 4+54i ... .. .. ..$ : cplx [1:3, 1:4] 15+65i 16+66i 17+67i ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : cplx 1+51i .. .. ..$ : cplx 14+64i .. ..$ field2:List of 2 .. .. ..$ : cplx [1:3, 1:4] 2+52i 3+53i 4+54i ... .. .. ..$ : cplx [1:3, 1:4] 15+65i 16+66i 17+67i ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : cplx 1+51i .. .. ..$ : cplx 14+64i .. ..$ field2:List of 2 .. .. ..$ : cplx [1:3, 1:4] 2+52i 3+53i 4+54i ... .. .. ..$ : cplx [1:3, 1:4] 15+65i 16+66i 17+67i ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : cplx 1+51i .. .. ..$ : cplx 14+64i .. ..$ field2:List of 2 .. .. ..$ : cplx [1:3, 1:4] 2+52i 3+53i 4+54i ... .. .. ..$ : cplx [1:3, 1:4] 15+65i 16+66i 17+67i ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : cplx 1+51i .. .. ..$ : cplx 14+64i .. ..$ field2:List of 2 .. .. ..$ : cplx [1:3, 1:4] 2+52i 3+53i 4+54i ... .. .. ..$ : cplx [1:3, 1:4] 15+65i 16+66i 17+67i ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : cplx 1+51i .. .. ..$ : cplx 14+64i .. ..$ field2:List of 2 .. .. ..$ : cplx [1:3, 1:4] 2+52i 3+53i 4+54i ... .. .. ..$ : cplx [1:3, 1:4] 15+65i 16+66i 17+67i ... > stopifnot(identical(a9_obs, a9_exp)) > > ## Run the same test with compression > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a9_exp), + filename = filename, + compression = TRUE, + version = "MAT5") > a9_zlib_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a9_zlib_obs) List of 2 $ :List of 10 ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : num 1 .. .. ..$ : num 14 .. ..$ field2:List of 2 .. .. ..$ : num [1:3, 1:4] 2 3 4 5 6 7 8 9 10 11 ... .. .. ..$ : num [1:3, 1:4] 15 16 17 18 19 20 21 22 23 24 ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : num 1 .. .. ..$ : num 14 .. ..$ field2:List of 2 .. .. ..$ : num [1:3, 1:4] 2 3 4 5 6 7 8 9 10 11 ... .. .. ..$ : num [1:3, 1:4] 15 16 17 18 19 20 21 22 23 24 ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : num 1 .. .. ..$ : num 14 .. ..$ field2:List of 2 .. .. ..$ : num [1:3, 1:4] 2 3 4 5 6 7 8 9 10 11 ... .. .. ..$ : num [1:3, 1:4] 15 16 17 18 19 20 21 22 23 24 ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : num 1 .. .. ..$ : num 14 .. ..$ field2:List of 2 .. .. ..$ : num [1:3, 1:4] 2 3 4 5 6 7 8 9 10 11 ... .. .. ..$ : num [1:3, 1:4] 15 16 17 18 19 20 21 22 23 24 ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : int 1 .. .. ..$ : int 14 .. ..$ field2:List of 2 .. .. ..$ : int [1:3, 1:4] 2 3 4 5 6 7 8 9 10 11 ... .. .. ..$ : int [1:3, 1:4] 15 16 17 18 19 20 21 22 23 24 ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : num 1 .. .. ..$ : num 14 .. ..$ field2:List of 2 .. .. ..$ : num [1:3, 1:4] 2 3 4 5 6 7 8 9 10 11 ... .. .. ..$ : num [1:3, 1:4] 15 16 17 18 19 20 21 22 23 24 ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : int 1 .. .. ..$ : int 14 .. ..$ field2:List of 2 .. .. ..$ : int [1:3, 1:4] 2 3 4 5 6 7 8 9 10 11 ... .. .. ..$ : int [1:3, 1:4] 15 16 17 18 19 20 21 22 23 24 ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : int 1 .. .. ..$ : int 14 .. ..$ field2:List of 2 .. .. ..$ : int [1:3, 1:4] 2 3 4 5 6 7 8 9 10 11 ... .. .. ..$ : int [1:3, 1:4] 15 16 17 18 19 20 21 22 23 24 ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : int 1 .. .. ..$ : int 14 .. ..$ field2:List of 2 .. .. ..$ : int [1:3, 1:4] 2 3 4 5 6 7 8 9 10 11 ... .. .. ..$ : int [1:3, 1:4] 15 16 17 18 19 20 21 22 23 24 ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : int 1 .. .. ..$ : int 14 .. ..$ field2:List of 2 .. .. ..$ : int [1:3, 1:4] 2 3 4 5 6 7 8 9 10 11 ... .. .. ..$ : int [1:3, 1:4] 15 16 17 18 19 20 21 22 23 24 ... $ :List of 10 ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : cplx 1+51i .. .. ..$ : cplx 14+64i .. ..$ field2:List of 2 .. .. ..$ : cplx [1:3, 1:4] 2+52i 3+53i 4+54i ... .. .. ..$ : cplx [1:3, 1:4] 15+65i 16+66i 17+67i ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : cplx 1+51i .. .. ..$ : cplx 14+64i .. ..$ field2:List of 2 .. .. ..$ : cplx [1:3, 1:4] 2+52i 3+53i 4+54i ... .. .. ..$ : cplx [1:3, 1:4] 15+65i 16+66i 17+67i ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : cplx 1+51i .. .. ..$ : cplx 14+64i .. ..$ field2:List of 2 .. .. ..$ : cplx [1:3, 1:4] 2+52i 3+53i 4+54i ... .. .. ..$ : cplx [1:3, 1:4] 15+65i 16+66i 17+67i ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : cplx 1+51i .. .. ..$ : cplx 14+64i .. ..$ field2:List of 2 .. .. ..$ : cplx [1:3, 1:4] 2+52i 3+53i 4+54i ... .. .. ..$ : cplx [1:3, 1:4] 15+65i 16+66i 17+67i ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : cplx 1+51i .. .. ..$ : cplx 14+64i .. ..$ field2:List of 2 .. .. ..$ : cplx [1:3, 1:4] 2+52i 3+53i 4+54i ... .. .. ..$ : cplx [1:3, 1:4] 15+65i 16+66i 17+67i ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : cplx 1+51i .. .. ..$ : cplx 14+64i .. ..$ field2:List of 2 .. .. ..$ : cplx [1:3, 1:4] 2+52i 3+53i 4+54i ... .. .. ..$ : cplx [1:3, 1:4] 15+65i 16+66i 17+67i ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : cplx 1+51i .. .. ..$ : cplx 14+64i .. ..$ field2:List of 2 .. .. ..$ : cplx [1:3, 1:4] 2+52i 3+53i 4+54i ... .. .. ..$ : cplx [1:3, 1:4] 15+65i 16+66i 17+67i ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : cplx 1+51i .. .. ..$ : cplx 14+64i .. ..$ field2:List of 2 .. .. ..$ : cplx [1:3, 1:4] 2+52i 3+53i 4+54i ... .. .. ..$ : cplx [1:3, 1:4] 15+65i 16+66i 17+67i ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : cplx 1+51i .. .. ..$ : cplx 14+64i .. ..$ field2:List of 2 .. .. ..$ : cplx [1:3, 1:4] 2+52i 3+53i 4+54i ... .. .. ..$ : cplx [1:3, 1:4] 15+65i 16+66i 17+67i ... ..$ :List of 2 .. ..$ field1:List of 2 .. .. ..$ : cplx 1+51i .. .. ..$ : cplx 14+64i .. ..$ field2:List of 2 .. .. ..$ : cplx [1:3, 1:4] 2+52i 3+53i 4+54i ... .. .. ..$ : cplx [1:3, 1:4] 15+65i 16+66i 17+67i ... > stopifnot(identical(a9_zlib_obs, a9_exp)) > > ## > ## cell: case-10 > ## > a10_exp <- list(list(field1 = list(triu(Matrix(1:20, nrow = 4, ncol = 5, + sparse = TRUE))), + field2 = list(tril(Matrix(1:20, nrow = 5, ncol = 4, + sparse = TRUE, + byrow = TRUE)))), + list(field1 = list(array(c(1 + 21i, 0 + 0i, 0 + 0i, + 0 + 0i, 5 + 25i, 6 + 26i, + 0 + 0i, 0 + 0i, 9 + 29i, + 10 + 30i, 11 + 31i, 0 + 0i, + 13 + 33i, 14 + 34i, + 15 + 35i, 16 + 36i, + 17 + 37i, 18 + 38i, + 19 + 39i, 20 + 40i), + c(4, 5))), + field2 = list(array(c(1 - 21i, 5 - 25i, 9 - 29i, + 13 - 33i, 17 - 37i, 0 + 0i, + 6 - 26i, 10 - 30i, + 14 - 34i, 18 - 38i, 0 + 0i, + 0 + 0i, 11 - 31i, 15 - 35i, + 19 - 39i, 0 + 0i, 0 + 0i, + 0 + 0i, 16 - 36i, 20 - 40i), + c(5, 4))))) > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a10_exp), + filename = filename, + compression = FALSE, + version = "MAT5") > a10_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a10_obs) List of 2 $ :List of 2 ..$ field1:List of 1 .. ..$ :Formal class 'dgCMatrix' [package "Matrix"] with 6 slots .. .. .. ..@ i : int [1:14] 0 0 1 0 1 2 0 1 2 3 ... .. .. .. ..@ p : int [1:6] 0 1 3 6 10 14 .. .. .. ..@ Dim : int [1:2] 4 5 .. .. .. ..@ Dimnames:List of 2 .. .. .. .. ..$ : NULL .. .. .. .. ..$ : NULL .. .. .. ..@ x : num [1:14] 1 5 6 9 10 11 13 14 15 16 ... .. .. .. ..@ factors : list() ..$ field2:List of 1 .. ..$ :Formal class 'dgCMatrix' [package "Matrix"] with 6 slots .. .. .. ..@ i : int [1:14] 0 1 2 3 4 1 2 3 4 2 ... .. .. .. ..@ p : int [1:5] 0 5 9 12 14 .. .. .. ..@ Dim : int [1:2] 5 4 .. .. .. ..@ Dimnames:List of 2 .. .. .. .. ..$ : NULL .. .. .. .. ..$ : NULL .. .. .. ..@ x : num [1:14] 1 5 9 13 17 6 10 14 18 11 ... .. .. .. ..@ factors : list() $ :List of 2 ..$ field1:List of 1 .. ..$ : cplx [1:4, 1:5] 1+21i 0+0i 0+0i ... ..$ field2:List of 1 .. ..$ : cplx [1:5, 1:4] 1-21i 5-25i 9-29i ... > stopifnot(identical(a10_obs, a10_exp)) > > ## Run the same test with compression > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a10_exp), + filename = filename, + compression = TRUE, + version = "MAT5") > a10_zlib_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a10_zlib_obs) List of 2 $ :List of 2 ..$ field1:List of 1 .. ..$ :Formal class 'dgCMatrix' [package "Matrix"] with 6 slots .. .. .. ..@ i : int [1:14] 0 0 1 0 1 2 0 1 2 3 ... .. .. .. ..@ p : int [1:6] 0 1 3 6 10 14 .. .. .. ..@ Dim : int [1:2] 4 5 .. .. .. ..@ Dimnames:List of 2 .. .. .. .. ..$ : NULL .. .. .. .. ..$ : NULL .. .. .. ..@ x : num [1:14] 1 5 6 9 10 11 13 14 15 16 ... .. .. .. ..@ factors : list() ..$ field2:List of 1 .. ..$ :Formal class 'dgCMatrix' [package "Matrix"] with 6 slots .. .. .. ..@ i : int [1:14] 0 1 2 3 4 1 2 3 4 2 ... .. .. .. ..@ p : int [1:5] 0 5 9 12 14 .. .. .. ..@ Dim : int [1:2] 5 4 .. .. .. ..@ Dimnames:List of 2 .. .. .. .. ..$ : NULL .. .. .. .. ..$ : NULL .. .. .. ..@ x : num [1:14] 1 5 9 13 17 6 10 14 18 11 ... .. .. .. ..@ factors : list() $ :List of 2 ..$ field1:List of 1 .. ..$ : cplx [1:4, 1:5] 1+21i 0+0i 0+0i ... ..$ field2:List of 1 .. ..$ : cplx [1:5, 1:4] 1-21i 5-25i 9-29i ... > stopifnot(identical(a10_zlib_obs, a10_exp)) > > ## > ## cell: case-11 > ## > a11_exp <- list(list(field1 = "abcdefghijklmnopqrstuvwxyz", + field2 = "ABCDEFGHIJKLMNOPQRSTUVWXYZ"), + list(field1 = "1234567890!@#$%^&*()-_=+`~", # + field2 = "[{]}\\|;:'\",<.>/? ")) > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a11_exp), + filename = filename, + compression = FALSE, + version = "MAT5") > a11_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a11_obs) List of 2 $ :List of 2 ..$ field1: chr "abcdefghijklmnopqrstuvwxyz" ..$ field2: chr "ABCDEFGHIJKLMNOPQRSTUVWXYZ" $ :List of 2 ..$ field1: chr "1234567890!@#$%^&*()-_=+`~" ..$ field2: chr "[{]}\\|;:'\",<.>/? " > stopifnot(identical(a11_obs, a11_exp)) > > ## Run the same test with compression > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a11_exp), + filename = filename, + compression = TRUE, + version = "MAT5") > a11_zlib_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a11_zlib_obs) List of 2 $ :List of 2 ..$ field1: chr "abcdefghijklmnopqrstuvwxyz" ..$ field2: chr "ABCDEFGHIJKLMNOPQRSTUVWXYZ" $ :List of 2 ..$ field1: chr "1234567890!@#$%^&*()-_=+`~" ..$ field2: chr "[{]}\\|;:'\",<.>/? " > stopifnot(identical(a11_zlib_obs, a11_exp)) > > ## > ## cell: case-12 > ## > a12_exp <- list(structure(c(FALSE, TRUE, FALSE, TRUE, FALSE, + TRUE, FALSE, TRUE, FALSE, TRUE, + FALSE, TRUE, FALSE, TRUE, FALSE, + TRUE, FALSE, TRUE, FALSE, TRUE), + .Dim = 4:5), + structure(c(TRUE, FALSE, TRUE, + FALSE, TRUE, FALSE, TRUE, FALSE, + TRUE, FALSE, TRUE, FALSE, TRUE, + FALSE, TRUE, FALSE, TRUE, FALSE, + TRUE, FALSE), + .Dim = 4:5), + structure(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), + .Dim = c(5L, 5L)), + structure(c(TRUE, FALSE, + FALSE, FALSE, FALSE, TRUE, TRUE, + FALSE, FALSE, FALSE, TRUE, TRUE, + TRUE, FALSE, FALSE, TRUE, TRUE, + TRUE, TRUE, FALSE, TRUE, TRUE, + TRUE, TRUE, TRUE), + .Dim = c(5L, 5L))) > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a12_exp), + filename = filename, + compression = FALSE, + version = "MAT5") > a12_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a12_obs) List of 4 $ : logi [1:4, 1:5] FALSE TRUE FALSE TRUE FALSE TRUE ... $ : logi [1:4, 1:5] TRUE FALSE TRUE FALSE TRUE FALSE ... $ : logi [1:5, 1:5] TRUE TRUE TRUE TRUE TRUE FALSE ... $ : logi [1:5, 1:5] TRUE FALSE FALSE FALSE FALSE TRUE ... > stopifnot(identical(a12_obs, a12_exp)) > > ## Run the same test with compression > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a12_exp), + filename = filename, + compression = TRUE, + version = "MAT5") > a12_zlib_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a12_zlib_obs) List of 4 $ : logi [1:4, 1:5] FALSE TRUE FALSE TRUE FALSE TRUE ... $ : logi [1:4, 1:5] TRUE FALSE TRUE FALSE TRUE FALSE ... $ : logi [1:5, 1:5] TRUE TRUE TRUE TRUE TRUE FALSE ... $ : logi [1:5, 1:5] TRUE FALSE FALSE FALSE FALSE TRUE ... > stopifnot(identical(a12_zlib_obs, a12_exp)) > > ## > ## cell: case-13 > ## > a13_exp <- list(structure(list(), + .Names = character(0)), + list(field1 = list(), field2 = list()), + structure(list(field1 = numeric(0), + field2 = character(0)), + .Names = c("field1", "field2"))) > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a13_exp), + filename = filename, + compression = FALSE, + version = "MAT5") > a13_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a13_obs) List of 3 $ : Named list() $ :List of 2 ..$ field1: list() ..$ field2: list() $ :List of 2 ..$ field1: num(0) ..$ field2: chr(0) > stopifnot(identical(a13_obs, a13_exp)) > > ## Run the same test with compression > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a13_exp), + filename = filename, + compression = TRUE, + version = "MAT5") > a13_zlib_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a13_zlib_obs) List of 3 $ : Named list() $ :List of 2 ..$ field1: list() ..$ field2: list() $ :List of 2 ..$ field1: num(0) ..$ field2: chr(0) > stopifnot(identical(a13_zlib_obs, a13_exp)) > > ## > ## cell: case-14 > ## > a14_in <- list(c("a", "bb"), c("c", "dd")) > a14_exp <- list(list("a", "bb"), list("c", "dd")) > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a14_in), + filename = filename, + compression = FALSE, + version = "MAT5") > a14_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a14_obs) List of 2 $ :List of 2 ..$ : chr "a" ..$ : chr "bb" $ :List of 2 ..$ : chr "c" ..$ : chr "dd" > stopifnot(identical(a14_obs, a14_exp)) > > ## Run the same test with compression > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a14_in), + filename = filename, + compression = TRUE, + version = "MAT5") > a14_zlib_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a14_zlib_obs) List of 2 $ :List of 2 ..$ : chr "a" ..$ : chr "bb" $ :List of 2 ..$ : chr "c" ..$ : chr "dd" > stopifnot(identical(a14_zlib_obs, a14_exp)) > > ## > ## cell: case-15 > ## > a15_in <- list(c("a", "bb"), list(c("d", "eee"))) > a15_exp <- list(list("a", "bb"), list(list(list("d", "eee")))) > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a15_in), + filename = filename, + compression = FALSE, + version = "MAT5") > a15_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a15_obs) List of 2 $ :List of 2 ..$ : chr "a" ..$ : chr "bb" $ :List of 1 ..$ :List of 1 .. ..$ :List of 2 .. .. ..$ : chr "d" .. .. ..$ : chr "eee" > stopifnot(identical(a15_obs, a15_exp)) > > ## Run the same test with compression > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a15_in), + filename = filename, + compression = TRUE, + version = "MAT5") > a15_zlib_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a15_zlib_obs) List of 2 $ :List of 2 ..$ : chr "a" ..$ : chr "bb" $ :List of 1 ..$ :List of 1 .. ..$ :List of 2 .. .. ..$ : chr "d" .. .. ..$ : chr "eee" > stopifnot(identical(a15_zlib_obs, a15_exp)) > > ## > ## cell: case-16 > ## > a16_in <- list(c("a", "bb"), + Matrix(c(0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, + 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1), + nrow = 3, + ncol = 9, + byrow = TRUE, + sparse = TRUE)) > a16_exp <- list(list("a", "bb"), list(a16_in[[2]])) > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a16_in), + filename = filename, + compression = FALSE, + version = "MAT5") > a16_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a16_obs) List of 2 $ :List of 2 ..$ : chr "a" ..$ : chr "bb" $ :List of 1 ..$ :Formal class 'dgCMatrix' [package "Matrix"] with 6 slots .. .. ..@ i : int [1:3] 0 1 2 .. .. ..@ p : int [1:10] 0 0 0 0 0 0 0 1 2 3 .. .. ..@ Dim : int [1:2] 3 9 .. .. ..@ Dimnames:List of 2 .. .. .. ..$ : NULL .. .. .. ..$ : NULL .. .. ..@ x : num [1:3] 1 1 1 .. .. ..@ factors : list() > stopifnot(identical(a16_obs, a16_exp)) > > ## Run the same test with compression > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a16_in), + filename = filename, + compression = TRUE, + version = "MAT5") > a16_zlib_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a16_zlib_obs) List of 2 $ :List of 2 ..$ : chr "a" ..$ : chr "bb" $ :List of 1 ..$ :Formal class 'dgCMatrix' [package "Matrix"] with 6 slots .. .. ..@ i : int [1:3] 0 1 2 .. .. ..@ p : int [1:10] 0 0 0 0 0 0 0 1 2 3 .. .. ..@ Dim : int [1:2] 3 9 .. .. ..@ Dimnames:List of 2 .. .. .. ..$ : NULL .. .. .. ..$ : NULL .. .. ..@ x : num [1:3] 1 1 1 .. .. ..@ factors : list() > stopifnot(identical(a16_zlib_obs, a16_exp)) > > ## > ## cell: case-17 > ## > a17_in <- list(list(c("a", "bb")), list()) > a17_exp <- list(list(list(list("a", "bb"))), list()) > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a17_in), + filename = filename, + compression = FALSE, + version = "MAT5") > a17_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a17_obs) List of 2 $ :List of 1 ..$ :List of 1 .. ..$ :List of 2 .. .. ..$ : chr "a" .. .. ..$ : chr "bb" $ : list() > stopifnot(identical(a17_obs, a17_exp)) > > ## Run the same test with compression > filename <- tempfile(fileext = ".mat") > write.mat(list(a = a17_in), + filename = filename, + compression = TRUE, + version = "MAT5") > a17_zlib_obs <- read.mat(filename)[["a"]] > unlink(filename) > str(a17_zlib_obs) List of 2 $ :List of 1 ..$ :List of 1 .. ..$ :List of 2 .. .. ..$ : chr "a" .. .. ..$ : chr "bb" $ : list() > stopifnot(identical(a17_zlib_obs, a17_exp)) > > proc.time() user system elapsed 1.12 0.17 1.29