src/libK/bindings/R/rlibkriging/tests/testthat/bench-TicToc.R

In rlibkriging: Kriging Models using the 'libKriging' Library

f <- function(X) apply(X, 1, function(x) prod(sin((x-.5)^2)))
n <- 100
set.seed(123)
X <- cbind(runif(n),runif(n),runif(n))
y <- f(X)

times.n = 10
x = .5

library(rlibkriging, lib.loc="bindings/R/Rlibs")
r <- Kriging(y, X, "gauss")
for (j in 1:times.n) {
  invisible(logLikelihoodFun(r,rep(x,3),           bench=TRUE ))
}