[논문 리뷰] Accurate calculation of oblate spheroidal wave functions

Van Buren, L Arnie|arXiv (Cornell University)|2017. 08. 26.

Geophysics and Gravity Measurements참고 문헌 10인용 수 4

한 줄 요약

이 논문은 64비트 산술을 사용하여 광범위한 매개변수 범위에서 타원 구면 함수 R1ml 및 R2ml를 정확하게 계산하는 안정적인 수치적 방법을 제시한다. 기존 방법의 실패를 극복한다. 구면 함수를 포함하는 급수 전개, 적분 표현, 큰 c에 대한 고유값 근접성, 그리고 수정된 계수 계산 방법을 활용하여, 최소 8자리 정밀도를 확보하는 포트란 프로그램(Oblfcn)을 개발하였다. 이는 m ≤ 1000 및 c ≤ 2000 범위에서 성능을 보이며, 128비트 산술을 통해 15자리 이상의 정밀도와 전체 xi 범위를 커버할 수 있다.

ABSTRACT

Alternative expressions for calculating the oblate spheroidal radial functions of both kinds R1ml and R2ml are shown to provide accurate values over very large parameter ranges using 64 bit arithmetic, even where the traditional expressions fail. First is the expansion of the product of a radial function and the angular function of the first kind in a series of products of the corresponding spherical functions, with the angular coordinate being a free parameter. Setting the angular coordinate equal to zero leads to accurate values for R2ml when the radial coordinate xi is larger than 0.01 and l is somewhat larger than m. Allowing it to vary with increasing l leads to highly accurate values for R1ml over all parameter ranges. Next is the calculation of R2ml as an integral of the product of S1ml and a spherical Neumann function kernel. This is useful for smaller values of xi. Also used is the near equality of pairs of low order eigenvalues when the size parameter c is large that leads to accurate values for R2ml using neighboring accurate values for R1ml. A modified method is described that provides accurate values for the necessary expansion coefficients when c is large and l is near m and traditional methods fail. A resulting Fortran computer program Oblfcn almost always provides radial function values with at least 8 accurate decimal digits using 64 bit arithmetic for m up to at least 1000 with c up to at least 2000 when xi is greater than 0.000001 and c up to at least 5000 when xi is greater than 0.01. Use of 128 bit arithmetic extends the accuracy to 15 or more digits and extends xi to all values other than zero. Oblfcn is freely available.

연구 동기 및 목표

• 큰 매개변수 범위에서 기존 방법이 실패하는 타원 구면 함수 R1ml 및 R2ml의 수치적 불안정성과 실패 문제를 해결하기 위해.
• 표준 64비트 부동소수점 산술을 사용함에도 불구하고, m, c, xi의 극단적 값에서도 높은 정밀도를 유지할 수 있는 방법을 개발하기 위해.
• 이전 방법이 실패하는 작은 xi 값(1e-6까지)과 높은 m, c 값(최대 5000까지)의 범위로 타원 구면 함수의 계산 가능 범위를 확장하기 위해.
• 파동 전파 및 산란 문제에서 고정밀도 결과가 필요한 연구자들이 사용할 수 있는 신뢰성 있고 무료로 이용 가능한 계산 도구(Oblfcn)를 제공하기 위해.

제안 방법

• 첫 번째 종류의 구면 함수의 극좌표 및 각도 함수의 곱을, 각도 좌표를 자유 매개변수로 삼아 구면 함수의 곱으로 이루어진 급수로 전개하여, 영각도에서 R2ml의 정확한 평가를 가능하게 한다.
• l의 증가에 따라 조정되는 가변 각도 좌표를 사용하여, 모든 매개변수 범위에서 R1ml의 정확한 값 확보.
• R2ml를 구면 뉴먼 함수 커널과 함께 타원 구면 함수 S1ml를 포함하는 적분으로 표현하여, 작은 xi 값에서 효과적으로 작동하도록 한다.
• 큰 크기 매개변수 c에서 저차수 고유값이 근접함을 이용하여, 인접한 R1ml 값으로부터 R2ml를 정확하게 계산한다.
• 표준 방법이 실패하는 경우, c가 크고 l이 m에 가까운 경우에 대해 계수 계산을 수정한 알고리즘을 도입한다.
• 64비트 산술을 기본으로 사용하고, 높은 정밀도를 위해 128비트 산술도 지원하는 포트란 프로그램(Oblfcn)으로 방법을 구현한다.

실험 결과

연구 질문

• RQ1 표준 64비트 부동소수점 산술을 사용하여 극도로 큰 매개변수 범위에서 타원 구면 함수를 정확하게 계산할 수 있는가?
• RQ2 특히 작은 xi 및 높은 m 또는 c 값에서 R1ml 및 R2ml의 기존 계산 방법의 수치적 불안정성을 어떻게 극복할 수 있는가?
• RQ3 큰 c에서 고유값 근접성을 어떻게 활용하여 R1ml 값으로부터 R2ml의 정밀도를 향상시킬 수 있는가?
• RQ4 c가 크고 l이 m에 가까울 때, 표준 방법이 실패하는 상황에서 계수 계산에 필요한 수정은 무엇인가?
• RQ5 일반 목적의 무료 사용 가능한 계산 도구가 m, c, xi 전체 범위에서 일관된 고정밀도(예: 8자리 이상)를 달성할 수 있는가?

주요 결과

• 제안된 방법은 m 최소 1000 및 c 최소 2000까지 64비트 산술을 사용하여 R1ml 및 R2ml를 정확하게 계산할 수 있으며, xi > 0.000001 조건을 만족한다.
• xi > 0.01일 경우, 동일한 정밀도로 c 최대 5000까지 지원하여 적용 범위가 확장됨을 보여준다.
• Oblfcn에서 128비트 산술을 사용하면 정밀도가 15자리 이상으로 향상되고, 0을 제외한 모든 xi 값에서 계산이 가능하다.
• Oblfcn 프로그램은 지정된 매개변수 범위에서 근본적으로 최소 8자리 정확한 소수 자릿수를 보장한다.
• 급수 전개 및 적분 기반 접근법은 기존 표현식이 실패하는 경우, 특히 극단적 매개변수 영역에서 안정적이고 정확한 결과를 제공한다.
• c가 크고 l이 m에 가까울 때 발생하던 수치적 실패를 성공적으로 해결하기 위해 수정된 계수 계산 방법이 효과를 보였다.