[논문 리뷰] Solving the Klein-Gordon equation using Fourier spectral methods: A benchmark test for computer performance

Samar Aseeri, Oleg Batrasev외 13명|arXiv (Cornell University)|2015. 01. 19.

Numerical Methods and Algorithms참고 문헌 27인용 수 6

한 줄 요약

이 논문은 13개의 다양한 초고성능 컴퓨팅 시스템을 대상으로 2DECOMP&FFT 라이브러리를 사용하여 3차원 입방형 클라인-고든 방정식을 해결하기 위한 푸리에 스펙트럴 방법의 성능을 벤치마킹한다. 강력한 스케일링 성능을 입증하고, 인터커넥트 지연과 프로세서 아키텍처의 영향으로 인한 주요 성능 변동을 규명하며, 복잡한 PDE 해법기에서의 실질적이고 실제적인 HPC 시스템 및 수치 해법에 대한 벤치마크를 제공한다.

ABSTRACT

The cubic Klein-Gordon equation is a simple but non-trivial partial differential equation whose numerical solution has the main building blocks required for the solution of many other partial differential equations. In this study, the library 2DECOMP&FFT is used in a Fourier spectral scheme to solve the Klein-Gordon equation and strong scaling of the code is examined on thirteen different machines for a problem size of 512^3. The results are useful in assessing likely performance of other parallel fast Fourier transform based programs for solving partial differential equations. The problem is chosen to be large enough to solve on a workstation, yet also of interest to solve quickly on a supercomputer, in particular for parametric studies. Unlike other high performance computing benchmarks, for this problem size, the time to solution will not be improved by simply building a bigger supercomputer.

연구 동기 및 목표

• 다양한 고성능 컴퓨팅 시스템에서 3차원 입방형 클라인-고든 방정식을 해결하기 위한 푸리에 스펙트럴 방법의 강력한 스케일링 성능을 평가하기 위해.
• 다양한 컴퓨터 아키텍처와 인터커넥트가 FFT 기반 PDE 해법기의 성능에 미치는 영향을 평가하기 위해, 특히 512³과 같은 중간 크기의 문제에 대해.
• 합성 또는 유사 합성 워크로드를 초월하는 실제 PDE 시뮬레이션의 계산 요구사항을 반영하는 실질적이고 실제적인 HPC 시스템용 벤치마크를 확립하기 위해.
• 향후 비선형 PDE 시뮬레이션, 특히 높은 처리량이 필요한 매개변수 연구에 최적의 하드웨어 및 소프트웨어 스택 선택을 안내하기 위해.
• 이 벤치마크가 기존의 HPC 순위 체계(예: Linpack 또는 NAS Parallel Benchmarks)의 보완 또는 업데이트로 활용될 수 있는지 탐색하기 위해.

제안 방법

• 입방형 클라인-고든 방정식은 512³ 직각좌표계 격자에서 푸리에 스펙트럴 공간 이산화를 사용한 의사스펙트럴 방법으로 해결된다.
• 2DECOMP&FFT 라이브러리를 활용하여 병렬 FFT를 수행하며, 2차원 블록 순환 분할 방식에서 MPI_Alltoallv 또는 MPI_Alltoall를 사용하여 프로세서 간 통신을 수행한다.
• 리더십 수준의 슈퍼컴퓨터와 워크스테이션을 포함한 13개의 이질적 HPC 시스템에서 강력한 스케일링 실험을 수행하여 해답 도달 시간과 효율성을 측정한다.
• 부동소수점 연산 수, 통신 대역폭, 인터커넥트 지연을 기준으로 성능를 분석하며, 다양한 프로세서 유형과 네트워크 구조 간의 결과를 비교한다.
• 초고성능 시스템에 적합한 크기이지만 워크스테이션에서 로컬로 후처리가 가능한 크기의 문제를 대상으로 하여 매개변수 연구가 가능하도록 한다.
• 벤치마크는 시간 적분, 비선형 항 평가, FFT 기반 스펙트럴 변환을 포함한 실제 PDE 시뮬레이션 워크로드를 반영하도록 설계되었다.

실험 결과

연구 질문

• RQ1 다양한 인터커넥트와 프로세서 유형을 가진 이질적 HPC 아키텍처에서 3차원 클라인-고든 방정식에 대한 푸리에 스펙트럴 방법의 강력한 스케일링 성능는 어떻게 변하는가?
• RQ2 현대 슈퍼컴퓨터에서 FFT 기반 PDE 해법기의 성능에 영향을 미치는 인터커넥트 지연과 메모리 대역폭의 차이 정도는 어느 정도인가?
• RQ3 이 벤치마크는 Linpack이나 NAS Parallel Benchmarks와 같은 기존 HPC 벤치마크의 의미 있는 대안 또는 보완으로서 실질적인 PDE 시뮬레이션 성능 평가에 활용될 수 있는가?
• RQ4 고정된 문제 크기에서 최신 인텔 프로세서는 기존 또는 다른 아키텍처(예: 파wrPC, AMD)와 비교해 FFT 기반 PDE 해법기 성능에서 어떻게 다를까?
• RQ5 이질적 HPC 시스템 간 FFT 기반 스펙트럴 방법의 이식성과 성능 이식성에 대해 어떤 통찰을 얻을 수 있는가?

주요 결과

• 유사한 최대 부동소수점 성능을 가진 시스템 간에도 인터커넥트 지연과 대역폭의 차이로 인해 푸리에 스펙트럴 해법기의 성능에 상당한 차이가 발생한다.
• 동일한 컴퓨팅 칩을 사용하고도 베두르(Vedur)는 헤토르(Hector)보다 성능이 현저히 열 劣하므로, 인터커넥트 지연이 주요 성능 저하 요인임을 시사한다.
• 비콘(Beacon), 호넷(Hornet), 마레노스트룸(Marenostrum), 스탬피드(Stampede) 등에서 신형 인텔 프로세서는 더 높은 부동소수점 처리량과 더 나은 메모리 컨트롤러 덕분에 이전 아키텍처(파워PC, AMD, 구형 인텔)를 초월한다.
• 관측된 성능 차이는 대역폭과 MPI_Alltoall 성능만으로는 완전히 예측할 수 없으며, FLOP/사이클 및 메모리 컨트롤러 효율성과 같은 마이크로아키텍처적 특성의 중요성을 강조한다.
• 이 벤치마크는 이 문제 크기에서는 더 큰 슈퍼컴퓨터를 구축한다고 해서 해답 도달 시간이 향상되지 않음을 보여주며, 알고리즘적 및 아키텍처적 선택이 지배하는 성능 정체점이 있음을 시사한다.
• 본 연구는 비선형 PDE의 매개변수 연구에 특히 중요한 스펙트럴 방법과 FFT를 활용하는 복잡한 PDE 시뮬레이션의 계산 요구사항을 반영하는 실질적이고 실제적인 HPC 시스템용 벤치마크를 제공한다.