# Alert! If use original factorlist as outer loop, after removal, the elem index in the list changed
# but the loop counter does not - will cause step over matching some factors
# For every factor, find a cluque which factor is a subset of
for factor in factorlist[:]:  # Every matching obtain a copy of factorlist
    fvar = factor.var
    for cid in range(len(cliques)):
        clique = cliques[cid]
        #print("Now matching set_fvar",set(fvar))
        #print("With set_clique", set(clique))
        if set(fvar).issubset(set(clique)):
            # The clique potential is the product of its assigned factors
            clique_factors[cid] = factor_product(clique_factors[cid], factor)
            factorlist.remove(factor)
            break