[논문 리뷰] Deflated BiCGStab for linear equations in QCD problems

Abdou Abdel-Rehim, Ronald B. Morgan외 1명|ArXiv.org|2007. 10. 10.

Matrix Theory and Algorithms참고 문헌 2인용 수 5

한 줄 요약

이 논문은 라티스 양성동역학(QCD)에서 다중 우변 및 다중 이동을 포함한 대규모 복소선형 연립방정식을 해결하기 위해 Deflated BiCGStab(D-BiCGStab)을 제안한다. Wilson 페르미온에서 γ₅ 대칭성을 통해 오른쪽 고유벡터로부터 유도된 왼쪽 고유벡터를 초기 추측에 투영함으로써 D-BiCGStab는 수렴 속도를 가속화한다. 이론적 분석과 실험 결과에 따르면, 비근접 κ 문제에서는 행렬-벡터 곱셈 횟수를 약 5배 감소시키며, 다중 질량을 가진 비틀림-질량 QCD에서는 20~50% 감소시킨다.

ABSTRACT

The large systems of complex linear equations that are generated in QCD problems often have multiple right-hand sides (for multiple sources) and multiple shifts (for multiple masses). Deflated GMRES methods have previously been developed for solving multiple right-hand sides. Eigenvectors are generated during solution of the first right-hand side and used to speed up convergence for the other right-hand sides. Here we discuss deflating non-restarted methods such as BiCGStab. For effective deflation, both left and right eigenvectors are needed. Fortunately, with the Wilson matrix, left eigenvectors can be derived from the right eigenvectors. We demonstrate for difficult problems with kappa near kappa_c that deflating eigenvalues can significantly improve BiCGStab. We also will look at improving solution of twisted mass problems with multiple shifts. Projecting over previous solutions is an easy way to reduce the work needed.

연구 동기 및 목표

• 라티스 QCD에서 다중 우변이 있는 대규모 복소선형 연립방정식을 풀이할 때 표준 BiCGStab의 느린 수렴 문제를 해결하기 위해.
• 이전에 재시작되지 않은 GMRES에 사용된 탈출 기법을 재시작되지 않은 방법인 BiCGStab에 확장하여 효율성을 향상시키기 위해.
• 다중 이동과 짝짓기 전처리를 포함한 비틀림-질량 페르미온 문제의 수렴을 향상시키기 위해.
• 특히 Wilson 페르미온 체계에서 왼쪽 고유벡터가 오른쪽 고유벡터로부터 유도될 수 있음을 활용하여, 양쪽 고유벡터를 효과적으로 탈출에 활용하기 위해.
• 이전에 계산된 해를 기반으로 최소 잔여항 투영을 통해 초기 추측을 개선하여 계산 비용을 절감하기 위해.

제안 방법

• 첫 번째 우변을 해결하기 위해 GMRES-DR을 사용하고, 작은 고유값에 해당하는 근사 오른쪽 고유벡터를 추출한다.
• Wilson Dirac 연산자에 대해 성립하는 γ₅ 대칭성 성질 γ₅A_Wγ₅ = A_W†를 이용해 오른쪽 고유벡터에서 왼쪽 고유벡터를 유도한다.
• 정규직교화된 왼쪽 및 오른쪽 고유벡터 행렬 U와 V를 사용해 k차원의 탈출 부분공간을 구성한다.
• 작은 k×k 선형계 U†AVy = U†r₀를 풀어 개선된 초기 추측 x₀⁺ = x₀ + Vy를 위한 계수 y를 계산한다.
• 탈출된 초기 추측을 사용해 표준 BiCGStab를 적용하여 이후의 우변을 효율적으로 해결한다.
• 비틀림-질량 문제의 경우, 각 새로운 이동에 대해 이전에 계산된 해를 기반으로 최소 잔여항 투영을 수행해 초기 추측을 개선한다.

실험 결과

연구 질문

• RQ1 비재시작 Krylov 해법기인 BiCGStab에 대해 다중 우변 문제에 탈출을 효과적으로 적용할 수 있는가?
• RQ2 오직 오른쪽 고유벡터를 사용하는 것보다 양쪽 고유벡터를 모두 사용할 경우 수렴 속도 향상이著명한가?
• RQ3 GMRES-DR에서 유도된 고유벡터의 정확도가 near-critical κ QCD 문제에서 D-BiCGStab의 성능에 미치는 영향은 무엇인가?
• RQ4 이전 해를 기반으로 한 투영을 통한 탈출이 다중 이동을 가진 비틀림-질량 페르미온 문제에서 행렬-벡터 곱셈 횟수를 줄일 수 있는가?
• RQ5 이동된 시스템을 해결하는 순서(가장 무거운 질량부터 또는 그 반대)가 투영 기반 초기 추측 개선에 영향을 미치는가?

주요 결과

• 16⁴ 및 20³×32 격자에서 near-critical κ 문제에 대해 D-BiCGStab는 매우 정확한 고유벡터(L3 케이스)를 사용할 경우 행렬-벡터 곱셈 횟수를 약 5배 감소시켰다.
• 고유벡터의 정확도가 높을수록 수렴 성능 향상이 뚜렷했으며, 이는 첫 번째 우변의 수렴 잔여항 노름(10⁻¹⁴, L3 케이스)으로 측정되었다.
• 11개의 질량 값이 있는 비틀림-질량 페르미온 문제에서 이전 해를 기반으로 한 투영은 총 행렬-벡터 곱셈 횟수를 20~50% 감소시켰다. 이는 8,540회에서 5,180~5,870회로 감소한 것이다. 해결 순서에 따라 다소 변동이 있었지만, 유사한 감소율을 보였다.
• 가장 무거운 질량 또는 가장 가벼운 질량을 먼저 해결하는지 여부에 관계없이 투영 방법은 효과적이었으며, 두 경우 모두 비슷한 감소율을 기록했다.
• Wilson 페르미온 체계에서는 왼쪽 고유벡터를 오른쪽 고유벡터에 γ₅ 곱하기만 하면 되므로 추가 작업이 필요하지 않았다.
• κ_c 근처의 어려운 문제에서는 오직 오른쪽 고유벡터만을 사용하는 탈출 기법으로는 작은 고유값 방향의 성분을 충분히 감소시킬 수 없었으며, 양쪽 고유벡터를 모두 사용하는 것이 필수적이었다.