From 3a312d50b5aba08e72039f1ebcde7b12c62a1e9f Mon Sep 17 00:00:00 2001 From: Michael Orlitzky Date: Sat, 30 May 2015 20:05:32 -0400 Subject: [PATCH] Add a few tests for LL(K). --- mjo/cone/cone.py | 66 ++++++++++++++++++++++++++++++++++++++++++++++-- 1 file changed, 64 insertions(+), 2 deletions(-) diff --git a/mjo/cone/cone.py b/mjo/cone/cone.py index 2e3dc8a..4b01936 100644 --- a/mjo/cone/cone.py +++ b/mjo/cone/cone.py @@ -87,8 +87,70 @@ def LL(K): OUTPUT: - A ``MatrixSpace`` object `M` such that every matrix `L \in M` is - Lyapunov-like on this cone. + A list of matrices forming a basis for the space of all + Lyapunov-like transformations on the given cone. + + EXAMPLES: + + The trivial cone has no Lyapunov-like transformations:: + + sage: L = ToricLattice(0) + sage: K = Cone([], lattice=L) + sage: LL(K) + [] + + The Lyapunov-like transformations on the nonnegative orthant are + simply diagonal matrices:: + + sage: K = Cone([(1,)]) + sage: LL(K) + [[1]] + + sage: K = Cone([(1,0),(0,1)]) + sage: LL(K) + [ + [1 0] [0 0] + [0 0], [0 1] + ] + + sage: K = Cone([(1,0,0),(0,1,0),(0,0,1)]) + sage: LL(K) + [ + [1 0 0] [0 0 0] [0 0 0] + [0 0 0] [0 1 0] [0 0 0] + [0 0 0], [0 0 0], [0 0 1] + ] + + Only the identity matrix is Lyapunov-like on the `L^{3}_{1}` and + `L^{3}_{\infty}` cones [Rudolf et al.]_:: + + sage: L31 = Cone([(1,0,1), (0,-1,1), (-1,0,1), (0,1,1)]) + sage: LL(L31) + [ + [1 0 0] + [0 1 0] + [0 0 1] + ] + + sage: L3infty = Cone([(0,1,1), (1,0,1), (0,-1,1), (-1,0,1)]) + sage: LL(L3infty) + [ + [1 0 0] + [0 1 0] + [0 0 1] + ] + + TESTS: + + The inner product `\left< L\left(x\right), s \right>` is zero for + every pair `\left( x,s \right)` in the discrete complementarity set + of the cone:: + + sage: K = random_cone(max_dim=8, max_rays=10) + sage: C_of_K = discrete_complementarity_set(K) + sage: l = [ (L*x).inner_product(s) for (x,s) in C_of_K for L in LL(K) ] + sage: sum(map(abs, l)) + 0 """ V = K.lattice().vector_space() -- 2.44.2