gp function

Gaussian process emulator construction