From 806b8e39118cadb56714b99589f77bc9fee5828d Mon Sep 17 00:00:00 2001 From: Michael Orlitzky Date: Sun, 10 Jan 2016 22:39:54 -0500 Subject: [PATCH] Add another test for motzkin_decomposition(). --- mjo/cone/cone.py | 8 ++++++++ 1 file changed, 8 insertions(+) diff --git a/mjo/cone/cone.py b/mjo/cone/cone.py index adba809..05c55ba 100644 --- a/mjo/cone/cone.py +++ b/mjo/cone/cone.py @@ -135,6 +135,14 @@ def motzkin_decomposition(K): sage: S.lineality() == S.dim() True + A strictly convex cone should be equal to its strictly convex component:: + + sage: set_random_seed() + sage: K = random_cone(max_ambient_dim=8, strictly_convex=True) + sage: (P,_) = motzkin_decomposition(K) + sage: K.is_equivalent(P) + True + The generators of the components are obtained from orthogonal projections of the original generators [Stoer-Witzgall]_:: -- 2.44.2