src/libK/bindings/R/rlibkriging/tests/testthat/test-SaveLoad.R

In rlibkriging: Kriging Models using the 'libKriging' Library

library(testthat)

#library(rlibkriging, lib.loc="bindings/R/Rlibs")
#library(rlibkriging)

# f <- function(X) apply(X, 1, function(x) prod(sin((x-.5)^2)))
f <- function(X) apply(X, 1, function(x)
  prod(sin(2*pi*( x * (seq(0,1,l=1+length(x))[-1])^2 )))
)
n <- 20
set.seed(123)
X <- cbind(runif(n),runif(n))
y <- f(X)
d = ncol(X)

x=seq(0,1,,51)
contour(x,x,matrix(f(as.matrix(expand.grid(x,x))),nrow=length(x)),nlevels = 30)
points(X)

k <- Kriging(y, X,"gauss",parameters = list(theta=matrix(runif(40),ncol=2)))
print(k)

rlibkriging::save(k, filename="k.h5")

k2 <- rlibkriging::load(filename="k.h5")
print(k2)

test_that("Save/Load NuggetKriging", expect_true( print(k) == print(k2)))

nuk <- NuggetKriging(y, X,"gauss",parameters = list(theta=matrix(runif(40),ncol=2)))
print(nuk)

rlibkriging::save(nuk, filename="nuk.h5")

nuk2 <- rlibkriging::load(filename="nuk.h5")
print(nuk2)

test_that("Save/Load NuggetKriging", expect_true( print(nuk) == print(nuk2)))

nok <- NoiseKriging(y, rep(0.1^2,nrow(X)), X,"gauss",parameters = list(theta=matrix(runif(40),ncol=2)))
print(nok)

rlibkriging::save(nok, filename="nok.h5")

nok2 <- rlibkriging::load(filename="nok.h5")
print(nok2)

test_that("Save/Load NoiseKriging", expect_true( print(nok) == print(nok2)))