Multigrid yöntemi - Multigrid method

İçinde Sayısal analiz, bir multigrid yöntemi (MG yöntemi) bir algoritma çözmek için diferansiyel denklemler kullanarak hiyerarşi nın-nin ayrılıklar. Bunlar, adı verilen bir teknik sınıfına örnektir. çoklu çözünürlük yöntemleri, ortaya çıkan problemlerde çok yararlı çoklu ölçekler davranış. Örneğin, birçok temel gevşeme yöntemleri kısa ve uzun dalga boylu bileşenler için farklı yakınsama oranları sergilerler, bu da bu farklı ölçeklerin, bir Fourier analizi multigrid yaklaşımı.[1] MG yöntemleri çözücü olarak kullanılabileceği gibi ön şartlandırıcılar.

Multigrid'in ana fikri, temel bir yinelemeli yöntemin (genellikle kısa dalga boyu hatasını azaltan gevşetme olarak bilinir) yakınsamasını bir küresel ince ızgara çözüm yaklaşımının zaman zaman düzeltilmesi, kaba problem. Kaba problem, çözümü daha ucuz olsa da, kısa ve uzun dalga boyu hatalarına sahip olması bakımından ince ızgara problemine benzer. Ayrıca, gevşeme ve daha kaba ızgaralara hitap etme kombinasyonu ile de çözülebilir. Bu yinelemeli süreç, doğrudan çözüm maliyetinin, ince ızgara üzerindeki bir gevşetme taramasının maliyetine kıyasla ihmal edilebilir olduğu bir ızgaraya ulaşılıncaya kadar tekrarlanır. Bu multigrid döngüsü, tipik olarak, tüm hata bileşenlerini, ince ızgara ağ boyutundan bağımsız olarak, sabit bir miktarın altında sınırlandırılmış olarak azaltır. Multigrid için tipik uygulama, sayısal çözümdür. eliptik kısmi diferansiyel denklemler iki veya daha fazla boyutta.[2]

Multigrid yöntemleri, yaygın ayrıklaştırma tekniklerinden herhangi biri ile kombinasyon halinde uygulanabilir. Örneğin, sonlu eleman yöntemi multigrid bir yöntem olarak yeniden biçimlendirilebilir.[3] Bu durumlarda multigrid yöntemler günümüzde bilinen en hızlı çözüm teknikleri arasındadır. Diğer yöntemlerin aksine, multigrid yöntemler, keyfi bölgeleri tedavi edebilmeleri ve sınır şartları. Bağımlı değiller denklemlerin ayrılabilirliği veya denklemin diğer özel özellikleri. Ayrıca, daha karmaşık simetrik olmayan ve doğrusal olmayan denklem sistemleri için yaygın olarak kullanılmıştır. Lamé denklemleri nın-nin esneklik ya da Navier-Stokes denklemleri.[4]

Algoritma

Hızlı O (n) yakınsaması için yinelemeli Multigrid algoritmasının görselleştirilmesi.

Multigrid algoritmalarının birçok çeşidi vardır, ancak ortak özellikler bir hiyerarşi olmasıdır. ayrılıklar (ızgaralar) dikkate alınır. Önemli adımlar:[5][6]

  • Yumuşatma - örneğin birkaçını kullanarak yüksek frekans hatalarını azaltmak yinelemeler of Gauss – Seidel yöntemi.
  • Artık Hesaplama - bilgi işlem Ölçüm hatası yumuşatma işlem (ler) inden sonra.
  • Kısıtlama - altörnekleme artık daha kaba bir ızgaraya hata.
  • İnterpolasyon veya uzatma - daha kaba bir ızgarada hesaplanan bir düzeltmenin daha ince bir ızgaraya enterpolasyonu.
  • Düzeltme - Daha ince ızgaraya daha kalın ızgara çözümünün eklenmesi.

