expand function

Fitting low-rank nonstationary spatial Gaussian process models through dimension expansion