[논문 리뷰] Runge-Kutta Theory and Constraint Programming

Julien Alexandre dit Sandretto|arXiv (Cornell University)|2018. 04. 13.

Numerical Methods and Algorithms참고 문헌 16인용 수 3

한 줄 요약

이 논문은 안정성, 해밀토니안 성질, 순서 정확도를 보장하는 조건을 만족시키는 데에 구간 분석을 활용한 새로운 제약 프로그래밍 접근법을 제안한다. 룬게-쿠타 계수를 계산하기 위해 브타르 표기법 조건을 제약 만족 문제로 공식화하고 전역 목적 함수를 최적화함으로써, 기존의 부동소수점 방법보다 더 좁은 해 구속을 제공하는 검증된 시뮬레이션 환경에서 뛰어난 성능을 보이는 새로운 구간 계수를 가진 룬게-쿠타 스킴을 도출한다.

ABSTRACT

There exist many Runge-Kutta methods (explicit or implicit), more or less adapted to specific problems. Some of them have interesting properties, such as stability for stiff problems or symplectic capability for problems with energy conservation. Defining a new method suitable to a given problem has become a challenge. The size, the complexity and the order do not stop growing. This informal challenge to implement the best method is interesting but an important unsolved problem persists. Indeed, the coefficients of Runge-Kutta methods are harder and harder to compute, and the result is often expressed in floating-point numbers, which may lead to erroneous integration schemes. Here, we propose to use interval analysis tools to compute Runge-Kutta coefficients. In particular, we use a solver based on guaranteed constraint programming. Moreover, with a global optimization process and a well chosen cost function, we propose a way to define some novel optimal Runge-Kutta methods.

연구 동기 및 목표

• 부동소수점 산술을 사용할 경우 수치적 통합에서 오류가 발생할 위험이 있는 정확한 룬게-쿠타 계수를 계산하는 데 점점 어려움이 증가하고 있는 데에 대응하기 위해.
• 부동소수점 계수를 구간 계수로 대체함으로써 순서, 안정성, 해밀토니안 성질과 같은 핵심 성질을 유지하기 위해.
• 제약 프로그래밍과 전역 최적화를 활용하여 새로운 룬게-쿠타 방법을 체계적이고 검증 가능한 방식으로 발견하기 위해.
• 검증된 시뮬레이션 환경에서 구간 계수 스킴의 효과성을 입증하기 위해.
• 이론적 성질을 엄격한 계산 하에 유지하면서 최적의 구간 계수를 가진 룬게-쿠타 방법을 생성하기 위한 프레임워크를 제공하기 위해.

제안 방법

• 순서 조건, 안정성, 해밀토니안 성질 등의 룬게-쿠타 방법 조건을 구간 변수에 대한 비선형 방정식 및 부등식의 체계로 공식화한다.
• 모든 해가 유효하고 포함관계가 정확한 구간 상자들을 찾기 위해 보장 가능한 제약 프로그래밍(CP) 솔버를 사용한다.
• 지역 절단 오차(LTE)를 다음 차수의 방법에 비해 최소화하도록 설계된 목적 함수를 사용하여 전역 최적화 절차를 적용한다. 이는 최적의 스킴을 선호한다.
• 건설 및 검증 과정 全 주기에 걸쳐 구간 산술을 적용하여 엄밀한 오차 경계와 성질 보장을 유지한다.
• 실험적 평가를 위해 결과 스킴을 DynIbex 검증된 시뮬레이션 툴에 구현한다.
• 기존의 RadauIIA, 가우스-레지온드르, RK4 등의 방법과의 비교를 통해 결과를 검증한다.

실험 결과

연구 질문

• RQ1 룬게-쿠타 방법에 구간 계수를 사용할 경우, 부동소수점 산술에서 자주 잃어버리는 순서, 안정성, 해밀토니안 성질 등의 본질적 성질을 유지할 수 있는가?
• RQ2 구간 계수를 가진 브타르 표기법 조건에서 발생하는 고차원적이고 과소정의된 방정식 체계를 효과적으로 해결하기 위해 제약 프로그래밍을 활용할 수 있는가?
• RQ3 다음 차수의 방법에 비해 지역 절단 오차(LTE)를 최소화하는 최적화 과정을 통해 경쟁력 있는 정확도와 효율성을 가진 룬게-쿠타 스킴를 도출할 수 있는가?
• RQ4 작은 계수 구간을 가진 구간 계수 스킴이 부동소수점 대비 더 좁은 해 구속을 제공할 수 있는가, 검증된 시뮬레이션 환경에서?
• RQ5 매개변수 불확실성 또는 강성 있는 역학 시스템이 존재할 경우, 구간 기반 방법이 수 contracitivity와 강건성 측면에서 얼마나 향상되는가?

주요 결과

• 구간 계수를 가진 3차 스킴인 S3O4는 5차 RadauIIA 방법과 유사한 성능을 보였지만 계산 속도가 더 빠르며, 높은 효율성을 입증했다.
• 구간 계수를 가진 3차 스킴인 S3O5는 6차 가우스-레지온드르 방법보다 5배 빠른 속도로 실행되었고, 최종 해 구속이 훨씬 좁았다.
• 구간 계수를 가진 명시적 3차 스킴인 ERK33은 RK4와 유사한 시간을 소비했지만, 해 구속 폭은 RK4의 1.2배에 불과하여 뛰어난 경쟁력을 보였다.
• 계수의 지름이 1×10⁻¹⁰에 불과한 구간 계수를 사용함에도 불구하고, 새로운 방법의 최종 해 구속이 종종 정확한 계수를 가진 고전적 방법보다 좁았으며, 이는 개선된 수축성(contractivity)을 시사한다.
• 제약 프로그래밍 접근법은 계수를 포함관계가 보장되는 구간으로 표현하더라도 순서와 안정성 성질을 유지하는 새로운 룬게-쿠타 스킴을 성공적으로 도출하였다.
• 최적화 과정은 LTE를 효과적으로 최소화하여, 더 높은 차수의 방법에 가까운 성능을 보이는 스킴를 도출하였으며, 목적 함수의 효과성을 검증하였다.