Dolgov, Sergey and Savostyanov, Dmitry (2020) Parallel cross interpolation for high-precision calculation of high-dimensional integrals. Computer Physics Communications, 246. p. 106869. DOI https://doi.org/10.1016/j.cpc.2019.106869
Dolgov, Sergey and Savostyanov, Dmitry (2020) Parallel cross interpolation for high-precision calculation of high-dimensional integrals. Computer Physics Communications, 246. p. 106869. DOI https://doi.org/10.1016/j.cpc.2019.106869
Dolgov, Sergey and Savostyanov, Dmitry (2020) Parallel cross interpolation for high-precision calculation of high-dimensional integrals. Computer Physics Communications, 246. p. 106869. DOI https://doi.org/10.1016/j.cpc.2019.106869
Abstract
We propose a parallel version of the cross interpolation algorithm and apply it to calculate high-dimensional integrals motivated by Ising model in quantum physics. In contrast to mainstream approaches, such as Monte Carlo and quasi Monte Carlo, the samples calculated by our algorithm are neither random nor form a regular lattice. Instead we calculate the given function along individual dimensions (modes) and use these values to reconstruct its behaviour in the whole domain. The positions of the calculated univariate fibres are chosen adaptively for the given function. The required evaluations can be executed in parallel along each mode (variable) and over all modes. To demonstrate the efficiency of the proposed method, we apply it to compute high-dimensional Ising susceptibility integrals, arising from asymptotic expansions for the spontaneous magnetisation in two-dimensional Ising model of ferromagnetism. We observe strong superlinear convergence of the proposed method, while the MC and qMC algorithms converge sublinearly. Using multiple precision arithmetic, we also observe exponential convergence of the proposed algorithm. Combining high-order convergence, almost perfect scalability up to hundreds of processes, and the same flexibility as MC and qMC, the proposed algorithm can be a new method of choice for problems involving high-dimensional integration, e.g. in statistics, probability, and quantum physics.
Item Type: | Article |
---|---|
Uncontrolled Keywords: | High-dimensional integration; High precision; Tensor train format; Cross interpolation; Ising integrals; Parallel algorithms |
Divisions: | Faculty of Science and Health Faculty of Science and Health > Mathematics, Statistics and Actuarial Science, School of |
SWORD Depositor: | Unnamed user with email elements@essex.ac.uk |
Depositing User: | Unnamed user with email elements@essex.ac.uk |
Date Deposited: | 07 Feb 2020 12:59 |
Last Modified: | 07 Aug 2024 19:58 |
URI: | http://repository.essex.ac.uk/id/eprint/26652 |
Available files
Filename: 1-s2.0-S0010465519302565-main.pdf
Licence: Creative Commons: Attribution 3.0