어떤 적용 사례에서 부가적인 사전 조정 체계가 곱셈 체계보다 우수합니까?


12

도메인 분해 (DD) 및 멀티 그리드 (MG) 방법 모두에서, 블록 업데이트 또는 거친 보정의 적용을 가산 적 또는 곱셈으로 구성 할 수 있습니다. 포인트 솔버의 경우 이것은 Jacobi와 Gauss-Seidel 반복의 차이점입니다. S ( x o l d , b ) = x n e w 로 작용하는 대한 곱셈 스무더가 다음 과 같이 적용됩니다.Ax=bS(xold,b)=xnew

xi+1=Sn(Sn1(...,S1(xi,b)...,b),b)

첨가제 스무더가 다음과 같이 적용됩니다.

xi+1=xi+=0nλ(S(xi,b)xi)

약간의 감쇠 . 일반적인 합의는 곱셈 스무더가 훨씬 더 빠른 수렴 속성을 갖는 것으로 보이지만 궁금합니다. 어떤 상황에서 이러한 알고리즘의 부가 변형의 성능이 더 좋습니까?λi

보다 구체적으로, 첨가 변이체가 곱셈 변이체보다 현저히 우수하고 /하거나 성능이 우수한 유스 케이스를 가지고있는 사람이 있습니까? 이것에 대한 이론적 인 이유가 있습니까? 멀티 그리드에 대한 대부분의 문헌 은 Additive 방법에 대해 상당히 비관적 이지만 DD 컨텍스트에서는 Add Schwarz와 같이 많이 사용됩니다. 이것은 선형 및 비선형 솔버를 구성하는 훨씬 일반적인 문제와 어떤 유형의 구조가 잘 수행되고 병렬로 잘 수행되는지에 이르기까지 확장됩니다.

답변:


6

추가 방법은 더 많은 동시성을 노출시킵니다. 해당 동시성을 사용할 수 있으면 일반적으로 곱셈 방법보다 빠릅니다. 예를 들어, 거친 수준의 멀티 그리드는 일반적으로 대기 시간이 제한됩니다. 거친 수준을 더 작은 하위 통신 자로 이동하면 더 작은 수준과 독립적으로 해결할 수 있습니다. 곱셈 체계를 사용하면 모든 수준의 프로세스가 거친 수준이 해결되는 동안 기다려야합니다.

또한 알고리즘이 모든 수준에서 감소가 필요한 경우, 곱셈법이 순차적으로 수행하도록 첨가제 변형이 알고리즘을 통합 할 수 있습니다.


이것이 내가 얻을 것이라고 생각한 답변이므로 질문을 계속 진행할 것입니다. DD 및 MG를 포함하여 부가 적으로 적용되는 방법뿐만 아니라 필드 분할 (DD 형으로 간주되지만 실제로는 다른 특성을 가질 수 있음) 또는 PDE 분할이 실제로 곱셈 변형보다 성능, 견고성 또는 안정성 측면에서 더 나은 상황이 있습니까?
Peter Brune

1
많은 알고리즘의 곱셈 버전은 더 많은 정보를 저장해야하지만 때로는 대략 빠른 속도로 수렴합니다. 때로는 추가 변형이 대칭이지만 곱셈 대칭을 만드는 것이 훨씬 더 효과적 일 수 있습니다. 필드 분할을 사용하면 추가 해결을 추가 할 때 전제 조건이 더 근사해질 수 있습니다. 원하는 경우 PETSc Stokes 예제로이를 설명 할 수 있습니다. 첨가제는 항상 벡터화 / 동시에 더 쉬워 지지만, 성능 및 성능은 문제 및 아키텍처에 따라 다릅니다.
Jed Brown

5

SPD 문제의 경우 이미 언급했듯이 몇 가지 이유로 부가적인 방법이 MG 스무딩에 더 좋습니다.

@Article{Adams-02, 
author = {Adams, M.~F. and Brezina, M. and Hu, J. J. and Tuminaro, R. S.}, 
title = {Parallel multigrid smoothing: polynomial versus {G}auss-{S}eidel}, 
journal = {J. Comp. Phys.}, 
year = {2003}, 
volume = {188}, 
number = {2}, 
pages = {593-610} }

그러나 곱셈법은 MG 매끄럽게하기 위해 정확한 스펙트럼 특성을 가지고 있습니다. 즉, 감쇠가 필요하지 않습니다. 이것은 다항식 평활화가 그리 좋지 않은 쌍곡선 문제에 큰 도움이 될 수 있습니다.


0

나는 @Jed가 말한 것을 다시 언급 할 것이다 : Multiplicative 메소드는 적어도 Additive 메소드 (무증상)뿐만 아니라 항상 수렴하기 때문에 동시성에만 의존하지만 아키텍처에 따라 다르다.


기술적으로 올바른 것은 아니지만 Gauss-Seidel에 대한 반복 매트릭스의 스펙트럼은 Jacobi보다 균일하지 않습니다 (예 : 하나의 고유 값이 하나의 Jacobi 반복으로 종료 됨). Mark (Jed로 로그 오프하는 방법 ...)
Jed Brown
당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.