[論文レビュー] Conquering the pre-computation in two-dimensional harmonic polynomial transforms

Richard Mikaël Slevinsky|arXiv (Cornell University)|Nov 21, 2017

Advanced Numerical Analysis Techniques被引用数 13

ひとこと要約

本稿では、二変数ジャコビ多項式構造を持つ球面やその他の幾何に適応可能な、超最適な前処理および記憶量複雑度を達成する2次元調和多項式変換のスケルトン化手法を提案する。この手法は、二進分割と対称正定値帯状一般固有値問題を活用することで、前処理および記憶量の複雑度を 𝒪(n^{3/2} log n) に低減する。このアプローチにより、数値的に安定した高速な変換が可能となり、実行時間は僅か 𝒪(log n) の増加にとどまる。

ABSTRACT

We describe a skeletonization of the spherical harmonic connection problem that reduces the storage and pre-computation to superoptimal complexities at the cost of increasing the execution time by the modest multiplicative factor of $\mathcal{O}(\log n)$. One advantage of accelerating the spherical harmonic connection problem over accelerating synthesis and analysis is that neighbouring layers (in steps of two) may be expanded in eachother's bases. The proposed skeletonization maximizes this interconnectivity by overlaying a dyadic partitioning on the connection problem. We derive the symmetric-definite banded generalized eigenvalue problem required to accelerate spherical harmonic transforms. We also include a full analysis of the weighted normalized Jacobi connection problem with applications to fast harmonic polynomial transforms on the disk, triangle, rectangle, deltoid, wedge, and any other geometry with a bivariate analogue of Jacobi polynomials.

研究の動機と目的

• 大規模な応用（例：数値的気象予測）において顕著に高まる球面および調和多項式変換における前処理および記憶量コストの低減。
• 漸近的に最適ではあるが実用的でない前処理を必要とする既存手法の限界を克服すること。
• 二変数ジャコビ多項式構造を活用することで、ディスク、三角形、長方形など多様な幾何に適用可能な統一的フレームワークの構築。
• 二進分割戦略を用いて隣接する調和多項式層間の相互接続性を最大限に活用し、数値的安定性とデータのスパarsityを向上させること。
• 構造化された固有値ソルバーを用いることで、前処理および記憶量の超最適な複雑度を達成しながら、後退安定性および高い相対精度を維持すること。

提案手法

• 球面調和関連問題にスケルトン化技術を適用し、対称正定値帯状一般固有値問題に還元する。
• 高速多重極法で加速された分割統治型固有値ソルバーを用い、最適な漸近的複雑度と実用的な効率性を達成する。
• 作用素の対称化およびコレスキー分解を用いて関連問題を再定式化し、固有問題における帯域と対称性を保持する。
• 重み付き正規化ジャコビ枠組みにおける乗法作用素の可換性を活用し、関連係数の構造的表現を維持する。
• 二進分割を用いて隣接層の展開間の作業をバランスさせ、関連行列におけるデータのスパarsityを向上させる。
• 特に度数が悪いコレスキー分解に対して、拡張精度算術を用いて高い相対精度を確保する。

実験結果

リサーチクエスチョン

• RQ1 調和多項式変換における前処理および記憶量の複雑度を、標準的な 𝒪(n² log n) の閾値未満に低下させつつ、数値的安定性を損なわずに実現できるか。
• RQ2 隣接する調和多項式層間の相互接続性をどの程度活用することで、関連問題の高速化が可能になるか。
• RQ3 二変数ジャコビ多項式設定において、関連係数を対称正定値帯状一般固有値問題の固有関数として表現できるか。
• RQ4 本手法は、数値的気象予測などの実世界の応用で関連する大規模帯域幅において、実際のスケーリング特性を示すか。
• RQ5 この文脈において、度数が悪い一般固有値問題を解く際、拡張精度算術が精度を保持するために果たす役割は何か。

主な発見

• 提案手法により、前処理および記憶量の複雑度が標準的な 𝒪(n² log n) より顕著に低い超最適な 𝒪(n^{3/2} log n) に低減された。
• 実行時間は僅か 𝒪(log n) の乗法的要因の増加にとどまり、大規模応用において実用的である。
• 対称正定値帯状一般固有値問題の定式化により、関連係数の表現における後退安定性および高い相対精度が達成された。
• 本手法は、ディスク、三角形、長方形、デルトイド、ワッペンなど、多様な幾何に、それらの二変数ジャコビ多項式類似物を介して適用可能である。
• 二進分割の活用により、隣接層間の相互接続性が最大化され、各層が互いの基底で相互に展開可能となり、数値的効率が向上した。
• 拡張精度算術は、特に度数が悪い計量作用素のコレスキー分解において高い相対精度を維持するために不可欠であり、前処理コストに定数因子しか追加しない。