Tek bir yinelemeyi çözme hızı ile söz konusu yinelemeyle yakınsama oranı arasında değişen ödünleşimlere sahip birçok multigrid yöntem seçeneği vardır. 3 ana tip V-Döngüsü, F Döngüsü ve W Döngüsüdür. Bir ayrık 2D problemi F-Cycle, V-Cycle yinelemesine göre% 83 daha fazla zaman alırken, W-Cycle yinelemesi% 125 daha fazla sürer. Sorun bir 3D etki alanında kuruluysa, bir F-Döngüsü yinelemesi ve bir W-Döngüsü yinelemesi, bir V-Döngüsü yinelemesinin göz ardı edilmesine göre sırasıyla yaklaşık% 64 ve% 75 daha fazla zaman alır. Genel giderler. Tipik olarak, W-Döngüsü, F-Döngüsüne benzer yakınsama üretir. Ancak, durumlarda konveksiyon-difüzyon yüksek problemler Péclet numaraları W-Cycle, F-Cycle'a göre yineleme başına yakınsama oranında üstünlük gösterebilir. Düzeltme operatörlerinin seçimi, içerdikleri için son derece çeşitlidir Krylov alt uzayı yöntemler ve olabilir önceden koşullandırılmış.

Herhangi bir geometrik multigrid döngü yinelemesi, bir ızgara hiyerarşisi üzerinde gerçekleştirilir ve bu nedenle özyineleme kullanılarak kodlanabilir. İşlev kendisini daha küçük boyutlu (daha kaba) parametrelerle çağırdığından, en kaba ızgara özyinelemenin durduğu yerdir. Sistemin yüksek olduğu durumlarda durum numarası, düzeltme prosedürü, daha ince ızgaraya uzatılmış iri ızgara çözeltisinin sadece bir kısmı eklenecek şekilde değiştirilir.

Bu adımlar, MATLAB tarzı sözde kodda gösterildiği gibi 1 yineleme için kullanılabilir. V-Cycle Multigrid:

işleviphi =V_Cycle(phi, f, h)	Tek tip bir h aralığı ızgarası üzerinde Poisson denklemini ( nabla ^ 2 phi = f) çözmek için% Yinelemeli V-Cycle Multigrid	Ön Düzeltme Yüzdesi	phi = yumuşatma(phi,f,h);		Kalan Hesaplama Hatalarının Yüzdesi	r = artık(phi,f,h);		% Kısıtlama	rhs = kısıtlama(r);	eps = sıfırlar(boyut(rhs));	En küçük ızgara boyutunda yinelemeyi% durdur, aksi takdirde özyinelemeye devam etEğer smallest_grid_size_is_achieved        	eps = yumuşatma(eps,rhs,2*h);Başka         	eps = V_Cycle (eps, rhs, 2 * h);son		% Uzatma ve Düzeltme	phi = phi + uzatma(eps);		% Sonrası Yumuşatma	phi = yumuşatma(phi,f,h);son

Aşağıdaki temsil eder F döngüsü multigrid. Bu multigrid döngüsü, yineleme başına V-Döngüsünden daha yavaştır ancak daha hızlı yakınsama ile sonuçlanır.

işleviphi =F_Cycle(phi, f, h)	Tek tip bir h aralığı ızgarası üzerinde Poisson denklemini ( nabla ^ 2 phi = f) çözmek için% Yinelemeli F-döngüsü multigrid	Ön yumuşatma	phi = yumuşatma(phi,f,h);		Kalan Hesaplama Hatalarının Yüzdesi	r = artık(phi,f,h);		% Kısıtlama	rhs = kısıtlama(r);	eps = sıfırlar(boyut(rhs));	En küçük ızgara boyutunda yinelemeyi% durdur, aksi takdirde özyinelemeye devam etEğer smallest_grid_size_is_achieved        	eps = yumuşatma(eps,rhs,2*h);Başka         	eps = F_Cycle (eps, rhs, 2 * h);son		% Uzatma ve Düzeltme	phi = phi + uzatma(eps);		% Yeniden yumuşatma	phi = yumuşatma(phi,f,h);	Kalan hataları hesapla	r = artık(phi,f,h);		% Kısıtlama	rhs = kısıtlama(r);	En küçük ızgara boyutunda yinelemeyi% durdur, aksi takdirde özyinelemeye devam etEğer smallest_grid_size_is_achieved        	eps = yumuşatma(eps,rhs,2*h);Başka         	eps = V_Cycle (eps, rhs, 2 * h);son		Uzatma ve Düzeltme Yüzdesi	phi = phi + uzatma(eps);		% Sonrası yumuşatma	phi = yumuşatma(phi,f,h);son

