[논문 리뷰] Local Exponential Methods: a domain decomposition approach to exponential time integration of PDEs

Luca Bonaventura|arXiv (Cornell University)|2015. 05. 09.

Numerical methods for differential equations참고 문헌 19인용 수 3

한 줄 요약

이 논문은 전역 행렬 지수를 국소적이고 병렬 처리 가능한 계산으로 대체하는 도메인 분할 접근법인 국소 지수 방법(Local Exponential Methods, LEM)을 제안한다. 이는 공간적 감쇠 성질을 활용하여 Iserles의 이론적 경계에 의해 뒷받침되는 행렬 지수의 공간적 감쇠 특성을 이용함으로써 계산 비용을 감소시키고, 작은 국소 지수 행렬의 저장을 가능하게 하여 고 쿠라우트 수에서 표준 방법과 유사한 정확도를 확보하면서도 확장성 향상을 이룬다.

ABSTRACT

A local approach to the time integration of PDEs by exponential methods is proposed, motivated by theoretical estimates by A.Iserles on the decay of off-diagonal terms in the exponentials of sparse matrices. An overlapping domain decomposition technique is outlined, that allows to replace the computation of a global exponential matrix by a number of independent and easily parallelizable local problems. Advantages and potential problems of the proposed technique are discussed. Numerical experiments on simple, yet relevant model problems show that the resulting method allows to increase computational efficiency with respect to standard implementations of exponential methods.

연구 동기 및 목표

• 전체 행렬 지수를 국소적이고 병렬 처리 가능한 계산으로 대체함으로써 PDE의 지수 시간 적분에 대한 계산 비용을 줄이기 위해.
• Iserles의 희소 행렬에 대한 감쇠 추정에 의해 이론적으로 뒷받침되는 행렬 지수의 공간 국소성 특성을 활용하기 위해.
• 반복 Krylov 방법을 피하기 위해 Padé 근사법을 통해 작은 국소 지수 행렬을 저장할 수 있도록 하기 위해.
• 특히 대기 및 해양 모델링에서 중요한 고차수 및 이방성 메esh에서의 확장성과 효율성을 향상시키기 위해.
• 고수직 쿠라우트 수를 가진 문제에 대해 방향 분할의 완전히 국소적인 대안을 제공하기 위해.

제안 방법

• 계산 도메인을 겹치는 부분을 가진 서브도메인으로 분할하기 위해 겹치는 도메인 분할을 적용한다. 이는 원거리 영향을 최소화하기 위한 버퍼 영역을 포함한다.
• 각 서브도메인에 대해 스케일링과 승법을 적용한 Padé 근사법을 사용하여 국소 행렬 지수를 계산함으로써 저장 및 재사용이 가능하게 한다.
• 국소 행렬의 지수를 사용하여 각 서브도메인 내에서 해를 전진시키며, 경계 노드는 공유된 것으로 간주한다.
• 물리적 전파 속도와 시간 간격에 기반하여 서브도메인을 겹쳐서 설정함으로써 원거리 노드의 영향이 무시 가능하도록 한다.
• 각 서브문제를 독립적으로 병렬로 풀어 통신 비용을 줄이고 확장성을 향상시킨다.
• 지수 로젠브룩 또는 기타 지수 시간 적분법과 통합하여 시간에 대한 고차 정확도를 유지한다.

실험 결과

연구 질문

• RQ1 지수 시간 적분에 있어 전역 행렬 지수를 근사하기 위해 국소 지수 행렬을 사용할 수 있는가? 이때 정확도는 유지되는가?
• RQ2 겹치는 서브도메인을 가진 도메인 분할이 기존 지수 방법에 비해 계산 비용을 줄이는가?
• RQ3 버퍼 영역 크기가 고 쿠라우트 수에서 효율성에 어느 정도 영향을 미치는가?
• RQ4 고차수 또는 이방성 메쉬에서 Padé 근사법을 통해 국소 지수 행렬을 효율적으로 저장하고 계산할 수 있는가?
• RQ5 강한 이방성 특성을 보이는 모델, 예를 들어 대기 및 해양 시뮬레이션에서 LEM은 반응형 또는 반라그랑주 방법과 경쟁력이 있는가?

주요 결과

• 선형 이동 문제에서 쿠라우트 수 7로 실행할 때, 국소 지수 방법(LEM)은 4차 수준의 명시적 룬게쿠타 방법과 유사한 2차 정확도를 확보한다.
• 이방성 메쉬에서의 비선형 2차원 버거스 방정식에 대해, LEM은 수직 쿠라우트 수 6, 수평 쿠라우트 수 1 이하에서 정확도를 유지하였으며, 열 방향 도메인 분할을 사용하였다.
• 겹치는 버퍼 영역으로 인한 오버헤드가 있음에도 불구하고, 효율적인 국소 행렬 계산 및 저장 덕분에 기존 지수 방법에 비해 계산 비용을 감소시켰다.
• Padé 근사법을 통해 국소 지수 행렬을 저장함으로써 반복 Krylov 방법의 필요성을 피할 수 있었고, 특히 서브도메인 크기가 작을 경우 단계당 비용을 크게 감소시켰다.
• 고차수 불연속 갈레르킨 방법 및 이방성 메쉬에서 LEM은 고쿠라우트 수가 높은 경우가 흔하고 전통적 방법이 비효율적인 상황에서 유망한 성능을 보였다.
• 강한 수직 이방성 특성을 가진 모델, 예를 들어 기상 예측 및 해양 모델링에서 방향 분할의 완전히 국소적인 대안을 제공한다.