- sage: K = random_cone(max_dim = 8, strictly_convex=True, solid=True)
- sage: K_S = rho(K)
- sage: P = rho(K_S.dual()).dual()
- sage: P.is_proper()
- True
- sage: P = rho(K_S, K_S.dual())
- sage: P.is_proper()
- True
-
- Test the proposition in our paper concerning the duals, where the
- subspace `W` is the span of `K^{*}`::
-
- sage: set_random_seed()
- sage: K = random_cone(max_dim = 8, solid=False, strictly_convex=False)
- sage: K_W = rho(K, K.dual())
- sage: K_star_W_star = rho(K.dual()).dual()
- sage: basically_the_same(K_W, K_star_W_star)
- True
-
- ::
-
- sage: set_random_seed()
- sage: K = random_cone(max_dim = 8, solid=True, strictly_convex=False)
- sage: K_W = rho(K, K.dual())
- sage: K_star_W_star = rho(K.dual()).dual()
- sage: basically_the_same(K_W, K_star_W_star)
- True
-
- ::
-
- sage: set_random_seed()
- sage: K = random_cone(max_dim = 8, solid=False, strictly_convex=True)
- sage: K_W = rho(K, K.dual())
- sage: K_star_W_star = rho(K.dual()).dual()
- sage: basically_the_same(K_W, K_star_W_star)
- True
-
- ::
-
- sage: set_random_seed()
- sage: K = random_cone(max_dim = 8, solid=True, strictly_convex=True)
- sage: K_W = rho(K, K.dual())
- sage: K_star_W_star = rho(K.dual()).dual()
- sage: basically_the_same(K_W, K_star_W_star)
+ sage: J = random_cone(max_ambient_dim = 8)
+ sage: K = Cone(random_sublist(J.rays(), 0.5), lattice=J.lattice())
+ sage: K_W_star = _restrict_to_space(K, J.span()).dual()
+ sage: K_star_W = _restrict_to_space(K.dual(), J.span())
+ sage: _basically_the_same(K_W_star, K_star_W)