[논문 리뷰] Zolotarev Quadrature Rules and Load Balancing for the FEAST Eigensolver

Stefan Guettel, Eric Polizzi외 2명|arXiv (Cornell University)|2014. 07. 30.

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

한 줄 요약

이 논문은 FEAST 고유값 해법기에서 전통적인 가우스 및 트라프레지오드 적분 규칙보다 우수한 대안으로 조로타레프 적분 규칙을 제안하며, 병렬 계산에서 수렴 속도와 로드 밸런싱을 크게 향상시킨다. 부호 함수에 대한 최적의 유리 근사치를 활용함으로써, 고유값이 잘 분리되어 있지 않거나 검색 공간 차원이 작을 경우에도 더 빠르고 예측 가능한 수렴을 달성한다.

ABSTRACT

The FEAST method for solving large sparse eigenproblems is equivalent to subspace iteration with an approximate spectral projector and implicit orthogonalization. This relation allows to characterize the convergence of this method in terms of the error of a certain rational approximant to an indicator function. We propose improved rational approximants leading to FEAST variants with faster convergence, in particular, when using rational approximants based on the work of Zolotarev. Numerical experiments demonstrate the possible computational savings especially for pencils whose eigenvalues are not well separated and when the dimension of the search space is only slightly larger than the number of wanted eigenvalues. The new approach improves both convergence robustness and load balancing when FEAST runs on multiple search intervals in parallel.

연구 동기 및 목표

• 대규모 희소 일반화 고유값 문제를 위한 FEAST 고유값 해법기의 수렴 안정성과 효율성을 향상시키기.
• FEAST이 병렬적으로 다수의 검색 간격에서 실행될 때 발생하는 로드 불균형 문제를 해결하기.
• 기존의 적분 규칙(가우스, 트라프레지오드)을 대체하여 더 나은 스펙트럼 필터링을 위한 조로타레프 기반 유리 근사치를 도입하기.
• 최적의 유리 근사치를 통해 이론적으로 탄탄하고 분석적으로 예측 가능한 수렴 인자를 FEAST에 제공하기.
• 특히 고유값이 잘 분리되어 있지 않은 경우에도 높은 정확도를 유지하면서 검색 공간 차원을 줄일 수 있도록 하기.

제안 방법

• FEAST 방법을 유리 행렬 함수 $ r_m(M) $를 사용한 암묵적 직교화 단계를 가진 부분공간 반복과 동치로 기술하며, 여기서 $ M = B^{-1}A $이다.
• 고유값이 원하는 영역에 투영되는 지표 함수 $ f(z) = 1/(2\bar{\rho}i) \int_{\Gamma} (\gamma - z)^{-1} d\gamma $에 대한 유리 함수 $ r_m(z) $의 근사 오차를 통해 FEAST의 수렴을 모델링한다.
• 부호 함수에 대한 최적의 근사치로 알려진 조로타레프 유리 함수를 도입하며, 이는 등분산성과 간격 경계에서 급격한 전이 특성을 지닌다.
• 원하는 고유값 영역와 원하지 않는 고유값 영역에서 이상적인 스텝 함수로부터 최대 편차를 최소화하는 조로타레프 이론 기반의 유리 필터를 구성한다.
• 기존의 표준 적분 기반 유리 함수 대신 조로타레프 기반 유리 함수를 사용하여 FEAST 알고리즘에 이러한 필터를 적용함으로써 스펙트럼 투영 정확도를 향상시킨다.
• 조로타레프 필터는 예측 가능한 분석적 수렴 인자를 제공하며, 무한대에서 감쇠되지 않는 특성과 날카운 전이 특성 덕분에 필요한 검색 공간 차원을 줄일 수 있음을 보여준다.

실험 결과

연구 질문

• RQ1 FEAST 고유값 해법기의 수렴 행동은 스펙트럼 프로젝터에 사용된 유리 근사치의 품질에 어떻게 의존하는가?
• RQ2 필터 설계에서 최적성으로 유명한 조로타레프 유리 함수는 기존의 가우스 및 트라프레지오드 적분 규칙보다 FEAST에서 더 빠르고 강건한 수렴을 이끌 수 있는가?
• RQ3 FEAST가 병렬적으로 다수의 검색 간격에 적용될 때, 조로타레프 기반 필터는 로드 밸런싱을 어느 정도 향상시키는가?
• RQ4 조로타레프 필터의 무한대에서 감쇠되지 않는 행동이 FEAST 내 선형 시스템 해법의 안정성과 정확성에 어떤 영향을 미치는가?
• RQ5 조로타레프 적분 규칙의 극점과 가중치는 실질적인 FEAST 구현에서 수치 조건과 수렴에 어떻게 영향을 미치는가?

주요 결과

• 조로타레프 기반 FEAST 변종은 고유값이 잘 분리되어 있지 않거나 검색 공간이 원하는 고유값 수와 약간만 큰 경우 가우스 및 트라프레지오드 적분 규칙보다 더 빠른 수렴을 달성한다.
• Caffeinep2 행렬에 대해, $ m=8 $일 때, 구간 $[-711, -4]$에서 가우스 방법은 39회의 반복을 필요로 했지만 조로타레프 방법은 8회로 줄었고, $[-4, 1.995]$에서는 5회에서 9회로 변화하여 더 나은 로드 밸런싱을 보였다.
• $ m=32 $일 때, 첫 번째 간격에서는 조로타레프 방법이 3회, 가우스 방법이 5회 반복을 필요로 했고, 두 번째 간격에서는 각각 2회와 3회로 나타나 일관되고 균형 잡힌 수렴을 보였다.
• 조로타레프 필터의 간격 경계에서 급격한 전이 특성 덕분에 더 작은 검색 공간 차원으로도 높은 정확도를 유지할 수 있어 계산 비용을 감소시켰다.
• 조로타레프 기반 방법의 수렴 인자는 분석적으로 예측 가능하며, 다양한 문제 사례에서 안정적으로 유지되며, 가우스 적분 규칙과 달리 기복적인 수렴 행동을 보이지 않는다.
• 수치 실험 결과, 극점들이 실수축에 가까워져도 안정성 문제 없이 작동했으며, 이는 극점의 허수부에 비례하는 가중치의 비율이 선형 시스템 해법의 오차를 감쇠시키기 때문인 것으로 나타났다.