def get_zonotope_vertices_and_bounds(active_subspace_eigvecs):
    """
    active_subspace_eigvecs : np.ndarray (num_vars,num_active_vars)
    """
    import active_subspaces as asub
    Y, X = asub.domains.zonotope_vertices(active_subspace_eigvecs)
    zonotope_max = Y.max(axis=0)
    zonotope_min = Y.min(axis=0)
    assert zonotope_min.shape[0] == active_subspace_eigvecs.shape[1]
    return Y.T, zonotope_min, zonotope_max
