From: Michael Orlitzky <michael@orlitzky.com>
Date: Thu, 12 Nov 2015 20:14:24 +0000 (-0500)
Subject: Add another test for the lineality space of Z(K) being LL(K).
X-Git-Url: https://gitweb.michael.orlitzky.com/?a=commitdiff_plain;h=c8c4a0927a7a83a5b9f0cd3c33ff450212b95b99;hp=b397966ded7306845c51c90efb559686f500eb04;p=sage.d.git

Add another test for the lineality space of Z(K) being LL(K).
---

diff --git a/mjo/cone/cone.py b/mjo/cone/cone.py
index 8f8f0d2..1ab6b97 100644
--- a/mjo/cone/cone.py
+++ b/mjo/cone/cone.py
@@ -301,6 +301,13 @@ def Z_transformation_gens(K):
         sage: z_cone.linear_subspace() == lls
         True
 
+    And thus, the lineality of Z is the Lyapunov rank::
+
+        sage: set_random_seed()
+        sage: K = random_cone(min_ambient_dim = 1, max_ambient_dim = 6)
+        sage: z_cone  = Cone([ z.list() for z in Z_transformation_gens(K) ])
+        sage: z_cone.lineality() == K.lyapunov_rank()
+        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