123

Kriging in Tensor Train data format

Abstract

Combination of low-tensor rank techniques and the Fast Fourier transform (FFT) based methods had turned out to be prominent in accelerating various statistical operations such as Kriging, computing conditional covariance, geostatistical optimal design, and others. However, the approximation of a full tensor by its low-rank format can be computationally formidable. In this work, we incorporate the robust Tensor Train (TT) approximation of covariance matrices and the efficient TT-Cross algorithm into the FFT-based Kriging. It is shown that here the computational complexity of Kriging is reduced to O(dr3n)\mathcal{O}(d r^3 n), where nn is the mode size of the estimation grid, dd is the number of variables (the dimension), and rr is the rank of the TT approximation of the covariance matrix. For many popular covariance functions the TT rank rr remains stable for increasing nn and dd. The advantages of this approach against those using plain FFT are demonstrated in synthetic and real data examples.

View on arXiv
Comments on this paper