Benzer şekilde prosedürler, 1 yineleme için MATLAB tarzı sözde kodda gösterildiği gibi değiştirilebilir. W döngüsü multigrid belirli durumlarda daha da yüksek bir yakınsama oranı için:

işleviphi =W_cycle(phi, f, h)	Tek tip bir h aralığı ızgarası üzerinde Poisson denklemini ( nabla ^ 2 phi = f) çözmek için% Yinelemeli W döngüsü multigrid	Ön yumuşatma	phi = yumuşatma(phi,f,h);		Kalan Hesaplama Hatalarının Yüzdesi	r = artık(phi,f,h);		% Kısıtlama	rhs = kısıtlama(r);	eps = sıfırlar(boyut(rhs));	En küçük ızgara boyutunda yinelemeyi% durdur, aksi takdirde özyinelemeye devam etEğer smallest_grid_size_is_achieved        	eps = yumuşatma(eps,rhs,2*h);Başka         	eps = W_döngü (eps, rhs, 2 * h);son		% Uzatma ve düzeltme	phi = phi + uzatma(eps);		% Yeniden yumuşatma	phi = yumuşatma(phi,f,h);	Kalan hataları hesapla	r = artık(phi,f,h);		% Kısıtlama	rhs = kısıtlama(r);	En küçük ızgara boyutunda yinelemeyi% durdur, aksi takdirde özyinelemeye devam etEğer smallest_grid_size_is_achieved        	eps = yumuşatma(eps,rhs,2*h);Başka         	eps = W_döngü (eps, rhs, 2 * h);son		% Uzatma ve düzeltme	phi = phi + uzatma(eps);		% Sonrası yumuşatma	phi = yumuşatma(phi,f,h);son

Hesaplamalı maliyet

2 boyutlu bir problem kurulumunu varsayarsak, hesaplama çeşitli multigrid döngüleri için grid hiyerarşisi boyunca farklı hareket eder.

Bu yaklaşım, diğer yöntemlere göre, genellikle kullanılan ayrık düğümlerin sayısıyla doğrusal olarak ölçeklenme avantajına sahiptir. Başka bir deyişle, bilinmeyenlerin sayısı ile orantılı olan bir dizi işlemde bu sorunları belirli bir doğrulukla çözebilir.

Bir ızgarada yaklaşık olarak (belirli bir doğrulukla) çözülebilen bir diferansiyel denklem olduğunu varsayalım. belirli bir ızgara noktası yoğunluğuyla . Ayrıca, herhangi bir şebekede bir çözüm olduğunu varsayalım. belirli bir çaba ile elde edilebilir daha kaba bir ızgaradaki bir çözümden . Buraya, "komşu" ızgaralardaki ızgara noktalarının oranıdır ve ızgara hiyerarşisi boyunca sabit olduğu varsayılır ve tek bir ızgara noktası için sonucu hesaplama çabasını sürekli olarak modellemektir.

Daha sonra çözümü grid üzerinde elde etme çabası için aşağıdaki tekrarlama ilişkisi elde edilir. :

