[논문 리뷰] Several new quadrature formulas for polynomial integration in the triangle

Mark A. Taylor, Beth Wingate외 1명|ArXiv.org|2005. 01. 27.

Numerical Methods and Algorithms참고 문헌 12인용 수 28

한 줄 요약

이 논문은 다각형에서 다항식을 정확히 차수 25까지 적분할 수 있는 최적의 구적공식을 계산하기 위한 새로운 카디널 함수 알고리즘을 제시한다. 다항식 공간의 차원과 일치하도록 구적점의 수를 제약함으로써, 뉴턴 유형 최적화를 통해 점의 위치와 가중치를 효율적으로 구하고, 7개의 새로운 공식을 도출하여 이전 결과를 향상시켰으며, 모든 가중치는 양수이고 모든 점은 삼각형 내부에 위치한다.

ABSTRACT

We present several new quadrature formulas in the triangle for exact integration of polynomials. The points were computed numerically with a cardinal function algorithm which imposes that the number of quadrature points $N$ be equal to the dimension of a lower dimensional polynomial space. Quadrature forumulas are presented for up to degree $d=25$, all which have positive weights and contain no points outside the triangle. Seven of these quadrature formulas improve on previously known results.

연구 동기 및 목표

• 고차수까지 정확히 적분할 수 있는 최적의 구적공식을 삼각형에서 수치적으로 효율적으로 계산하는 방법을 개발하기 위해.
• 이전 방법의 한계를 극복하기 위해, 미지수의 수를 줄이면서도 최적화를 위한 해석적 그래디언트 표현을 유지하는 카디널 함수 알고리즘을 사용하기 위해.
• 모든 가중치가 양수이고 모든 점이 삼각형 내부에 엄격히 위치하는 구적규칙을 도출하여 문헌에 알려진 결과를 향상시키기 위해.
• 차수 25까지의 다항식 적분을 위한 최적 또는 근접 최적의 점 수를 확보하기 위해, 카디널 함수 제약 조건과 자유도 경계 조건을 모두 만족시키기 위해.

제안 방법

• 이 방법은 다항식 공간 $ \mathcal{P}_d $ 의 차원과 일치하도록 구적점의 수 $ N $ 을 설정하는 카디널 함수 알고리즘을 사용한다. 즉, $ N = \frac{1}{2}(d+1)(d+2) $ 이다.
• 직교된 Kornwinder-Dubiner 다항식 기저 함수를 포함하는 선형 시스템을 풀어 뉴턴-코테스 가중치를 도출함으로써 $ \mathcal{P}_d $ 의 정확한 적분을 보장한다.
• 뉴턴-코테스 가중치의 다변수 일반화를 사용하여 가중치의 변화와 점의 위치 변화 간의 해석적 관계를 설정함으로써 최적화를 위한 효율적인 그래디언트 계산을 가능하게 한다.
• 기하학적 대칭성 제약 없이도, 기울기 강하 또는 뉴턴 반복을 적용하여 점의 위치를 최적화함으로써 더 고차수의 다항식 $ \mathcal{P}_{d+e} $ 를 적분하도록 한다.
• 자유도 경계 조건 $ \dim\mathcal{P}_{d+e} \leq 3N $ 을 사용하여 이론적으로 달성 가능한 적분 차수 $ d+e $ 의 한계를 결정한다.
• 해결책은 정규화된 $ L^2 $-노름이 2인 정규직교 기저 함수들에 대한 최대 구적오차를 계산하여 검증된다.

실험 결과

연구 질문

• RQ1 카디널 함수 알고리즘이 점의 수가 최적이고 차수 25까지 정확히 적분 가능한 삼각형 내의 구적공식을 생성할 수 있는가?
• RQ2 이러한 방법이 대칭성을 부여하지 않더라도 모든 가중치가 양수이고 삼각형 외부에 점이 존재하지 않는 구적규칙을 도출할 수 있는가?
• RQ3 계산된 구적공식이 삼각형 내 다항식 적분 분야에서 기존 문헌에 알려진 최고의 결과를 향상시키는가?
• RQ4 카디널 함수 제약 조건 하에서 $ N = \dim\mathcal{P}_d $ 개의 점을 사용할 때, 정확히 적분할 수 있는 다항식의 최대 차수는 얼마인가?

주요 결과

• 차수 $ d+e = 16 $ 에서 $ 25 $ 까지의 다항식 적분을 위한 7개의 새로운 구적공식이 이전에 발표된 결과를 개선하였으며, 점의 수를 줄였다.
• 차수 $ d+e = 13 $ 에서는 36개의 점이 삼각형 내부에 위치하고 가중치가 모두 양수인 새로운 공식이 도출되었으며, 이는 이전 결과보다 대칭적인 점이 도메인 외부에 존재하는 경우를 향상시켰다.
• 모든 계산된 구적규칙은 가중치가 모두 양수이며, 모든 점이 삼각형 내부에 엄격히 위치해 있으며, 기하학적 대칭성 제약 조건이 없음에도 불구하고 이를 만족한다.
• 모든 정규직교 기저 함수에 대한 최대 구적오차는 $ 4.3 \times 10^{-14} $ 이하였으며, 이는 높은 수치 정확도를 의미한다.
• 이 방법은 카디널 함수 제약 조건과 자유도 경계 조건 $ \dim\mathcal{P}_{d+e} \leq 3N $ 을 모두 만족하는 차수 25까지의 최적 또는 근접 최적의 공식을 성공적으로 계산하였다.
• 차수 $ d+e = 11, 13, 14, 16, 17, 18, 25 $ 에 대해 비대칭 해를 도출하였으며, 'asym'으로 표기함으로써 대칭성이 최적 성능을 위해 필수적이지 않음을 시사한다.