[論文レビュー] Computing the Mertens function on a GPU

Eugene Kuznetsov|arXiv (Cornell University)|Jul 31, 2011

Neural Networks and Applications参考文献 2被引用数 7

ひとこと要約

この論文は、再帰的総和恒等式と並列配列更新を活用して、GPUを用いた高速化アルゴリズムを提示し、$ O(x^{2/3+\theta}) $ 時間でメルテンス関数 $ M(x) $ を計算する。$ x \approx 1.161 \times 10^{19} $ において、$ M(x)/x^{1/2} = -0.585768 $ の新たな極値を報告し、$ \exp(10^{15}) $ まで探索を拡張したが、メルテンス予想の反例は見つからなかった。

ABSTRACT

A GPU implementation of an algorithm to compute the Mertens function in O(x2/3+{\ko}) time is discussed. Results for x up to $10^{22}$, and a new extreme value for $M(x)/x^{1/2}$, -0.585768 ($M(x) \approx -1.996 \ast 10^9$ at $x \approx 1.161 \ast 10^{19}$), are reported.An approximate algorithm is used to examine values of M(x) for x up to $\exp{(10^{15})}$.

研究の動機と目的

• GPU並列処理を活用して、メビウス関数 $ \mu(n) $ による $ M(x) = \sum_{n=1}^x \mu(n) $ のメルテンス関数の計算を高速化すること。
• 特にメルテンス予想の検証を目的として、これまでの限界を超えて $ x $ の計算範囲を拡大すること。
• 総和恒等式に内在する再帰的構造を活用し、複数の値を同時に計算できるスケーラブルなアルゴリズムの開発。
• $ x \gg \exp(10^{10}) $ の場合に必要な高精度三角関数評価における数値精度の問題を解決すること。

提案手法

• 項を $ \lfloor n/x \rfloor $ の値ごとにグループ化することで、異なる項の数を $ O(\sqrt{n}) $ に削減できる恒等式 $ M(n) = 1 - \sum_{x=2}^n M(\lfloor n/x \rfloor) $ を使用。
• 二段階アプローチを採用：(1) $ y \leq u $ に対してCPUベースのエラトステネスのふるいにより $ M(y) $ を計算し、(2) $ u = O(n^{2/3+\theta}) $ として、$ x > 1 $ に対してGPU並列更新により $ M(\lfloor n/x \rfloor) $ を計算。
• 部分的素因数積を8ビット整数に格納する対数近似を適用し、ふるい処理中のメモリ帯域とキャッシュ効率を向上。
• ビット操作により平方因子を持つ数（MSBがセット）をフラグし、蓄積された対数素因数積の符号と大きさから $ \mu(y) $ を計算。
• 大きな $ n $ に対してブロック単位の処理を実装し、各ブロックごとに高精度定数テーブルを更新することで、大規模な $ x $ における余弦評価の精度を維持。
• 各ブロックに対して、$ b_i' = (b_i + z_i x_0) \mod 2\pi $ のようにずらした参照テーブルを用いて、近似解析に用いる三角関数項の精度を保持。

実験結果

リサーチクエスチョン

• RQ1 GPU並列処理を活用して、スケールの大きなメルテンス関数を効率的に計算できるか？
• RQ2 $ |M(x)/x^{1/2}| $ の最大値は何か？また、メルテンス予想で提唱された境界を超えるか？
• RQ3 $ x \approx \exp(10^{15}) $ までスケーリング可能であり、かつ数値精度を維持できるか？
• RQ4 計算範囲内でメルテンス予想 $ |M(n)| < n^{1/2} $ の反例は存在するか？

主な発見

• アルゴリズムは $ O(x^{2/3+\theta}) $ 時間計算量を達成し、$ x $ が $ 10^{22} $ に達する範囲まで $ M(x) $ を正しく計算できた。
• $ x \approx 1.161 \times 10^{19} $ において、$ M(x)/x^{1/2} = -0.585768 $ の新たな極値が発見され、$ M(x) \approx -1.996 \times 10^9 $ であった。
• $ \exp(10^{15}) $ までの範囲で、メルテンス予想の反例は一切発見されず、前回の探索範囲を2桁分拡大した。
• $ M(x)/x^{1/2} $ の代理指標である $ |q_{10^6}(x)| $ の最大観測値は、$ x \approx \exp(5.0586 \times 10^{14}) $ で $ 0.95652 $ であった。
• $ q_{10^6}(x) $ の極値の分布は対称的であり、大きな正または負の値に符号の非対称性の兆候は認められなかった。
• 高精度のブロック単位の定数更新により、$ \exp(10^{15}) $ まで正確な余弦評価が可能となり、GPUの単精度浮動小数点精度の損失を緩和した。