[論文レビュー] Determinant Approximations

Ilse C. F. Ipsen, Dean Lee|arXiv (Cornell University)|May 2, 2011

Matrix Theory and Algorithms参考文献 15被引用数 31

ひとこと要約

本稿では、行列式のトレース-ログ展開 det(M) = exp(trace(log(M))) を用いて、複雑で非エルミート行列に対する行列式の逐次近似を提示する。誤差境界は M_D^(-1)M_off の固有値半径に基づく。この手法は1〜3反復で2〜3桁の有効数字の精度に達し、ガウスの消去法と比較してメモリ使用量を顕著に削減する。また、ブロック対角近似を用いて、フィッシャーの不等式およびハダマールの不等式を非エルミートの場合に拡張する。

ABSTRACT

A sequence of approximations for the determinant and its logarithm of a complex matrixis derived, along with relative error bounds. The determinant approximations are derived from expansions of det(X)=exp(trace(log(X))), and they apply to non-Hermitian matrices. Examples illustrate that these determinant approximations are efficient for lattice simulations of finite temperature nuclear matter, and that they use significantly less space than Gaussian elimination. The first approximation in the sequence is a block diagonal approximation; it represents an extension of Fischer's and Hadamard's inequalities to non-Hermitian matrices. In the special case of Hermitian positive-definite matrices, block diagonal approximations can be competitive with sparse inverse approximations. At last, a different representation of sparse inverse approximations is given and it is shown that their accuracy increases as more matrix elements are included.

研究の動機と目的

• 有限温度の核物質格子シミュレーションに現れる非エルミート的でスパースな行列に対する、効率的で低メモリの行列式近似法の開発。
• ブロック対角近似を用いて、古典的行列式不等式（フィッシャーの不等式およびハダマールの不等式）を非エルミート行列に拡張すること。
• 行列式対数の級数展開に基づく、次第に高精度になる行列式近似の相対誤差境界の提供。
• ブロック対角近似が、非エルミート行列においてスパース逆行列法を上回る精度と効率を達成できることの実証。
• 本手法が、完全ピボット選択を伴うガウスの消去法と比較して、顕著に少ないメモリ使用量で済むことの示唆。

提案手法

• 行列 M を M = M_D + M_off に分解する。ここで M_D はブロック対角行列、M_off はブロック対角成分以外の要素からなる。
• 恒等式 det(M) = det(M_D) * exp(trace(log(I + M_D^(-1)M_off))) を用いて、行列式の対数の級数展開を導出する。
• テイラー展開 log(I + X) = X - X^2/2 + X^3/3 - ... を適用し、log(I + M_D^(-1)M_off) を近似する。誤差境界は M_D^(-1)M_off の固有値半径 ρ に依存する。
• 行列式の対数の近似として δ_j = trace(∑_{i=1}^j (-1)^{i+1} (M_D^(-1)M_off)^i / i) を計算し、行列式自体の近似として Δ_j = exp(δ_j) を得る。
• チェス盤型スパーシティなどのゾーン分割構造を活用し、格子シミュレーションにおける計算コストを低減する。
• 二部ゾーン構造の場合、trace((M_D^(-1)M_off)^p) = 0（p が奇数のとき）であることを確認し、級数を偶数次項のみに簡略化できる。

実験結果

リサーチクエスチョン

• RQ1 ブロック対角近似は、非エルミート行列に対するフィッシャーの不等式およびハダマールの不等式を拡張できるか？
• RQ2 行列式近似の誤差境界は、M_D^(-1)M_off の固有値半径にどのように依存するか？
• RQ3 格子シミュレーションにおいて、低精度の行列式近似はガウスの消去法よりも速く、かつより少ないメモリで計算可能か？
• RQ4 非エルミート行列において、ブロック対角近似の精度はスパース逆行列近似と比較してどうか？
• RQ5 級数近似の収束性は、M_D^(-1)M_off の固有値分布に依存するか？

主な発見

• 核物質シミュレーションの例では、ブロック対角近似 δ_0 = ln(det(M_D)) が、行列式の対数に対して2桁の精度を達成する。
• 2反復目（δ_2）では、対数の近似が3桁の精度に達し、実部の絶対誤差 ≈ 0.48、虚部の絶対誤差 ≈ 0.0025 となる。
• δ_j の絶対誤差は ρ^j にほぼ比例して減少し、ρ ≈ 0.6613 が M_D^(-1)M_off の固有値半径であるため、幾何的収束が示唆される。
• 本手法は最初の3つの近似で162n個の非ゼロ要素を要するガウスの消去法と比較して、非ゼロ要素を49n個にまで削減し、メモリ使用量を三分の二以上削減する。
• 虚部の対数は実部よりも速く収束しており、成分ごとの収束挙動の違いが示唆される。
• 誤差境界に現れる定数 c ≈ 554 は、M_D^(-1)M_off の多くの固有値が ρ よりもはるかに小さい大きさであるため、過剰に悲観的である。