[논문 리뷰] A fast minimal residual solver for overlap fermions

Artan Boriçi, Alban Allkoci|ArXiv.org|2006. 02. 09.

Particle physics theoretical and experimental studies참고 문헌 12인용 수 4

한 줄 요약

이 논문은 라티스 QCD에서 이동 단위원 시스템을 위한 새로운 세항 재귀를 제안하여 아르놀디 과정에 적합하게 하고, 오버랩 페르미온 전파함수를 위한 두 가지 새로운 크릴로프 부분공간 해법기 SUOM과 SHUMR를 구축한다. SHUMR 알고리즘은 모든 시험된 격자 크기와 쿼크 질량에서 이전에 가장 빠른 방법인 SUMR보다 10–15% 더 빠른 수렴 속도를 보이며, 대규모 오버랩 페르미온 선형 시스템을 해결하는 데 있어 뛰어난 효율성을 입증한다.

ABSTRACT

Computing quark propagators with overlap fermions requires the solution of a shifted unitary linear system. Jagels and Reichel have shown that for such systems it is possible to construct a minimal residual algorithm by short recurrences. The Jülich-Wuppertal group have found this algorithm to be the fastest among overlap solvers. In this paper we present a three-term recurrence for the Arnoldi unitary process. Using the new recurrence we construct a minimal residual solver which is the fastest among all Krylov subspace algorithms considered so far for the overlap inversion.

연구 동기 및 목표

• 라티스 QCD에서 오버랩 페르미온을 사용한 쿼크 전파함수 계산을 위한 더 효율적인 반복 해법기를 개발하기 위해.
• 표준 GMRES에서 이동 단위원 시스템에 대해 긴 재귀로 인한 계산 블로킹 문제를 해결하기 위해.
• 오버랩 페르미온 시스템에서 잔여항 최소화의 최적성을 유지하면서도 짧은 재귀를 갖는 크릴로프 부분공간 방법을 구축하기 위해.
• 기존의 SUMR 알고리즘을 개선하기 위해 새로운 단위원 아르놀디 과정과 세항 재귀를 도입하기 위해.
• 실제 라티스 구조에서 기존의 방법들인 CGNE, CR, SUMR과의 비교를 통해 신규 해법기를 검증하기 위해.

제안 방법

• 표준 아르놀디 과정 대비 계산 비용과 메모리 사용을 줄이는 새로운 세항 재귀를 단위원 아르놀디 과정에 제안한다.
• 크릴로프 부분공간에 수직인 잔여항 벡터를 구성하는 이동 단위원 직교 방법(SUOM)을 유도한다.
• 새로운 재귀를 사용하여 크릴로프 부분공간에서 잔여항 노름을 최소화하는 이동 단위원 최소잔여법(SHUMR)을 도입한다.
• 과거의 벡터를 모두 축적하지 않고도 기저 벡터의 정규직교성을 유지하는 짧은 재귀를 통해 수정된 아르놀디 과정을 구현한다.
• 네우버거 오버랩 디랙 연산자(이동 단위원 행렬)의 구조를 활용하여 단위원성을 이용하고 짧은 재귀를 가능하게 한다.
• MATLAB에서 알고리즘을 구현하고, 고성능 계산을 위한 C++ 및 uBLAS 라이브러리로의 이식을 위한 코드를 제공한다.

실험 결과

연구 질문

• RQ1 이동 단위원 아르놀디 과정에 대해 세항 재귀를 유도할 수 있는가? 이를 통해 이동 단위원 시스템에 대해 짧은 재귀 크릴로프 해법기를 가능하게 할 수 있는가?
• RQ2 새로운 SHUMR 알고리즘이 오버랩 페르미온 시스템에서 이전에 최적의 성능을 보였던 SUMR 해법기보다 더 빠른 수렴 속도를 보이는가?
• RQ3 다양한 쿼크 질량과 격자 결합 상수에서 새로운 해법기(SUOM 및 SHUMR)의 수렴 속도와 안정성은 어떻게 비교되는가?
• RQ4 크릴로프 부분공간 방법에서 오버랩 페르미온에 대해 기하학적 최적성과 대수적 최적성이 수렴 행동에 미치는 영향은 무엇인가?
• RQ5 새로운 재귀 방법이 대규모 라티스 QCD 시뮬레이션에서 수치적 안정성과 효율성을 향상시킬 수 있는가?

주요 결과

• 새로운 세항 재귀는 이동 단위원 아르놀디 과정을 통해 처음으로 짧은 재귀 크릴로프 부분공간 방법을 가능하게 하여 계산 및 메모리 오버헤드를 감소시킨다.
• 모든 시험된 구성에서 SHUMR 알고리즘이 SUMR보다 10–15% 더 빠르게 수렴하며, 경량 쿼크 질량(m=0.005)과 다양한 β 값(β=5.7, 6.0)에서도 동일하게 성능을 발휘한다.
• SHUMR는 경량 쿼크 질량에서 초선형 수렴을 보이며, 잔여항 노름 약 10⁻³(β=6.0)과 10⁻⁵(β=5.7) 근처에서 수렴 속도가 점점 증가한다. 이는 CR의 선형 수렴과는 대조적이다.
• SUOM과 SUMR은 수렴 속도에서 근소한 격차를 보이지만, SHUMR은 항상 둘을 앞서며, 크릴로프 부분공간 구성 방식이 수렴 효율성에 영향을 준다는 것을 시사한다.
• SHUMR 및 SUOM과 같은 새로운 해법기는 CGNE 및 CG-CHI보다 더 효율적이며, 더 무거운 쿼크 질량에서도 성능이 열등하다.
• SHUMR 알고리즘은 기하학적으로 최적이며, 잔여항 노름을 최소화하는 데서도 SUMR보다 더 빠른 수렴 속도를 보이므로, SHUMR에서 크릴로프 부분공간 탐색이 더 효율적임을 시사한다.