[論文レビュー] $3$-dimensional Continued Fraction Algorithms Cheat Sheets

Sebastien Labbé|arXiv (Cornell University)|Nov 26, 2015

Mathematical Dynamics and Fractals参考文献 6被引用数 8

ひとこと要約

本論文は、7種類の3次元連分数アルゴリズム（Brun, Selmer, Poincaré, Fully Subtractive, Arnoux-Rauzy-Poincaré, Reverse, Cassaigne）の包括的な「 Cheat Sheets 」を提示する。各アルゴリズムの力学的性質、不変測度、リャプノフ指数、記号的力学系、およびデジタル幾何学について詳細に記述している。Sageとslabbeパッケージを用いて実装されており、不変測度の密度、シリンダ集合、S-Adic語の複雑度、不均一性、双対的置換といった再現可能な計算結果が得られており、特にリャプノフ指数のアルゴリズム間比較が重要な分析として含まれる。

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 $\mathbb{R}^d$. We consider multidimensional continued fraction algorithms that acts symmetrically on the positive cone $\mathbb{R}^d_+$ for $d=3$. 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 with the optional Sage package exttt{slabbe-0.2.spkg}. The information includes the $n$-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 $S$-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.

研究の動機と目的

• 3次元連分数アルゴリズムの力学的および数論的性質を体系的かつ文書化し、研究および教育的利用を目的とする。
• オープンソースソフトウェア（Sageとslabbeパッケージ）を用いた統一的な計算フレームワークを提供し、多次元連分数の分析を可能にする。
• リャプノフ指数と不均一性の指標を用いて、さまざまな3次元アルゴリズムの効率性および収束特性を比較する。
• アルゴリズム的置換から生成されるS-Adic語、因子複雑度、およびデジタル平面を生成・分析する。
• 公開されたSageコードと実験データを用いて、結果の再現性を確保する。

提案手法

• 本研究では、Sageとslabbe-0.2.spkgパッケージを用いた計算フレームワークを採用し、2-単体上での行列コキューブおよび射影写像を用いて各アルゴリズムを実装・分析する。
• 各アルゴリズムについて、1億回の反復を30軌道で行った上で、n-シリンダ、不変測度の密度、自然拡張、リャプノフ指数を計算・可視化する。
• 記号的力学系は置換とS-Adic語を用いて生成され、代表的な方向（例：(1,e,π)）について因子複雑度と不均一性が計算される。
• 双対的置換は元の置換から導出され、3次元空間におけるデジタル平面生成を可視化するためのE₁*パッチが計算される。
• リャプノフ指数は、ビンゴン・デレクロワ氏のCコードをCythonに移植したバージョンを用いて計算され、高精度な統計的推定が保証される。
• 図や表を含むすべての結果は再現可能なSageコードによって生成されており、図やデータの完全な再現が可能である。

実験結果

リサーチクエスチョン

• RQ1 異なる3次元連分数アルゴリズムのリャプノフ指数は、収束速度および安定性の観点でどのように比較できるか？
• RQ2 各アルゴリズムが2-単体上に持つ不変測度の構造とその密度関数はどのようなものか？
• RQ3 各アルゴリズムが生成するS-Adic語の記号的力学系および因子複雑度は、一般の方向に対してどのように振る舞うか？
• RQ4 生成されたデジタル平面の不均一性は何か？ また、アルゴリズム間でどのように変化するか？
• RQ5 双対的置換とE₁*パッチは、アルゴリズムの幾何学的および算術的性質をどのように反映しているか？

主な発見

• Arnoux-Rauzy-Poincaréアルゴリズムは、最大の第一リャプノフ指数（θ₁ ≈ 0.44290）と、1−θ₂/θ₁ ≈ 1.38879 の比を示し、強い収束性を示している。
• Brunアルゴリズムは、θ₁ ≈ 0.30449 および 1−θ₂/θ₁ ≈ 1.36833 を示し、他のアルゴリズムと比較して相対的に高い収束速度を示している。
• Reverseアルゴリズムは、θ₁ ≈ 0.40489 および 1−θ₂/θ₁ ≈ 1.25489 を示し、Brunより速い収束を示すが、AR-Pほど安定ではない。
• Fully SubtractiveおよびPoincaréアルゴリズムは極めて小さなリャプノフ指数（θ₁ ≈ 2.5e−9 および 6.9e−9）を示し、ほぼゼロまたは発散的挙動を示している。
• CassaigneおよびSelmerアルゴリズムのS-Adic語の因子複雑度は、n ≤ 20 において p_w(n) ≈ 2n+1 として線形に増加し、低複雑度および高い規則性を示している。
• 合計が200である19,701の方向についてのS-Adic語の不均一性は一貫した分布パターンを示しており、CassaigneおよびSelmerアルゴリズムは低い不均一性値を達成している。