def pi_cone(K1, K2=None):
if K2 is None:
K2 = K1
gens = K1.positive_operators_gens(K2)
L = ToricLattice(K1.lattice_dim()*K2.lattice_dim())
def pi_cone(K1, K2=None):
if K2 is None:
K2 = K1
gens = K1.positive_operators_gens(K2)
L = ToricLattice(K1.lattice_dim()*K2.lattice_dim())