Hangi uygulama örneklerinde, ek ön koşullandırma şemaları çarpımsal olanlardan daha üstündür?


12

Hem alan ayrışması (DD) hem de çoklu-ızgara (MG) yöntemlerinde, blok güncellemelerinin veya kaba düzeltmelerin uygulanmasını katkı veya çoklayıcı olarak oluşturabilir . Noktasal çözücüler için, Jacobi ve Gauss-Seidel iterasyonları arasındaki fark budur. için gibi çarpma özelliği daha yumuşak olarak uygulanır.Ax=bS(xold,b)=xnew

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

ve katkı maddesi daha pürüzsüz olarak uygulanır.

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

bazı sönümleme için . Genel fikir birliği, çarpma düzleştiricilerinin çok daha hızlı yakınsama özelliklerine sahip olduğu görülüyor, ama merak ediyorum: Hangi durumlarda bu algoritmaların katkı varyantlarının performansı daha iyi?λi

Daha spesifik olarak, Katkı varyantının çarpma varyantından önemli ölçüde daha iyi performans göstermesi gereken ve / veya daha iyi performans gösteren herhangi bir kullanım durumu var mı? Bunun teorik nedenleri var mı? Multigrid hakkındaki çoğu literatür , Additive yöntemi hakkında oldukça kötümserdir , ancak DD bağlamında Schwarz katkı maddesi olarak çok kullanılır. Bu aynı zamanda doğrusal ve doğrusal olmayan çözücülerin oluşturulmasına ilişkin daha genel bir konuya ve hangi tür yapıların iyi performans göstereceğine ve paralel olarak iyi performans göstereceğine de uzanmaktadır.

Yanıtlar:


6

Katkı yöntemleri daha fazla eşzamanlılık ortaya çıkarır. Bu eşzamanlılığı kullanabiliyorsanız genellikle çoğullayıcı yöntemlerden daha hızlıdırlar. Örneğin, kaba multigrid seviyeleri tipik olarak gecikme ile sınırlıdır. Kaba seviyeleri daha küçük alt iletişimcilere taşırsanız, daha ince seviyelerden bağımsız olarak çözülebilirler. Çarpımsal bir şemada, kaba seviyeler çözülürken tüm procları beklemek zorundadır.

Ayrıca, algoritmanın her seviyede azalmaya ihtiyacı varsa, katkı varyantı, çarpımsal metot bunları sırayla gerçekleştirmek zorunda kaldığı yerlerde birleştirebilir.


Bu alacağımı düşündüğüm cevap, bu yüzden sanırım soru ile daha da ileri gideceğim. DD ve MG dahil olmak üzere ilave olarak uygulanan yöntemlerin yanı sıra alan ayırma (DD benzeri kabul edilebilir, ancak uygulamada farklı özelliklere sahip olabilir) veya PDE bölünmesinin performans, sağlamlık veya kararlılık açısından çarpma varyantından daha iyi olduğu durumlar var mı?
Peter Brune

1
Birçok algoritmanın çarpımsal versiyonlarının daha fazla bilgi depolaması gerekir, ancak bazen kabaca hızlı bir şekilde birleşirler. Bazen katkı varyantları simetriktir, ancak çarpımsal simetrik yapmak çok daha fazla iş olabilir. Fieldsplit ile bu ekstra çözücüleri eklediğinizde önkoşul daha yaklaşık hale gelebilir. İsterseniz bunu bir PETSc Stokes örneği ile gösterebiliriz. Katkı maddesinin vektörelleştirilmesi / daha eşzamanlı olması her zaman daha kolaydır, ancak bundan kazanılan herhangi bir performans probleme ve mimariye özgüdür.
Jed Brown

5

SPD problemleri için, daha önce de belirtildiği gibi birkaç nedenden dolayı MG yumuşatma için ek yöntemler daha iyidir ve birkaç tane daha:

@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} }

Bununla birlikte, çoğullayıcı yöntemler bir MG pürüzsüzleştirici için kutudan çıkar çıkmaz doğru spektral özelliklere sahiptir, yani sönümlemeye ihtiyaç duymazlar. Bu, polinom düzleştirmenin çok hoş olmadığı hiperbolik problemler için büyük bir kazanç olabilir.


0

@Jed'in söylediklerini yeniden ifade edeceğim: Multiplicative yöntem her zaman en azından Additive yöntemini (asimptotik olarak) birleştirir, böylece sadece eşzamanlılığa göre kazanırsınız, ancak bu mimariye bağlıdır.


Teknik olarak doğru değil, örneğin Gauss-Seidel'in yineleme matrisinin spektrumları Jacobi'den eşit olarak üstün değildir (örneğin, bir özdeğer bir Jacobi yinelemesi ile öldürülür). Mark (Jed olarak nasıl oturumu kaparım ...)
Jed Brown
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.