Paralelní metody pro výpočet spektra matic

školitel: doc. Ing. Tomáš Oberhuber, Ph.D.
e-mail: zobrazit e-mail
typ práce: bakalářská práce, diplomová práce
zaměření: MI_MM, MI_AMSM, MINF
klíčová slova: spektrum matric, numerická matematika, GPU, paralelní algoritmy
popis: Výpočet spektra matic – tedy jejich vlastních čísel a vlastních vektorů – je zásadní součástí mnoha algoritmů a vědeckých aplikací. Tato úloha má mnoho aplikací např. v: 1. Kvantové a jaderné fyzice (řešení operátorů, spektrální rozklady), 2. analýze stability dynamických systémů (např. lineární aproximace), 3. strojovém učení (např. PCA, kernelové metody, spektrální klasifikace), 4. strukturální mechanice (vlastní módy vibrací), 5. grafové analýze (např. Laplaceovo spektrum, Fiedlerův vektor). V oblasti semidefinitního programování (SDP) je klíčové ověřování, zda matice splňuje podmínku pozitivní semidefinitnosti, což znamená, že všechna její vlastní čísla musí být nezáporná. Pro velké řídké matice je výpočet spektra výpočetně velmi náročný a bez efektivní paralelizace prakticky neřešitelný. To platí zejména pro úlohy v oblasti SDP, kde se pozitivně semidefinitní podmínky uplatňují nad maticemi o rozměrech stovky až tisíce – často v rámci iterativního algoritmu. Cíle práce: 1. Prostudovat základní a pokročilé metody pro výpočet spektra: Jacobiho metoda, QR algoritmus, Lanczosova metoda, divide-and-conquer, Jacobi–Davidson apod. 2. Zaměřit se na metody vhodné pro kompletní spektrum středně velkých hermitovských/symetrických matic. 3. Implementovat paralelní variantu vybrané metody (např. na CPU pomocí OpenMP nebo na GPU pomocí CUDA/TNL). 4. Ověřit výkonnost na testovacích maticích ze SDP úloh (např. ze solverů nebo testovacích sad). 5. Porovnat s referenčními knihovnami (např. LAPACK, Eigen, cuSolver, SLEPc). Přínosy pro studenta: 1. Získá hluboké znalosti v oblasti lineární algebry a numerických metod. 2. Naučí se vytvářet efektivní kód v jazyce C++ včetně využití moderních vlastností tohoto jazyka, včetně programování s pomocí knihovny TNL s možností využití počítání na GPU. 3. Vyzkouší si práci se spektrálními daty v reálných aplikacích (SDP, ML, výpočetní fyzika).
literatura: 1. Saad, Yousef. Numerical methods for large eigenvalue problems: revised edition. Society for Industrial and Applied Mathematics, 2011. 2. Golub, Gene H., and Charles F. Van Loan. Matrix computations. JHU press, 2013. 3. Gu, Ming, and Stanley C. Eisenstat. \"A divide-and-conquer algorithm for the symmetric tridiagonal eigenproblem.\" SIAM Journal on Matrix Analysis and Applications 16.1 (1995): 172-191.
naposledy změněno: 20.05.2025 09:07:57

za obsah této stránky zodpovídá: Pavel Strachota | naposledy změněno: 9.9.2021
Trojanova 13, 120 00 Praha 2, tel. +420 770 127 494
České vysoké učení technické v Praze | Fakulta jaderná a fyzikálně inženýrská | Katedra matematiky