Krylov 부분 공간 방법을 멀티 그리드에보다 매끄럽게 사용할 수 있습니까?


15

내가 아는 한, 멀티 그리드 솔버는 Jacobi, Gauss-Seidel 및 SOR과 같은 반복적 인 스무더를 사용하여 다양한 주파수에서 오류를 완화합니다. Krylov subspace 방법 (conjugate gradient, GMRES 등)을 대신 사용할 수 있습니까? 나는 그들이 "스무더 너 (smoothener)"로 분류되지는 않는다고 생각하지만, 거친 그리드 솔루션을 근사화하는데 사용될 수 있습니다. 표준 멀티 그리드 방법에서와 마찬가지로 솔루션과 유사한 수렴을 볼 수 있습니까? 아니면 문제에 의존합니까?

답변:


18

그렇습니다. 그러나 Krylov 방법에는 일반적으로 큰 평활 특성이 없습니다. 이는 잔류 또는 적합한 오류 표준을 최소화하는 적응 방식으로 전체 스펙트럼을 대상으로하기 때문입니다. 여기에는 일반적으로 거친 격자가 잘 처리 할 수있는 일부 저주파 (장파장) 모드가 포함됩니다. Krylov 스무더는 또한 멀티 그리드주기를 비선형으로 만들므로 멀티 그리드를 외부 Krylov 방법의 전제 조건으로 사용하는 경우 외부 방법은 "유연성"이어야합니다 (예 : GCR 또는 FGMRES).

Krylov 스무더를 사용하면 계산해야하는 도트 제품 수가 크게 증가하여 병목 현상이 크게 줄어 듭니다. 그러나 이러한 매력적이지 않더라도 Krylov 스무더는 때로는 유용한 보간 연산자를 사용할 수없는 어려운 문제에 유용합니다.

λ최대11(0.1λ최대,1.1λ최대)15510)는 를 추정하는 데 GMRES 또는 CG의)를 사용 하므로 사용자가 직접 계산할 필요는 없습니다. 의 추정치는 일부 대수 멀티 그리드 방법에서도 전략을 선택하는 데 사용됩니다.λ최대λ최대

Adams, Brezina, Hu 및 Tuminaro (2003) 는 다항식 스무더의 병렬 및 알고리즘 성능에 대한 훌륭한 논문입니다. 다항식 스무더는 비대칭 문제에 대해 효과가 떨어지거나 공식화하기 어려운 경향이 있으므로 Gauss-Seidel 또는보다 정교한 (블록 / 분산) 이완 기법을 사용하는 것이 좋습니다.


다항식 및 / 또는 크리 로프 스무더에 대한 유용한 자료를 제안 할 수 있습니까? 나는 실제로 한 번도 들어 본 적이 없다 :)
Paul

@JedBrown : PDE 또는 이중선 형태 의미에서 "타원형"을 의미합니까 (즉, 연산자 또는 주요 기호의 모든 고유 값이 양수임을 의미합니까?) 나는 당신이 포인트 블록 Jacobi에 대해 이야기하고 있기 때문에 후자를 가정하고 있습니다.
Jack Poulson

폴 참조를 추가했습니다. @Jack 엄밀히 말하면, 이산 연산자는 SPD 여야하지만 실제로는 스펙트럼이 너무 잘 분포되지 않는 한 방법이 작동하는 경향이 있습니다.
Jed Brown
당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.