[논문 리뷰] Speeding up Domain Wall Fermion Algorithms using QCDLAB

Artan Boriçi|ArXiv.org|2007. 03. 21.

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

한 줄 요약

이 논문은 도메인 월 페르미온에 대한 격자 QCD 알고리즘의 프로토타이핑을 가속화하기 위해 저수준 선형 대수 연산을 추상화한 고수준 MATLAB/OCTAVE 기반 계산 프레임워크인 QCDLAB 1.1을 소개한다. 페르미온 행렬-벡터 연산과 카일로프 부분공간 해법기의 빠른 구현 및 테스트를 가능하게 함으로써, QCDLAB 1.1은 복잡한 격자 양자장 이론 시뮬레이션의 개발 시간을 크게 단축시키면서도 표준 QCD 알고리즘과의 호환성을 유지한다.

ABSTRACT

Simulating lattice QCD with chiral fermions and indeed using Domain Wall Fermions continues to be challenging project however large are concurrent computers. One obvious bottleneck is the slow pace of prototyping using the low level coding which prevails in most, if not all, lattice projects. Recently, we came up with a new proposal, namely QCDLAB, a high level language interface, which we believe will boost our endeavours to rapidly code lattice prototype applications in lattice QCD using MATLAB/OCTAVE language and environment. The first version of the software, QCDLAB 1.0 offers the general framework on how to achieve this goal by simulating set of the lattice Schwinger model { t http://phys.fshn.edu.al/qcdlab.html}. In this talk we introduce QCDLAB 1.1, which extends QCDLAB 1.0 capabilities for real world lattice computations with Wilson and Domain Wall fermions.

연구 동기 및 목표

• 도메인 월 페르미온과 같은 카이랄 페르미온 형식에 대해 격자 QCD 시뮬레이션의 느린 프로토타이핑 주기를 해결하기 위해.
• 격자 알고리즘 테스트를 위한 고수준이고 접근 가능한 환경을 제공함으로써 연구자들이 진입 장벽을 낮추기 위해.
• 계산 비용이 많이 들고 수렴에 민감한 5차원 도메인 월 페르미온 행렬을 위한 효율적인 해법기의 개발 및 벤치마킹을 가속화하기 위해.
• 고성능 생산 코드로의 이식 전에 알고리즘 아이디어를 신속하게 테스트할 수 있도록 하기 위해.
• 최소한의, 빠르고 상호작용 가능한 계산 환경을 통해 새로운 페르미온 형식과 해법기의 탐색을 지원하기 위해.

제안 방법

• QCDLAB 1.1은 MATLAB/OCTAVE 기반의 고수준 인터페이스를 제공하며, 그 광범위한 선형 대수 및 수치 계산 기능을 활용한다.
• 특수 함수인 Mult_Dirac_W와 Mult_DWF를 통해 디랙-윌슨 및 도메인 월 페르미온 행렬-벡터 곱셈과 같은 핵심 격자 연산을 구현한다.
• 희소 행렬 연산을 지원하며, SU(3) 연결 변수에 대해 9N×2 행렬로 최적화된 게이지 장 구성 설정을 포함한다.
• BiCGStab(bicg5)와 같은 표준 카일로프 부분공간 해법기와 헤르미트 문제에 적합한 고유값 해법기인 InversePower 및 PowerMethod를 통합한다.
• 페르미온 이론에 필수적인 감마행렬 연산(g. e., mult_gamma5)과 카이랄 투영(P5plus, P5minus)을 위한 유틸리티 함수를 포함한다.
• 코드는 최소한의 구조, 이식 가능성, 디버깅 용이성을 고려하여 설계되었으며, 생산 스케일 성능보다는 상호작용 기반 개발에 중점을 둔다.

실험 결과

연구 질문

• RQ1 MATLAB/OCTAVE와 같은 고수준 언어 환경이 도메인 월 페르미온에 대한 격자 QCD 알고리즘의 프로토타이핑에 소요되는 시간을 상당히 줄일 수 있는가?
• RQ2 QCDLAB 1.1은 알고리즘 정밀도를 유지하면서도 저수준 격자 코드의 동작을 얼마나 효과적으로 재현할 수 있는가?
• RQ3 고수준 추상화의 사용이 복잡한 페르미온 행렬 연산과 반복 해법기의 개발 및 테스트를 얼마나 향상시킬 수 있는가?
• RQ4 QCDLAB 1.1은 저수준 C/C++ 프로그래밍 없이도 5차원 도메인 월 페르미온 행렬을 해결하기 위한 새로운 알고리즘 접근법의 신속한 탐색을 지원할 수 있는가?
• RQ5 알고리즘 검증 및 수렴 테스트 측면에서 QCDLAB 1.1의 성능와 정확도는 생산 수준의 격자 코드와 비교해 어떻게 되는가?

주요 결과

• QCDLAB 1.1은 고수준 MATLAB/OCTAVE 함수를 사용하여 도메인 월 페르미온에 대한 격자 QCD 알고리즘의 빠른 프로토타이핑을 가능하게 하여, 새로운 페르미온 해법기 개발 시간을 단축시켰다.
• Mult_DWF를 통한 5차원 도메인 월 페르미온 행렬-벡터 곱셈이 성공적으로 구현되었으며, 표준 및 절삭된 형식 모두를 지원한다.
• bicg5 및 InversePower와 같은 카일로프 부분공간 해법기의 포함으로 인해 격자 양자장 이론에서 발생하는 큰 희소 선형 연립방정식의 효율적 해법이 가능해졌다.
• 감마행렬 연산, 특히 카이랄 투영 및 감마행렬 곱셈을 포함한 디랙 행렬 연산의 전 범위를 지원하여, 카이랄 페르미온 시뮬레이션에 필수적인 요소를 구현하였다.
• QCDLAB 1.1은 연구자가 알고리즘 아이디어를 일주일이 아닌 수 분 내로 테스트할 수 있도록 최소한의, 상호작용 가능하고 디버깅이 쉬운 환경을 제공한다.
• 프레임워크는 확장 가능하며, 모듈식 함수 설계 덕분에 미래에 무비우스 또는 절삭 오버랩 페르미온과 같은 고급 페르미온 형식의 통합을 지원한다.