From a0fc026db4351a7d74bb066b3c79a64a981b9c2b Mon Sep 17 00:00:00 2001 From: Michael Orlitzky Date: Fri, 23 Oct 2015 10:46:06 -0400 Subject: [PATCH] Add a test for the lineality of the cone of positive operators. --- mjo/cone/cone.py | 12 ++++++++++++ 1 file changed, 12 insertions(+) diff --git a/mjo/cone/cone.py b/mjo/cone/cone.py index e9d70e4..8f8f0d2 100644 --- a/mjo/cone/cone.py +++ b/mjo/cone/cone.py @@ -197,6 +197,18 @@ def positive_operator_gens(K): sage: actual == expected True + The lineality of the cone of positive operators is given by the + corollary in my paper:: + + sage: set_random_seed() + sage: K = random_cone(max_ambient_dim = 5) + sage: n = K.lattice_dim() + sage: pi_of_K = positive_operator_gens(K) + sage: L = ToricLattice(n**2) + sage: actual = Cone([p.list() for p in pi_of_K], lattice=L).lineality() + sage: expected = n**2 - K.dim()*K.dual().dim() + sage: actual == expected + 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 -- 2.44.2