We study a preconditioner for a Hermitian positive definite linear system, which is obtained as the solution of a matrix nearness problem based on the Bregman \emph{log determinant} divergence. The preconditioner is on the form of a Hermitian positive definite matrix plus a low-rank matrix. For this choice of structure, the generalised eigenvalues of the preconditioned system are easily calculated, and we show that the preconditioner is optimal in the sense that it minimises the $\ell_2$ condition number of the preconditioned matrix. We develop practical numerical approximations of the preconditioner based on the randomised singular value decomposition (SVD) and the Nystr\"om approximation and provide corresponding approximation results. Furthermore, we prove that the Nystr\"om approximation is in fact also a matrix approximation in a range-restricted Bregman divergence and establish several connections between this divergence and matrix nearness problems in different measures. Numerical examples are provided to support the theoretical results.
翻译:我们研究了一种用于埃尔米特正定线性系统的预处理器,该预处理器通过基于布雷格曼对数行列式散度的矩阵逼近问题求解得到。该预处理器采用埃尔米特正定矩阵加低秩矩阵的结构形式。针对这种结构选择,预处理系统的广义特征值可被简便计算,并且我们证明了该预处理器在最小化预处理矩阵的ℓ2条件数意义上是最优的。我们基于随机奇异值分解(SVD)和Nyström近似开发了预处理器的实用数值近似方法,并给出了相应的近似结果。此外,我们证明Nyström近似实际上也是范围受限布雷格曼散度意义下的矩阵近似,并建立了该散度与不同测度下矩阵逼近问题之间的若干联系。最后通过数值算例验证了理论结果。