def remove_subcomplexes(cxs, max_size=3):
    """
    Given a list of complexes, return a list without any complexes sized
    max_size or lower that are a subset of another complex.
    """
    cxs_copy = list(cxs)
    return [c for c in cxs if len(c) > max_size or len([c2 for c2 in cxs_copy
        if c.issubset(c2) and not c==c2]) == 0]