R/diagABA.R

diagABA <- function(A, B){
  D <- B
  D[lower.tri(D, diag=FALSE)] <- 0
  D[upper.tri(D, diag=FALSE)] <- 2*D[upper.tri(D, diag=FALSE)]
  D <- (A%*%D)*A
  rowSums(D)
}
IRSN/DiceKriging documentation built on May 8, 2019, 1:25 p.m.