[논문 리뷰] Analysis of a Time Multigrid Algorithm for DG-Discretizations in Time

Martin J. Gander, Martin Neumüller|arXiv (Cornell University)|2014. 09. 18.

Advanced Numerical Methods in Computational Mathematics참고 문헌 21인용 수 7

한 줄 요약

이 논문은 시간에 대한 불연속 갈레르킨(DG) 이산화에 대해 시간 다중격자 알고리즘을 제시하고 분석한다. 블록 자코비 스무스너와 시간에 대한 표준 제약/확장 연산을 사용한다. 점근적으로 최적의 스무스닝 파rameter를 유도하고, 두 격자 수렴 추정을 엄밀히 증명하여, 대규모 시스템에서 파라볼릭 PDE의 효율적인 병렬 시간 해법을 가능하게 한다. 강한 및 약한 스케일링 성능이 뛰어나다.

ABSTRACT

We present and analyze for a scalar linear evolution model problem a time multigrid algorithm for DG-discretizations in time. We derive asymptotically optimized parameters for the smoother, and also an asymptotically sharp convergence estimate for the two grid cycle. Our results hold for any A-stable time stepping scheme and represent the core component for space-time multigrid methods for parabolic partial differential equations. Our time multigrid method has excellent strong and weak scaling properties for parallelization in time, which we show with numerical experiments.

연구 동기 및 목표

• 시간에 대한 DG를 사용한 진화 문제의 시간 이산화에 대해 강력한 다중격자 방법을 개발하고 분석한다.
• 기존의 공간만을 고려한 다중격자 방법의 한계를 극복하여 시간 방향으로의 효율적 병렬화를 가능하게 한다.
• 시간 성분을 독립적으로 분석함으로써 공간-시간 다중격자 방법에 대한 이론적 기반을 마련한다.
• 두 격자 순환에 대해 점근적으로 최적의 스무스닝 파rameter와 수렴 한계를 확립한다.
• 대규모 슈퍼컴퓨터에서 제안된 시간 다중격자 알고리즘의 강한 및 약한 스케일링 성능을 시현한다.

제안 방법

• 시간에 대한 불연속 갈레르킨(DG) 이산화를 사용하며, 이는 시간 경계에서 연속성의 약한 강제화에 기반한다. 이의 차수는 $ p_t $이다.
• DG 시간 이산화에서 유도된 선형 연립방정식에 블록 자코비 스무스너를 적용하며, 점근적으로 최적화된 스무스닝 파rameter를 사용한다.
• 시간 방향으로 표준 제약 및 확장 연산자를 적용하여 시간 차원에서의 완전한 해상도를 확보한다.
• 두 격자 방법의 최적 스무스닝 파arameter와 수렴 인자 결정을 위해 푸리에 모드 분석을 수행한다.
• 두 격자 수렴 분석을 통해 시간 간격 크기와 무관한 점근적으로 날카로운 수렴 추정을 확립한다.
• Monte Rosa 및 Vienna Scientific Cluster 슈퍼컴퓨터에서의 수치 실험을 통해 스케일링 성능을 검증한다.

실험 결과

연구 질문

• RQ1 DG 이산화를 사용하는 시간 다중격자 방법에서 블록 자코비 스무스너에 대한 최적의 스무스닝 파aram터는 무엇인가?
• RQ2 시간 다중격자 방법이 DG 시간 이산화에 대해 메esh-독립적인 수렴 성능를 달성할 수 있는가?
• RQ3 두 격자 방법의 수렴 행동은 시간 이산화 차수 $ p_t $와 스무스닝 파aram터에 어떻게 의존하는가?
• RQ4 대규모 병렬 아키텍처에서 제안된 시간 다중격자 알고리즘의 강한 및 약한 스케일링 성능은 어떠한가?
• RQ5 이 방법은 공간 해법이 효율적인 파라볼릭 PDE에 대해 전체 공간-시간 다중격자로 확장될 수 있는가?

주요 결과

• 블록 자코비 스무스너에 대한 점근적으로 최적의 스무스닝 파aram터가 도출되었으며, 이는 두 격자 방법의 최적 수렴 속도를 보장한다.
• 엄밀한 두 격자 수렴 추정이 증명되었으며, 시간 간격 크기와 무관하고 임의의 A-안정적인 시간 적분 방법에 대해 유효하다.
• 메서드는 뛰어난 강한 스케일링 성능를 보이며, 코어 수가 두 배로 증가할수록 계산 시간이 약 절반으로 줄어들며, 지역 문제 크기에 의해 한계에 도달한다.
• 약한 스케일링 성능 또한 뛰어나며, 문제 크기와 코어 수가 비례적으로 증가함에 따라 총 계산 시간이 거의 일정하게 유지된다.
• 131,120,896개의 자유도를 가진 대규모 문제에서, 2048개의 코어에서 병렬 실행 시 순차적 전진 대입 방법 대비 1000배의 속도 향상을 달성한다.
• 다양한 다항식 차수 $ p_t $에 걸쳐도 안정적인 성능를 유지하며, $ p_t $ 증가에 따라 계산 시간이 유리하게 스케일링된다.