[論文レビュー] A reduced fast component-by-component construction of lattice point sets with small weighted star discrepancy

Ralph Kritzinger, Helene Laimer|arXiv (Cornell University)|Jan 28, 2015

Mathematical Approximation and Integration参考文献 16被引用数 3

ひとこと要約

本稿では、重みが急速に減少する場合に、重み付きスター分散が小さいラティス点集合を構築するための高速化された削減版成分別最適化（CBC）アルゴリズムを提示する。重みの構造的生成ベクトルと成分選択の最適化により、重みとシフトパラメータの和分条件のもとで強多項式トレーキタビリティを達成し、高次元における準モンテカルロ積分を効率的に行える。

ABSTRACT

The weighted star discrepancy of point sets appears in the weighted Koksma-Hlawka inequality and thus is a measure for the quality of point sets with respect to their performance in quasi-Monte Carlo algorithms. A special choice of point sets are lattice point sets whose generating vector can be obtained one component at a time such that the resulting lattice point set has a small weighted star discrepancy. In this paper we consider a reduced fast component-by-component algorithm which significantly reduces the construction cost for such generating vectors provided that the weights decrease fast enough.

研究の動機と目的

• 準モンテカルロ積分のため、重み付きスター分散が小さいラティス点集合の構築コストを低減すること。
• 構造的生成ベクトルを活用することで、次元とサンプルサイズに応じてスケーラブルに動作する成分別最適化アルゴリズムの開発。
• 得られるラティスルールが高次元積分問題において強多項式トレーキタビリティを達成する条件の特定。
• 既存のCBC構成を改善し、操作回数を削減しながらも最適な分散バウンドを維持すること。

提案手法

• 本手法は、重み列とシフトパラメータに基づいて関連する成分にのみ焦点を当てた、成分別に生成ベクトルを構築する削減版高速CBCアルゴリズムを用いる。
• 基本$b$の桁とシフトパラメータ$w_j$を用いた生成ベクトルの変換を導入し、重み付きスター分散の計算をより効率的に行う。
• 識別子$ R_{N,\boldsymbol{\nu}}^s(\textbf{z}) = \frac{1}{N} \sum_{k=0}^{N-1} \eta_s(k) $を活用して、分散関連項を効率的に計算する。
• 調和級数のバウンディング$ S_N \leq 4\log N $を適用し、分散推定における誤差項を制御する。
• 構築コストは標準CBC法の$ O(sN\log N) $に対して、$ O(N\log N + \min\{s,t\}N + N\sum_{d=1}^{\min\{s,t\}} (m - w_d)b^{-w_d}) $にまで低減される。
• 積型重み$ \gamma_{\mathfrak{u}} = \prod_{j\in\mathfrak{u}} \gamma_j $を用い、$ \gamma_j \leq 1 $、非増加であり、$ \sum_j \gamma_j b^{w_j} < \infty $を仮定する。

実験結果

リサーチクエスチョン

• RQ1 重み付きスター分散が小さいラティス点集合の構築コストを、分散の低さを保ちながら低減できるか。
• RQ2 重みとシフトパラメータにどのような条件下で、削減版CBCアルゴリズムが強多項式トレーキタビリティを達成するか。
• RQ3 計算複雑性の観点から、削減版CBCアルゴリズムは標準CBCアルゴリズムに比べてどの程度効率的か。
• RQ4 シフトパラメータ$ w_j $は、分散と構築コストにどのような影響を及ぼすか。
• RQ5 変換された生成ベクトルと調和級数の推定値を用いて、分散バウンディングをきめ細かく制御できるか。

主な発見

• 削減版高速CBCアルゴリズムにより、構築コストが$ O(N\log N + \min\{s,t\}N + N\sum_{d=1}^{\min\{s,t\}} (m - w_d)b^{-w_d}) $にまで低減され、標準CBC法の$ O(sN\log N) $に比べて著しく安価になる。
• 条件$ \sum_{j=1}^\infty \gamma_j b^{w_j} < \infty $のもとで、強多項式トレーキタビリティが達成され、任意の$ \delta > 0 $に対して$ R_{N,\boldsymbol{\gamma}}^s(\mathbf{z}) \leq c_{\boldsymbol{\gamma},\delta} N^{\delta - 1} $が成り立つ。
• 重み付きスター分散は$ D_{N,\boldsymbol{\gamma}}^*(\mathbf{z}) \leq O(1/N) + \frac{1}{2} R_{N,\boldsymbol{\gamma}}^s(\mathbf{z}) $を満たし、和分条件のもとで$ 1/N $に比例する収束を保証する。
• 条件$ N \geq (c_{\boldsymbol{\gamma},\delta} \varepsilon^{-1})^{1/(1 - \delta)} $を満たすとき、$ R_{N,\boldsymbol{\gamma}}^s(\mathbf{z}) \leq \varepsilon $が保証され、強多項式トレーキタビリティが確認される。
• バウンディング$ \prod_{j=1}^s (1 + 13\gamma_j b^{w_j} \log N) \leq (1 + \sigma_d^{-1})^d N^{(\sigma_0 + 1)\sigma_d} $により、和分可能な重みのもとでトレーキタビリティ結果の導出が可能となる。
• 理論的枠組みにより、削減版CBCアルゴリズムが計算負荷を著しく低減しながらも、最適な分散バウンディングを維持することが確認された。特に重みが急速に減少する場合には顕著な効果を示す。