[논문 리뷰] Optimal preconditioning for problems of negative order

Rob Stevenson, Raymond van Venetië|arXiv (Cornell University)|2018. 03. 14.

Advanced Numerical Methods in Computational Mathematics참고 문헌 10인용 수 5

한 줄 요약

이 논문은 연속적인 조각별 선형 함수로 이산화된 반대 차수의 유계로 가역성 있는 연산자를 사용하여 음의 차수 연산자에 대한 최적의 조절행렬을 구성한다. 이 방법은 메시 셀 수에 비례하여 선형적으로 확장 가능하며, 비대각 행렬의 역행렬 계산이 필요 없고, 메시의 균형 조절이나 바리센트 리파인먼트 요구 조건 없이 작동하여 고차수 불연속 및 연속 갈레르킨 이산화에 대해 강건하고 효율적인 해결책을 제공한다.

ABSTRACT

Optimal preconditioners for operators of negative order discretized by (dis)continuous piecewise polynomials of any order are constructed from a boundedly invertible operator of opposite order discretized by continuous piecewise linears. Besides the cost of the application of the latter discretized operator, the other cost of the preconditioner scales linearly with the number of mesh cells. Compared to earlier proposals, the preconditioner has the following advantages: It does not require the inverse of a non-diagonal matrix; it applies without any mildly grading assumption on the mesh; and it does not require a barycentric refinement of the mesh underlying the trial space.

연구 동기 및 목표

• 유한요소 이산화의 맥락에서 음의 차수 연산자에 대한 최적의 조절행렬을 개발하는 것.
• 계산적으로 비용이 많이 들 수 있는 비대각 행렬의 역행렬 계산이 필요 없도록 하는 것.
• 메시의 균형 조절이나 시험 공간 메시의 바리센트 리파인먼트 요구 조건 없이도 강건한 성능을 확보하는 것.
• 메시 셀 수에 비례하여 계산 비용이 선형적으로 증가하는 것을 보장하는 것.
• 모든 차수의 조각별 다항식 이산화에 대해 연속 및 불연속 갈레르킨 방법 모두에서 효율성과 안정성을 유지하는 것.

제안 방법

• 반대 차수의 유계로 가역성 있는 연산자를 기반으로 조절행렬을 구성하며, 특히 연속적인 조각별 선형 함수를 사용한다.
• 반대 차수 연산자의 이산화를 선택함으로써 우수한 조건수와 역행렬 계산의 용이성을 확보한다.
• 변분 공식화에서 양의 차수와 음의 차수 연산자 간의 쌍대성 덕분에 조건행렬이 최적성을 확보한다.
• 조절행렬을 적용하는 데 소요되는 비용이 메시 셀 수에 비례하여 선형적으로 증가함을 보장한다.
• 매우 부드러운 메시의 균형 조절이나 바리센트 리파인먼트와 같은 메시에 의존하는 가정에 의존하지 않는다.
• 이 방법은 일반적이며 고차수의 임의의 다항식 차수를 가진 불연속 및 연속 갈레르킨 방법에 모두 적용 가능하다.

실험 결과

연구 질문

• RQ1 비대각 행렬의 역행렬 계산이 필요 없이 음의 차수 연산자에 대한 최적의 조절행렬을 구성할 수 있는가?
• RQ2 매우 부드러운 메시의 균형 조절을 가정하지 않더라도 조건행렬이 효과적으로 유지되는가?
• RQ3 시험 공간 메시의 바리센트 리파인먼트가 필요 없도록 할 수 있는가?
• RQ4 조절행렬의 계산 비용은 메시 셀 수에 따라 어떻게 변화하는가?
• RQ5 모든 차수의 조각별 다항식 이산화에 대해 연속 및 불연속 갈레르킨 방법 모두에서 조건행렬이 강건한가?

주요 결과

• 제안된 조건행렬은 비대각 행렬의 역행렬 계산이 필요 없이 음의 차수 문제에 대해 최적의 수렴 속도를 달성한다.
• 메시에 대한 균형 조절 가정 없이도 비정규 메시에서도 여전히 효과적인 성능을 유지한다.
• 조건행렬을 적용하는 데 소요되는 계산 비용은 메시 셀 수에 비례하여 선형적으로 증가한다.
• 이 방법은 모든 다항식 차수의 불연속 및 연속 갈레르킨 이산화에 모두 적용 가능하다.
• 메시의 바리센트 리파인먼트가 필요 없어 구현이 단순해지고 오버헤드가 감소한다.
• 구성은 반대 차수의 유계로 가역성 있는 연산자를 기반으로 하여 안정성과 최적성을 보장한다.