Gaussian Process Modeling for Heterogeneous Functions