Pytanie: Jakie metody są dostępne w celu dokładnego i wydajnego obliczenia struktury rzadkości matrycy elementów skończonych? Informacje: Pracuję nad solwerem Poissona Równania Ciśnienia, stosując metodę Galerkina z kwadratową podstawą Lagrange'a, napisaną w C, i używam PETSc do rzadkiego przechowywania macierzy i procedur KSP. Aby efektywnie korzystać z PETSc, muszę wstępnie …
Liczby zmiennoprzecinkowe pojedynczej precyzji zajmują połowę pamięci, a na nowoczesnych komputerach (nawet na GPU) operacje można wykonywać z nimi prawie dwukrotnie szybciej niż podwójna precyzja. Wiele kodów FDTD, które znalazłem, używa wyłącznie arytmetyki i pamięci o pojedynczej precyzji. Czy istnieje ogólna zasada, kiedy dopuszczalne jest stosowanie pojedynczej precyzji do rozwiązywania …
Korzystam z MATLAB, aby rozwiązać problem polegający na rozwiązywaniu za każdym razem, gdy b zmienia się z czasem. Obecnie realizuję to za pomocą MATLAB-a :Ax=bAx=b\mathbf{A} \mathbf{x}=\mathbf{b}bb\mathbf{b}mldivide x = A\b Mam elastyczność, aby wykonać tyle wstępnych obliczeń, ile potrzeba, więc zastanawiam się, czy istnieje szybsza i / lub dokładniejsza metoda niż …
Szukam biblioteki, która wykonuje operacje macierzowe na dużych macierzach rzadkich bez poświęcania stabilności numerycznej. Macierze będą miały wartości 1000+ na 1000+, a wartości macierzy będą zawierać się w przedziale od 0 do 1000. Będę wykonywać algorytm rachunku indeksu, więc będę generował (rzadkie) wektory rzędowe macierzy szeregowo. Gdy rozwijam każdy rząd, …
Pracuję nad biblioteką macierzy zawierającą tylko nagłówki, aby zapewnić pewien rozsądny stopień możliwości algebry liniowej w tak prostym pakiecie, jak to możliwe, i próbuję zbadać, jaki jest obecny stan techniki: obliczanie SVD złożona macierz. Robię dwufazowy rozkład, dwukieragonalizację, a następnie obliczanie wartości osobliwych. W tej chwili używam metody gospodarstwa domowego …
Biorąc pod uwagę dwie macierze i , to, że, aby znaleźć kierunków i takie, że W postaci macierzowej próbuję zminimalizować normę Frobeniusa A - \ mbox {diag} (x) \ cdot B \ cdot \ mbox {diag} (y) = A - B \ circ (xy ^ \ top) .ZAAAbBBxxxyyyA - diag …
Jak na przykład MATLAB oblicza SVD danej macierzy? Zakładam, że odpowiedź prawdopodobnie obejmuje obliczenie wektorów własnych i wartości własnych A*A'. Jeśli tak jest, chciałbym również wiedzieć, jak to oblicza?
Dano mi macierz , która jest symetryczna, odwracalna, dodatnia określona i gęsta. Muszę sprawdzić, czy gdzie J jest macierzą wszystkich jedynek.Q det ( Q ) = det ( 12 I - Q - J )12×1212×1212 \times 12QQQJdet ( Q ) = det ( 12 I- Q - J)( 1 )det(Q)=det(12I−Q−J)(1)\det(Q) …
Poniższe równanie macierzowe w dla danych macierzy i pojawia się w mojej pracy jako charakterystyka macierzy kowariancji. Dowiedziałem się, że to równanie jest znane, zwłaszcza w teorii ciągłego sterowania czasem, jako równanie Lapunowa i że istnieją różne dobrze znane algorytmy jego rozwiązywania, które wykorzystują szczególną naturę tego równania liniowego. Σ …
Biorąc pod uwagę system gdzie , przeczytałem, że w przypadku gdy iteracja Jacobiego jest używana jako solver, metoda nie zbiegnie się, jeśli ma wartość niezerową składnik zerowa przestrzeni . Jak zatem można formalnie stwierdzić, że pod warunkiem, że ma niezerowy składnik obejmujący zerową przestrzeń , metoda Jacobiego jest niespójna? Zastanawiam …
Na przykład biblioteki rzadkich macierzy C ++, których użyłem - Eigen i SuiteSparse, wydają się nie mieć żadnej funkcji SVD dla rzadkiej macierzy. Ciekawe, czy SVD jest trudniejsze niż rzadka matryca QR / LU?
Pracuję nad jakimś źle uwarunkowanym dużym rzadkim liniowym układem równań. Chcę użyć podwójnej podwójnej arytmetyki lub czterokrotnej podwójnej arytmetyki, aby je rozwiązać. Wiem, że istnieje pakiet o nazwie MPACK opracowany przez Nakata, Maho, który może wykonywać numeryczne obliczenia liniowe algebraiczne w oparciu o podwójną arytmetykę. Jest on jednak przeznaczony do …
Jestem nowy w pakiecie PETSc. Mam macierz A ~ 4000x4000 w formacie matrycowym i chcę, aby PETSc rozwiązał ten problem za pomocą wielu procesorów. Wiem, jak rozwiązać system na jednym procesorze, ale nie wiem, jak rozdzielić macierz i wektory między różne procesory. Czy istnieje prosty zestaw instrukcji, aby to zrobić?
Używamy plików cookie i innych technologii śledzenia w celu poprawy komfortu przeglądania naszej witryny, aby wyświetlać spersonalizowane treści i ukierunkowane reklamy, analizować ruch w naszej witrynie, i zrozumieć, skąd pochodzą nasi goście.
Kontynuując, wyrażasz zgodę na korzystanie z plików cookie i innych technologii śledzenia oraz potwierdzasz, że masz co najmniej 16 lat lub zgodę rodzica lub opiekuna.