[논문 리뷰] A viscosity solution approach to the Monge-Ampere formulation of the Optimal Transportation Problem

Jean‐David Benamou, Brittany D. Froese외 1명|arXiv (Cornell University)|2012. 08. 23.

Geometric Analysis and Curvature Flows참고 문헌 31인용 수 12

한 줄 요약

이 논문은 상태 제약 조건이 있는 몽체-암페르 방정식을 통한 최적 운반 문제를 해결하기 위한 수렴성 보장된 유한차분 방법을 제시한다. 표적 영역의 부호거리 함수를 기반으로 한 해밀턴-자코비 방정식을 사용하여 제2의 경계값 문제를 재구성한다. 이 방법은 단조성 스킴과 필터링 스킴 이론을 활용하여 유일한 점성해로의 수렴을 보장하며, 일반적인 볼록 표적을 가진 운반 문제의 안정적인 수치적 해를 가능하게 한다.

ABSTRACT

In this work we present a numerical method for the Optimal Mass Transportation problem. Optimal Mass Transportation (OT) is an active research field in mathematics.It has recently led to significant theoretical results as well as applications in diverse areas. Numerical solution techniques for the OT problem remain underdeveloped. The solution is obtained by solving the second boundary value problem for the MA equation, a fully nonlinear elliptic partial differential equation (PDE). Instead of standard boundary conditions the problem has global state constraints. These are reformulated as a tractable local PDE. We give a proof of convergence of the numerical method, using the theory of viscosity solutions. Details of the implementation and a fast solution method are provided in the companion paper arXiv:1208.4870.

연구 동기 및 목표

• 일반적인 볼록 표적 영역을 가진 최적 운반 문제를 수치적으로 수렴하는 방법을 개발한다.
• 표준 경계 조건이 없는 몽체-암페르 방정식의 상태 제약 조건 문제를 다룬다.
• 표적 영역의 부호거리 함수를 이용하여 제2의 경계값 문제를 재구성함으로써 국소적이고 암묵적인 경계 처리를 가능하게 한다.
• 점성해 이론과 필터링 스킴을 활용하여 유한차분 이산화의 수렴성을 증명한다.
• 비볼록 및 퇴화된 경우를 포함한 최적 운반 문제의 효율적 수치적 해를 위한 기초를 제공한다.

제안 방법

• 최적 운반 경계 조건 ∇u(X) = Y를 표적 영역의 부호거리 함수 d(∇u(x), Y) = 0 을 포함하는 해밀턴-자코비 방정식으로 재구성한다. 여기서 x ∈ ∂X 이다.
• 거리 함수의 볼록성을 이용하여 해밀턴-자코비 경계 조건의 볼록성을 보장함으로써 해 u의 볼록성을 유지한다.
• 몽체-암페르 PDE와 경계 해밀턴-자코비 방정식을 단조성 유한차분 스킴을 사용하여 이산화하며, 기울기 조건을 활용하여 안정한 스킴을 구성한다.
• PDE와 경계 조건을 조합하여 거의 타원형 유한차분 스킴을 구성함으로써, 필터링 스킴 이론을 통해 수렴성을 보장한다.
• 점성해 이론과 바르레-수오가니디 프레임워크를 활용하여 이산해가 유일한 점성해로 수렴함을 증명한다.
• 함께 제공된 논문 [6] 에서 설명된 바와 같이 빠른 뉴턴 해법을 구현하여 효율적인 계산을 수행한다.

실험 결과

연구 질문

• RQ1 몽체-암페르 방정식의 제2의 경계값 문제는 국소적 PDE 시스템으로 재구성될 수 있는가? 이는 안정적인 수치적 해를 가능하게 하는가?
• RQ2 부호거리 함수와 해밀턴-자코비 경계 조건을 기반으로 한 유한차분 스킴은 최적 운반 문제의 점성해로 수렴하는가?
• RQ3 단조성 및 필터링 스킴은 상태 제약 조건이 있는 완전 비선형 타원형 PDE에 어떻게 적용되어 수렴성을 보장하는가?
• RQ4 부드럽고 비부드러운 표적 영역을 가진 최적 운반 문제에 대해 이 수치 방법의 수렴 속도는 어떠한가?
• RQ5 비볼록 표적과 비정상적인 경우(예: 분리된 소스 영역)를 처리할 수 있는가?

주요 결과

• 타원 영역에 대한 수치 결과로, 소스 및 표적 이산화 크기에서 일계 수렴을 달성함을 입증하였다.
• 타원에서 타원으로의 매핑의 경우, 32×32 그리드 해상도에서 최대 오차 0.1163에서 362×362에서 0.0056으로 감소하여 수렴성을 확인하였다.
• 분리된 소스 영역(예: 두 개의 반원이 원으로 매핑되는 경우)에 대해서도 최적 매핑을 성공적으로 계산하였으며, 고해상도에서 최대 오차가 약 0.0039로 안정화됨을 보였다.
• 비볼록 표적과 비부드러운 밀도에 대해서도 스킴은 안정적이고 수렴성을 유지하였으며, 직선과 곡선 경계를 혼합한 복합 기하 표적에 대한 계산된 매핑을 통해 이를 입증하였다.
• 필터링 스킴 이론을 활용하여 수렴성 증명를 수립하였으며, 해밀턴-자코비 재구성에 의해 상태 제약 조건을 다룰 수 있도록 바르레-수오가니디 프레임워크를 확장하였다.
• 이 방법은 일반적인 볼록 지지 영역과 밀도를 가진 최적 운반 문제의 정확한 해를 가능하게 하며, 퇴화 또는 불연속성으로 인해 기존 PDE 방법이 실패하는 경우에도 적용 가능하다.