Diğer yumuşatma operatörlerine kıyasla Multigrid Döngülerin Yakınsama Oranı. Multigrid, tipik düzleştirme operatörlerinden daha hızlı birleşir. F-Döngüsü ve W-Döngüsü neredeyse eşit sağlamlıkta çalışır.
Diğer yumuşatma operatörleri ile karşılaştırmalı olarak Multigrid Döngülerin Yakınsama Oranlarının Örneği.

Ve özellikle, en iyi ızgarayı buluyoruz o

Bu iki ifadeyi birleştirmek (ve ) verir

Kullanmak Geometrik seriler, sonra buluruz (sonlu için )

yani, bir çözüm elde edilebilir zaman. Bir istisna olduğu belirtilmelidir. ör. 1B problemde kullanılan W döngüsü multigrid; sonuçlanır karmaşıklık.

Multigrid ön koşullandırma

Kasıtlı olarak azaltılmış toleransa sahip bir multigrid yöntemi, verimli bir ön koşullayıcı harici bir yinelemeli çözücü için, ör.[7] Çözüm hala şurada elde edilebilir: çözücü olarak multigrid yönteminin kullanıldığı durumda olduğu gibi zaman. Multigrid ön koşullandırma, pratikte doğrusal sistemler için bile kullanılır, tipik olarak yineleme başına bir döngü ile, örn. Hypre. Tamamen multigrid çözücüye göre ana avantajı, özellikle doğrusal olmayan problemler için açıktır, örn. özdeğer sorunlar.

Orijinal denklemin matrisi veya bir özdeğer problemi simetrik pozitif tanımlı (SPD) ise, ön koşullayıcı genellikle SPD olacak şekilde yapılandırılır, böylece standart eşlenik gradyan (CG) yinelemeli yöntemler hala kullanılabilir. Bu tür empoze edilen SPD kısıtlamaları, ön koşullandırıcının yapımını karmaşıklaştırabilir, örneğin koordineli ön ve son düzleştirme gerektirebilir. Ancak, önceden koşullandırılmış en dik iniş ve esnek CG yöntemleri SPD lineer sistemler için ve LOBPCG simetrik özdeğer problemlerinin hepsi gösterilmiştir[8] ön koşullandırıcı SPD değilse sağlam olacaktır.

Bramble – Pasciak – Xu ön koşullandırıcı

Başlangıçta Xu’nun Ph.D. tez [9]ve daha sonra Bramble-Pasciak-Xu'da yayınlandı,[10] BPX-ön koşullandırıcı, kısmi diferansiyel denklemler tarafından tanımlanan bilim ve mühendislikte modellerin ayrıklaştırılmasından ortaya çıkan büyük ölçekli cebir sistemlerini çözmek için iki ana çoklu-basamaklı yaklaşımdan biridir (diğeri, V-döngüsü gibi klasik multigrid algoritmasıdır). Alt uzay düzeltme çerçevesi göz önüne alındığında,[11] BPX ön koşullandırıcı, klasik V döngüsünün ardışık bir alt uzay düzeltme yöntemi olduğu paralel bir alt uzay düzeltme yöntemidir. BPX ön koşullandırıcının doğal olarak daha paralel olduğu ve bazı uygulamalarda klasik V-döngüsü multigrid yönteminden daha sağlam olduğu bilinmektedir. Yöntem, 1990'dan beri araştırmacılar ve uygulayıcılar tarafından yaygın olarak kullanılmaktadır.

Genelleştirilmiş multigrid yöntemler

Multigrid yöntemleri birçok farklı şekilde genelleştirilebilir. Zaman aşamalı bir çözümde doğal olarak uygulanabilirler. parabolik kısmi diferansiyel denklemler veya doğrudan zamana bağlı olarak uygulanabilir. kısmi diferansiyel denklemler.[12] İçin çok düzeyli teknikler üzerine araştırma hiperbolik kısmi diferansiyel denklemler devam ediyor.[13] Multigrid yöntemleri ayrıca şunlara da uygulanabilir: integral denklemler veya içindeki sorunlar için istatistiksel fizik.[14]

