[論文レビュー] Double Precision Computation of the Logistic Map Depends on Computational Modes of the Floating-point Processing Unit

Michiro Yabuki, Takashi Tsuchiya|arXiv (Cornell University)|May 14, 2013

Chaos control and synchronization参考文献 4被引用数 9

ひとこと要約

本稿は、浮動小数点演算装置（FPU）の計算モード（[53:53] 対 [53:64]）およびコード内で使用されるロジスティック写像の代数的表現の違いに応じて、二重精度計算におけるロジスティック写像に顕著な仮想周期性および分岐図の差異が生じることを示している。主な貢献は、これらの実装固有の要因がカオス的システムのシミュレーションにおける再現可能性に顕著に影響することを明らかにしたことであり、信頼できる数値的研究所のためには、計算モードおよび式の表現を標準化して報告する必要があることを示唆している。

ABSTRACT

Today's most popular CPU can operate in two different computational modes for double precision computations. This fact is not very widely recognized among scientific computer users. The present paper reports the differences the modes bring about using the most thoroughly studied system in chaos theory, the logistic map. Distinct virtual periods due to finite precision come about depending on the computational modes for the parameter value corresponding to fully developed chaos. For other chaotic regime various virtual periods emerge depending on the computational modes and the mathematical expressions of the map. Differences in the bifurcation diagrams due to the modes and the expressions are surveyed exhaustively. A quantity to measure those differences is defined and calculated.

研究の動機と目的

• 浮動小数点演算装置（FPU）の異なる計算モードが、カオス的システムの二重精度シミュレーションに与える影響を調査すること。
• ロジスティック写像の代数的表現のバリエーションが、カオス的領域における数値的結果に与える影響を検討すること。
• FPUモードおよび写像の定式化に起因する仮想周期性および分岐図の差異を定量化すること。
• 再現可能性を確保するため、数値的カオス研究においてFPU計算モードおよび数学的表現の標準化を提唱すること。

提案手法

• 本研究では、2つのFPU計算モードを比較した：[53:53]（両側とも53ビット仮数）および[53:64]（左側は53ビット、右側は64ビット）。
• 数式的に同等のロジスティック写像の6通りの表現形（L1–L6）を実装し、式依存の数値的アーティファクトをテストした。
• 仮想周期性の分析は、パラメータ値 a=4.0 および a=0x400fffffffffffff（4.0に非常に近い値）で、1,000個の均等に分布した初期条件を用いてロジスティック写像を反復して行った。
• 分岐図は、異なるFPUモードおよび写像表現の下で生成され、構造的差異を比較した。
• FPUモードおよび表現形ごとの分岐図の差異を区別するための定量的指標を定義し、計算した。
• 有限精度算術における仮数部の演算（m(x)）を用いた理論的分析により、特定の写像表現の代数的同等性を証明した。

実験結果

リサーチクエスチョン

• RQ1 [53:53] および [53:64] のFPU計算モードが、二重精度算術におけるロジスティック写像の仮想周期性にどのように影響するか？
• RQ2 ロジスティック写像の異なる代数的表現（例：ax(1−x) 対 4.0x−4.0xx）が、カオス的領域でどのように異なる数値的結果を生じるか？
• RQ3 FPU計算モードおよび写像表現が、二重精度シミュレーションにおける分岐図の構造にどのように影響するか？
• RQ4 FPUモードおよび表現形に起因する分岐図の差異を、定量的指標が信頼性を持って区別できるか？
• RQ5 なぜFPUモードおよび写像表現の標準化が、数値的カオス研究における再現可能性にとって不可欠なのか？

主な発見

• a=4.0 の場合、[53:53] モードとグループAの表現（L1–L3）を用いた場合、周波数678で主な仮想周期5,638,349が観測され、WangらおよびWagnerの先行研究と一致した。
• [53:64] モードとグループAの表現を用いた場合、最も頻度の高い仮想周期は588（周波数945,815）であり、FPUモードの変更に起因する周期性の顕著なシフトが示された。
• a=0x400fffffffffffff（4.0に非常に近い値）の場合、FPUモードと写像表現の12通りの組み合わせすべてが異なる仮想周期を示し、カテゴリ間で重複は一切なかった。
• 理論的分析により、表現 (L1)=(L2)=(L3) および (L4)=(L5)=(L6) が有限精度算術下で代数的に同等であることが確認されたが、FPUモードの影響により異なる数値的結果が得られた。
• 分岐図はFPUモードおよび表現グループごとに顕著に異なり、提案された定量的指標によりこれら差異が明確に区別された。
• 本研究は、仮想周期性および分岐構造がロジスティック写像そのものに内在するものではなく、FPUモードおよびコード定式化の結果であると結論づけた。これらが報告されない場合、再現性が損なわれる。