[논문 리뷰] A FEAST Algorithm with oblique projection for generalized eigenvalue problems

Guojian Yin, Raymond H. Chan외 1명|arXiv (Cornell University)|2014. 04. 07.

Matrix Theory and Algorithms참고 문헌 27인용 수 7

한 줄 요약

이 논문은 비헤르미트형 일반화된 고유값 문제에 대해 기존 FEAST 알고리즘을 확장하여 기울기 투영과 둘레 적분을 사용함으로써 복소 평면 내 사용자가 정의한 영역 내에 있는 고유값을 안정적이고 정확하게 계산할 수 있도록 한다. 이 방법은 특히 악조건이거나 비헤르미트인 경우에 전통적인 솔버보다 정확도와 효율성 면에서 뛰어나다.

ABSTRACT

The contour-integral based eigensolvers are the recent efforts for computing the eigenvalues inside a given region in the complex plane. The best-known members are the Sakurai-Sugiura (SS) method, its stable version CIRR, and the FEAST algorithm. An attractive computational advantage of these methods is that they are easily parallelizable. The FEAST algorithm was developed for the generalized Hermitian eigenvalue problems. It is stable and accurate. However, it may fail when applied to non-Hermitian problems. In this paper, we extend the FEAST algorithm to non-Hermitian problems. The approach can be summarized as follows: (i) to construct a particular contour integral to form a subspace containing the desired eigenspace, and (ii) to use the oblique projection technique to extract desired eigenpairs with appropriately chosen test subspace. The related mathematical framework is established. We also address some implementation issues such as how to choose a suitable starting matrix and design good stopping criteria. Numerical experiments are provided to illustrate that our method is stable and efficient.

연구 동기 및 목표

• 기존 FEAST 알고리즘이 비헤르미트형 일반화된 고유값 문제에서 실패하는 문제를 해결하기 위해.
• 사용자가 지정한 영역 내 고유값을 계산하기 위한 안정적이고 정확한 고유값 해법을 개발하기 위해.
• 비헤르미트 시스템에서 둘레 적분 기반 고유값 해법을 위한 수학적 프레임워크를 구축하기 위해.
• 모든 원하는 고유값이 포착됨을 보장하는 정지 기준을 제공하기 위해.
• 특정 고유값 계산에서 CIRR 및 MATLAB의 eig와 같은 기존 방법보다 수치적 성능을 향상시키기 위해.

제안 방법

• 주어진 영역 내에 있는 고유공간을 포함하는 부분공간을 형성하기 위해 둘레 적분을 구성한다.
• 원래 FEAST에서 사용하는 레일리-리츠 절차를 대체하기 위해 신중히 선택한 시험 부분공간을 사용한 기울기 투영을 적용하여 원하는 고유쌍을 추출한다.
• 추출된 부분공간 내에서 일반화된 고유값 해법을 사용하여 원하는 고유쌍을 정확하게 계산한다.
• 알고리즘 매개변수를 안내하기 위해 둘레 내부에 있는 고유값의 수에 대한 상한을 추정하는 방법을 유도한다.
• 수렴성과 완전성을 보장하기 위해 상대 잔여 norm 기반의 적응형 정지 기준을 구현한다.
• 수렴성과 안정성을 향상시키기 위해 적절한 시작 행렬을 선택하는 전략을 도입한다.

실험 결과

연구 질문

• RQ1 기울기 투영과 둘레 적분을 사용하여 FEAST 알고리즘이 비헤르미트형 일반화된 고유값 문제에 일반화될 수 있는가, 안정성과 정확성을 유지할 수 있는가?
• RQ2 비헤르미트 시스템에서 레일리-리츠 절차에 비해 기울기 투영이 고유쌍 추출을 어떻게 향상시킬 수 있는가?
• RQ3 모든 둘레 내 고유값이 포착됨을 보장하는 효과적인 정지 기준은 무엇인가?
• RQ4 정확도와 계산 비용 측면에서 블록_CIRR 및 MATLAB의 eig와 비교할 때 이 방법의 성능는 어떠한가?
• RQ5 강건한 성능을 위해 고려해야 할 주요 구현 고려사항은 무엇인가 (예: 시작 행렬 선택 및 조건 수치 등)?

주요 결과

• 제안된 방법은 비헤르미트형 일반화된 고유값 문제에 대해 주어진 둔환 내 모든 고유값을 성공적으로 계산한다.
• 수치 실험 결과, MATLAB의 eig 함수보다 높은 정확도를 확보하였으며, 소형에서 중형 문제에서는 상대 오차가 약 10^-15 수준이다.
• 대규모 문제(예: 크기 4800)에서는 영역 내 169개의 고유값을 상대 오차 5.10×10^-8로 계산하여 강건성을 입증하였다.
• 특정 고유값만 필요할 경우, eig보다 훨씬 빠르며, 전체 스펙트럼을 계산하는 eig조차도 성능 향상을 보였다.
• Block_CIRR와 비교해 본 결과, 특히 악조건이거나 특이 행렬 B를 가진 경우 더 정확하고 안정적이지만, 반복마다 오른쪽 항의 수가 많아져 CPU 시간이 더 소요된다.
• 알고리즘은 쉽게 병렬 처리가 가능하며, 모든 테스트 문제에서 10회 이내의 반복 수로 수렴함을 보여, 빠른 수렴성을 입증하였다.