Eftekhari, Aryan Pasadakis, Dimosthenis Bollhöfer, Matthias Scheidegger, Simon Schenk, Olaf 2021-05-12 The ℓ1-regularized Gaussian maximum likelihood method is a common approach for sparse precision matrix estimation, but one that poses a computational challenge for high-dimensional datasets. We present a novel ℓ1- regularized maximum likelihood method for performant large-scale sparse precision matrix estimation utilizing the block structures in the underlying computations. We identify the computational bottlenecks and contribute a block coordinate descent update as well as a block approximate matrix inversion routine, which is then parallelized using a shared-memory scheme. We demonstrate the effectiveness, accuracy, and performance of these algorithms. Our numerical examples and comparative results with various modern open-source packages reveal that these precision matrix estimation methods can accelerate the computation of covariance matrices by two to three orders of magnitude, while keeping memory requirements modest. Furthermore, we conduct large-scale case studies for applications from finance and medicine with several thousand random variables to demonstrate applicability for real-world datasets. application/pdf https://susi.usi.ch/global/documents/319335 https://n2t.net/ark:/12658/srd1319335 https://susi.usi.ch/documents/319335/files/Schenk_jocs_2021.pdf eng info:eu-repo/semantics/altIdentifier/doi/10.1016/j.jocs.2021.101389 info:eu-repo/semantics/altIdentifier/ark/12658/srd1319335 info:eu-repo/semantics/openAccess CC BY-NC Journal of computational science. - Elsevier. - 2021, vol. 53, p. 13 Covariance matrices Graphical model Optimization Gaussian Markov random field Machine learning application info:eu-repo/classification/udc/004 Block-enhanced precision matrix estimation for large-scale datasets http://purl.org/coar/resource_type/c_6501