vectors = [ W(tp.list()) for tp in tensor_products ]
# Create the *dual* cone of the positive operators, expressed as
- # long vectors. WARNING: takes forever unless we pass check=False
- # to Cone().
+ # long vectors. WARNING: check=True is necessary even though it
+ # makes Cone() take forever. For an example take
+ # K = Cone([(1,0,0),(0,0,1),(0,0,-1)]).
pi_dual = Cone(vectors, ToricLattice(W.dimension()))
# Now compute the desired cone from its dual...
vectors = [ W(m.list()) for m in tensor_products ]
# Create the *dual* cone of the cross-positive operators,
- # expressed as long vectors. WARNING: takes forever unless we pass
- # check=False to Cone().
+ # expressed as long vectors. WARNING: check=True is necessary
+ # even though it makes Cone() take forever. For an example take
+ # K = Cone([(1,0,0),(0,0,1),(0,0,-1)]).
Sigma_dual = Cone(vectors, lattice=ToricLattice(W.dimension()))
# Now compute the desired cone from its dual...