[論文レビュー] H-Matrix and Block Error Tolerances

Andrew Bradley|arXiv (Cornell University)|Oct 12, 2011

Electromagnetic Scattering and Analysis参考文献 5被引用数 11

ひとこと要約

本稿では、階層行列（H行列）の近似を向上させるために、新たな手法である修正相対誤差マッピング（MREM）を提案する。MREMは、全誤差許容誤差を各ブロック固有の許容誤差にマッピングする方法により、過剰な精度を回避する。特異性の順序が1より大きいカーネルに対して、標準的なブロック相対誤差マッピング（BREM）と比較して1.5〜5倍の高い圧縮率を達成し、計算コストは低く抑えられ、誤差制御もきめ細やかになる。特に、対角優勢な行列において顕著である。

ABSTRACT

We describe a new method to map the requested error tolerance on an H-matrix approximation to the block error tolerances. Numerical experiments show that the method produces more efficient approximations than the standard method for kernels having singularity order greater than one, often by factors of 1.5 to 5 and at a lower computational cost.

研究の動機と目的

• 標準的なH行列構築法ではブロック許容誤差が一様に設定されるため、過剰な精度が生じ、圧縮率が低下するという非効率性を是正すること。
• 全誤差許容誤差を、必要な精度よりもわずかに高い近似が得られるように、各ブロックの許容誤差にマッピングする手法の開発。
• 不要な精度を低減させつつ、全誤差境界を満たすことを保証する計算効率の向上。
• さまざまな特異性の順序と幾何形状を有する問題に対して、特に高めの対角優勢性を示すものに対して、手法の評価。
• 特異性の順序が1より大きいカーネルにおいて、MREMが標準的手法を一貫して圧縮率と効率の点で上回ることの実証。

提案手法

• 全誤差フロベニウスノルムへの各ブロックの寄与度に応じて許容誤差を割り当てるための導出不等式（式7）に基づく、新たなブロック許容誤差マッピング戦略を提案。
• 各ブロックのノルムに比例するように、ブロック許容誤差を正規化要因として全行列のフロベニウスノルム $\|B\|_{\text{F}}$ を用いる。
• 特異値分解（SVD）または類似手法による各ブロックの低ランク近似（LRA）を、導出されたブロック許容誤差で制約し、全誤差境界を満たすようにする。
• 全行列の $\|B\|_{\text{F}}$ を、総構築時間の約1%のコストで効率的に推定することで、追加コストを最小限に抑えつつ手法を実装可能にする。
• さまざまな問題サイズと幾何形状において、MREMと標準的な BREM 法を比較し、達成誤差、圧縮率、改善係数を測定。
• 行列-ベクトル積（MVP）の誤差境界を用いて、ブロック許容誤差と全精度の関係を定式化し、$\|Bx - \bar{B}x\|_2 / \|B\|_{\text{F}}\|x\|_2 \leq \varepsilon$ を保証する。

実験結果

リサーチクエスチョン

• RQ1 より効率的なブロック許容誤差マッピング戦略により、H行列近似における過剰な精度を低減させつつ、所定の全誤差許容誤差を維持できるか？
• RQ2 さまざまなカーネルタイプと幾何形状において、MREMと標準的な BREM 法との間で、圧縮率および計算コストの改善はどの程度か？
• RQ3 MREMの性能向上は、カーネルの特異性の順序に依存するか。特に、順序が1より大きい場合に顕著か？
• RQ4 行列の対角優勢性が、BREM と比較して MREM の有効性にどの程度影響を与えるか？
• RQ5 特に $\|B\|_{\text{F}}$ の推定において、MREM は最小限の計算コストで実装可能か？

主な発見

• 特異性の順序が1より大きいカーネルに対して、MREM は BREM と比較して1.5〜5倍の高い圧縮率を達成し、近似の効率が顕著に向上する。
• 156 × 402 パッチを有する地球物理学的断層すべり問題において、MREM は圧縮率75.74を達成したのに対し、BREM は16.89にとどまり、改善係数は4.48に達した。
• MREM における達成誤差は $1.71 \times 10^{-6}$ であり、要求された $10^{-5}$ 許容誤差に近く、BREM は $1.28 \times 10^{-8}$ を達成しており、顕著な過剰精度が確認された。
• MREM は BREM よりも劣ることがほとんどなく、劣る場合でも差は小さく、特に特異性の順序が1のカーネルでは顕著である。
• 改善係数は、行列の対角優勢性および特異性の順序が高くなるに従い増加し、高次の特異性に対して最も顕著な向上が得られた。
• 全行列の $\|B\|_{\text{F}}$ の推定コストは、全体の構築時間の約1%にとどまり、MREM が実世界の応用に実用的であることを示している。