-# Sage doesn't load ~/.sage/init.sage during testing (sage -t), so we
-# have to explicitly mangle our sitedir here so that "mjo.cone"
-# resolves.
-from os.path import abspath
-from site import addsitedir
-addsitedir(abspath('../../'))
-
from sage.all import *
-from mjo.cone.cone import lyapunov_rank, random_element
def rearrangement_cone(p,n):
r"""
The Lyapunov rank of the rearrangement cone of order ``p`` in ``n``
dimensions is ``n`` for ``p == 1`` or ``p == n`` and one otherwise::
- sage: all([ lyapunov_rank(rearrangement_cone(p,n)) == n
+ sage: all([ rearrangement_cone(p,n).lyapunov_rank() == n
....: for n in range(2, 10)
....: for p in [1, n-1] ])
True
- sage: all([ lyapunov_rank(rearrangement_cone(p,n)) == 1
+ sage: all([ rearrangement_cone(p,n).lyapunov_rank() == 1
....: for n in range(3, 10)
....: for p in range(2, n-1) ])
True
sage: for n in range(2,10):
....: for p in range(1, n-1):
....: K = rearrangement_cone(p,n)
- ....: v = random_element(K)
+ ....: v = K.random_element()
....: if not has_rearrangement_property(v,p): print v
"""