code: def get_link_domain(link, dist):
      
      domain = np.array([-np.inf, -1, 0, 1, np.inf])
      domain = domain[~np.isnan(link.link(domain, dist))]
      return [domain[0], domain[-1]]
docs: tool to identify the domain of a given monotonic link function
  
      Parameters
      ----------
      link : Link object
      dist : Distribution object
  
      Returns
      -------
      domain : list of length 2, representing the interval of the domain.