[논문 리뷰] Low discrepancy constructions in the triangle

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

Mathematical Approximation and Integration참고 문헌 16인용 수 6

한 줄 요약

이 논문은 단위 삼각형에서의 수치 적분을 위한 두 가지 새로운 준몬테카를로 구성법을 제안하며, 이는 실수의 비율이 0으로 수렴한다. 첫 번째 구성법은 $12/\sqrt{N}$ 이하의 비율을 갖는 확장 가능한 디지털 반-더 코르푸트 유형 수열이며, 두 번째 구성법은 2차 무리수 각도로 2차원 정수 격자를 회전시켜 삼각형 내부에 있는 점들만 유지함으로써 최적의 $O(\log N / N)$ 비율을 달성한다. 두 구성법 모두 유계 변동성을 요구하지 않으며, 리만 적분 가능한 함수에 대해 수렴 가능성을 보장한다.

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.

연구 동기 및 목표

• 삼각형 내에서 실수의 비율이 0으로 수렴하는 준몬테카를로 점 집합을 개발함으로써, 리만 적분 가능한 함수에 대해 수렴 가능성을 보장한다.
• 단체(삼각형)로의 코크스마-플라우카 부등식 프레임워크를 확장하여 삼각형 적분의 오차 한계에 대한 이론적 기초를 제공한다.
• 실수의 비율이 낮다는 것이 증명된 명시적이고 결정적인 수열을 구성함으로써 이전 연구에서 부족했던 이러한 구성법의 결여를 보완한다.
• 컴퓨터 그래픽스 및 수치 적분 응용 분야에서 실용적이며, 확장 가능하고 스캐러빙 지원 기능을 갖춘 구성법을 확보한다.

제안 방법

• 영역을 하위 삼각형으로 재귀적으로 이진 분할하는 방식으로 반-더 코르푸트 수열을 삼각형에 적응시킨다.
• 기수 4의 자릿수를 삼각형 하위 영역에 매핑하여 디지털 수열을 구성함으로써 확장 가능성을 확보하고 디지털 스캐러빙을 가능하게 한다.
• 2차 무리수 탄젠트를 갖는 각도로 스케일된 정수 격자를 회전시키고, 삼각형 내부에 있는 점들만 유지한다.
• 첸과 트라바글리니(2007) 및 브란돌리니 등(2013)의 정리들을 활용하여, 회전된 격자가 $O(\log N / N)$ 비율을 달성함을 증명한다.
• 필라즈와 쿨스의 비율과 브란돌리니 등이 정의한 평행사변형 비율 간의 동치성을 확립함으로써, 두 비율 모두 $N \to \infty$일 때 0으로 수렴함을 보인다.
• 부호가 있는 비율 함수와 평행사변형으로의 기하학적 분해를 활용하여 하위 삼각형의 비율을 근사한다.

실험 결과

연구 질문

• RQ1 단위 초입방에서의 사상이 아닌, 삼각형 내부에서 직접 저비율 수열을 구성할 수 있는가?
• RQ2 삼각형 내에서 결정적인 점 집합이 달성할 수 있는 최적의 비율 속도는 무엇인가?
• RQ3 삼각형에서 리만 적분 가능한 함수에 대해 유계 변동성을 요구하지 않고 준몬테카를로 적분이 수렴할 수 있는가?
• RQ4 삼각형 기반 QMC 수열에 디지털 스캐러빙을 적용하여 부드러운 함수에 대해 수렴 속도를 향상시킬 수 있는가?
• RQ5 최적의 $O(\log N / N)$ 비율을 갖는 격자 유사 수열을 삼각형 내에서 직접 구성할 수 있는가?

주요 결과

• 삼각형 내에서의 디지털 반-더 코르푸트 구성법은 최대 $12/\sqrt{N}$의 평행사변형 비율을 갖으며, $N = 4^k$일 때 정확한 값은 $2/(3\sqrt{N}) - 1/(9N)$이다.
• 회전된 격자 구성법은 최적의 비율 속도인 $O(\log N / N)$을 달성하며, 삼각형 내에서 달성 가능한 최고의 비율 속도와 일치한다.
• 디지털 구성법의 자릿수를 스캐러빙하면 부드러운 함수에 대해 적분 정확도가 향상되며, $N = 4^k$일 때 루트 평균 제곱 오차가 $O(N^{-1})$이 된다.
• 두 구성법 모두 삼각형에서 리만 적분 가능한 함수에 대해 QMC 추정치가 진짜 적분값으로 수렴함을 보장한다. 이는 유계 변동성이 없어도 가능하다.
• 필라즈와 쿨스의 측도 비율은 브란돌리니 등이 정의한 평행사변형 비율의 두 배 이내이므로, 둘 다 $N \to \infty$일 때 0으로 수렴한다.
• 이론적 프레임워크는 표준 단위 삼각형 외에도 임의의 비퇴화 삼각형에 대해 수렴 가능성을 지원한다.