English

An efficient and fast sparse grid algorithm for high-dimensional numerical integration

Numerical Analysis 2022-10-27 v1 Numerical Analysis

Abstract

This paper is concerned with developing an efficient numerical algorithm for fast implementation of the sparse grid method for computing the dd-dimensional integral of a given function. The new algorithm, called the MDI-SG ({\em multilevel dimension iteration sparse grid}) method, implements the sparse grid method based on a dimension iteration/reduction procedure, it does not need to store the integration points, neither does it compute the function values independently at each integration point, instead, it re-uses the computation for function evaluations as much as possible by performing the function evaluations at all integration points in a cluster and iteratively along coordinate directions. It is shown numerically that the computational complexity (in terms of CPU time) of the proposed MDI-SG method is of polynomial order O(Nd3)O(Nd^3 ) or better, compared to the exponential order O(N(logN)d1)O(N(\log N)^{d-1}) for the standard sparse grid method, where NN denotes the maximum number of integration points in each coordinate direction. As a result, the proposed MDI-SG method effectively circumvents the curse of dimensionality suffered by the standard sparse grid method for high-dimensional numerical integration.

Keywords

Cite

@article{arxiv.2210.14313,
  title  = {An efficient and fast sparse grid algorithm for high-dimensional numerical integration},
  author = {Huicong Zhong and Xiaobing Feng},
  journal= {arXiv preprint arXiv:2210.14313},
  year   = {2022}
}

Comments

28 pages, 12 tables, 8 figures. arXiv admin note: text overlap with arXiv:2210.13658