[논문 리뷰] A High-Performance Numerical Library for Solving Eigenvalue Problems: FEAST Solver v2.0 User's Guide

Eric Polizzi|arXiv (Cornell University)|2012. 03. 19.

Matrix Theory and Algorithms인용 수 8

한 줄 요약

FEAST v2.0는 복소 평면 내 지정된 구역 또는 둘레 내에서 모든 고유값과 고유벡터를 효율적으로 계산할 수 있도록 설계된 고성능 수치 라이브러리로, 에르미트 및 비에르미트 고유값 문제를 해결하는 데 사용됩니다. 이 라이브러리는 둘레 통합과 유리식 필터링을 조합한 새로운 알고리즘적 접근을 활용하여 공유 메모리 및 분산 메모리 아키텍처에서 견고하고 확장 가능한 성능을 제공합니다.

ABSTRACT

The FEAST eigensolver package is a free high-performance numerical library for solving the Hermitian and non-Hermitian eigenvalue problems, and obtaining all the eigenvalues and (right/left) eigenvectors within a given search interval or arbitrary contour in the complex plane. Its originality lies with a new transformative numerical approach to the traditional eigenvalue algorithm design - the FEAST algorithm. The FEAST eigensolver combines simplicity and efficiency and it offers many important capabilities for achieving high performance, robustness, accuracy, and scalability on parallel architectures. FEAST is both a comprehensive library package, and an easy to use software. It includes flexible reverse communication interfaces and ready to use predefined interfaces for dense, banded and sparse systems. The current version v3.0 of the FEAST package can address both Hermitian and non-Hermitian eigenvalue problems (real symmetric, real non-symmetric, complex Hermitian, complex symmetric, or complex general systems) on both shared-memory and distributed memory architectures (i.e contains both FEAST-SMP and FEAST-MPI packages). This User's guide provides instructions for installation setup, a detailed description of the FEAST interfaces and a large number of examples.

연구 동기 및 목표

• 다양한 아키텍처에서 대규모 고유값 문제를 효율적으로 해결하는 데 발생하는 계산적 과제를 해결하기 위해.
• 특정 구간 또는 둘레를 대상으로 설정할 수 있는 확장 가능하고 견고한 고유값 해법 제공를 위해.
• 실대칭, 복소 에르미트, 일반적인 복소 시스템을 포함한 에르미트 및 비에르미트 시스템 모두를 지원하기 위해.
• 공유 메모리(SMP) 및 분산 메모리(MPI) 시스템에서 최적화된 병렬 처리를 통해 고성능을 제공하기 위해.
• 과학 계산 응용 프로그램에의 통합을 용이하게 하기 위해 역방향 통신 및 사전 정의된 인터페이스를 포함한 유연한 인터페이스 제공을 위해.

제안 방법

• FEAST 알고리즘은 사용자가 정의한 영역 내 고유값을 분리하기 위해 복소 평면에서의 둘레 통합을 사용합니다.
• 수렴성과 정확도를 향상시키기 위해 고유값 문제를 축소된 크리로프 부분공간에 투영하기 위해 유리식 필터링을 적용합니다.
• 맞춤형 인터페이스를 통해 밀도, 밴드, 희소 행렬 형식을 모두 지원합니다.
• 다중 프로세서에 걸친 분할 및 정복 전략을 사용하여 확장성 향상을 도모합니다.
• 기존의 해법기나 응용 프로그램에 유연하게 통합할 수 있도록 역방향 통신 인터페이스를 제공합니다.
• 공유 메모리 및 분산 메모리 병렬 컴퓨팅 환경을 지원하기 위해 FEAST-SMP 및 FEAST-MPI 버전이 포함되어 있습니다.

실험 결과

연구 질문

• RQ1 어떻게 복소 평면 내 특정 구간 또는 둘레 내에서 고유값 문제를 효율적으로 해결할 수 있는가?
• RQ2 어떤 알고리즘적 접근이 에르미트 및 비에르미트 고유값 문제 모두에 대해 고성능과 견고성을 보장하는가?
• RQ3 FEAST 알고리즘이 공유 메모리 및 분산 메모리 아키텍처에서 효과적으로 병렬화될 수 있는가?
• RQ4 어떤 인터페이스 설계 선택이 과학 계산 워크플로우에서의 사용성과 통합의灵活性를 극대화하는가?
• RQ5 FEAST 방법이 다양한 행렬 유형과 시스템 구성에서 정확성과 확장성을 얼마나 잘 유지하는가?

주요 결과

• FEAST v2.0은 실대칭, 복소 에르미트, 일반적인 복소 시스템을 포함한 다양한 행렬 유형에서 에르미트 및 비에르미트 고유값 문제를 성공적으로 해결합니다.
• FEAST-SMP와 FEAST-MPI 아키텍처에서 모두 높은 성능과 확장성을 입증합니다.
• 둘레 기반 유리식 필터링을 통해 사용자가 지정한 관심 영역 내 고유값과 고유벡터를 정확하게 추출할 수 있습니다.
• 역방향 통신 및 사전 정의된 인터페이스의 유연성 덕분에 기존 과학 소프트웨어 스택에 간편하게 통합할 수 있습니다.
• 투영 기반 접근 덕분에 악조건이거나 대규모 문제에서도 견고성과 정확성이 유지됩니다.
• 빠른 구현을 위한 광범위한 문서화와 예제 코드를 포함한 포괄적이고 사용자 友好的 인터페이스를 제공합니다.