R/combinepmp.R

In mcrucifix/gp: Gaussian process calibration for emulation

combinepmp <-
function(Xp, xmp)
 {
  # combine a subspace with one element belonging to the remaining space 
  if (length(Xp) ==0) return (xmp) else
  {
  xmp <- as.matrix(xmp)
  Xmp <- matrix( xmp, nrow(Xp), length(xmp), byrow=TRUE, dimnames=list(NULL, colnames(xmp)))
  X <- cbind(Xp, Xmp)
  X <- X[, order(colnames(X))]
  X
  }
}