def METHOD_NAME(node_status_at_link):
    """Find links that are inactive.
    A link is inactive if it connects two boundary nodes or one of
    its nodes is closed.
    Parameters
    ----------
    node_status_at_link : ndarray of int, shape `(n_links, 2)`
        Node status a link tail and head.
    Returns
    -------
    ndarray of bool, shape `(n_links, )`
        True if link is isactive.
    Examples
    --------
    >>> from landlab.grid.linkstatus import is_inactive_link
    >>> from landlab import NodeStatus
    >>> is_inactive_link([NodeStatus.CORE, NodeStatus.CLOSED])
    array([ True], dtype=bool)
    >>> is_inactive_link([NodeStatus.FIXED_GRADIENT, NodeStatus.FIXED_VALUE])
    array([ True], dtype=bool)
    >>> is_inactive_link([[NodeStatus.FIXED_GRADIENT, NodeStatus.CLOSED],
    ...                   [NodeStatus.CORE, NodeStatus.CORE]])
    array([ True, False], dtype=bool)
    """
    node_status_at_link = np.asarray(node_status_at_link).reshape((-1, 2))
    is_core = node_status_at_link == NodeStatus.CORE
    is_fixed_value = node_status_at_link == NodeStatus.FIXED_VALUE
    is_fixed_gradient = node_status_at_link == NodeStatus.FIXED_GRADIENT
    is_closed = node_status_at_link == NodeStatus.CLOSED
    is_boundary_node = is_fixed_value | is_fixed_gradient | is_closed
    return (
        (is_boundary_node[:, 0] & is_boundary_node[:, 1])
        | (is_closed[:, 0] & is_core[:, 1])
        | (is_core[:, 0] & is_closed[:, 1])
    )