Abstract
This paper explores the problem of finding a rank-structured approximation to a given Hermitian positive definite matrix. We first study a fundamental matrix nearness problem which seeks an approximation whose off-diagonal block must have low rank. Our investigation focuses on three objective functions: the spectral norm, Frobenius divergence, and log-determinant divergence. We derive closed-form solutions for the fundamental nearness problem, extending our analysis to include problems with a range restriction on the off-diagonal block. Building upon these results, we propose a practical, scalable algorithm to compute a positive definite hierarchical off-diagonal low-rank approximation to a Hermitian positive definite matrix. Through numerical experiments, we explore the use of our approximation as a general-purpose preconditioner for the conjugate gradient method. Our results demonstrate that the log-determinant divergence-based approximation exhibits good performance across a diverse set of test problems.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.