R/correlatedRNG.R

Defines functions corrnorm2 corrnorm

corrnorm <- function(C, n) {
  U <- chol(C)
  nC <- dim(C)[1]
  N <- (n %% nC) + n
  result <- array(rnorm(N), dim=c(N/nC,nC)) %*% U
  return(invisible(result[1:n]))
}

corrnorm2 <- function(C, n) {
  U <- chol(C)
  nC <- dim(C)[1]
  N <- (n %% nC) + n
  result <- array(rnorm(N), dim=c(N/nC,nC)) %*% U
  return(invisible(result))
}
etmc/hadron documentation built on Sept. 17, 2022, 7:33 p.m.