[논문 리뷰] Two-dimensional generalization of the Muller root-finding algorithm and its applications

Plamen Fiziev, Denitsa Staicova|arXiv (Cornell University)|2010. 05. 28.

Model Reduction and Neural Networks참고 문헌 18인용 수 7

한 줄 요약

이 논문은 복소수 변수를 가진 두 개의 비선형 초월방정식 시스템을 해결하기 위한 Müller 수치해법의 두 차원 일반화를 제안한다. 특히, 공통 Heun 함수를 포함하는 문제에 매우 효과적이며, 뉴턴의 방법과 브로이덴의 방법보다 수렴 속도가 더 빠르다. 이는 유사 정규 모드 계산과 같은 도전적인 경우—스카우즈실드 및 킬러 블랙홀의 경우에서 도함수 평가에 실패하거나 수렴 속도가 느린 경우—에서 특히 유리하다.

ABSTRACT

We propose a new algorithm for solving a system of two nonlinear transcendental equations with two complex variables based on the Muller algorithm. The two-dimensional Muller algorithm is tested on systems of different type and is found to work comparably to Newton's method and Broyden's method in many cases. The new algorithm is particularly useful in systems featuring the Heun functions whose complexity may make the already known algorithms not efficient enough or not working at all. In those specific cases, the new algorithm gives distinctly better results than the other two methods. As an example for its application in physics, the new algorithm was used to find the quasi-normal modes (QNM) of Schwarzschild black hole described by the Regge-Wheeler equation. The numerical results obtained by our method are compared with the already published QNM frequencies and are found to coincide to a great extent with them. Also discussed are the QNM of the Kerr black hole, described by the Teukolsky Master equation.

연구 동기 및 목표

• 특수 함수인 Heun 함수를 포함하는 복잡한 비선형 시스템을 해결할 때 뉴턴의 방법과 브로이덴의 방법의 한계를 해결한다.
• 공통 Heun 함수를 포함하는 시스템에서 도함수 평가의 계산적 과제와 열악한 수렴 성능에 대응한다.
• 이론 물리학의 스펙트럼 문제, 예를 들어 블랙홀의 준정규 모드에 적합한 도함수 없이도 신속하게 수렴하는 알고리즘을 개발한다.
• 기본 반복 방법이 실패하거나 너무 느리게 수렴하는 시스템에 대해 강력한 수치적 대안을 제공한다.
• 실제 물리 문제, 즉 레이지-웨일러 및 툴콜스키 방정식을 통한 블랙홀 준정규 모드에 대해 알고리즘의 유효성을 입증한다.

제안 방법

• 원래 1차원 Müller 알고리즘—초기 수렴 속도가 초월선형적(~1.84)인 2차 보간법—을 두 복소수 변수로 일반화한다.
• 세 개의 초기 점을 기반으로 한 2차원 보간을 사용하여 도함수 계산을 피함으로써 근을 근사화한다.
• 1차원 Müller 서브루틴을 반복적으로 두 차원 프레임워크 내에서 사용하여 근의 추정치를 정밀화한다.
• 1차원 서브루틴에서 최대 반복 수 $ P $ 를 조정하여 수렴을 튜닝함으로써 전체 성능과 수렴 속도에 영향을 미친다.
• 함수 순서와 선형 조합을 특별히 고려하여 수렴 신뢰도를 향상시키는 방식으로 두 복소수 방정식 시스템에 알고리즘을 적용한다.
• Heun 함수를 포함하는 시스템을 사용한 수치적 테스트에서 뉴턴의 방법과 브로이덴의 방법과의 결과를 비교하여 알고리즘을 검증한다.

실험 결과

연구 질문

• RQ1 도함수 없이도 두 차원 비선형 시스템에서 뉴턴의 방법과 브로이덴의 방법보다 성능이 뛰어난 알고리즘이 가능할 수 있는가?
• RQ2 도함수 평가가 어려운 공통 Heun 함수를 포함하는 시스템에서 두 차원 Müller 알고리즘이 어떻게 작동하는가?
• RQ3 초기 추정치와 함수 순서가 새로운 알고리즘의 수렴 행동과 근 선택에 어떤 영향을 미치는가?
• RQ4 1차원 Müller 서브루틴에서 반복 수 $ P $ 를 조정하면 두 차원 방법의 수렴 속도와 정확도에 상당한 영향을 미치는가?
• RQ5 새로운 알고리즘이 레지-웨일러 및 툴콜스키 방정식으로 기술되는 스키아우즈실드 및 킬러 블랙홀의 준정규 모드 주파수를 신뢰성 있게 계산할 수 있는가?

주요 결과

• 공통 Heun 함수를 포함하는 시스템에서 도함수 평가에 실패하거나 정확하지 않을 경우, 두 차원 Müller 알고리즘이 뉴턴의 방법과 브로이덴의 방법보다 더 빠른 수렴 속도를 보였다.
• 스카우즈실드 블랙홀 준정규 모드의 경우, 새로운 방법으로 계산된 주파수는 기존에 발표된 결과와 매우 유사하게 일치하여 정확도가 검증되었다.
• 어떤 경우에는 1차원 서브루틴에서 최대 반복 수 $ P $ 를 증가시켜도 총 계산 시간이 감소하는 경우가 있었으며, 이는 $ P $ 에 대한 비단조화적 의존성임을 시사한다.
• 알고리즘의 성능은 함수 순서와 초기 조건에 민감하다. 함수 순서를 재정렬하거나 방정식의 선형 조합을 사용하면 발산 또는 원하지 않는 근으로의 수렴을 방지할 수 있다.
• 브로이덴의 방법이 일반적으로 강건하더라도, Heun 함수를 포함하는 시스템에서는 두 차원 Müller 방법이 더 빠른 수렴 속도를 보였다.
• 실근 또는 단순한 비선형성을 포함하는 시스템에서는 두 차원 Müller 방법이 뉴턴의 방법과 브로이덴의 방법과 유사하게 성능을 보였지만, 특수 함수를 포함하는 복잡하고 특이한 경우에서 두드러진 성능을 발휘했다.