3-dimensional Continued Fraction Algorithms Cheat Sheets

-dimensional Continued Fraction Algorithms Cheat Sheets

Sébastien Labbé111 Université de Liège, Bât. B37 Institut de Mathématiques, Grande Traverse 12, 4000 Liège, Belgium, slabbe@ulg.ac.be.
Abstract

Multidimensional Continued Fraction Algorithms are generalizations of the Euclid algorithm and find iteratively the gcd of two or more numbers. They are defined as linear applications on some subcone of . We consider multidimensional continued fraction algorithms that acts symmetrically on the positive cone for . We include well-known and old ones (Poincaré, Brun, Selmer, Fully Subtractive) and new ones (Arnoux-Rauzy-Poincaré, Reverse, Cassaigne).

For each algorithm, one page (called cheat sheet) gathers a handful of informations most of them generated with the open source software Sage [sage] with the optional Sage package slabbe-0.2.spkg [labbe_slabbe_2015]. The information includes the -cylinders, density function of an absolutely continuous invariant measure, domain of the natural extension, lyapunov exponents as well as data regarding combinatorics on words, symbolic dynamics and digital geometry, that is, associated substitutions, generated -adic systems, factor complexity, discrepancy, dual substitutions and generation of digital planes.

The document ends with a table of comparison of Lyapunov exponents and gives the code allowing to reproduce any of the results or figures appearing in these cheat sheets.

Brun algorithm

Definition

On , the map

is defined by

where is the permutation of such that [MR0111735].

Matrix Definition

The partition of the cone is where

The matrices are given by the rule

The map on and the projective map on are:

Matrices

Cylinders

Density function

The density function of the invariant measure of for the Brun algorithm is [arnoux_symmetric_2015]:

on the part .

Invariant measure

Natural extension

Two sequences and defined such that

Lyapunov exponents

(using 30 orbits of 100000000 iterations each)
30 succesfull orbits min mean max std

Substitutions

-adic word example

Using vector :

Factor Complexity of is

Discrepancy

Discrepancy [MR593979] for all 19701 -adic words with directions such that :

Dual substitutions

E one star

Using vector , the 9-th iteration on the unit cube is:

Selmer algorithm

Definition

On , the map

is defined by

where is the permutation of such that [MR0130852].

Matrix Definition

The partition of the cone is where

The matrices are given by the rule

The map on and the projective map on are:

Matrices

Cylinders

Density function

The sorted version of admits a -finite invariant measure which is absolutely continuous with respect to Lebesgue measure on the central part and its density is known [schweiger].

Invariant measure

Natural extension

Two sequences and defined such that

Lyapunov exponents

(using 30 orbits of 100000000 iterations each)
30 succesfull orbits min mean max std

Substitutions

-adic word example

Using vector :

Factor Complexity of is

Discrepancy

ValueError: On input=[198, 1, 1], algorithm Selmer loops on (1.0, 1.0, 0.0)

Dual substitutions

E one star

Using vector , the 13-th iteration on the unit cube is:

Poincaré algorithm

Definition

On , the map

is defined by

where is the permutation of such that [MR1336331].

Matrix Definition

The partition of the cone is where

The matrices are given by the rule

The map on and the projective map on are:

Matrices

Cylinders

Density function

The sorted version of admits a -finite invariant measure which is absolutely continuous with respect to Lebesgue measure and its density is known [schweiger, MR1336331].

Invariant measure

Natural extension

Two sequences and defined such that

Lyapunov exponents

(using 30 orbits of 100000000 iterations each)
20 succesfull orbits min mean max std

Substitutions

-adic word example

Using vector :

Factor Complexity of is

Discrepancy

Discrepancy [MR593979] for all 19701 -adic words with directions such that :

Dual substitutions

E one star

Using vector , the 5-th iteration on the unit cube is:

Fully Subtractive algorithm

Definition

On , the map

is defined by

where is the permutation of such that [schweiger].

Matrix Definition

The partition of the cone is where

The matrices are given by the rule

The map on and the projective map on are:

Matrices

Cylinders

Density function

The sorted version of admits a -finite invariant measure which is absolutely continuous with respect to Lebesgue measure and its density is known [schweiger].

Invariant measure

Natural extension

Two sequences and defined such that

Lyapunov exponents

(using 30 orbits of 100000000 iterations each)
24 succesfull orbits min mean max std

Substitutions

-adic word example

Using vector :

Factor Complexity of is

Discrepancy

ValueError: On input=[198, 1, 1], algorithm Fully Subtractive loops on (197.0, 1.0, 0.0)

Dual substitutions

E one star

Using vector , the 7-th iteration on the unit cube is:

Arnoux-Rauzy-Poincaré algorithm

Definition

On , the map

is defined by

where is the permutation of such that [2015_berthe_factor].

Matrix Definition

The subcones are

The matrices are given by the rule

The map on and the projective map on are:

Matrices

Cylinders

Density function

The density of the absolutely continuous invariant measure is unknown [arnoux_symmetric_2015].

Invariant measure

Natural extension

Two sequences and defined such that

Lyapunov exponents

(using 30 orbits of 100000000 iterations each)
30 succesfull orbits min mean max std

Substitutions

-adic word example

Using vector :

Factor Complexity of is

Discrepancy

Discrepancy [MR593979] for all 19701 -adic words with directions such that :

Dual substitutions

E one star

Using vector , the 5-th iteration on the unit cube is:

Reverse algorithm

Definition

On , the map

is defined by