[논문 리뷰] New Recurrence Relationships between Orthogonal Polynomials which Lead to New Lanczos-type Algorithms

Muhammad Farooq, Abdellah Salhi|arXiv (Cornell University)|2014. 03. 03.

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

한 줄 요약

이 논문은 선형 시스템을 해결하기 위한 새로운 랑츠 유형 알고리즘을 개발하기 위해 형식적 직교다항식(FOPs) 간의 새로운 재귀 관계를 제안한다. 전통적인 차수 차이 제약을 초월하여, 이전에 탐색되지 않은 재귀 구조—특히 $A_{12}, A_{13}, A_{16}, A_{19}, B_{13}, B_{15}, B_{16}$—를 식별한다. 이러한 구조는 수학적으로 타당하며 알고리즘 구현에 적합하여 안정성과 수렴성 향상이 가능한 새로운 랑츠 방법의 변종을 가능하게 한다.

ABSTRACT

Lanczos methods for solving $ extit{A} extbf{x}= extbf{b}$ consist in constructing a sequence of vectors $( extbf{x}_k), k=1,...$ such that $ extbf{r}_{k}= extbf{b}- extit{A} extbf{x}_{k}= extit{P}_{k}( extit{A}) extbf{r}_{0}$,, where $ extit{P}_{k}$ is the orthogonal polynomial of degree at most $k$ with respect to the linear functional $c$ defined as $c(ξ^i)=( extbf{y}, extit{A}^i extbf{r}_{0})$. Let $ extit{P}^{(1)}_{k}$ be the regular monic polynomial of degree $k$ belonging to the family of formal orthogonal polynomials (FOP) with respect to $c^{(1)}$ defined as $c^{(1)}(ξ^{i})=c(ξ^{i+1})$. All Lanczos-type algorithms are characterized by the choice of one or two recurrence relationships, one for $ extit{P}_{k}$ and one for $ extit{P}^{(1)}_{k}$. We shall study some new recurrence relations involving $ extit{P}_{k}$ and $ extit{P}^{(1)}_{k}$ and their possible combination to obtain new Lanczos-type algorithms. We will show that some recurrence relations exist, but cannot be used to derive Lanczos-type algorithms, while others do not exist at all.

연구 동기 및 목표

• 표준 차수 차이 제약을 초월하여 형식적 직교다항식(FOPs) $P_k$와 $P_k^{(1)}$ 간의 새로운 재귀 관계를 탐색한다.
• 이러한 새로운 재귀 관계 중에서 수학적으로 타당하며, 랑츠 유형 알고리즘을 구성하는 데 사용될 수 있는 관계를 규명한다.
• 직교성과 정규화 조건을 만족하더라도 잔여 벡터 $\mathbf{r}_k$와 반복값 $\mathbf{x}_k$를 계산하는 데 부적합한 관계를 식별하고 제외한다.
• 새로운 재귀 관계를 존재하지만 사용 불가능한 것, 전혀 존재하지 않는 것, 알고리즘 설계에 유용한 것의 세 부류로 분류한다.
• 유효한 $A_i/B_j$ 재귀 쌍의 조합을 통해 새로운 랑츠 유형 알고리즘을 체계적으로 생성할 수 있는 프레임워크를 제공한다.

제안 방법

• 기존의 차수 차이가 1 또는 2인 것에 국한되지 않고, 양변의 다항식 차수 차이가 최대 3까지 허용하는 형식적 직교다항식(FOP) 프레임워크를 확장하여 재귀 관계를 분석한다.
• 직교성과 정규화 조건을 사용하여 $P_k(x)$에 대한 재귀 관계 $A_i$와 $P_k^{(1)}(x)$에 대한 $B_j$를 유도하고 평가한다.
• $P_k^{(1)}$의 내적 구조를 정의하기 위해 선형 함수 $c^{(1)}(\xi^i) = c(\xi^{i+1})$를 적용하여 재귀적 구성이 가능하게 한다.
• 계수 일치와 직교성 제약 조건을 활용하여 고차수 재귀에서 계수($C_k, D_k, F_k, G_k$)에 대한 연립방정식을 유도한다.
• 계수 시스템의 가역성을 평가하기 위해 행렬식 분석($\Delta'_k$)을 활용하여, 타당한 재귀 구조를 위한 비특이성 보장을 확보한다.
• 유효한 재귀 쌍에서 유도된 알고리즘 조합, 예를 들어 $A_{13}/B_{13}$, $A_{16}/B_{16}$, $A_{19}/B_{15}$ 등을 제안하여 새로운 랑츠 유형 해법기를 생성한다.

실험 결과

연구 질문

• RQ1 표준 차수 차이 제약을 초월하여 형식적 직교다항식 $P_k$와 $P_k^{(1)}$ 간의 새로운 재귀 관계는 무엇이 존재하는가?
• RQ2 최근에 식별된 이러한 재귀 관계 중에서 수학적으로 타당하나, 랑츠 유형 알고리즘 구현에 부적합한 관계는 무엇인가?
• RQ3 직교성 또는 정규화 조건 위반으로 인해 존재하지 않는 재귀 관계는 무엇인가?
• RQ4 유용한 재귀 쌍($A_i/B_j$)을 체계적으로 조합하여 새로운 랑츠 유형 알고리즘을 생성할 수 있는가?
• RQ5 FOP 기반 랑츠 방법에서 고차수 계수 시스템의 가역성과 안정성을 보장하는 조건은 무엇인가?

주요 결과

• 관계 $A_{11}, A_{17}, B_{11}, B_{12}$는 직교성과 정규화의 근본적인 모순으로 인해 존재하지 않는다.
• 관계 $A_{14}, A_{15}, A_{18}, B_{14}$는 수학적으로 존재하지만, 유효한 $\mathbf{r}_k$ 및 $\mathbf{x}_k$ 벡터를 생성하지 못하므로 랑츠 유형 알고리즘에 사용될 수 없다.
• 관계 $A_{12}, A_{13}, A_{16}, A_{19}, B_{13}, B_{15}, B_{16}$는 수학적으로 타당하며 알고리즘 구현에 적합하다.
• 관계 $A_{12}$는 독립적으로도 새로운 랑츠 유형 알고리즘을 생성할 수 있으며, [10:Salhi]에서 이를 보여주었다.
• 유효한 알고리즘 조합으로는 $A_{13}/B_{13}$, $A_{13}/B_{15}$, $A_{13}/B_{16}$, $A_{16}/B_{13}$, $A_{16}/B_{15}$, $A_{16}/B_{16}$, $A_{19}/B_{13}$, $A_{19}/B_{15}$, $A_{19}/B_{16}$ 등이 있으며, 모두 새로운 랑츠 유형 해법기를 도출한다.
• 재귀 계수 시스템의 가역성은 행렬식 $\Delta'_k$의 영이 아님으로써 결정되며, 이는 계수 $D_k, F_k, G_k$에 대한 유일한 해를 보장한다.