[논문 리뷰] A two-grid accelerated sweeping preconditioner for the Helmholtz equation

Christiaan C. Stolk|arXiv (Cornell University)|2014. 12. 01.

Electromagnetic Scattering and Analysis참고 문헌 13인용 수 2

한 줄 요약

이 논문은 헬름홀츠 방정식에 대해 이중 그리드 가속화 스위핑 조건부 조건자(preconditioner)를 제안하며, 계산 비용을 줄이기 위해 스위핑 조건부 조건자를 해상도가 1/8인 굴곡된 문제에만 적용한다. 이 방법은 새로운 굴곡 스케일 연산자와 병렬화 전략을 결합하여 상태의 최고 수준 성능을 달성하고, 10^8 자유도까지의 상당한 비용 절감을 이룬다. 이는 파장 당 점 수가 적은 경우에도 빠른 수렴을 가능하게 한다.

ABSTRACT

Helmholtz solvers based on sweeping preconditioners have been applied successfully in large 3-D examples. However, such solvers are relatively expensive, due to the high cost per iteration. We propose to reduce this cost by using the sweeping preconditioner at the coarse level of a two-grid method, so that it is applied to an eight times smaller problem. For this to be attractive, the two-grid method must converge rapidly, even when the discretization is done using few points per wavelength. new coarse scale operator was developed to accomplish this, using the ideas of [C.C.Stolk, M. Ahmed and S.K. Bhowmik, A multigrid method for the Helmholtz equation with optimized coarse grid corrections, DOI:10.1137/13092349X]. Several strategies to parallellize the algorithm are proposed. Experiments with a 3-D, MPI implementation show state-of-the-art performance, and a large cost reduction compared to pure sweeping methods for problems with up to 10^8 degrees of freedom.

연구 동기 및 목표

• 대규모 3차원 문제에서 스위핑 조건부 조건자의 높은 반복당 비용을 줄이기 위해.
• 파장 당 점 수가 적은 경우에도 이중 그리드 방법의 수렴을 빠르게 하기 위해.
• 굴곡된 문제에서 정확성과 효율성을 유지하는 새로운 굴곡 스케일 연산자를 개발하기 위해.
• 이중 그리드 알고리즘의 확장성을 향상시키기 위해 병렬화 전략을 설계하기 위해.
• 기본 스위핑 조건부 조건자와 비교해 상당한 계산 비용 절감을 달성하기 위해.

제안 방법

• 반복당 계산 부담을 줄이기 위해 스위핑 조건부 조건자를 자유도가 1/8인 굴곡된 문제에만 적용한다.
• 다중 격자 방법에서 영감을 얻은 새로운 굴곡 스케일 연산자를 도입하여 최적화된 굴곡 격자 보정을 통해 굴곡 격자에서의 수렴을 유지한다.
• 이중 그리드 프레임워크를 사용하여 굴곡 격자 보정을 축소된 문제에 스위핑 조건부 조건자를 적용함으로써 계산한다.
• 다중 격자에서 영감을 얻은 굴곡 연산자를 활용해 파장 당 점 수가 적은 경우에도 강건하고 빠른 수렴을 보장한다.
• MPI를 사용해 병렬 실행을 구현하여 대규모 3차원 문제에서의 확장성을 향상시킨다.
• 병렬화 전략을 통합하여 굴곡 격자 해법 및 조건부 조건화 연산을 프로세서 간에 효율적으로 분배한다.

실험 결과

연구 질문

• RQ1 스위핑 조건부 조건자를 이중 그리드 방법의 굴곡 수준에서만 효과적으로 적용하여 전체 계산 비용을 줄일 수 있는가?
• RQ2 파장 당 점 수가 적은 헬름홀츠 문제에서 이중 그리드 방법의 수렴을 빠르게 하기 위해 새로운 굴곡 스케일 연산자를 어떻게 설계할 수 있는가?
• RQ3 대규모 3차원 헬름홀츠 문제에서 기존 스위핑 조건부 조건자와 비교해 성능 향상과 비용 절감 수준은 어느 정도 달성할 수 있는가?
• RQ4 대규모 문제에서 MPI를 사용해 이중 그리드 스위핑 조건부 조건자를 얼마나 효율적으로 병렬화할 수 있는가?
• RQ5 메esh가 굴곡된 경우(파장 당 점 수가 적은 경우)에도 제안된 방법이 정확성과 강건성을 유지하는가?

주요 결과

• 이중 그리드 가속화 스위핑 조건부 조건자는 자유도가 최대 10^8에 이르는 3차원 헬름홀츠 문제에서 최고 수준의 성능을 달성한다.
• 비용이 높은 조건부 조건자를 굴곡된 문제에만 적용함으로써 기존 스위핑 조건부 조건자와 비교해 계산 비용을 상당히 절감한다.
• 새로운 굴곡 스케일 연산자는 파장 당 점 수가 적은 경우에도 빠른 수렴을 가능하게 하며, 강건성을 유지한다.
• MPI 구현은 효과적인 병렬화를 보이며 대규모 시뮬레이션을 지원한다.
• 정확성과 확장성은 유지하면서도 대규모 3차원 문제에서 상당한 비용 절감을 달성한다.
• 특히 고주파 및 대규모 문제에서 순수한 스위핑 조건부 조건자보다 계산 효율성이 뛰어나다.