- 1. M.S. Gowda and J. Tao. On the bilinearity rank of a proper cone
- and Lyapunov-like transformations, Mathematical Programming, 147
+ .. [Gowda/Tao] M.S. Gowda and J. Tao. On the bilinearity rank of a proper
+ cone and Lyapunov-like transformations, Mathematical Programming, 147
- 2. G. Rudolf, N. Noyan, D. Papp, and F. Alizadeh, Bilinear
+ .. [Rudolf et al.] G. Rudolf, N. Noyan, D. Papp, and F. Alizadeh, Bilinear
optimality constraints for the cone of positive polynomials,
Mathematical Programming, Series B, 129 (2011) 5-31.
EXAMPLES:
optimality constraints for the cone of positive polynomials,
Mathematical Programming, Series B, 129 (2011) 5-31.
EXAMPLES:
sage: positives = Cone([(1,)])
sage: lyapunov_rank(positives)
sage: positives = Cone([(1,)])
sage: lyapunov_rank(positives)
sage: quadrant = Cone([(1,0), (0,1)])
sage: lyapunov_rank(quadrant)
2
sage: quadrant = Cone([(1,0), (0,1)])
sage: lyapunov_rank(quadrant)
2
- sage: octant = Cone([(1,0,0), (0,1,0), (0,0,1)])
+ sage: octant = Cone([(1,0,0), (0,1,0), (0,0,1)])
sage: L31 = Cone([(1,0,1), (0,-1,1), (-1,0,1), (0,1,1)])
sage: lyapunov_rank(L31)
1
sage: L31 = Cone([(1,0,1), (0,-1,1), (-1,0,1), (0,1,1)])
sage: lyapunov_rank(L31)
1
sage: L3infty = Cone([(0,1,1), (1,0,1), (0,-1,1), (-1,0,1)])
sage: lyapunov_rank(L3infty)
1
sage: L3infty = Cone([(0,1,1), (1,0,1), (0,-1,1), (-1,0,1)])
sage: lyapunov_rank(L3infty)
1
sage: L31 = Cone([(1,0,1), (0,-1,1), (-1,0,1), (0,1,1)])
sage: octant = Cone([(1,0,0), (0,1,0), (0,0,1)])
sage: L31 = Cone([(1,0,1), (0,-1,1), (-1,0,1), (0,1,1)])
sage: octant = Cone([(1,0,0), (0,1,0), (0,0,1)])
- Two isomorphic cones should have the same Lyapunov rank. The cone
- ``K`` in the following example is isomorphic to the nonnegative
+ Two isomorphic cones should have the same Lyapunov rank [Rudolf et al.]_.
+ The cone ``K`` in the following example is isomorphic to the nonnegative
octant in `\mathbb{R}^{3}`::
sage: K = Cone([(1,2,3), (-1,1,0), (1,0,6)])
octant in `\mathbb{R}^{3}`::
sage: K = Cone([(1,2,3), (-1,1,0), (1,0,6)])
sage: K = Cone([(2,2,4), (-1,9,0), (2,0,6)])
sage: lyapunov_rank(K) == lyapunov_rank(K.dual())
sage: K = Cone([(2,2,4), (-1,9,0), (2,0,6)])
sage: lyapunov_rank(K) == lyapunov_rank(K.dual())
sage: K1 = random_cone(max_dim=10, max_rays=10)
sage: K2 = random_cone(max_dim=10, max_rays=10)
sage: K1 = random_cone(max_dim=10, max_rays=10)
sage: K2 = random_cone(max_dim=10, max_rays=10)
sage: K = random_cone(max_dim=10, max_rays=10)
sage: lyapunov_rank(K) == lyapunov_rank(K.dual())
True
sage: K = random_cone(max_dim=10, max_rays=10)
sage: lyapunov_rank(K) == lyapunov_rank(K.dual())
True
+ The Lyapunov rank of a proper polyhedral cone in `n` dimensions can
+ be any number between `1` and `n` inclusive, excluding `n-1`
+ [Gowda/Tao]_ (by accident, this holds for the trivial cone in a
+ trivial space as well)::
+
+ sage: K = random_cone(max_dim=10, strictly_convex=True, solid=True)
+ sage: b = lyapunov_rank(K)
+ sage: n = K.lattice_dim()
+ sage: 1 <= b and b <= n
+ True
+ sage: b == n-1
+ False
+