+ ::
+
+ sage: set_random_seed()
+ sage: K = random_cone(max_ambient_dim = 8,
+ ....: strictly_convex=False,
+ ....: solid=True)
+ sage: K_S = _rho(K)
+ sage: K_SP = _rho(K_S.dual()).dual()
+ sage: K_SP.is_proper()
+ True
+ sage: K_SP = _rho(K_S, K_S.dual())
+ sage: K_SP.is_proper()
+ True
+
+ ::
+
+ sage: set_random_seed()
+ sage: K = random_cone(max_ambient_dim = 8,
+ ....: strictly_convex=True,
+ ....: solid=True)
+ sage: K_S = _rho(K)
+ sage: K_SP = _rho(K_S.dual()).dual()
+ sage: K_SP.is_proper()
+ True
+ sage: K_SP = _rho(K_S, K_S.dual())
+ sage: K_SP.is_proper()
+ True
+
+ Test Proposition 7 in our paper concerning the duals and
+ restrictions. Generate a random cone, then create a subcone of
+ it. The operation of dual-taking should then commute with rho::
+
+ sage: set_random_seed()
+ sage: J = random_cone(max_ambient_dim = 8,
+ ....: solid=False,
+ ....: strictly_convex=False)
+ sage: K = Cone(random_sublist(J.rays(), 0.5), lattice=J.lattice())
+ sage: K_W_star = _rho(K, J).dual()
+ sage: K_star_W = _rho(K.dual(), J)
+ sage: _basically_the_same(K_W_star, K_star_W)
+ True
+
+ ::
+
+ sage: set_random_seed()
+ sage: J = random_cone(max_ambient_dim = 8,
+ ....: solid=True,
+ ....: strictly_convex=False)
+ sage: K = Cone(random_sublist(J.rays(), 0.5), lattice=J.lattice())
+ sage: K_W_star = _rho(K, J).dual()
+ sage: K_star_W = _rho(K.dual(), J)
+ sage: _basically_the_same(K_W_star, K_star_W)
+ True
+
+ ::
+
+ sage: set_random_seed()
+ sage: J = random_cone(max_ambient_dim = 8,
+ ....: solid=False,
+ ....: strictly_convex=True)
+ sage: K = Cone(random_sublist(J.rays(), 0.5), lattice=J.lattice())
+ sage: K_W_star = _rho(K, J).dual()
+ sage: K_star_W = _rho(K.dual(), J)
+ sage: _basically_the_same(K_W_star, K_star_W)
+ True
+
+ ::
+
+ sage: set_random_seed()
+ sage: J = random_cone(max_ambient_dim = 8,
+ ....: solid=True,
+ ....: strictly_convex=True)
+ sage: K = Cone(random_sublist(J.rays(), 0.5), lattice=J.lattice())
+ sage: K_W_star = _rho(K, J).dual()
+ sage: K_star_W = _rho(K.dual(), J)
+ sage: _basically_the_same(K_W_star, K_star_W)
+ True
+
+ """
+ if K2 is None:
+ K2 = K