CoCalc Shared Filessage-surface_dynamics.sagews
Author: Harald Schilly
Views : 24

# Sage Surface Dynamics Module

https://pypi.org/project/surface_dynamics/

version()
'SageMath version 8.2, Release Date: 2018-05-05'
from surface_dynamics.all import Origami o = Origami('(1,2)','(1,3)') o
(1,2)(3) (1,3)(2)
o.sum_of_lyapunov_exponents()
Error in lines 1-1 Traceback (most recent call last): File "/cocalc/lib/python2.7/site-packages/smc_sagews/sage_server.py", line 1044, in execute exec compile(block+'\n', '', 'single', flags=compile_flags) in namespace, locals File "", line 1, in <module> File "surface_dynamics/flat_surfaces/origamis/origami_dense.pyx", line 3215, in surface_dynamics.flat_surfaces.origamis.origami_dense.Origami_dense_pyx.sum_of_lyapunov_exponents return self.teichmueller_curve().sum_of_lyapunov_exponents() File "/ext/sage/sage-8.2_1604/local/lib/python2.7/site-packages/surface_dynamics/flat_surfaces/origamis/teichmueller_curve.py", line 367, in sum_of_lyapunov_exponents K = Integer(1)/Integer(12) * sum(m*(m+2)/(m+1) for m in self.stratum().zeros()) File "/ext/sage/sage-8.2_1604/local/lib/python2.7/site-packages/surface_dynamics/flat_surfaces/origamis/teichmueller_curve.py", line 226, in stratum return self[0].stratum() File "surface_dynamics/flat_surfaces/origamis/origami_dense.pyx", line 3063, in surface_dynamics.flat_surfaces.origamis.origami_dense.Origami_dense_pyx.stratum degrees = self.vertex_degrees() File "surface_dynamics/flat_surfaces/origamis/origami_dense.pyx", line 3258, in genexpr return sorted((v.degree() for v in self.vertices()), reverse=True) File "surface_dynamics/flat_surfaces/origamis/origami_dense.pyx", line 3258, in genexpr return sorted((v.degree() for v in self.vertices()), reverse=True) File "surface_dynamics/flat_surfaces/origamis/origami_dense.pyx", line 3240, in surface_dynamics.flat_surfaces.origamis.origami_dense.Origami_dense_pyx.vertices return [v for v in self._vertices() if v.degree()] File "sage/misc/cachefunc.pyx", line 2909, in sage.misc.cachefunc.CachedMethod.__get__ (build/cythonized/sage/misc/cachefunc.c:15723) Caller = CachedMethodCaller(self, inst, File "sage/misc/cachefunc.pyx", line 1813, in sage.misc.cachefunc.CachedMethodCaller.__init__ (build/cythonized/sage/misc/cachefunc.c:9837) cachedmethod._cachedfunc.argfix_init() File "sage/misc/cachefunc.pyx", line 851, in sage.misc.cachefunc.CachedFunction.argfix_init (build/cythonized/sage/misc/cachefunc.c:4504) self._argument_fixer = ArgumentFixer(self.f, File "sage/misc/function_mangling.pyx", line 120, in sage.misc.function_mangling.ArgumentFixer.__init__ (build/cythonized/sage/misc/function_mangling.c:1434) arg_names, varargs, varkw, defaults = sage_getargspec(f) File "/ext/sage/sage-8.2_1604/local/lib/python2.7/site-packages/sage/misc/sageinspect.py", line 1503, in sage_getargspec source = sage_getsource(obj) File "/ext/sage/sage-8.2_1604/local/lib/python2.7/site-packages/sage/misc/sageinspect.py", line 1848, in sage_getsource t = sage_getsourcelines(obj) File "/ext/sage/sage-8.2_1604/local/lib/python2.7/site-packages/sage/misc/sageinspect.py", line 2170, in sage_getsourcelines return sage_getsourcelines(obj.__class__) File "/ext/sage/sage-8.2_1604/local/lib/python2.7/site-packages/sage/misc/sageinspect.py", line 2168, in sage_getsourcelines return sage_getsourcelines(B) File "/ext/sage/sage-8.2_1604/local/lib/python2.7/site-packages/sage/misc/sageinspect.py", line 2171, in sage_getsourcelines raise err IOError: source code not available
o.lyapunov_exponents_approx() # abs tol 0.05
Error in lines 1-1 Traceback (most recent call last): File "/cocalc/lib/python2.7/site-packages/smc_sagews/sage_server.py", line 1044, in execute exec compile(block+'\n', '', 'single', flags=compile_flags) in namespace, locals File "", line 1, in <module> File "surface_dynamics/flat_surfaces/origamis/origami_dense.pyx", line 1611, in surface_dynamics.flat_surfaces.origamis.origami_dense.Origami_dense_pyx.lyapunov_exponents_approx nb_vectors = self.genus()-1 File "surface_dynamics/flat_surfaces/origamis/origami_dense.pyx", line 3082, in surface_dynamics.flat_surfaces.origamis.origami_dense.Origami_dense_pyx.genus return self.stratum().genus() File "surface_dynamics/flat_surfaces/origamis/origami_dense.pyx", line 3063, in surface_dynamics.flat_surfaces.origamis.origami_dense.Origami_dense_pyx.stratum degrees = self.vertex_degrees() File "surface_dynamics/flat_surfaces/origamis/origami_dense.pyx", line 3258, in genexpr return sorted((v.degree() for v in self.vertices()), reverse=True) File "surface_dynamics/flat_surfaces/origamis/origami_dense.pyx", line 3258, in genexpr return sorted((v.degree() for v in self.vertices()), reverse=True) File "surface_dynamics/flat_surfaces/origamis/origami_dense.pyx", line 3240, in surface_dynamics.flat_surfaces.origamis.origami_dense.Origami_dense_pyx.vertices return [v for v in self._vertices() if v.degree()] File "sage/misc/cachefunc.pyx", line 2909, in sage.misc.cachefunc.CachedMethod.__get__ (build/cythonized/sage/misc/cachefunc.c:15723) Caller = CachedMethodCaller(self, inst, File "sage/misc/cachefunc.pyx", line 1813, in sage.misc.cachefunc.CachedMethodCaller.__init__ (build/cythonized/sage/misc/cachefunc.c:9837) cachedmethod._cachedfunc.argfix_init() File "sage/misc/cachefunc.pyx", line 851, in sage.misc.cachefunc.CachedFunction.argfix_init (build/cythonized/sage/misc/cachefunc.c:4504) self._argument_fixer = ArgumentFixer(self.f, File "sage/misc/function_mangling.pyx", line 120, in sage.misc.function_mangling.ArgumentFixer.__init__ (build/cythonized/sage/misc/function_mangling.c:1434) arg_names, varargs, varkw, defaults = sage_getargspec(f) File "/ext/sage/sage-8.2_1604/local/lib/python2.7/site-packages/sage/misc/sageinspect.py", line 1503, in sage_getargspec source = sage_getsource(obj) File "/ext/sage/sage-8.2_1604/local/lib/python2.7/site-packages/sage/misc/sageinspect.py", line 1848, in sage_getsource t = sage_getsourcelines(obj) File "/ext/sage/sage-8.2_1604/local/lib/python2.7/site-packages/sage/misc/sageinspect.py", line 2170, in sage_getsourcelines return sage_getsourcelines(obj.__class__) File "/ext/sage/sage-8.2_1604/local/lib/python2.7/site-packages/sage/misc/sageinspect.py", line 2168, in sage_getsourcelines return sage_getsourcelines(B) File "/ext/sage/sage-8.2_1604/local/lib/python2.7/site-packages/sage/misc/sageinspect.py", line 2171, in sage_getsourcelines raise err IOError: source code not available
o.veech_group()
Arithmetic subgroup with permutations of right cosets S2=(2,3) S3=(1,2,3) L=(1,2) R=(1,3)
H_5(2^4)
q.components()
[Q_2(1^4)^hyp]
c = q.components()[0] c
Q_2(1^4)^hyp
c.orientation_cover_component()
H_5(2^4)^odd
from surface_dynamics.all import AbelianStrata AbelianStrata(genus=3).list()
[H_3(4), H_3(3, 1), H_3(2^2), H_3(2, 1^2), H_3(1^4)]
from surface_dynamics.all import OrigamiDatabase, AbelianStratum O = OrigamiDatabase() q = O.query(("stratum","=",AbelianStratum(2)), ("nb_squares","=",5)) q.number_of()
2
for o in q: print("%s\n- - - - - - - -" % o)
(1)(2)(3)(4,5) (1,2,3,4)(5) - - - - - - - - (1)(2)(3,4,5) (1,2,3)(4)(5) - - - - - - - -
Q12_reg = QuadraticStratum(12).regular_component() Q12_reg.lyapunov_exponents_H_plus(nb_iterations=2**20) # abs tol 0.05
[0.6619595333890421, 0.44903083817447104, 0.22962202171849683, 0.08650300278009455]
Q12_reg.lyapunov_exponents_H_minus(nb_iterations=2**20) # abs tol 0.05
Error in lines 1-1 Traceback (most recent call last): File "/cocalc/lib/python2.7/site-packages/smc_sagews/sage_server.py", line 1044, in execute exec compile(block+'\n', '', 'single', flags=compile_flags) in namespace, locals File "", line 1, in <module> File "/ext/sage/sage-8.2_1604/local/lib/python2.7/site-packages/surface_dynamics/flat_surfaces/quadratic_strata.py", line 924, in lyapunov_exponents_H_minus return perm.lyapunov_exponents_H_plus(**kargs) File "/ext/sage/sage-8.2_1604/local/lib/python2.7/site-packages/surface_dynamics/interval_exchanges/cover.py", line 864, in lyapunov_exponents_H_plus dimensions, projections, None, verbose) File "src/cysignals/signals.pyx", line 263, in cysignals.signals.python_check_interrupt File "src/cysignals/signals.pyx", line 96, in cysignals.signals.sig_raise_exception KeyboardInterrupt