Başka bir çoklu çözünürlük yöntemleri seti temel alır dalgacıklar. Bu dalgacık yöntemleri multigrid yöntemlerle birleştirilebilir.[15][16] Örneğin, dalgacıkların bir kullanımı, sonlu eleman yaklaşımını çok düzeyli bir yöntem açısından yeniden formüle etmektir.[17]

Uyarlanabilir multigrid sergiler uyarlanabilir ağ iyileştirme yani, hesaplamanın kendisine bağlı bir şekilde hesaplama ilerledikçe ızgarayı ayarlar.[18] Buradaki fikir, şebekenin çözünürlüğünü yalnızca ihtiyaç duyulan çözüm bölgelerinde artırmaktır.

Cebirsel multigrid (AMG)

Multigrid yöntemlerin pratik olarak önemli uzantıları, çok düzeyli hiyerarşiyi oluşturmak için hiçbir kısmi diferansiyel denklemin veya geometrik problem arka planının kullanılmadığı teknikleri içerir.[19] Böyle cebirsel multigrid yöntemler (AMG) operatör hiyerarşisini doğrudan sistem matrisinden oluşturur. Klasik AMG'de hiyerarşi seviyeleri, herhangi bir geometrik yorumlama olmaksızın basitçe bilinmeyenlerin alt kümeleridir. (Daha genel olarak, kaba ızgara bilinmeyenleri, ince ızgara bilinmeyenlerinin belirli doğrusal kombinasyonları olabilir.) Bu nedenle, AMG yöntemleri belirli sınıflar için kara kutu çözücüler haline gelir. seyrek matrisler. AMG, esas olarak geometrik multigrid'in uygulanmasının çok zor olduğu durumlarda avantajlı olarak kabul edilir,[20] ancak genellikle sadece gerçek bir multigrid uygulaması için gerekli olan kodlamayı engellediği için kullanılır. İlk olarak klasik AMG geliştirilirken, ilgili bir cebirsel yöntem, düzleştirilmiş toplama (SA) olarak bilinir.

Yeni bir genel bakış makalesinde [21] Jinchao Xu ve Ludmil Zikatanov tarafından "cebirsel multigrid" yöntemler soyut bir bakış açısıyla anlaşılmıştır. Birleşik bir çerçeve geliştirdiler ve mevcut cebirsel multigrid yöntemler tutarlı bir şekilde türetilebilir. Optimal kaba uzayın ve yarı-optimal uzayların nasıl inşa edileceğine dair soyut teori türetildi. Ayrıca, uygun varsayımlar altında, soyut iki seviyeli AMG yönteminin doğrusal sistemin boyutuna, katsayı değişimine ve anizotropiye göre homojen bir şekilde birleştiğini kanıtladılar. Soyut çerçeveleri, klasik AMG, enerji minimizasyon AMG, düzleştirilmemiş ve düzleştirilmiş toplama AMG ve spektral AMG gibi mevcut AMG yöntemlerinin çoğunu kapsar.

Zaman yöntemlerinde multigrid

Multigrid yöntemleri de çözüm için benimsenmiştir. ilk değer problemleri.[22]Burada özellikle ilgi çekici olan, zaman içinde paralel multigrid yöntemleridir:[23]klasik aksine Runge-Kutta veya doğrusal çok adımlı yöntemler sunabilirler eşzamanlılık zamansal yönde. iyi bilinen Parareal Zaman içinde paralel entegrasyon yöntemi de zaman içinde iki seviyeli bir multigrid olarak yeniden formüle edilebilir.

Neredeyse tekil problemler için multigrid

