[論文レビュー] Computational linear algebra over finite fields

Jean‐Guillaume Dumas, Clément Pernet|arXiv (Cornell University)|Apr 17, 2012

Coding theory and cryptography参考文献 57被引用数 3

ひとこと要約

この論文は、ビットパッキング、ビットスライcing、およびクロネッカー代入技術を用いて、有限体上の密行列線形代数のための最適化アルゴリズムを提示する。特に、F₂などの小さな体における行列乗算、逆行列計算、分解の効率化に注力している。F₂上での行列乗算において、O(n³/log n)の時間計算量にまで低減され、サブキュービックな複雑性を達成しており、M4RI、Magma、Sageなどの高性能ライブラリに統合されている。

ABSTRACT

We present here algorithms for efficient computation of linear algebra problems over finite fields.

研究の動機と目的

• 密行列線形代数の計算性能ボトルネックを、低レベルの算術演算とデータ表現の最適化によって解消すること。
• 特に小さな有限体に対して、行列乗算の計算複雑性を標準のO(n³)の境界を下回るようにアルゴリズムを開発すること。
• コンパクトな表現（ビットパッキング、ビットスライcing、クロネッカー代入）と再帰的ブロックアルゴリズムを用いて、高性能実装を可能にすること。
• Ple分解を用いて、標準的な線形代数演算（LU分解、逆行列計算、上三角行列の解法）を任意のランクプロファイルに対応可能に一般化すること。
• 提案されたアルゴリズムを、暗号理論や符号理論における実用的応用を想定した、既存の高性能コンピューティングシステムおよびライブラリに統合すること。

提案手法

• F₂および小さな拡張体に対して、1つのマシンワードに複数の体要素を格納するビットパッキングを用い、64ビットワード上で並列演算を可能にする。
• F₃、F₅、F₇、およびF_{2³}、F_{3²}の要素を、複数のビットを同時にブール演算で表現することでビットスライcingを適用し、算術コストを最小限に抑える。
• 拡張体の多項式要素を整数として表現するためのクロネッカー代入を用い、REDQアルゴリズムによる効率的な還元を可能にする。
• 行列乗算に帰着するように設計された再帰的ブロックアルゴリズム（例：PLE、trsm、trmm、trtri）を用い、最高の既知の行列乗算アルゴリズムと同等の複雑性上限を達成する。
• F₂に対しては、四人の男の方法（Greasing）を用い、事前計算で2ᵏ個の行の線形結合を生成することで、O(n³/log n)の複雑性を達成する。
• 適合分割を用いた再帰的ブロック分解を活用し、上三角行列および行列乗算演算を小さな再帰呼び出しとGEMMサブルーチンに表現する。

実験結果

リサーチクエスチョン

• RQ1 事前計算とビットレベル並列処理を用いて、F₂上での行列乗算をO(n³)を超える速度で高速化できるか？
• RQ2 コンパクトな表現（ビットパッキング、ビットスライcing）は、密行列線形代数における有限体算術のコストをどのように低減できるか？
• RQ3 高度なデータ表現および再帰的アルゴリズムを用いることで、小さな有限体上での行列乗算の達成可能な複雑性はどの程度か？
• RQ4 任意のランクプロファイルに対して、PLE分解を行列乗算の時間で計算可能か？また、PLU分解への一般化はどのように実現されるか？
• RQ5 再帰的ブロックアルゴリズムを用いた上三角行列の解法および逆行列計算は、数値的安定性を保ちつつ、有限体上での体演算回数をどの程度削減できるか？

主な発見

• 四人の男の方法（Greasing）を用いることで、F₂上での行列乗算はO(n³/log n)の複雑性を達成し、グレイコード順序によりビット演算回数をk·2ᵏ·nから2ᵏ·nに削減した。
• ビットスライcingにより、F₃、F₅、F₇ではそれぞれ6、20、17回のブール演算で加算が効率的に行える。
• REDQアルゴリズムにより、64ビット整数算術を用いて素数モジュロの係数還元を効率的に行い、体要素のコンパクトな表現を可能にした。
• PLE分解アルゴリズムにより、ランクを明らかにするガウスの消去法の複雑性が行列乗算の複雑性にまで低下し、M(n)時間で実行可能となった。ここでM(n)は2つのn×n行列の乗算に要する時間である。
• ブロック分解とGEMMサブルーチンを用いて、再帰的上三角行列の解法（trsm）、行列乗算（trmm）、逆行列計算（trtri）アルゴリズムを実装し、最適な漸近的複雑性を達成した。
• 提案されたアルゴリズムは、M4RI、Magma、Sage、NTL、FLINT、IMLなどの高性能ライブラリに統合され、計算代数および暗号理論分野での実用的応用を可能にした。