[논문 리뷰] A preconditioner based on the shift-splitting method for generalized saddle point problems

Davod Khojasteh Salkuyeh, Mohsen Masoudi외 1명|arXiv (Cornell University)|2015. 06. 15.

Matrix Theory and Algorithms참고 문헌 4인용 수 4

한 줄 요약

이 논문은 비대칭 정부호(1,1)-블록과 대칭 준정부호(2,2)-블록을 가진 일반화된 매끄러운 점 문제를 위한 시프트-스플릿팅 조절자를 제안한다. 이 방법은 조건부 수렴하지 않는 반복적 방법에 기반하며, 나비에-스토크스 문제에 대한 수치 실험 결과, 비조절된 GMRES에 비해 GMRES 반복 횟수와 CPU 시간을 크게 감소시킨다.

ABSTRACT

In this paper, we propose a preconditioner based on the shift-splitting method for generalized saddle point problems with nonsymmetric positive definite (1,1)-block and symmetric positive semidefinite $(2,2)$-block. The proposed preconditioner is obtained from an basic iterative method which is unconditionally convergent. We also present a relaxed version of the proposed method. Some numerical experiments are presented to show the effectiveness of the method.

연구 동기 및 목표

• 유한요소 이산화에 의해 유도되는 큰 희소 일반화된 매끄러운 점 문제를 위한 효과적인 조절자를 개발하기 위해.
• 기존에 대칭 블록에 국한된 시프트-스플릿팅 방법을 (1,1)-블록이 비대칭 정부호인 경우로 확장하기 위해.
• 제안된 반복 방법의 조건부 수렴 없이 수렴성을 보장하고, GMRES와 같은 카일로프 부분공간 해법기용으로 안정적인 조절자를 유도하기 위해.
• GMRES(30)를 사용하여 벤치마크 나비에-스토크스 문제에서 조절자의 성능을 평가하기 위해.

제안 방법

• 방법은 (1,1)-블록에 대해 이동 파rameter α와 β를 기반으로 매끄러운 점 시스템을 M 및 N 행렬로 분할한다.
• 모든 α, β > 0에 대해 반복 행렬 Γα,β = M⁻¹N의 스펙트럼 반경이 1 미만임을 증명하여 조건부 수렴 없이 수렴함을 보장한다.
• 조절자 P_MGSS는 Mα,β = ½[(αI + A) Bᵀ; -B (βI + C)]로 정의되며, 이는 대칭 정부호 행렬이면서 GMRES에 적합하다.
• 완화된 형태인 P_RMGSS는 (1,1)-블록의 이동을 생략하여 구현을 단순화하면서도 효과를 유지한다.
• 내부 시스템은 잔차 감소 허용오차 10⁻², 최대 40회 반복으로, GMRES(10)을 사용하여 풀어낸다.
• 방법은 정적 상태 나비에-스토크스 방정식의 Q1-P0 유한요소 이산화에 대해 IFISS 패키지를 사용하여 테스트되었다.

실험 결과

연구 질문

• RQ1 비대칭 정부호(1,1)-블록을 가진 일반화된 매끄러운 점 문제에 대해 시프트-스플릿팅 방법을 확장할 수 있는가?
• RQ2 제안된 조절자가 이러한 문제에 대해 조건부 수렴 없이 수렴성을 보장하는가?
• RQ3 나비에-스토크스 문제에서 GMRES 반복 횟수와 CPU 시간을 줄이는 데 조절자가 얼마나 효과적인가?
• RQ4 α = 0.01 및 β = 0.001의 파라미터 선택이 수렴 행동에 어떤 영향을 미치는가?

주요 결과

• 제안된 MGSS 방법은 조건부 수렴 없이 수렴하며, 모든 α, β > 0에 대해 반복 행렬의 스펙트럼 반경이 1보다 작다.
• 64×64 격자에서, 조절된 GMRES(30)는 28회의 반복과 21.48초가 소요되었고, 비조절된 경우는 3554회의 반복과 110.1초가 소요되었다.
• 32×32 격자에서, 조절자는 반복 횟수를 비조절된 경우 608회에서 12회로 줄였으며, CPU 시간도 4.95초에서 3.58초로 단축시켰다.
• 완화된 조절자 RMGSS 역시 뛰어난 성능을 보였으며, 16×16 격자에서는 반복 횟수를 8회로 줄였고, 32×32 격자에서는 12회로 줄였다.
• 모든 격자 크기에서 비조절된 GMRES에 비해 조절자가 뚜렷이 뛰어난 성능을 보였으며, 특히 더 세밀한 격자에서 두드러졌다.
• 이론적 분석은 A가 비대칭 정부호일 경우에도 방법이 여전히 수렴함을 확인하였으며, 이는 이전 연구가 대칭 A에 국한된 바를 확장한 것이다.