powered by
gaussian.condition(y, .A, .B, .C, mu.x = rep(0, nrow(.A)), mu.y = rep(0, nrow(.B)), U = chol(.B))