[論文レビュー] Complexity Estimates for Two Uncoupling Algorithms

Alin Bostan, Frédéric Chyzak+1 more|Jan 23, 2013

Polynomial and algebraic computation参考文献 22被引用数 5

ひとこと要約

本稿は、2つの分離化アルゴリズム—サイクリック・ベクタ法（CVM）およびダニレフスキー＝バルカトゥー＝ツルゼル（DBZ）—に対する最初の厳密な複雑性解析を提供し、一般の場合に両者とも多項式時間の複雑性を持つことを示している。本稿では、出力サイズに対して準最適な複雑性を達成する高速版CVMを提案し、DBZと比較してほぼ2桁の速度向上を実現した。さらに、両アルゴリズムの深い構造的関係を解明し、CVMの出力複雑性に関する長年の誤解を解消した。

ABSTRACT

Uncoupling algorithms transform a linear differential system of first order into one or several scalar differential equations. We examine two approaches to uncoupling: the cyclic-vector method (CVM) and the Danilevski-Barkatou-Zürcher algorithm (DBZ). We give tight size bounds on the scalar equations produced by CVM, and design a fast variant of CVM whose complexity is quasi-optimal with respect to the output size. We exhibit a strong structural link between CVM and DBZ enabling to show that, in the generic case, DBZ has polynomial complexity and that it produces a single equation, strongly related to the output of CVM. We prove that algorithm CVM is faster than DBZ by almost two orders of magnitude, and provide experimental results that validate the theoretical complexity analyses.

研究の動機と目的

• 分離化アルゴリズム、特にサイクリック・ベクタ法（CVM）に対して理論的複雑性解析が不足している問題に対処すること。CVMは、出力複雑性が高いため、長年にわたり無視されてきた。
• CVMおよびDBZの計算複雑性を厳密に分析し、CVMが過剰に複雑な出力を生成するという一般的な誤解に反論すること。
• 出力サイズに対して準最適な複雑性を達成するCVMの高速版を設計すること。
• CVMとDBZの間の構造的およびアルゴリズム的関係を確立し、一般の場合に両者が同一の結果を生成することを示すこと。
• Magmaを用いた大規模行列を用いた広範な実験的評価を通じて、理論的複雑性境界の妥当性を検証すること。

提案手法

• 本稿は、DBZの操作を基底変換の観点から再解釈し、CVMで用いられるサイクリック・ベクタ基底と関連付ける新しい代数的解釈を導入する。
• CVMが生成するスカラー微分方程式の次数にきつい上限を導出し、出力次数が指数的ではなく、二次的に増加することを示した。これは、かつて懸念されていた指数的増加とは対照的である。
• 高速行列および多項式演算を活用し、高速行列乗算を応用することで、複雑性が ≈ n^{θ+1}d（2 ≤ θ ≤ 3）となるCVMの高速版を設計した。
• 著者らは、一般の場合にDBZとCVMが同一のスカラー方程式を生成することを証明した。これは、DBZが複数の式に分解されているように見えるにもかかわらず成立する。
• Magmaのタイミングデータを用いた線形回帰分析により、理論的複雑性境界を実験的に検証した。アルゴリズムの性能をさまざまな行列サイズおよび次数で分析した。
• 実装では、Magmaの最適化された多項式および行列演算を用い、スケーラビリティおよび実用的効率性を評価するために、メモリ使用量とパフォーマンスの監視を徹底的に行った。

実験結果

リサーチクエスチョン

• RQ1 サイクリック・ベクタ法（CVM）は、一般的に信じられているように、指数的複雑性の出力を本当に生成するのか？
• RQ2 ダニレフスキー＝バルカトゥー＝ツルゼル（DBZ）アルゴリズムの実際の計算複雑性は何か？ そしてCVMと比較してどうなるか？
• RQ3 出力サイズに対して準最適な複雑性を達成するCVMの高速版を設計できるか？
• RQ4 CVMとDBZの間に、一般の場合に同一の出力を生じる理由を説明する構造的関係が存在するか？
• RQ5 理論的複雑性境界は、大規模行列における実際のパフォーマンスと一致するか？

主な発見

• サイクリック・ベクタ法（CVM）は、一般の場合に多項式時間の複雑性を持つ。出力次数は指数的ではなく、二次的に増加する。これは、長年にわたり信じられていた仮定に反する。
• 出力サイズが ≈ n^3d であるのに対し、複雑性が ≈ n^{θ+1}d（2 ≤ θ ≤ 3）である、出力サイズに対して準最適なCVMの高速版を設計した。
• DBZも一般の場合に多項式時間の複雑性を持つ。理論的複雑性は ≈ n^5d であり、一般条件下ではCVMと同一のスカラー方程式を生成する。
• 実験的結果により、高速CVMバージョンが大規模行列に対してDBZの約2桁の速度で上回ることが確認され、理論的予測と一致した。
• CVMのパフォーマンスを支配するのは、ストルジョアンの行列逆行列計算アルゴリズムの実装であり、そのメモリ制限により、現在のMagmaベースのコードでは入力サイズが n ≤ 130 に制限される。
• DBZの実験的指数は予測された5ではなく6.01であった。これは、多項式ではなく有理関数係数の取り扱いに起因し、数値安定性とパフォーマンスに影響を与えていると考えられる。