Próbuję dowiedzieć się, czy istnieje szybszy sposób obliczenia wszystkich wartości własnych i wektorów własnych bardzo dużej i rzadkiej macierzy przylegania niż przy użyciu scipy.sparse.linalg.eigsh O ile mi wiadomo, metody te wykorzystują tylko rzadkość i atrybuty symetrii macierzy. Macierz...
Zarówno w przypadku metod dekompozycji domen (DD), jak i metod wielosiatkowych (MG), można skomponować zastosowanie aktualizacji bloków lub zgrubnych korekt jako addytywne lub multiplikatywne . W przypadku solverów punktowych jest to różnica między iteracjami Jacobiego i Gaussa-Seidela....
Załóżmy, że mam dwie macierze Nx2, Mx2 reprezentujące odpowiednio wektory N, M 2d. Czy istnieje prosty i dobry sposób obliczenia odległości między każdą parą wektorów (n, m)? Łatwy, ale nieefektywny sposób to oczywiście: d = zeros(N, M); for i = 1:N, for j = 1:M, d(i,j) = norm(n(i,:) -...
Obliczenia równoległe są często modelowane przy użyciu deterministycznej lokalnej szybkości obliczeń, obciążenia związanego z opóźnieniami i przepustowości sieci. W rzeczywistości są one zmienne przestrzennie i niedeterministyczne. Badania takie jak Skinner i Kramer (2005) obserwują rozkłady...
Jest to prosty problem z wyliczaniem symetrii. Podaję tutaj pełne tło, ale nie jest wymagana znajomość chemii kwantowej. Integralną dwóch cząstek jest: ⟨ ja j | K L ⟩ = ∫ ψ * I ( x ) ψ * j ( x ' ) ψ k ( x ) ψ L ( x ' )⟨ Ja j | K L ⟩⟨ij|kl⟩\langle ij|kl\rangle I ma następujące 4 symetrie: ⟨ í j...
Rozważ dwa komputery z różnymi konfiguracjami sprzętu i oprogramowania. Podczas uruchamiania dokładnie tego samego seryjnego kodu Naviera-Stokesa na każdej platformie wykonanie x ity odpowiednio dla komputera 1 i 2 wymaga czasu xiy. W tym przypadku jest różnicą czasu iteracji między komputerem 1 a...
Szukam najważniejszej literatury i slajdów do modelowania bieżących i przyszłych kosztów energii operacji zmiennoprzecinkowych i transferu danych w procesorze, pamięci, sieci i pamięci. Oznacziłem to pytanie jako wiki społeczności i wolę ograniczyć każdą odpowiedź do następującego formatu: Tytuł,...
Czy mnożenie macierzy (zarówno Mat * Mat, jak i Mat * Vec) skaluje się z liczbą niezerowych lub z rozmiarem macierzy? Lub jakąś kombinację tych dwóch. Co z kształtem. Na przykład mam macierz 100 x 100 ze 100 wartościami lub macierz 1000 x 1000 ze 100 wartościami. Czy podczas kwadratowania tych...
C ++ 11 wprowadza semantykę przenoszenia, która może na przykład poprawić wydajność kodu w sytuacjach, w których C ++ 03 musiałby wykonać konstrukcję kopii lub przypisanie kopii. W tym artykule opisano, że następujący kod przyspiesza 5-krotnie po skompilowaniu z C + 11: vector<vector<int>...
Muszę dokonać transformacji współrzędnych między dwoma układami odniesienia (osiami). W tym celu należy pomnożyć trzy macierze ( ) ze względu na zastosowanie niektórych osi pośrednich. Pomyślałem o dwóch podejściach do rozwiązania tego:3×33×33\times3 Metoda nr 1 : Wykonanie mnożenia bezpośrednio,...
Próbuję zrozumieć różnicę między tymi dwiema kartami graficznymi dla komputerów akademickich, szczególnie dla komponentu DGEMM. Jeśli spojrzymy na surowe statystyki, oba mają ten sam układ GK110, mają porównywalne statystyki w praktycznie każdej kategorii i, moim zdaniem, mają tę samą...
Mam problem, w którym muszę znaleźć wszystkie pozytywne (jak w wartości własnej dodatniej) pary własne małej (zwykle mniejszej niż 60 x 60) macierzy niesymetrycznej. Mogę przestać obliczać, kiedy wartość własna jest mniejsza niż pewien próg. Wiem, że wartości własne są prawdziwe. Jakieś sugestie...
double trap(double func(double), double b, double a, double N) { double j; double s; double h = (b-a)/(N-1.0); //Width of trapezia double func1 = func(a); double func2; for (s=0,j=a;j<b;j+=h){ func2 = func(j+h); s = s + 0.5*(func1+func2)*h; func1 = func2; } return s; } Powyżej...