[논문 리뷰] Optimal Point Sets for Total Degree Polynomial Interpolation in Moderate Dimensions

Max Gunzburger, Aretha L. Teckentrup|arXiv (Cornell University)|2014. 07. 11.

Advanced Numerical Analysis Techniques참고 문헌 17인용 수 7

한 줄 요약

이 논문은 중간 차원(최대 10차원)에서 총 차수 다항식 근사에 대한 최적의 보간 점을 계산하는 새로운 알고리즘을 제안한다. 이 알고리즘은 바르멘드 행렬 계산을 요구하지 않으며, 리프시츠 상수를 최소화한다. 이 방법은 지금까지 알려진 바 중에서 가장 낮은 리프시츠 상수를 도출하며, 차원과 차수에 대해 비선형 성장을 보이며, 단위 입방체와 단위 구 내에서 점들이 경계 근처에 집중되어 있음을 드러낸다.

ABSTRACT

This paper is concerned with Lagrange interpolation by total degree polynomials in moderate dimensions. In particular, we are interested in characterising the optimal choice of points for the interpolation problem, where we define the optimal interpolation points as those which minimise the Lebesgue constant. We give a novel algorithm for numerically computing the location of the optimal points, which is independent of the shape of the domain and does not require computations with Vandermonde matrices. We perform a numerical study of the growth of the minimal Lebesgue constant with respect to the degree of the polynomials and the dimension, and report the lowest values known as yet of the Lebesgue constant in the unit cube and the unit ball in up to 10 dimensions.

연구 동기 및 목표

• 중간 차원에서 총 차수 다항식 보간에 대해 리프시츠 상수를 최소화하는 보간 점 집합을 규명하는 것.
• 영역 기하구조에 의존하지 않으며, 명시적인 바르멘드 행렬 계산이 없는 수치 알고리즘을 개발하는 것.
• 최대 10차원에서 단위 입방체와 단위 구에서 다항식 보간의 최소 리프시츠 상수의 가장 낮은 알려진 값을 계산하고 보고하는 것.
• 최적 점 집합의 구조와 분포를 분석하여, 특히 영역 경계 근처 집중 현상과의 관계를 규명하는 것.
• 이론적 추측과 일치하는 다항식 차수와 차원에 따른 리프시츠 상수의 경험적 성장률을 제공하는 것.

제안 방법

• 바르멘드 행렬의 행렬식이나 선형 시스템을 계산하지 않고도, 대규모 비선형 최적화 문제를 수치적으로 해결하여 리프시츠 상수를 최소화하는 알고리즘을 사용한다.
• 신뢰 영역 접근법을 활용하여 점 구성의 반복적 개선을 수행하며, 리프시츠 상수가 보간 점에 대한 부드러운 함수임을 활용한다.
• 영역에 종속되지 않는 방법으로, 단위 입방체와 단위 구와 같은 임의의 컴acts 영역에 적용 가능하다.
• 라그랑주 기본 다항식의 절댓값 합의 최댓값을 통해 리프시츠 상수의 기저 독립적 표현을 사용한다.
• 명시적인 바르멘드 행렬을 형성하지 않음으로써 고차수 다항식 기저와 관련된 악조건 문제를 피한다.
• 나쁜 국소 최솟값을 피하기 위해 대칭적이고 잘 분포된 초기 점 구성 방식을 선택한다.

실험 결과

연구 질문

• RQ1 중간 차원에서 총 차수 다항식 보간에 대해 리프시츠 상수를 최소화하는 최적의 보간 점 구성은 무엇인가?
• RQ2 최소 리프시츠 상수는 다항식 차수와 공간 차원이 증가함에 따라 어떻게 증가하는가?
• RQ3 최적 점 집합은 어떤 구조적 성질을 보이며, 특히 영역 기하구조와 대칭성과의 관계는 어떠한가?
• RQ4 바르멘드 행렬 계산에 의존하지 않는 수치 알고리즘을 설계할 수 있는가?
• RQ5 계산된 리프시츠 상수는 문헌에서 이전에 보고된 값과 비교해 어떻게 되는가?

주요 결과

• 단위 구와 단위 입방체에 대한 계산된 리프시츠 상수는 지금까지 알려진 바 중 가장 낮은 값이며, 특히 n ≥ 5인 경우 이전 작업에 비해 상당한 향상이 있었다.
• d=2인 경우 최소 리프시츠 상수는 약 √(n+1) 비례로 증가하며, d=3인 경우는 n에 대해 선형 증가하며, 이는 이론적 기대와 일치한다.
• 최적 점 집합은 영역 경계 근처에서 높은 밀도를 보이며, 반사 및 회전 대칭성을 가지며 강한 기하학적 구조를 나타낸다.
• d=2와 d=3인 경우 점 집합은 각각 공심 원과 구를 따라 군집하는 경향이 있으나, 정확히 가우스-레지온드르-로바토 점들과 일치하지는 않는다.
• 알고리즘은 바르멘드 행렬과 관련된 수치적 문제를 성공적으로 피하고 10차원까지 효과적으로 확장된다.
• 브리아니 등 이전 결과와의 비교에서, 모든 테스트 차수에서 리프시츠 상수가 일관되고 측정 가능한 감소를 보였다.