[论文解读] 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_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)。
• 利用区域分区结构（如棋盘式稀疏性）以利用对称性，降低格点模拟中的计算成本。
• 验证对于双分区域结构，当 p 为奇数时，trace((M_D^(-1)M_off)^p) = 0，从而将级数简化为仅含偶数阶项。

实验结果

研究问题

• 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 个非零元，内存使用量减少逾三分之二。
• 对数的虚部收敛速度比实部快，表明各分量收敛行为不同。
• 误差界中的常数 c ≈ 554 过于保守，因为 M_D^(-1)M_off 的许多特征值模长远小于 ρ。