Büyük bir seyrek matrisin (~ 2 milyon element) Lanczos köşegenleştirmesini yapıyorum. Lanzcos algoritmasındaki adımların neredeyse tamamı, yakınsama olup olmadığını kontrol etmek için Lanczos matrisinin köşegenleştirilmesi dışında GPU'da paralel olarak yapılır. Bunun için Sayısal Tarifler TQLI algoritması kullanıyorum. Üçgen bir matrisin paralel veya kolayca paralelleştirilebilen öz sistemi bulmak için yöntemler var mı? TQLI'nın paralel bir versiyonu var mı?
Yanıtlar:
Seri veya paralel eigensistemleri çözmek için birçok farklı yönteme arayüz içeren SLEPc gibi bir kütüphane kullanmanızı öneririm . Kullanım kılavuzu özdeğer problemleri çözmek için birçok farklı yöntemlerle de söz edilmektedir.
TQL paralelleştirilemez.
Standart paralel algoritma Cuppen algoritmasıdır:
JJM Cuppen, Simetrik üçgensel özdeğer için bir bölme ve fethetme yöntemi, 1980.
http://www.springerlink.com/content/t21365q2gh702714/
Ayrıca bakınız:
F. Tisseur, Dağıtılmış bellek mimarilerindeki simetrik özdeğer problemi için paralel bir bölme ve fethetme algoritması, 1999
http://eprints.ma.man.ac.uk/981/01/covered/MIMS_ep2007_225.pdf
http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.75.4109&rep=rep1&type=pdf
http://www14.in.tum.de/konferenzen/Jass09/courses/2/Kleine_Albers_paper.pdf