[논문 리뷰] Multilevel Picard approximation algorithm for semilinear partial integro-differential equations and its complexity analysis

Ariel Neufeld, Sizhou Wu|arXiv (Cornell University)|2022. 05. 19.

Fractional Differential Equations Solutions인용 수 4

한 줄 요약

이 논문은 고차원 비선형 포아송-적분미분방정식(PIDEs)을 해결하기 위한 다수준 피카르(MLP) 근사 알고리즘을 제안한다. 점성해에 대한 파인만-카크 표현을 유도하고 시간 이산화와 몬테카를로 샘플링을 활용함으로써, 차원의 저주 없이 수렴성을 확보하며, 계산 복잡도는 차원 $d$와 정밀도의 역수 $
^{-1}$에 대해 다항식적으로 유 bounds된다. 이 방법은 최대 10,000 차원에서 수치적으로 검증되었다.

ABSTRACT

In this paper we introduce a multilevel Picard approximation algorithm for semilinear parabolic partial integro-differential equations (PIDEs). We prove that the numerical approximation scheme converges to the unique viscosity solution of the PIDE under consideration. To that end, we derive a Feynman-Kac representation for the unique viscosity solution of the semilinear PIDE, extending the classical Feynman-Kac representation for linear PIDEs. Furthermore, we show that the algorithm does not suffer from the curse of dimensionality, i.e. the computational complexity of the algorithm is bounded polynomially in the dimension $d$ and the reciprocal of the prescribed accuracy $\varepsilon$. We also provide a numerical example in up to 10'000 dimensions to demonstrate its applicability.

연구 동기 및 목표

• 고차원 비선형 포화 파라볼릭 PIDEs를 해결하기 위한 수치적 방법을 개발하여 차원의 저주를 극복한다.
• 연구 대상 PIDE의 유일한 점성해에 대한 파인만-카크 표현을 수립한다.
• 다수준 피카르 근사가 점성해로 수렴하는 것을 증명한다.
• 알고리즘의 계산 복잡도를 분석하고, 차원 $d$와 정밀도 $
  ^{-1}$에 대해 다항식적으로 척도가 되는지 보여준다.
• 최대 10,000 차원에서의 수치 예제를 통해 방법의 실용적 적용 가능성을 보여준다.

제안 방법

• 이 방법은 PIDE의 점성해에 대한 파인만-카크 표현에서 유도된 다수준 피카르 반복 체계에 기반한다.
• 해는 비틀림, 분산 및 레비 점프 성분을 포함하는 점프-확산 확률미분방정식(SDE)의 기댓값으로 표현된다.
• SDE에 시간 이산화를 적용하고, 기댓값은 몬테카를로 샘플링을 통해 근사한다.
• 알고리즘은 각 수준에서 이산화 및 샘플링 파rameter의 계층적 구조를 사용하여 근사를 정밀화하는 내재된 구조를 갖는다.
• 점프의 보상 적분은 이전 연구에서 선형 PIDEs에 대해 사용된 기법을 따르며, 레비 측도의 몬테카를로 샘플링을 통해 근사된다.
• 수렴성 및 복잡도 분석은 확률적 고정점 추론, 그론발의 부등식 및 이산화된 과정에 대한 모멘트 추정에 기반한다.

실험 결과

연구 질문

• RQ1 비선형 PIDEs에 대해 차원의 저주를 피할 수 있는 다수준 피카르 체계를 구성할 수 있는가?
• RQ2 계수의 적절한 정규성 및 성장 조건 하에서 제안된 알고리즘이 유일한 점성해로 수렴하는가?
• RQ3 비국소적(점프) 항을 포함하는 비선형 PIDEs에 대해 고전적인 선형 경우에서 확장된 파인만-카크 표현을 수립할 수 있는가?
• RQ4 알고리즘의 계산 복잡도는 차원 $d$와 정밀도 $
  $에 대해 어떻게 표현되는가?
• RQ5 10,000 차원과 같은 고차원 환경에서 이 방법은 얼마나 잘 작동하는가?

주요 결과

• 다수준 피카르 알고리즘이 계수의 온건한 정규성 및 성장 조건 하에서 비선형 PIDE의 유일한 점성해로 수렴한다.
• 비선형 PIDEs에 대한 새로운 파인만-카크 표현이 유도되었으며, 고전적인 선형 경우에서 비선형 PIDEs로의 확장을 가능하게 하였다.
• 알고리즘의 계산 복잡도는 차원 $d$와 주어진 정밀도의 역수 $
  ^{-1}$에 대해 다항식적으로 유 bounds되며, 이로써 차원의 저주를 극복하였다.
• 최대 10,000 개의 공간 차원을 가진 테스트 케이스에서 수치적으로 검증되었으며, 고차원에서의 실용 가능성 확인되었다.
• 수렴성 분석은 모멘트 추정과 확률적 고정점 추론에 기반하며, 이산화 및 샘플링된 근사의 오차에 엄밀한 경계를 제공한다.
• 오차 경계는 차원 $d$, 시간 영역 $T$, 그리고 점프 강도에 따라 달라지며, $M$ (점프 샘플 수) 및 $N$ (시간 단계 수)와 같은 매개변수를 통해 명시적으로 제어된다.