Pytania otagowane jako linear-algebra

Pytania dotyczące algorytmicznych / obliczeniowych aspektów algebry liniowej, w tym rozwiązania układów liniowych, problemów z najmniejszymi kwadratami, problemów własnych i innych podobnych kwestii.

4
Kiedy transformacje ortogonalne przewyższają eliminację Gaussa?
Jak wiemy, metody transformacji ortogonalnej (obroty Givensa i odbicia Housholder'a) dla układów równań liniowych są droższe niż eliminacja Gaussa, ale teoretycznie mają lepsze właściwości stabilności w tym sensie, że nie zmieniają numeru stanu układu. Chociaż znam tylko jeden akademicki przykład matrycy, która psuje Gaussowska eliminacja z częściowym przestawieniem. Istnieje powszechna …
2
Biblioteki do rozwiązywania rzadkich układów liniowych
Istnieje wiele różnych bibliotek, które rozwiązują rzadki liniowy układ równań, jednak trudno mi zrozumieć, jakie są różnice. O ile wiem, istnieją trzy główne pakiety: Trilinos , PETSc i Intel MKL . Wszyscy mogą wykonywać rzadkie rozwiązania macierzy, wszystkie są szybkie (o ile mogę stwierdzić, nie byłem w stanie znaleźć solidnych …
3
Czy można rozwiązywać ukośne i stałe symetryczne układy liniowe w czasie kwadratowym po wstępnym obliczeniu?
Czy istnieje metoda do rozwiązania układów liniowych formy gdzie jest stałą macierzą SPD, a są dodatnimi macierzami diagonalnymi?O(n3+n2k)O(n3+n2k)O(n^3+n^2 k)kkk(Di+A)xi=bi(Di+A)xi=bi(D_i + A) x_i = b_iAAADiDiD_i Na przykład, jeżeli każdy jest skalarem, wystarczy obliczyć SVD . Jest to jednak podział na ogólne powodu braku przemienności.DiDiD_iAAADDD Aktualizacja : Jak dotąd odpowiedzi są „nie”. …
1
Aktualizacja diagonalna symetrycznej dodatniej określonej macierzy
AAA jest symetryczną macierzą dodatnią z oznaczeniem dodatnim (SPD). jest rzadką macierzą diagonalną. jest duże ( > 10000), a liczba niezerowych w wynosi zwykle 100 ~ 1000.n×nn×nn \times nGGGnnnnnnGGG AAAZostał factorized w Cholesky'iego postaci jako .LDLTLDLTLDL^T Jak skutecznie zaktualizować i gdy staje się ?LLLDDDAAAA+GA+solA+G
3
Jakich wskazówek powinienem użyć, szukając dobrych metod przygotowania do określonego problemu?
W przypadku rozwiązania dużych układów liniowych metodami iteracyjnymi często interesujące jest wprowadzenie wstępnego kondycjonowania, np. Zamiast tego rozwiąż , gdzie jest tutaj stosowane do lewego wstępnego kondycjonowania układu . Zazwyczaj powinniśmy mieć ten i zapewnić podstawę (znacznie bardziej) wydajnego rozwiązania lub zmniejszenia zasobów obliczeniowych (np. Pamięci) w porównaniu z rozwiązaniem …
5
20% spadek wydajności za ładny projekt oprogramowania
Piszę małą bibliotekę do rzadkich obliczeń macierzowych, aby nauczyć się jak najlepiej korzystać z programowania obiektowego. Bardzo ciężko pracowałem nad stworzeniem ładnego modelu obiektowego, w którym części (rzadkie macierze i wykresy opisujące ich strukturę połączeń) są bardzo luźno połączone. Moim zdaniem kod jest znacznie bardziej rozszerzalny i łatwiejszy w utrzymaniu. …
2
Przestrzeń zerowa prostokątnej gęstej matrycy
Biorąc pod uwagę gęstą macierz jaki jest najlepszy sposób na znalezienie podstawy zerowej w granicach tolerancji ?A∈Rm×n,m>>n;max(m)≈100000A∈Rm×n,m>>n;max(m)≈100000A \in R^{m \times n}, m >> n; max(m) \approx 100000 ϵϵ\epsilon Czy na tej podstawie mogę powiedzieć, że niektóre cols są liniowo zależne od ? Innymi słowy, po obliczeniu podstawy zerowej przestrzeni, jakie …
4
Dlaczego odbicia Householdera nie mogą diagonalizować matrycy?
Przy obliczaniu faktoryzacji QR w praktyce używa się odbić Householdera, aby wyzerować dolną część matrycy. Wiem, że w przypadku obliczania wartości własnych macierzy symetrycznych najlepsze, co można zrobić z odbiciami Householdera, to doprowadzić go do postaci tridiagonal. Czy istnieje oczywisty sposób, aby zobaczyć, dlaczego nie można w pełni przekątnej w …
2
Czy jest jakiś sposób na „podwójne przygotowanie wstępne”
Pytanie: Załóżmy, że masz dwa różne (uwzględnione) warunki wstępne dla symetrycznej dodatniej macierzy : i gdzie odwrotność czynników wynosi łatwy do zastosowania.ZAZAAA ≈ BT.bZA≈bT.bA \approx B^TBA ≈ CT.do,ZA≈doT.do,A \approx C^TC,B , BT., C., C.T.b,bT.,do,doT.B, B^T, C, C^T Gdy jest to możliwe, aby wykorzystać informacje z zarówno i do budowania lepszego …
3
Wydajne obliczanie odwrotności pierwiastka kwadratowego z macierzy
Częstym problemem w statystyce jest obliczanie pierwiastka kwadratowego odwrotnego symetrycznej dodatniej macierzy określonej. Jaki byłby najbardziej efektywny sposób obliczenia tego? Natknąłem pewnym literaturze (które nie zostały jeszcze przeczytane) oraz jakiegoś przypadkowego kodu R tutaj , które będę tutaj dla wygody odtworzenia # function to compute the inverse square root of …
3
metoda wielosiatkowa do rozwiązania PDE
Potrzebuję prostego wyjaśnienia metody wielosiatkowej lub literatury na ten temat. Znam metody iteracyjne, w tym BiCGStab, CG, GS, Jacobi i kondycjonowanie wstępne, ale jestem początkującym w metodzie wielosiatkowej. Czy ktoś może wyjaśnić to szczegółowo lub przynajmniej podać wyraźnie pseudokod lub kod źródłowy, nawet z dobrą literaturą dla początkujących? Dzięki!
1
Czy można zastosować metodę podprzestrzeni Kryłowa jako wygładzającą dla wielu sieci?
O ile mi wiadomo, solwery wielosieciowe używają iteracyjnych wygładzaczy, takich jak Jacobi, Gauss-Seidel i SOR, aby tłumić błąd przy różnych częstotliwościach. Czy zamiast tego można zastosować metodę podprzestrzeni Kryłowa (jak gradient sprzężony, GMRES itp.)? Nie sądzę, że są one klasyfikowane jako „wygładzacze”, ale można ich użyć do przybliżenia rozwiązania zgrubnej …
Korzystając z naszej strony potwierdzasz, że przeczytałeś(-aś) i rozumiesz nasze zasady używania plików cookie i zasady ochrony prywatności.
Licensed under cc by-sa 3.0 with attribution required.