code: def get_moments(model,x):
      '''
      Moments (mean and sdev.) of a GP model at x
  
      '''
      input_dim = model.X.shape[1]
      x = reshape(x,input_dim)
      fmin = min(model.predict(model.X)[0])
      m, v = model.predict(x)
      s = np.sqrt(np.clip(v, 0, np.inf))
      return (m,s, fminf get_moments(model,x):
      '''
      Moments (mean and sdev.) of a GP model at x
  
      '''
      input_dim = model.X.shape[1]
      x = reshape(x,input_dim)
      fmin = min(model.predict(model.X)[0])
      m, v = model.predict(x)
      s = np.sqrt(np.clip(v, 0, np.inf))
      return (m,s, fmin)
docs: Moments (mean and sdev.) of a GP model at x