code: def EI_gaussian(mean, var, thresh):
      """Expected Improvement of Gaussian over threshold
  
      (estimated analytically)
      """
      sigma = np.sqrt(var)
      score = (mean - thresh) / sigma
      n = scipy.stats.norm
      return sigma * (score * n.cdf(score) + n.pdf(score))
  
imports: import numpy
  import scipy.stats
label_desc: Use image node_id 3 for calling a global function with example usage: EI_gaussian(mean, var, thresh) and returns: unknown