[논문 리뷰] PETSc/TS: A Modern Scalable ODE/DAE Solver Library

Shrirang Abhyankar, Jed Brown외 4명|arXiv (Cornell University)|2018. 06. 04.

Numerical methods for differential equations참고 문헌 25인용 수 147

한 줄 요약

논문은 PETSc 프레임워크 위에 구축된 PETSc/TS로, forward/adjoint sensitivities, 다양한 타임스텝 스킴, 적응 오차 제어, 유연한 해석기 옵션을 지원하는 ODE/DAE 해결을 위한 확장 가능하고 일반-purpose 라이브러리이다.

ABSTRACT

High-quality ordinary differential equation (ODE) solver libraries have a long history, going back to the 1970s. Over the past several years we have implemented, on top of the PETSc linear and nonlinear solver package, a new general-purpose, extensive, extensible library for solving ODEs and differential algebraic equations (DAEs). Package includes support for both forward and adjoint sensitivities that can be easily utilized by the TAO optimization package, which is also part of PETSc. The ODE/DAE integrator library strives to be highly scalable but also to deliver high efficiency for modest-sized problems. The library includes explicit solvers, implicit solvers, and a collection of implicit-explicit solvers, all with a common user interface and runtime selection of solver types, adaptive error control, and monitoring of solution progress. The library also offers enormous flexibility in selection of nonlinear and linear solvers, including the entire suite of PETSc iterative solvers, as well as several parallel direct solvers.

연구 동기 및 목표

• PETSc와 통합된 확장 가능하고 일반-purpose ODE/DAE 해석기 라이브러리 제공.
• 적응 오차 제어와 함께 명시적, 암시적, IMEX 타임스테핑 스킴에 대한 통합 인터페이스 제공.
• TS 프레임워크 내에서 순방향 및 어드조인트 민감도 계산을 가능하게 하며 TAO 최적화와의 통합.
• 유연한 설계와 성능 향상을 위한 선형 및 비선형 해석기와 야코비안 처리 도구의 폭넓은 지원.
• 이벤트 처리, 모니터링 및 시각화 도구를 통해 해 해결 과정 돕기

제안 방법

• PETSc 아키텍처와 그 핵심 클래스(Vec, Mat, KSP, SNES, DM) 및 이들이 TS 객체를 어떻게 뒷받침하는지 설명.
• ODE/DAE 문제를 위한 TS API(F(t,u,du/dt), G(t,u))와 그 야코비안, 그리고 시프트 기반 야코비안의 사용 방법 설명.
• 사용 가능한 타임스테핑 스킴(명시적, 암시적, 분할 IMEX/ARK)과 이들의 안정성/정확도 특성 상세화.
• 사용자 제공 야코비안을 돕기 위한 야코비안 계산 도구(매트릭스 프리, 컬러링, 자동 테스트) 상세화.
• DM 객체가 이산화/메시 관리와 해석기를 연결해 공간 이산화를 포함하는 PDE를 처리하는 방법 개요.
• 설정 및 실행 예시를 C와 Python(petsc4py)에서 TS 해 solver의 설정 및 실행 예시를 제공

실험 결과

연구 질문

• RQ1 분산 메모리 아키텍처에서 현대적인 ODE/DAE 해석기를 매우 확장 가능하도록 설계하면서도 다양한 문제 클래스(ODEs, DAEs, 강성/비강성, 질량 행렬)를 위한 유연성을 유지하려면 어떻게 해야 하는가?
• RQ2 야코비안을 사용자로부터 공급하고 검증하는 실용적 메커니즘은 무엇이며, 행렬 프리 및 컬러링 기법을 어떻게 통합해 이 부담을 완화할 수 있는가?
• RQ3 순방향 및 어드조인트 민감도가 TS에 얼마나 효과적으로 통합되고 TAO 최적화 패키지에서 얼마나 활용될 수 있는가?
• RQ4 한 API가 명시적, 암시적 및 IMEX 방법(분할 Runge-Kutta 스킴 포함)과 적응 오차 제어 및 이벤트 처리를 어떻게 수용할 수 있는가?

주요 결과

• TS 인터페이스는 명시적, 암시적, IMEX 방법을 포함한 광범위한 시간 적분 스킴을 지원하며, 적응 오차 제어와 런타임 구성 가능성을 갖추고 있다.
• 사용자는 F와 G 및 그들의 야코비안을 유연하게 조합할 수 있으며, 매트릭스 프리 야코비안 적용, 명시적 야코비안의 컬러링, 계산 검증을 위한 자동 테스트 도구를 제공한다.
• PETSc 기반 구현은 최소 공개 API와 런타임 옵션 시스템을 통해 반복 및 직접 해를 포함한 선형 및 비선형 해석기 생태계를 활용한다.
• 프레임워크는 이벤트 처리, 해 과정 모니터링 및 시각화, 네트워크 및 화학 반응에 대한 고수준 인터페이스를 제공하는 해석 과정 관리 메커니즘을 제공한다.
• 코드 예시를 통해 C와 Python에서 TS를 끝까지 사용하는 방법을 보여주며 문제 정의, 야코비안 설정 및 해 해결 루프를 포함한다.