Abstract:We study the approximation of integrals $\int_D f(\boldsymbol{x}^\top A) \mathrm{d} \mu(\boldsymbol{x})$, where $A$ is a matrix, by quasi-Monte Carlo (QMC) rules $N^{-1} \sum_{k=0}^{N-1} f(\boldsymbol{x}_k^\top A)$. We are interested in cases where the main cost arises from calculating the products $\boldsymbol{x}_k^\top A$. We design QMC rules for which the computation of $\boldsymbol{x}_k^\top A$, $k = 0, 1, \ldots, N-1$, can be done fast, and for which the error of the QMC rule is similar to the standard QMC error. We do not require that $A$ has any particular structure.
For instance, this approach can be used when approximating the expected value of a function with a multivariate normal random variable with a given covariance matrix, or when approximating the expected value of the solution of a PDE with random coefficients.
The speed-up of the computation time is sometimes better and sometimes worse than the fast QMC matrix-vector product from [Dick, Kuo, Le Gia, and Schwab, Fast QMC Matrix-Vector Multiplication, SIAM J. Sci. Comput. 37 (2015)]. As in that paper, our approach applies to (polynomial) lattice point sets, but also to digital nets (we are currently not aware of any approach which allows one to apply the fast method from the aforementioned paper of Dick, Kuo, Le Gia, and Schwab to digital nets).
Our method does not use FFT, instead we use repeated values in the quadrature points to derive a reduction in the computation time. This arises from the reduced CBC construction of lattice rules and polynomial lattice rules. The reduced CBC construction has been shown to reduce the computation time for the CBC construction. Here we show that it can also be used to also reduce the computation time of the QMC rule.
What problem does this paper attempt to address?
### What problem does this paper attempt to solve?
This paper aims to solve the problem of efficient approximate calculation of high - dimensional integrals. In particular, when using the Quasi - Monte Carlo (QMC) method, it focuses on how to reduce the computational cost of the matrix - vector multiplication $\mathbf{x}_k^\top A$. Specifically, the paper is concerned with the approximate calculation of integrals in the following form:
\[
\int_D f(\mathbf{x}^\top A) \, d\mu(\mathbf{x}),
\]
where $A$ is a matrix, $f$ is a function, $D\subseteq \mathbb{R}^s$, and $\mu$ is a measure. The QMC method approximates the above integral by a rule in the following form:
\[
Q_N(f)=\frac{1}{N}\sum_{k = 0}^{N-1}f(\mathbf{x}_k^\top A),
\]
where $\mathbf{x}_0,\mathbf{x}_1,\ldots,\mathbf{x}_{N - 1}\in D$ are deterministic quadrature points.
The main computational cost comes from calculating the matrix - vector multiplications $\mathbf{x}_k^\top A$ for $N$ points, which requires $O(Nst)$ operations, where $s$ is the dimension and $t$ is the number of columns of matrix $A$.
### The method proposed in the paper
To reduce the computational cost, the paper proposes a fast QMC matrix - vector multiplication algorithm, which takes advantage of certain repetitive structures of the point set. Specifically, the authors introduce the concept of "simplified CBC construction", which makes the calculation of $\mathbf{x}_k^\top A$ more efficient. This method is applicable not only to lattice rules and polynomial lattice rules, but also to digital nets, and can be used in combination with random shifts, which previous methods could not do.
### Application scenarios
This computational problem is very common in statistics and uncertainty quantification. For example, similar problems are encountered when calculating the expected value of a multivariate normal distribution or the expected value of the solution of a partial differential equation with random coefficients. In addition, in financial mathematics, for example, when calculating the returns of financial products depending on multiple assets, such problems are also encountered.
### Main contributions
1. **Improve computational efficiency**: By taking advantage of the repetitive structure of the point set, the time complexity of calculating $\mathbf{x}_k^\top A$ is significantly reduced.
2. **Wide range of applicability**: It is applicable not only to lattice rules and polynomial lattice rules, but also to digital nets, and can be used in combination with random shifts.
3. **Theoretical support**: A detailed algorithm description and computational complexity analysis are provided, proving the effectiveness of this method.
In conclusion, this paper proposes a new fast QMC matrix - vector multiplication algorithm, which can significantly reduce the computational cost while maintaining the accuracy of the QMC rule, and is applicable to a variety of high - dimensional integral problems.