Deep Intrinsic Coregionalization Gaussian Processes for Multiple Outputs

Fits two-layer deep intrinsic coregionalization Gaussian process surrogate models for multiple outputs observed at shared input locations. Uses Markov chain Monte Carlo sampling for kernel parameters and latent coordinates, with plug-in coregionalization estimates. Provides burn-in removal and thinning, posterior predictive means and within-location output covariances, and sequential experimental design using an active learning criterion. The method is described by Chang and Sung (2026) . Adapts code from the 'deepgp' package, with underlying deep Gaussian process methodology described by Sauer, Gramacy, and Higdon (2023) .