code: def log_relative_gauss(z):
      
      if z < -6:
          return 1, -1.0e12, -1
      if z > 6:
          return 0, 0, 1
      else:
          logphi = -0.5 * (z * z + l2p)
          logPhi = np.log(.5 * special.erfc(-z / sq2))
          e = np.exp(logphi - logPhi)
      return e, logPhi, 0
docs: log_relative_gauss