Sum of Kronecker products representation and its Cholesky factorization for spatial covariance matrices from large grids

Jian Cao, Marc G. Genton, David E. Keyes, George M. Turkiyyah

Research output: Contribution to journalArticlepeer-review

Abstract

The sum of Kronecker products (SKP) representation for spatial covariance matrices from gridded observations and a corresponding adaptive-cross-approximation-based framework for building the Kronecker factors are investigated. The time cost for constructing an -dimensional covariance matrix is and the total memory footprint is , where is the number of Kronecker factors. The memory footprint under the SKP representation is compared with that under the hierarchical representation and found to be one order of magnitude smaller. A Cholesky factorization algorithm under the SKP representation is proposed and shown to factorize a one-million dimensional covariance matrix in under 600 seconds on a standard scientific workstation. With the computed Cholesky factor, simulations of Gaussian random fields in one million dimensions can be achieved at a low cost for a wide range of spatial covariance functions.
Original languageEnglish (US)
Pages (from-to)107165
JournalComputational Statistics & Data Analysis
DOIs
StatePublished - Jan 2021

Fingerprint

Dive into the research topics of 'Sum of Kronecker products representation and its Cholesky factorization for spatial covariance matrices from large grids'. Together they form a unique fingerprint.

Cite this