[논문 리뷰] Overcoming the curse of dimensionality in the numerical approximation of high-dimensional semilinear elliptic partial differential equations

Christian Beck, Lukas Gonon외 1명|arXiv (Cornell University)|2020. 03. 01.

Advanced Numerical Methods in Computational Mathematics참고 문헌 59인용 수 29

한 줄 요약

이 논문은 리프시츠 비선형성을 가진 고차원 반선형 타원형 편미분방정식(PDE)을 위한 새로운 전역역사 재귀 다수준 피카르(MLP) 근사 방법을 제안하고 분석한다. 이 방법이 주어진 정밀도 $\varepsilon$를 달성하기 위한 계산 비용이 차원 $d$와 $\varepsilon^{-1}$에 대해 최대 다항식적으로 증가함을 증명함으로써, 이 맥락에서 차원의 저주를 처음으로 극복한다.

ABSTRACT

Recently, so-called full-history recursive multilevel Picard (MLP) approximation schemes have been introduced and shown to overcome the curse of dimensionality in the numerical approximation of semilinear parabolic partial differential equations (PDEs) with Lipschitz nonlinearities. The key contribution of this article is to introduce and analyze a new variant of MLP approximation schemes for certain semilinear elliptic PDEs with Lipschitz nonlinearities and to prove that the proposed approximation schemes overcome the curse of dimensionality in the numerical approximation of such semilinear elliptic PDEs.

연구 동기 및 목표

• 고차원 반선형 타원형 PDE에서 차원의 저주를 극복하는 엄밀한 방법의 부재를 해결하기 위해.
• 이전에 타원형 PDE가 아닌 쌍곡형 PDE에만 적용된 성공적인 MLP 프레임워크를 리프시츠 비선형성을 가진 타원형 PDE로 확장하기 위해.
• 이러한 타원형 PDE의 수치 근사가 차원에 따라 지수적 비용 증가 없이 효율적으로 계산될 수 있음을 보여주는 이론적 기반을 마련하기 위해.
• 타원형 설정에서 제안된 MLP 방법에 대한 완전한 오차 및 복잡도 분석을 제공하기 위해.

제안 방법

• 해의 스토크래틱 표현을 피카르 공식을 통해 유도한 스토크래틱 고정점 방정식(SFPE)을 통해 사용한다.
• 브라운 운동 경로에 대한 기대값을 반복적으로 정밀화함으로써 해를 근사하는 재귀적 다수준 피카르 반복을 도입한다.
• 각 수준에서 독립적인 몬테카를로 샘플을 생성하기 위해 i.i.d. 브라운 운동과 i.i.d. 정지 시간 $R^\theta$를 사용한다.
• 알고리즘은 $U^{d,\theta}_n(x)$ 근사를 재귀적으로 계산하며, 전역역사 재귀 구조를 통해 선형항과 비선형항을 결합한다.
• 해의 차원에 대한 의존성의 증가를 억제하기 위해, 차원 $d$에 따라 스케일링되는 선형 연산자 $B_d$를 도입한다.
• 함수 평가 횟수에 대한 재귀적 bound $\mathfrak{C}_{d,n}$을 통해 복잡도를 분석하며, 이는 $d$와 $\varepsilon^{-1}$에 대해 다항식적으로 증가함을 보였다.

실험 결과

연구 질문

• RQ1 효율적인 파라볼릭 PDE에 적용된 MLP 근사 방법이 반선형 타원형 PDE로 확장될 수 있는가?
• RQ2 제안된 타원형 PDE에 대한 MLP 방법이 차원의 저주를 피할 수 있는가?
• RQ3 계산 비용이 차원 $d$와 목표 정밀도 $\varepsilon$에 대해 정확히 어떻게 의존하는가?
• RQ4 비선형성과 확산 연산자에 대한 가정이 수렴성과 복잡도에 어떻게 영향을 미치는가?

주요 결과

• 정밀도 $\varepsilon$를 달성하기 위한 계산 비용이 차원 $d$와 역정밀도 $\varepsilon^{-1}$에 대해 최대 다항식적으로 증가함을 증명함으로써, 차원의 저주가 존재하지 않음을 입증한다.
• $[-cd^c, cd^c]^d$ 영역에서 $L^2$ 기준 오차가 $\varepsilon$ 이하로 유계이며, 이 유계성은 $d$에 대해 균일하게 유지된다.
• 복잡도 유계 $\mathfrak{C}_{d,\mathfrak{N}_{\varepsilon,d}}$는 어떤 $\kappa > 0$에 대해 최대 $\kappa d^\kappa \varepsilon^{-\kappa}$로 유계이며, 이는 다항식 증가를 의미한다.
• 비선형성 $f_d$와 확산 연산자 $B_d$에 대한 약한 가정 하에, 방법이 잘 정의되고 적분 가능하다는 것이 분석에 의해 입증된다.
• 다른 분기 브라운 운동 방법이 큰 시간 수평에서 실패하는 것과 달리, 이 방법은 모든 시간 수평과 모든 차원에서 수렴성을 확보한다.
• 이론적 프레임워크는 비선형성이 리프시츠 조건을 만족하고 영역이 $d$에 따라 증가하더라도, MLP 방법이 차원 증가에 대해 강건함을 확인한다.