[논문 리뷰] Two hierarchies of spline interpolations. Practical algorithms for multivariate higher order splines

Cristian C. Lalescu|ArXiv.org|2009. 05. 21.

Advanced Numerical Analysis Techniques참고 문헌 3인용 수 12

한 줄 요약

이 논문은 다변수 고차 스플라인 보간법의 두 계층—헤르미트 스플라인과 격자 스플라인—을 정규 직사각형 격자 위에서 다항식 기저 함수를 사용하여 제시한다. 이 기저 함수들은 테일러 전개와 중심 유한차분에서 유도된다. 격자 스플라인을 구현하기 위한 명시적 알고리즘을 제공하며, 셀 간의 도함수 연속성을 보장하고, 충분히 규칙적인 함수에 대해 부드러운 근사로 수렴함을 보여준다.

ABSTRACT

A systematic construction of higher order splines using two hierarchies of polynomials is presented. Explicit instructions on how to implement one of these hierarchies are given. The results are limited to interpolations on regular, rectangular grids, but an approach to other types of grids is also discussed.

연구 동기 및 목표

• 정규 직사각형 격자에서 고차 스플라인 보간을 위한 체계적이고 실용적인 알고리즘을 개발한다.
• 모든 도함수를 저장하지 않고도 격자 노드에서만 알려진 함수에 대해 부드러운 보간 함수를 구성한다.
• 다항식 기저 함수를 사용하여 셀 간에 지정된 차수까지 도함수 연속성을 확보한다.
• 국소적 계산과 최소한의 프로세서 간 통신을 통해 효율적인 병렬 구현을 가능하게 한다.
• 공유 노드에서 부드러움을 유지하면서 비직사각형 격자로 접근을 일반화한다.

제안 방법

• 두 계층을 사용한다: 격자 노드에서의 함수값과 도함수를 기반으로 하는 헤르미트 스플라인과, 중심 유한차분을 통해 오직 함수값만을 기반으로 하는 격자 스플라인.
• 격자 점에서 특정 도함수 조건을 만족하는 다항식 기저 함수 α(n,l)i(x)를 사용하여 Cm 연속성을 보장한다.
• 이항계수와 내림階승의 곱을 사용하여 기저 다항식 α(n,l)0(x)와 α(n,l)1(x)의 명시적 닫힌 표현식을 유도한다.
• 스플라인 평가를 함수값과 그 도함수의 선형 조합으로 재구성하여, 큰 선형 시스템을 풀지 않고도 효율적인 계산을 가능하게 한다.
• 1차원 기저 함수의 텐서곱을 사용하여 D차원 격자에 적용하고, 도메인 분할을 통해 병렬 계산에 적용한다.
• 중심 유한차분을 사용하여 도함수를 근사함으로써 메모리 요구량을 감소시키고 셀 경계에서의 도함수 연속성을 보장한다.

실험 결과

연구 질문

• RQ1 정규 격자에서 다변수 함수에 대해 고차 스플라인 보간을 어떻게 체계적으로 구성할 수 있는가?
• RQ2 메모리 제약 조건이 있을 때 모든 도함수를 저장하지 않고도 고차 스플라인을 어떻게 최적으로 계산할 수 있는가?
• RQ3 스플라인 평가가 오직 함수값과 유한차분에 의존하도록 할 수 있는 다항식 기저를 유도할 수 있는가?
• RQ4 결과로 도출된 스플라인 계층이 셀 간에 Cm 연속성을 유지하면서도 고차 정확도를 유지할 수 있는가?
• RQ5 함수의 유한 푸리에 표현이 있을 경우, 스플라인 근사의 수렴 행동은 어떻게 되는가?

주요 결과

• 격자 스플라인 계층은 중심 유한차분을 사용하여 도함수를 근사함으로써 셀 간에 Cm 연속성을 달성하며, 모든 도함수를 명시적으로 저장할 필요 없이도 된다.
• 이항계수와 x 및 (1−x)의 거듭제곱을 포함하는 명시적 닫힌 표현식이 기저 다항식 α(n,l)0(x)와 α(n,l)1(x)에 대해 유도되었다.
• 기저 다항식은 x=0과 x=1에서 요구되는 도함수 조건을 만족하여, 스플라인이 함수값과 도함수를 차수 m=(n−1)/2까지 정확히 일치함을 보장한다.
• 오차 ||s(n+1)−s(n)||는 n이 증가함에 따라 감소하며, 충분히 부드러운 함수에 대해 0으로 수렴함을 보여주어 진짜 함수로의 수렴을 시사한다.
• 모듈러 연산을 사용한 인덱싱을 통해 각 프로세서가 오직 두 개의 인접 도메인과만 통신하도록 제한함으로써 효율적인 병렬화를 지원한다.
• 타원형 격자로의 일반화가 가능하며, 공유 노드에서 테일러 전개가 일관되게 유지되면 부드러움이 유지된다.