[논문 리뷰] SYM-ILDL: Incomplete $LDL^{T}$ Factorization of Symmetric Indefinite and Skew-Symmetric Matrices

Chen Greif, Shiwen He외 1명|arXiv (Cornell University)|2015. 05. 28.

Matrix Theory and Algorithms참고 문헌 15인용 수 6

한 줄 요약

SYM-ILDL는 대칭 부정정 및 실수 반대칭 행렬의 완전하지 않은 LDLᵀ 분해를 위한 새로운 C++ 템플릿 소프트웨어 패키지를 소개한다. Crout 기반 알고리즘을 사용하며, 고급 피벗 전략(Bunch-Kaufman, 룩)과 재정렬(AMD, RCM)을 통합하였다. 반복적 해법기에서 뛰어난 수렴성을 보이며, 특히 룩 피벗 전략을 사용할 경우 수렴성이 뛰어나며, ILUTP 대비 약 2/3의 메모리만을 사용하여 대류-확산 문제에서 뚜렷한 성능 향상을 보였다.

ABSTRACT

SYM-ILDL is a numerical software package that computes incomplete $LDL^{T}$ (or `ILDL') factorizations of symmetric indefinite and real skew-symmetric matrices. The core of the algorithm is a Crout variant of incomplete LU (ILU), originally introduced and implemented for symmetric matrices by [Li and Saad, Crout versions of ILU factorization with pivoting for sparse symmetric matrices, Transactions on Numerical Analysis 20, pp. 75--85, 2005]. Our code is economical in terms of storage and it deals with real skew-symmetric matrices as well, in addition to symmetric ones. The package is written in C++ and it is templated, open source, and includes a MATLAB interface. The code includes built-in RCM and AMD reordering, two equilibration strategies, threshold Bunch-Kaufman pivoting and rook pivoting, as well as a wrapper to MC64, a popular matching based equilibration and reordering algorithm. We also include two built-in iterative solvers: SQMR preconditioned with ILDL, or MINRES preconditioned with a symmetric positive definite preconditioner based on the ILDL factorization.

연구 동기 및 목표

• PDE 이산화 및 안장점 문제에서 흔한 대칭 부정정 및 실수 반대칭 행렬을 대상으로 안정적이고 메모리 효율적인 완전하지 않은 LDLᵀ 분해 방법을 개발한다.
• 기존의 ILU 기반 Crout 알고리즘을 확장하여 부정정 및 반대칭 시스템을 대칭을 유지하는 피벗 전략으로 처리할 수 있도록 한다.
• 다양한 재정렬(AMD, RCM), 균형 조정(1-, 2-, ∞-노름), 피벗 전략(Bunch-Kaufman, 룩) 기법을 통합하여 수치적 안정성과 수렴성을 향상시킨다.
• MINRES에 대한 양의 정정 조건부 개선자(정규화된 ILDL에서 유도된)와 ILDL에 직접적으로 조건부 개선된 SQMR를 포함한 두 가지 내장된 조건부 개선 반복 해법기를 제공한다.
• 고성능 MEX 인터페이스를 갖춘 완전히 템플릿화되고 확장 가능하며 오픈소스인 C++ 라이브러리를 제공하여 과학 계산 분야에서 넓은 활용성을 확보한다.

제안 방법

• 대칭 및 반대칭 행렬에 적합한 왼쪽 봉우리 기반 Crout 변형을 사용하여, 분해 과정에서 대칭성 $PAP^T = LDL^T$ 를 유지한다.
• 임계값 기반 요소 제거 및 필로우 인 제어를 통해 희소성 유지와 함께 수치적 안정성 확보.
• 부정정 및 반대칭 시스템을 처리하기 위해 Bunch-Kaufman 및 룩 피벗 전략을 적용하여 붕괴를 방지하고 조건수 개선.
• AMD, RCM 등의 재정렬 기법과 최대노름, 1-노름, 2-노름 기반 균형 조정을 통합하여 필로우 인 감소 및 수렴성 향상.
• HSL_MC64를 래핑하여 고급 매칭 기반 재정렬 및 균형 조정을 통합하여 어려운 행렬에 대한 강건성 향상.
• 반대칭 행렬의 경우 $2\times2$ 블록 구조를 활용하여 MINRES용 양의 정정 조건부 개선자를 구성하고, 대칭 행렬의 경우 Gill 등(1992)의 기법을 응용한다.

실험 결과

연구 질문

• RQ1 실수 반대칭 행렬은 대각성분이 0이며 $2\times2$ 피벗이 필요하므로, Crout 기반 완전하지 않은 LDLᵀ 분해가 효과적으로 확장될 수 있는가?
• RQ2 대류-확산 문제에서 SYM-ILDL의 성능는 ILUTP 및 기타 최신 조건부 개선자(예: Scott & Tuma, 2014)와 비교해 수렴성 및 메모리 사용 측면에서 어떻게 나타나는가?
• RQ3 다양한 피벗 전략(룩 대비 부분 피벗)과 재정렬(AMD 대비 RCM)이 조건부 개선 반복 해법기의 수렴성에 미치는 영향은 무엇인가?
• RQ4 균형 조정은 반대칭 행렬의 ILDL에 대해 조건수 개선 및 수렴성 향상에 기여하는가, 아니면 삼중대각 블록의 구조적 지배성으로 인해 성능 저하를 초래하는가?
• RQ5 자기적으로 부정정 또는 반대칭인 ILDL 분해에서 유도된 양의 정정 조건부 개선자가 MINRES에 대해 얼마나 효과적인가?

주요 결과

• SYM-ILDL는 ILUTP와 유사한 필로우 인을 유지하면서도 GMRES 수렴성이 뛰어나다: $n=125000$일 때, ILDL-룩은 6회 반복만으로 수렴했지만 ILUTP는 58회가 소요되었으며, 메모리 사용은 약 2/3 수준이었다.
• 룩 피벗 전략은 부분 피벗 전략보다 일관되게 뛰어난 성능을 보이며, 대규모 대류-확산 문제에서 반복 횟수를 최대 80%까지 감소시켰다.
• 반대칭 행렬의 경우, 삼중대각 블록의 지배적 특성으로 인해 균형 조정이 효과적이지 않거나 오히려 해로운 것으로 나타났다. 따라서 최종 해법기 설정에서는 균형 조정을 비활성화하였다.
• 시스템의 조건부 고유값은 1 주변으로 강하게 군집되어 있었으며, 이는 특히 룩 피벗 전략을 사용할 경우 반복 해법기의 빠른 수렴을 의미한다.
• $70^3$ 대류-확산 행렬에서 ILDL-룩은 106.81초 만에 5회 반복으로 수렴했고, ILUTP는 876초 이내에 수렴하지 못했다. 이는 SYM-ILDL의 강건성을 입증한다.
• MEX 인터페이스는 원본 C++ 버전과 동등한 성능를 제공하여 MATLAB 워크플로우에 손실 없이 원활하게 통합 가능하다.