Bir dizi önemli fiziksel ve mühendislik uygulamasında neredeyse tekil problemler ortaya çıkar. Neredeyse tekil problemlerin basit ama önemli bir örneği, yer değiştirme formülasyonunda bulunabilir. doğrusal esneklik neredeyse sıkıştırılamaz malzemeler için. Tipik olarak, bu tür neredeyse tekil sistemleri çözmenin en büyük problemi, aşağıdaki gibi verilen neredeyse tekil operatörü tedavi etmek için kaynar. pozitif, ancak küçük parametreye göre sağlam bir şekilde . Buraya simetrik yarı belirsiz büyük operatör boş alan, süre simetrik pozitif tanımlı Şebeke. Bu tür neredeyse tekil problemler için sağlam ve hızlı bir multigrid yöntemi tasarlamaya çalışan birçok çalışma vardı. Parametrelere ulaşmak için tasarım ilkesi olarak genel bir kılavuz sağlanmıştır (örneğin, ağ boyutu ve Poisson oranı Neredeyse tekil operatörde görülen) bu tür neredeyse tekil sistemlere uygulanan multigrid yönteminin bağımsız yakınsama oranı,[24] yani, her ızgarada, düzleştirmenin uygulandığı bir uzay ayrıştırması, neredeyse tekil operatörün tekil kısmının sıfır uzayının yerel boş uzayların toplamına, kesişimine dahil edilecek şekilde inşa edilmelidir. boş uzay ve uzay ayrışmalarından kaynaklanan yerel uzaylar.

