R Under development (unstable) (2025-09-17 r88852 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2025 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.

> # Various tests about object creation and manipulation
> require(houba)
Loading required package: houba
> 
> # float mmatrix ------------------------
> A <- mmatrix("float", 10, 20)
> a <- as.matrix(A)
> stopifnot(typeof(a) == "double")
> 
> # set max size
> houba(max.size = 10)
$max.size
[1] 10

> stopifnot( typeof(A[1:10]) == "double" )
> stopifnot( typeof(A[1:11]) == "S4" )
> 
> # assignement and subsetting
> A[1,1] <- 12
> A[3,] <- 1.34
> stopifnot( all(abs(as.matrix(A[1:3, 1:2] - c(12, 0, 1.34, 0, 0, 1.34))) < 1e-6)  )
> 
> # int mmatrix ---------------------------
> B <- mmatrix("int", 10, 20)
> b <- as.matrix(B)
> stopifnot(typeof(b) == "integer")
> 
> # assignement and subsetting
> B[1,1] <- 12
> B[3,] <- 1.34
> stopifnot( all(as.matrix(B[1:3, 1:2] == c(12L, 0L, 1L, 0L, 0L, 1L))) )
> 
> # assignement with other mmatrix values
> houba(max.size = 0) # force non conversion to R 
$max.size
[1] 0

> B[2,] <- A[3,]
> stopifnot( all(as.vector(B[2, 1:4]) == c(1L, 1L, 1L, 1L)) )
> 
> # double mvector -------------------------
> V <- mvector("double", 10)
> v <- as.vector(V)
> V[1:4] <- pi
> stopifnot( all(as.vector(V[3:6]) == c(pi, pi, 0, 0)) ) 
> 
> 
> # int16 mmatrix ----------------------
> C <- mmatrix("short", 10, 20)
> C[] <- sample.int(200)
> 
> # create descriptor file 
> dsc <- descriptor.file(C)
Created descriptor file D:\temp\2025_09_18_11_20_17_21797\RtmpyAyUA3\mmatrix29e1440875b6d.desc 
> 
> # linking it to other object
> D <- read.descriptor(dsc)
> stopifnot( all(as.matrix(C) == as.matrix(D)))
> 
> # descriptor for mvector ----------------
> dsc <- descriptor.file(V)
Created descriptor file D:\temp\2025_09_18_11_20_17_21797\RtmpyAyUA3\mmatrix29e14f3520e4.desc 
> 
> # reading it
> Vbis <- read.descriptor(dsc, FALSE) #so NOT read-only
> 
> # modified V through Vbis
> Vbis[,] <- pi
> flush(Vbis)
> stopifnot( all(as.vector(V) == pi) )
> 
> proc.time()
   user  system elapsed 
   0.34    0.10    0.43