[論文レビュー] Low discrepancy constructions in the triangle

Kinjal Basu, Art B. Owen|arXiv (Cornell University)|Mar 11, 2014

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

ひとこと要約

本稿では、単位三角形上の数値積分のための2つの新しい準モンテカルロ構成を提示し、不均一性が消えることを達成している。最初のものは、不均一性が $12/\sqrt{N}$ で有界な拡張可能デジタルvan der Corput型列であり、2番目は2次無理数の角度で回転された2次元整数格子を用い、最適な $O(\log N / N)$ の不均一性を達成する。両者とも、有界 variation を要件としないリーマン可積分関数に対しても収束を保証する。

ABSTRACT

Most quasi-Monte Carlo research focuses on sampling from the unit cube. Many problems, especially in computer graphics, are defined via quadrature over the unit triangle. Quasi-Monte Carlo methods for the triangle have been developed by Pillands and Cools (2005) and by Brandolini et al. (2013). This paper presents two QMC constructions in the triangle with a vanishing discrepancy. The first is a version of the van der Corput sequence customized to the unit triangle. It is an extensible digital construction that attains a discrepancy below 12/sqrt(N). The second construction rotates an integer lattice through an angle whose tangent is a quadratic irrational number. It attains a discrepancy of O(log(N)/N) which is the best possible rate. Previous work strongly indicated that such a discrepancy was possible, but no constructions were available. Scrambling the digits of the first construction improves its accuracy for integration of smooth functions. Both constructions also yield convergent estimates for integrands that are Riemann integrable on the triangle without requiring bounded variation.

研究の動機と目的

• 三角形における準モンテカルロ点集合を、不均一性が消えるように開発すること。これにより、リーマン可積分関数の収束が保証される。
• 単体（三角形）へのKoksma-Hlawka不等式フレームワークの拡張を行い、三角形積分における誤差評価の理論的基盤を提供すること。
• 明示的で決定的かつ不均一性が理論的に低い点集合を構築することにより、先行研究におけるこのような構成の欠如を補う。
• コンピュータグラフィックスや数値積分への応用を想定し、拡張可能性とスクラッチングのサポートを備えた実用的な構成を確保すること。

提案手法

• ドメインを再帰的な2進分割により部分三角形に分割することで、van der Corput列を三角形に適応する。
• 4進数の桁を三角形の部分領域にマッピングするデジタル列を構築し、拡張可能性を保証するとともに、デジタルスクラッチングを可能にする。
• 2次無理数の正接を持つ角度でスケーリングされた整数格子を回転し、三角形内にある点のみを保持する。
• ChenとTravaglini（2007）およびBrandoliniら（2013）の定理を用いて、回転格子が $O(\log N / N)$ の不均一性を達成することを証明する。
• PillardsとCoolsの不均一性とBrandoliniらの平行四辺形不均一性の等価性を確立し、両者とも $N \to \infty$ のとき消えることを示す。
• 符号付き不均一性関数と平行四辺形への幾何的分解を用いて、部分三角形の不均一性を評価する。

実験結果

リサーチクエスチョン

• RQ1 単位立方体への写像ではなく、直接的に三角形内に低不均一性列を構築することは可能か？
• RQ2 決定的点集合が三角形で達成可能な最適な不均一性レートは何か？
• RQ3 有界 variation を要件としないリーマン可積分関数に対しても、準モンテカルロ積分が収束するか？
• RQ4 三角形に基づくQMC列にデジタルスクラッチングを適用することで、滑らかな関数のための収束性が向上するか？
• RQ5 最適な $O(\log N / N)$ の不均一性を達成する、格子に類似した直接的な構成が三角形に可能か？

主な発見

• 三角形におけるデジタルvan der Corput構成は、平行四辺形不均一性が最大で $12/\sqrt{N}$ である。$N = 4^k$ のとき、正確な値は $2/(3\sqrt{N}) - 1/(9N)$ である。
• 回転格子構成は、最適な不均一性レート $O(\log N / N)$ を達成し、三角形における不均一性の最良レートに一致する。
• デジタル構成の桁をスクラッチングすることで、滑らかな関数の積分精度が向上し、$N = 4^k$ のとき、平均二乗誤差が $O(N^{-1})$ に達する。
• 両構成とも、有界 variation を要件としない任意のリーマン可積分関数について、QMC推定値が真の積分値に収束することを保証する。
• PillardsとCoolsの測度の不均一性は、Brandoliniらの平行四辺形不均一性の2倍以下であるため、両者とも $N \to \infty$ のとき消える。
• 理論的枠組みは、標準単位三角形に限らず、任意の非退化三角形に対しても収束を支持する。