Notlar

  1. ^ Roman Wienands; Wolfgang Joppich (2005). Multigrid yöntemler için pratik Fourier analizi. CRC Basın. s. 17. ISBN  978-1-58488-492-7.
  2. ^ U. Trottenberg; C. W. Oosterlee; A. Schüller (2001). Multigrid. Akademik Basın. ISBN  978-0-12-701070-0.
  3. ^ Yu Zhu; Andreas C. Cangellaris (2006). Elektromanyetik alan modellemesi için çok çizgili sonlu eleman yöntemleri. Wiley. s. 132 ff. ISBN  978-0-471-74110-7.
  4. ^ Şah, Tasneem Mohammad (1989). Multigrid yönteminin analizi (Tez). Oxford Üniversitesi. Bibcode:1989STIN ... 9123418S.
  5. ^ M. T. Heath (2002). "Bölüm 11.5.7 Multigrid Yöntemleri". Bilimsel Hesaplama: Bir Giriş Araştırması. McGraw-Hill Yüksek Öğrenim. s. 478 ff. ISBN  978-0-07-112229-0.
  6. ^ P. Wesseling (1992). Multigrid Yöntemlerine Giriş. Wiley. ISBN  978-0-471-93083-9.
  7. ^ Andrew V Knyazev, Klaus Neymeyr. Lokal olarak optimal blok eşlenik gradyan yönteminde multigrid ön koşullayıcılar kullanarak simetrik özdeğer problemlerinin etkin çözümü. Sayısal Analiz Üzerine Elektronik İşlemler, 15, 38–55, 2003. http://emis.ams.org/journals/ETNA/vol.15.2003/pp38-55.dir/pp38-55.pdf
  8. ^ Henricus Bouwmeester, Andrew Dougherty, Andrew V Knyazev. Eşlenik Gradyan ve En Dik İniş Yöntemleri için Simetrik Olmayan Ön Koşullandırma. Procedia Computer Science, Cilt 51, Sayfa 276–285, Elsevier, 2015. https://doi.org/10.1016/j.procs.2015.05.241
  9. ^ Xu, Jinchao. Çok düzeyli yöntem teorisi. Cilt 8924558. Ithaca, NY: Cornell Üniversitesi, 1989.
  10. ^ Bramble, James H., Joseph E. Pasciak ve Jinchao Xu. "Paralel çok düzeyli ön koşullandırıcılar." Hesaplamanın Matematiği 55, no. 191 (1990): 1–22.
  11. ^ Xu, Jinchao. "Uzay ayrıştırma ve alt uzay düzeltmesiyle yinelemeli yöntemler." SIAM inceleme 34, no. 4 (1992): 581-613.
  12. ^ F. Hülsemann; M. Kowarschik; M. Mohr; U. Rüde (2006). "Paralel geometrik çoklu çizgili". Are Magnus Bruaset'de; Aslak Tveito (editörler). Paralel bilgisayarlarda kısmi diferansiyel denklemlerin sayısal çözümü. Birkhäuser. s. 165. ISBN  978-3-540-29076-6.
  13. ^ Örneğin, J. Blaz̆ek (2001). Hesaplamalı akışkanlar dinamiği: ilkeler ve uygulamalar. Elsevier. s. 305. ISBN  978-0-08-043009-6. ve Achi Brandt ve Rima Gandlin (2003). "Atmosferik Veri Asimilasyonu için Multigrid: Analiz". Thomas Y. Hou'da; Eitan Tadmor (eds.). Hiperbolik problemler: teori, sayısal, uygulamalar: Dokuzuncu Uluslararası Hiperbolik Sorunlar Konferansı'nın 2002 Bildirileri. Springer. s. 369. ISBN  978-3-540-44333-9.
  14. ^ Achi Brandt (2002). "Çok ölçekli bilimsel hesaplama: inceleme". Timothy J. Barth'da; Tony Chan; Robert Haimes (editörler). Çok ölçekli ve çok çözünürlüklü yöntemler: teori ve uygulamalar. Springer. s. 53. ISBN  978-3-540-42420-8.
  15. ^ Björn Engquist; Olof Runborg (2002). "Dalgacık tabanlı sayısal homojenizasyon uygulamaları ile". Timothy J. Barth'da; Tony Chan; Robert Haimes (editörler). Çok Ölçekli ve Çoklu Çözünürlük Yöntemleri. Cilt Hesaplamalı Bilim ve Mühendislikte Ders Notlarının 20'si. Springer. s. 140 ff. ISBN  978-3-540-42420-8.
  16. ^ U. Trottenberg; C. W. Oosterlee; A. Schüller (2001). Multigrid. ISBN  978-0-12-701070-0.
  17. ^ Albert Cohen (2003). Dalgacık Yöntemlerinin Sayısal Analizi. Elsevier. s. 44. ISBN  978-0-444-51124-9.
  18. ^ U. Trottenberg; C. W. Oosterlee; A. Schüller (2001). "Bölüm 9: Adaptive Multigrid". Multigrid. s. 356. ISBN  978-0-12-701070-0.
  19. ^ Yair Shapira (2003). "Cebirsel çoklu çizgili". Matris tabanlı multigrid: teori ve uygulamalar. Springer. s. 66. ISBN  978-1-4020-7485-1.
  20. ^ U. Trottenberg; C. W. Oosterlee; A. Schüller (2001). Multigrid. s. 417. ISBN  978-0-12-701070-0.
  21. ^ Xu, J. ve Zikatanov, L., 2017. Cebirsel multigrid yöntemler. Açta Numerica, 26, s. 591-721.
  22. ^ Hackbusch, Wolfgang (1985). "Parabolik çoklu şebeke yöntemleri". Uygulamalı Bilimler ve Mühendislikte Hesaplama Yöntemleri, VI: 189–197. Alındı 1 Ağustos 2015.
  23. ^ Horton Graham (1992). "Zamana paralel multigrid yöntemi". Uygulamalı Sayısal Yöntemlerde İletişim. 8 (9): 585–595. doi:10.1002 / cnm.1630080906.
  24. ^ Young-Ju Lee, Jinbiao Wu, Jinchao Xu ve Ludmil Zikatanov, Neredeyse Tekil Sistemler için Güçlü Altuzay Düzeltme Yöntemleri, Uygulamalı Bilimlerdeki Matematiksel Modeller ve Yöntemler, Cilt. 17, Sayı 11, s. 1937-1963 (2007)

Referanslar

Dış bağlantılar