[논문 리뷰] Preconditioning rectangular spectral collocation

Kui Du|arXiv (Cornell University)|2015. 10. 01.

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

한 줄 요약

이 논문은 일반 선형 제약 조건을 가진 고차 미분방정식에서 발생하는 선형 연립방정식의 심각한 불안정성 문제를 해결하기 위해 직사각형 스펙트럴 콜로케이션(RSC) 방법에 대한 가짜스펙트럴 통합 조정 방법을 제안한다. Birkhoff형 보간 기반 조정 행렬을 구성함으로써 시스템의 조건수를 콜로케이션 점 수에 영향을 받지 않도록 만들었으며, 이는 GMRES와 같은 반복 해법을 사용할 때 조건수에 관계없이 몇 번의 반복만으로도 빠른 수렴을 가능하게 한다. 이는 큰 N 값에서도 성립한다.

ABSTRACT

Rectangular spectral collocation (RSC) methods have recently been proposed to solve linear and nonlinear differential equations with general boundary conditions and/or other constraints. The involved linear systems in RSC become extremely ill-conditioned as the number of collocation points increases. By introducing suitable Birkhoff-type interpolation problems, we present pseudospectral integration preconditioning matrices for the ill-conditioned linear systems in RSC. The condition numbers of the preconditioned linear systems are independent of the number of collocation points. Numerical examples are given.

연구 동기 및 목표

• 콜로케이션 점 수 N이 증가함에 따라 직사각형 스펙트럴 콜로케이션(RSC) 방법에서 발생하는 선형 연립방정식의 극심한 불안정성 문제를 해결하기 위해.
• 결과로 얻어지는 선형 연립방정식의 조건수가 N에 독립적이게 하는 조정 행렬을 개발하여 안정적이고 효율적인 반복적 해법을 보장하기 위해.
• 이전에 표준 스펙트럴 콜로케이션에 사용된 가짜스펙트럴 통합 조정 프레임워크를 일반적인 선형 제약 조건이 있는 더 일반적인 직사각형 스펙트럴 콜로케이션 설정으로 확장하기 위해.
• 수치적으로 조건수가 N에 관계없이 빠르게 수렴하는 고차 미분방정식의 해법을 유지함을 보여주기 위해, 심한 불안정성이나 고주파 진동 문제에서도 성능이 우수함을 입증하기 위해.

제안 방법

• m개의 일반 선형 제약 조건이 있는 m차 미분 연산자의 이산화된 형태에 대해 정확한 역행렬이 되는 가짜스펙트럴 통합 행렬을 정의하기 위해 Birkhoff형 보간 문제를 도입한다.
• 도함수 평가를 위한 첫 번째 종류의 체비셰프 점들과 콜로케이션 점들을 위한 두 번째 종류의 체비셰프 점들 사이의 바리센트 재샘플링 행렬을 사용하여 오른쪽 조정 행렬을 구성함으로써 안정성과 정확성을 확보한다.
• 함수 값이 콜로케이션 점에 위치할 때 도함수 값을 별개의 점들에서 계산할 수 있도록 하는 직사각형 미분 행렬 D^{(m)}_{x→y}를 활용하여 m차 미분 연산자의 안정적인 이산화를 가능하게 한다.
• 가짜스펙트럴 통합 행렬을 RSC 시스템의 오른쪽 조정 행렬로 사용하여 원래의 불안정한 시스템을 N에 독립적인 조건수를 가지는 시스템으로 변환한다.
• 상대 오차 기준 1e-10과 재시작 수 40을 사용하여 GMRES 반복 해법을 적용하여 조정된 시스템을 풀며, 몇 번의 반복만으로도 수렴을 보장한다.
• 콜로케이션 점으로서 두 번째 종류의 체비셰프 점들을, 도함수 평가 점으로서 첫 번째 종류의 체비셰프 점들을 사용하여 스펙트럴 정확도와 조건수 최적화를 달성한다.

실험 결과

연구 질문

• RQ1 콜로케이션 점 수 N이 증가함에 따라 RSC 선형 연립방정식의 조건수가 유한하게 유지되는 조정 행렬을 구성할 수 있는가?
• RQ2 일반 선형 제약 조건이 있는 강성 또는 고주파 진동 문제에 대해 제안된 조정 전략은 어떻게 성능을 보이는가?
• RQ3 표준 스펙트럴 콜로케이션에서의 가짜스펙트럴 통합 조정 프레임워크를 임의의 선형 제약 조건이 있는 직사각형 스펙트럴 콜로케이션 설정으로 일반화할 수 있는가?
• RQ4 GMRES와 같은 반복 해법의 수렴 속도에 조정 행렬이 고차 미분방정식 문제에서 미치는 영향은 무엇인가?

주요 결과

• a_0(x)=2x와 제약 조건 (5.2)를 가진 RSC의 경우, 조건수는 N에 관계없이 약 3.19로 거의 일정하며, 비조정 시스템은 O(N^{2.5}) 비율로 증가한다.
• ε=0.01인 2차 문제에서 비조정 RSC의 조건수는 O(N^{4.5}) 비율로 증가하여 N=1024일 때 4.84e+11에 이를 뿐만 아니라, 조정된 RSC 시스템은 N에 관계없이 조건수 3.70e+05를 유지한다.
• ε=1인 P-RSC의 경우, N=1024일 때도 GMRES 해법이 단 8회 반복만으로 수렴하지만, 비조정 RSC는 1000회 반복 이내에 수렴하지 못한다.
• ε=0.01일 때 P-RSC는 N=1024에서 67회 반복만으로 수렴하지만, 비조정 RSC는 1000회 반복 이내에 수렴하지 못함을 보여주며, 조정의 필요성을 입증한다.
• 모든 N 및 ε 값에서 P-RSC의 최대 점별 오차는 1e-13에서 1e-15 수준을 유지하여 높은 수치적 정확도를 보여준다.
• N=2200일 때 정확한 해와 수치적 해의 비교를 통해, 고주파 진동 해에 대해서도 이 방법이 높은 정확도와 안정성을 유지함을 입증한다.