[논문 리뷰] Determinant Approximations

Ilse C. F. Ipsen, Dean Lee|arXiv (Cornell University)|2011. 05. 02.

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_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)로 계산한다.
• 체스보드 스parser 구조와 같은 존 분할 구조를 활용하여 대칭성을 이용하고 계산 비용을 줄인다.
• 이중 존 구조의 경우, 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) 이후, 로그의 정확도는 3자리가 되며, 실수부의 절대 오차는 약 0.48, 허수부는 약 0.0025이다.
• δ_j의 절대 오차는 약 ρ^j로 감소하며, 여기서 ρ ≈ 0.6613은 M_D^(-1)M_off의 스펙트럼 반경이다. 이는 기하급수적 수렴을 시사한다.
• 이 방법은 첫 세 근사값에 대해 49n개의 비제로 요소만 필요로 하며, 완전 피봇팅을 사용한 가우스 소거법의 162n개 비제로 요소보다 메모리 사용을 2/3 이상 줄였다.
• 로그의 허수부는 실수부보다 더 빨리 수렴함을 확인하여, 성분 간 수렴 행동이 다름을 시사한다.
• 오차 한계에 나타나는 상수 c ≈ 554는 지나치게 낙관적이며, M_D^(-1)M_off의 많은 고유값이 ρ보다 훨씬 작기 때문이다.