[논문 리뷰] Algebraic distance for anisotropic diffusion problems: multilevel results

Achi Brandt, James Brannick외 2명|arXiv (Cornell University)|2014. 09. 16.

Advanced Numerical Methods in Computational Mathematics참고 문헌 21인용 수 11

한 줄 요약

이 논문은 비정렬도가 있는 확산 문제를 위한 대수적 다중격자(AMG) 방법에서 연결 강도를 정의하기 위해 최소 제곱 최소화를 통한 유연한 테스트 벡터 기반의 대수적 거리 측정법을 제안한다. 이 방법은 매개변수 없는 안정적인 다운샘플링과 보간을 가능하게 하며, 격자에 따라 정렬되지 않은 비정렬도를 적응적으로 포착하여 조정 없이 다양한 비정렬도 방향에서 스케일러블한 다중격자 수렴을 달성한다.

ABSTRACT

In this paper we motivate, discuss the implementation and present the resulting numerics for a new definition of strength of connection which is based on the notion of algebraic distance. This algebraic distance measure, combined with compatible relaxation, is used to choose suitable coarse grids and accurate interpolation operators for algebraic multigrid algorithms. The main tool of the proposed measure is the least squares functional defined using a set of relaxed test vectors. The motivating application is the anisotropic diffusion problem, in particular problems with non-grid aligned anisotropy. We demonstrate numerically that the measure yields a robust technique for determining strength of connectivity among variables, for both two-grid and multigrid solvers. %We illustrate the use of the measure to construct, in addition, an adaptive aggregation form of interpolation for the targeted anisotropic problems. %Our approach is not a two-level approach -- we provide preliminary results that show its extendability to multigrid. The proposed algebraic distance measure can also be used in any other coarsening procedure, assuming a rich enough set of the near-kernel components of the matrix for the targeted system is known or computed.

연구 동기 및 목표

• 기존의 M행렬 히وري스틱에 의존하는 전통적 AMG가 비정렬도가 있는 확산 문제에서 실패하는 이유를 해결한다.
• 행렬의 부호 구조나 격자 정렬에 의존하지 않는 강력한 연결 강도 측정법을 개발한다.
• 조정 없이도 비정렬도 문제에 대해 효과적인 다중격자 선택과 정확한 보간을 가능하게 한다.
• 최소 제곱 보간과 호환성 있는 리콜리베이션을 사용하여 일정 계수 비정렬도 확산 문제에서 스케일러블한 다중격자 수렴을 달성한다.
• 근사 핵 성분이 알려져 있거나 계산 가능한 경우 일반적인 다운샘플링 절차를 지원할 수 있도록 프레임워크를 확장한다.

제안 방법

• 유연한 테스트 벡터 집합 위에서 최소 제곱 기능을 정의하여 변수 연결성을 정량화함으로써 대수적 거리 정의.
• 대수적 거리를 사용하여 다중격자 변수와 보간 스텐실 선택을 위한 연결 강도 그래프 구축.
• 대수적 거리와 호환성 있는 리콜리베이션을 조합하여 연결 강도 정의의 타당성과 정밀도를 검증 및 개선.
• 대수적 거리 측정법과 유연한 벡터를 사용하여 최소 제곱 보간 연산자 구축.
• 다중 스무딩 단계를 포함한 W사이클 부트스트랩 설정을 적용하여 보간과 다운샘플링을 위한 유연한 및 부트스트랩 벡터 계산.
• 유연한 CG와 함께 비선형 AMLI 조건화를 사용하여 해법 단계에서 수렴 성능 평가.

실험 결과

연구 질문

• RQ1 최소 제곱 최소화 기반의 대수적 거리가 비정렬도 확산 문제에 대해 강력한 연결 강도 측정법을 제공할 수 있는가?
• RQ2 제안된 방법이 매개변수 조정 없이도 비정렬도가 있는 확산 문제에서 스케일러블한 수렴을 달성하는가?
• RQ3 미세격자 이산화가 비정렬도 방향과 일치하지 않을 경우 방법의 성능은 어떠한가?
• RQ4 호환성 있는 리콜리베이션과 최소 제곱 보간이 어려운 비정렬도 케이스에서 효과적인 다중격자 근사치를 보장할 수 있는가?
• RQ5 변수 계수 비정렬도 문제에서 수렴 스케일러비리티에 대해 보간의 정밀도(c)가 미치는 영향은 무엇인가?

주요 결과

• 대수적 거리 측정법은 비정렬도가 있는 확산 문제에서 강한 연결성을 성공적으로 식별하여 효과적인 다중격자 선택을 가능하게 한다.
• ε = 0.1인 일정 계수 비정렬도 확산 문제에서, 모든 비정렬도 각도에서 4–7개의 다중격자 사이클 동안 사이클당 1.4–2.0회의 반복을 달성한다.
• 미세격자 이산화가 잘못 정렬된 경우(α = −π/4)에도, 방법은 실제 비정렬도 방향에 맞게 적응하는 보간을 구성한다.
• ε = 0.0001일 때에도 수렴이 유지되며, 사이클당 5–9회의 반복을 기록하여 고비정렬도에 대한 강건성을 입증한다.
• c = 2인 최소 제곱 보간을 사용할 경우 일정 계수 문제에서 거의 최적의 다수준 성능을 달성한다.
• 문제 크기에 따른 반복 횟수의 경미한 의존성은 변수 계수 케이스에서 완전한 스케일러비리티를 확보하기 위해 더 높은 정밀도 보간(c > 2)이 필요함을 시사한다.