From: Michael Orlitzky Date: Mon, 11 Jan 2016 03:39:54 +0000 (-0500) Subject: Add another test for motzkin_decomposition(). X-Git-Url: https://gitweb.michael.orlitzky.com/?a=commitdiff_plain;h=806b8e39118cadb56714b99589f77bc9fee5828d;p=sage.d.git Add another test for motzkin_decomposition(). --- 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]_::