[논문 리뷰] Generalization of the complex shifted Laplacian: on the class of expansion preconditioners for Helmholtz problems

Siegfried Cools, Wim Vanroose|arXiv (Cornell University)|2015. 01. 19.

Matrix Theory and Algorithms참고 문헌 25인용 수 3

한 줄 요약

이 논문은 헬름홀츠 문제를 위한 복소 이동 라플라시안(CSL) 전치조건자에 대한 일반화로, 역 헬름홀츠 연산자의 유한 테일러 급수 전개를 기반으로 한 전치조건자 클래스 EX(m)를 소개한다. EX(m)에 포함된 고차항은 케리로프 해법기의 수렴을 향상시키지만, 적용 비용이 증가하므로 EX(1), 즉 CSL과 동일한 경우가 실용적으로 가장 효율적인 선택이다.

ABSTRACT

In this note we propose the class of expansion preconditioners which forms a direct generalization of the CSL preconditioner. The construction of the EX(m) preconditioner is based upon a truncated Taylor series expansion of the original Helmholtz operator inverse. The EX(m) preconditioner is shown to significantly improve Krylov solver convergence rates for the Helmholtz problem for growing values of m. However, the addition of multiple terms in the expansion also increases the computational cost of applying the preconditioner. A thorough cost-benefit analysis of the addition of extra terms in the EX(m) preconditioner proves the CSL or EX(1) preconditioner to be the practically most efficient member of the EX(m) class. We suggest possible extensions to the expansion preconditioner class which are particularly interesting forma theoretical viewpoint.

연구 동기 및 목표

• 헬름홀츠 문제를 위한 복소 이동 라플라시안(CSL) 전치조건자를 더 넓은 전치조건자 클래스로 일반화하는 것.
• 전치조건자 전개에 포함된 항의 수를 늘릴 경우 수렴성 향상과 계산 비용 간의 상충 관계를 분석하는 것.
• 헬름홀츠 문제를 해결하는 데 있어서 실용적 효율성을 기준으로 EX(m) 클래스의 최적 성능을 갖는 성분을 규명하는 것.
• 향후 연구를 위한 전개 전치조건자 프레임워크의 이론적 확장 가능성을 탐색하는 것.

제안 방법

• 역 헬름홀츠 연산자의 유한 테일러 급수 전개를 통해 EX(m) 전치조건자를 구성하는 것.
• 역 연산자의 전개에서 유도된 유한 급수 형태로 전치조건자를 공식적으로 정의하는 것.
• 케리로프 부분공간 해법기에 EX(m) 전치조건자를 적용하여 헬름홀츠 문제의 수렴 속도를 가속화하는 것.
• 다양한 m 값에 대해 해법기의 수렴 속도와 전치조건자 적용 비용을 비교하는 비용-편익 분석을 수행하는 것.
• CSL 전치조건자를 기준선으로 삼으며, 이는 EX(m) 가족에서 m=1에 해당한다.
• 향후 개발을 위한 전개 전치조건자 클래스의 이론적 확장 가능성을 탐색하는 것.

실험 결과

연구 질문

• RQ1 유한 테일러 급수 전개의 항 수를 늘일수록 헬름홀츠 문제에 대한 케리로프 해법기의 수렴 속도에 어떤 영향을 미치는가?
• RQ2 고차항을 포함한 더 높은 순서의 EX(m) 전치조건자를 사용할 경우, 수렴성 향상과 계산 비용 증가 사이의 상충 관계는 어떻게 되는가?
• RQ3 실제로 EX(m) 전치조건자의 효율성을 최대화하는 최적의 m 값이 존재하는가?
• RQ4 EX(m) 전치조건자는 기존의 CSL 전치조건자를 어떻게 일반화하는가?
• RQ5 전개 전치조건자 프레임워크의 적용 범위를 향상시키기 위해 어떤 이론적 확장이 가능할 수 있는가?

주요 결과

• EX(m) 전치조건자 가족은 역 헬름홀츠 연산자의 테일러 급수 전개에서 유도된 고차항을 포함함으로써 CSL 전치조건자를 일반화한다.
• m의 값이 증가할수록 더 나은 전치조건 효과로 인해 헬름홀츠 문제에 대한 케리로프 해법기의 수렴 속도가 빨라진다.
• 그러나 전치조건자 적용 비용은 m이 증가함에 따라 증가하며, 특히 다수의 항을 평가해야 하는 데 기인한다.
• 수렴성 향상에도 불구하고 고차항을 포함한 더 높은 순서의 EX(m) 전치조건자는 실용적 효율성이 떨어진다.
• EX(1) 전치조건자, 즉 CSL 방법과 동일한 경우는 실용적으로 EX(m) 클래스에서 가장 효율적인 성분으로 확인되었다.
• 본 연구는 향후 연구에 유용할 수 있는 전개 전치조건자 클래스의 이론적 확장 가능성을 제안한다.