[논문 리뷰] On a weighted quasi-residual minimization strategy of the QMR method for solving complex symmetric shifted linear systems

Tomohiro Sogabe, Takeo Hoshi외 2명|ArXiv.org|2009. 02. 16.

Matrix Theory and Algorithms인용 수 4

한 줄 요약

이 논문은 전자 구조 시뮬레이션에서 발생하는 대규모 복소 대칭 이동 선형계를 해결하기 위해 QMR_SYM 방법에 대한 가중 카오풍선 최소잔차법(WQMR) 전략을 제안한다. 최적화된 가중 행렬을 도입함으로써, 근사 해 갱신에 소요되는 계산 비용을 크게 감소시켜, 시드 시스템 선택이 필요 없음에도 불구하고 이동 COCG 방법과 경쟁 가능한 성능을 달성한다.

ABSTRACT

We consider the solution of complex symmetric shifted linear systems. Such systems arise in large-scale electronic structure simulations and there is a strong need for the fast solution of the systems. With the aim of solving the systems efficiently, we consider a special case of the QMR method for non-Hermitian shifted linear systems and propose its weighted quasi-minimal residual approach. A numerical algorithm, referred to as shifted QMR\_SYM($B$), is given by the choice of a particularly cost-effective weight. Numerical examples are presented to show the performance of the shifted QMR\_SYM($B$) method.

연구 동기 및 목표

• 큰 수의 이동 선형계를 해결할 때, 이동 QMR_SYM 방법에서 근사 해 갱신에 소요되는 높은 계산 비용을 해결한다.
• 수렴 실패를 방지하기 위해 시드 시스템 선택이 철저히 필요로 하는 이동 COCG 방법의 한계를 극복한다.
• 카오리브 부분공간의 이동 불변성 성질을 활용하여 여러 이동 시스템 간에 기저 벡터를 재사용한다.
• 잔차 갱신 비용을 최소화하면서 수렴 행동에 영향을 주지 않는 비용 효율적인 가중 전략을 개발한다.
• 전자 구조 계산에서 복소 대칭 이동 선형계에 대한 기존 카오리브 방법의 강력하고 효율적인 대안을 제공한다.

제안 방법

• 비에르미트 이동 시스템에 대한 QMR 방법을 복소 대칭 케이스에 특화시켜, 이동 QMR_SYM 방법을 도출한다.
• 이동 수가 많을 경우, 이동 QMR_SYM 방법에서 근사 해 갱신이 가장 시간이 오래 걸리는 연산임을 규명한다.
• 특정 가중 행렬을 도입함으로써 근사 해 갱신 비용을 줄이는 가중 카오풍선 최소잔차법(WQMR) 전략을 제안한다.
• 잔차 및 해 갱신에서 실수 산술 연산을 가능하게 하는 비용 효율적인 가중 행렬 B를 선택하여 부동소수점 연산 수를 감소시킨다.
• WQMR 전략을 QMR_SYM 프레임워크에 통합하여, 수렴 성질을 유지하는 새로운 알고리즘인 shifted QMR_SYM(B)를 구성한다.
• 카오리브 부분공간의 이동 불변성을 활용하여, 이동 COCG 방법과 달리 시드 시스템 선택이 필요 없는 방법을 확보한다.

실험 결과

연구 질문

• RQ1 큰 수의 이동 선형계를 해결할 때, 이동 QMR_SYM 방법의 주요 계산 블로킹 요소는 무엇인가?
• RQ2 가중 카오풍선 최소잔차 전략은 수렴 성능에 영향을 주지 않으면서 근사 해 갱신 비용을 줄일 수 있는가?
• RQ3 계산 시간과 강건성 측면에서 제안된 shifted QMR_SYM(B) 방법은 이동 COCG 방법과 비교해 어떻게 성능을 내는가?
• RQ4 새로운 방법은 이동 COCG 방법의 핵심적인 한계인 시드 시스템 선택이 불필요한가?
• RQ5 제안된 방법은 이동 COCG 방법과 비슷한 효율성을 달성하면서도 더 넓은 범위의 복소 대칭 이동 선형계에 적용 가능한가?

주요 결과

• 이동 QMR_SYM 방법은 근사 해 갱신에서 복소수 산술 연산을 수행하기 때문에, 이동 COCG 방법보다 반복당 계산 비용이 크게 높다.
• 이동 시스템의 수가 임계값을 초과할 경우(예: m > 200), 이동 QMR_SYM 방법에서 근사 해 갱신이 가장 오래 걸리는 부분이 된다.
• 실수 가중 전략을 사용함으로써 제안된 shifted QMR_SYM(B) 방법은 이동 COCG 방법 수준의 계산 비용으로 감소시킨다.
• m = 1501일 때, shifted QMR_SYM(B)는 이동 COCG 방법과 비교해 약 1.0 수준의 계산 시간 비율을 달성하여 경쟁 가능한 성능을 보였다.
• 이동 COCG 방법과 달리 시드 시스템 선택이 필요 없음에도 불구하고 강건성을 유지한다.
• 대규모 전자 구조 문제에 대한 수치 실험을 통해, shifted QMR_SYM(B)가 많은 수의 이동이 있는 시스템에서 실현 가능하고 효율적인 대안임을 확인했다.