tests/testthat/optimCheck-testfunctions.R

In mlysy/optimCheck: Graphical and Numerical Checks for Mode-Finding Routines

#--- package test functions ----------------------------------------------------

# generate an n x p matrix of iid N(0,1)
rMnorm <- function(n, p) {
  if(missing(p)) p <- n
  matrix(rnorm(n*p), n, p)
}

# max of min of abs and rel error
max.xdiff <- function(x) {
  xdiff <- abs(diff(x))
  max(pmin(xdiff[,1], xdiff[,2]))
}