Gaussian Process Modeling of Heterogeneity and Discontinuities Using Voronoi Tessellations. Issue 1 (25th January 2021)