+ sage: set_random_seed()
+ sage: K = random_cone(max_ambient_dim=5)
+ sage: pi_of_K = positive_operator_gens(K)
+ sage: Z_of_K = Z_transformation_gens(K)
+ sage: L = ToricLattice(K.lattice_dim()**2)
+ sage: pi_star = Cone([p.list() for p in pi_of_K], lattice=L).dual()
+ sage: z_star = Cone([ z.list() for z in Z_of_K], lattice=L).dual()
+ sage: pi_star.linear_subspace() == z_star.linear_subspace()
+ True
+ """
+ # Matrices are not vectors in Sage, so we have to convert them
+ # to vectors explicitly before we can find a basis. We need these
+ # two values to construct the appropriate "long vector" space.
+ F = K.lattice().base_field()
+ n = K.lattice_dim()
+
+ # These tensor products contain generators for the dual cone of
+ # the cross-positive transformations.
+ tensor_products = [ s.tensor_product(x)
+ for (x,s) in K.discrete_complementarity_set() ]
+
+ # Turn our matrices into long vectors...
+ W = VectorSpace(F, n**2)
+ vectors = [ W(m.list()) for m in tensor_products ]
+
+ # Create the *dual* cone of the cross-positive operators,
+ # expressed as long vectors..
+ Sigma_dual = Cone(vectors, lattice=ToricLattice(W.dimension()))
+
+ # Now compute the desired cone from its dual...
+ Sigma_cone = Sigma_dual.dual()
+
+ # And finally convert its rays back to matrix representations.
+ # But first, make them negative, so we get Z-transformations and
+ # not cross-positive ones.
+ M = MatrixSpace(F, n)
+ return [ -M(v.list()) for v in Sigma_cone.rays() ]
+
+
+def Z_cone(K):
+ gens = Z_transformation_gens(K)
+ L = None
+ if len(gens) == 0:
+ L = ToricLattice(0)
+ return Cone([ g.list() for g in gens ], lattice=L)
+
+def pi_cone(K):
+ gens = positive_operator_gens(K)
+ L = None
+ if len(gens) == 0:
+ L = ToricLattice(0)
+ return Cone([ g.list() for g in gens ], lattice=L)