[논문 리뷰] Linear scaling algorithms for solving high-dimensional nonlinear parabolic differential equations

E Weinan, Martin Hutzenthaler외 2명|arXiv (Cornell University)|2016. 07. 12.

Differential Equations and Numerical Methods인용 수 3

한 줄 요약

이 논문은 파동-켈러 및 비즈무트-엘워티-리 공식을 조합하고 다중수준 피카르 반복 분해를 적용하여 고차원 비선형 포아송 편미분방정식을 해결하는 새로운 선형 스케일링 알고리즘을 제안한다. 임의의 $\delta > 0$ 에 대해 $O(d\varepsilon^{-(4+\delta)})$ 의 계산 복잡도를 달성하며, 기울기 독립 비선형성을 갖는 반선형 열 방정식에 대해 효율적인 수렴성을 입증한다.

ABSTRACT

We introduce a new family of numerical algorithms for approximating solutions of general high-dimensional semilinear parabolic partial differential equations at single space-time points. The algorithm is obtained through a delicate combination of the Feynman-Kac and the Bismut-Elworthy-Li formulas, and an approximate decomposition of the Picard fixed-point iteration with multilevel accuracy. The algorithm has been tested on a variety of semilinear partial differential equations that arise in physics and finance, with very satisfactory results. Analytical tools needed for the analysis of such algorithms, including a semilinear Feynman-Kac formula, a new class of semi-norms and their recursive inequalities, are also introduced. They allow us to prove for semilinear heat equations with gradient-independent nonlinearity that the computational complexity of the proposed algorithm is bounded by $O(d\,\varepsilon^{-(4+\delta)})$ for any $\delta \in (0,\infty)$ under suitable assumptions, where $d\in \mathbb{N}$ is the dimensionality of the problem and $\varepsilon\in(0,\infty)$ is the prescribed accuracy.

연구 동기 및 목표

• 물리학 및 금융 분야에서 고차원 비선형 포아송 편미분방정식의 계산 비가역성 문제를 해결하기 위해.
• 해당 편미분방정식을 효율적으로 해결하기 위해 차원에 대해 선형 스케일링을 보장하는 수치 알고리즘을 개발하기 위해.
• 일반적인 가정 하에 제안된 방법의 엄밀한 복잡도 경계를 설정하기 위해.
• 반선형 파인만-카프 공식과 재귀적 반노름 부등식을 포함한 새로운 분석 도구를 도입하기 위해.
• 알고리즘이 정확도 $\varepsilon$ 를 달성할 때 $O(d\varepsilon^{-(4+\delta)})$ 의 복잡도를 갖는다는 것을 증명하기 위해.

제안 방법

• 해결 방법의 확률적 표현을 위한 파인만-카프 공식과 도함수 추정을 위한 비즈무트-엘워티-리 공식을 조합한다.
• 오차를 수준 간에 제어하기 위해 피카르 고정점 반복의 다중수준 근사화를 시행한다.
• 오차 전파를 분석하기 위해 새로운 종류의 반노름과 그 재귀적 부등식을 도입한다.
• 해결 과정을 계층적인 수준으로 분해하여 효율적인 계산을 가능하게 한다.
• 물리학 및 금융 분야의 반선형 편미분방정식에 대해 고차원 시험 케이스를 대상으로 방법을 검증한다.
• 이론적 분석은 확률적 표현과 재귀적 오차 경계를 기반으로 복잡도 추정을 유도한다.

실험 결과

연구 질문

• RQ1 고차원 비선형 포아송 편미분방정식을 위한 수치 알고리즘이 차원에 대해 선형 스케일링을 달성할 수 있는가?
• RQ2 기울기 독립 비선형성을 갖는 반선형 열 방정식을 풀기 위해 확률적 고정점 반복을 사용할 경우 계산 복잡도는 어떻게 되는가?
• RQ3 파인만-카프 공식과 비즈무트-엘워티-리 공식을 어떻게 조합하여 다중수준 근사화를 가능하게 할 수 있는가?
• RQ4 고차원 확률적 편미분방정식 해법에서 오차를 제한하기 위해 어떤 새로운 분석 도구가 필요한가?
• RQ5 적절한 가정 하에 제안된 방법이 임의의 $\delta > 0$ 에 대해 $O(d\varepsilon^{-(4+\delta)})$ 의 복잡도를 달성할 수 있는가?

주요 결과

• 제안된 알고리즘은 기울기 독립 비선형성을 갖는 반선형 열 방정식을 해결할 때 임의의 $\delta > 0$ 에 대해 $O(d\varepsilon^{-(4+\delta)})$ 의 계산 복잡도를 달성한다.
• 계수 및 비선형성에 대한 적절한 정규성 및 유계성 가정 하에 복잡도 경계가 유지된다.
• 오차를 유지하면서 계산 비용을 감소시키기 위해 피카르 반복의 다중수준 분해를 활용한다.
• 일반적인 반선형 포아송 편미분방정식의 해를 표현하기 위해 새로운 반선형 파인만-카프 공식을 유도한다.
• 오차 전파를 제어하기 위해 새로운 종류의 반노름과 그 재귀적 부등식을 도입하고 이를 활용한다.
• 물리학 및 금융 분야의 벤치마크 문제에서 강력한 경험적 성능을 보여준다.