Pytania otagowane jako lapack

LAPACK (Linear Algebra PACKage) jest powszechnie używaną biblioteką podprogramów dla zadań numerycznej algebry liniowej, w tym rozwiązań liniowych układów równań, liniowych najmniejszych kwadratów, problemów z wartościami własnymi i rozkładu wartości osobliwych. Procedury LAPACK mogą być używane z językami Fortran, C i krewnymi oraz wieloma innymi językami.


3
Dlaczego moje skalowanie mnożenia macierzy-wektora?
Przepraszam za długi post, ale chciałem załączyć wszystko, co uważałem za istotne za pierwszym razem. Czego chcę Wdrażam równoległą wersję Krystalicznych metod podprzestrzeni dla gęstych matryc. Głównie GMRES, QMR i CG. Zdałem sobie sprawę (po profilowaniu), że moja procedura DGEMV była żałosna. Postanowiłem więc skoncentrować się na tym, izolując to. …



3
Zrozumienie, jak Numpy robi SVD
Stosowałem różne metody do obliczania zarówno rangi macierzy, jak i rozwiązania układu równań macierzowych. Natknąłem się na funkcję linalg.svd. Porównując to do własnych wysiłków związanych z rozwiązaniem systemu z eliminacją Gaussa, wydaje się on zarówno szybszy, jak i bardziej precyzyjny. Próbuję zrozumieć, jak to jest możliwe. O ile mi wiadomo, …

2
rozwiąż
Przenoszę istniejący kod z MATLAB do C ++ i mam liniowy system do rozwiązania (zamiast bardziej typowej formy A x = b )x A = bxZA=bxA=bA x = bZAx=bAx=b Matryca jest gęsta i ma ogólną formę, ale nie jest większa niż 1000 x 1000. Tak więc w MATLAB rozwiązaniem jest …

2
Jaka jest odpowiednia funkcja LAPACK za Matlabem [Q, R, E] = qr (A)?
Obecnie próbuje tanio obliczyć szacunkową dobry rangi dla macierzy . Dlatego obliczam kolumnę obrotową dekompostację QR za pomocąZAAA [Q,R,E]=qr(A) w Matlabie. Oceniam rangę za pomocąZAAA tol = size(A,n)*eps*norm(A,'fro'); r = sum(abs(diag(R))>tol) Działa to dobrze, a wykres wszystkich przekątnych wpisów R wygląda następująco: W przypadku przeniesienia całego algorytmu do C / …
12 matlab  lapack 


4
Wykładnicza macierz skośno-pustelniczej macierzy z fortranem 95 i LAPACK
Właśnie wciągam się w fortran 95 dla niektórych symulacji mechaniki kwantowej. Szczerze mówiąc, zostałem zepsuty przez Octave, więc wziąłem potęgowanie macierzy za pewnik. Biorąc pod uwagę (małe, ) pochylenie -Hermycka macierz o rozmiarze n × n , jaki jest najbardziej efektywny sposób użycia LAPACK do rozwiązania tego problemu? Nie używam …
11 fortran  lapack 



3
Macierz wykładniczy prawdziwej macierzy asymetrycznej z Fortran 95 i LAPACK
Niedawno zadałem pytanie w tym samym stylu dla matryc skośno-hermitowskich. Zainspirowany sukcesem tego pytania i po kilku godzinach uderzenia głową o ścianę patrzę na wykładniczą macierz prawdziwych matryc asymetrycznych. Droga do znalezienia wartości własnych i wektorów własnych wydaje się dość skomplikowana i obawiam się, że się zgubiłem. Tło: Jakiś czas …

3
Jak zacząć korzystać z LAPACK w c ++?
Jestem nowy w informatyce i nauczyłem się już podstawowych metod integracji, interpolacji, metod takich jak RK4, Numerov itp. Na c ++, ale ostatnio mój profesor poprosił mnie o nauczenie się, jak używać LAPACK do rozwiązywania problemów związanych z macierzami. Jak na przykład znajdowanie wartości własnych złożonej macierzy. Nigdy nie korzystałem …
10 lapack 

3
Z czego korzysta LAPACK
Procedura QR LAPACK przechowuje Q jako reflektory Householdera. Skaluje wektor odbicia pomocą , więc pierwszy element wyniku staje się , więc nie trzeba go przechowywać. I przechowuje osobny wektor , który zawiera potrzebne współczynniki skali. Zatem macierz reflektorów jest taka:vvv1/v11/v11/v_1111ττ\tauH=I−τvvT,H=I−τvvT,H=I-\tau v v^T, gdzie nie jest znormalizowane. Natomiast w podręcznikach matryca …

1
Algorytm równoważenia macierzy
Pisałem przybornik systemu sterowania od zera i wyłącznie w Python3 (bezwstydna wtyczka:) harold. Na podstawie moich wcześniejszych badań zawsze narzekałem na solver Riccati care.mz powodów technicznych / nieistotnych. Dlatego piszę własny zestaw procedur. Jednej rzeczy, której nie mogę znaleźć, jest uzyskanie algorytmu równoważenia o wysokiej wydajności, przynajmniej tak dobrego balance.m. …

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.