Pytania otagowane jako sparse-matrix

Pytania związane z przechowywaniem, montażem, operacjami i innymi aspektami radzenia sobie z rzadkimi macierzami, dla których przechowywane są tylko niezerowe elementy. Pytania, które nie dotyczą bezpośrednio macierzy rzadkich, ale inne sposoby użycia rzadkości należy oznaczyć tagiem [operator-rzadki].

17
Czy istnieje wysokiej jakości nieliniowy solver programowania dla Pythona?
Mam kilka trudnych, niewypukłych problemów globalnej optymalizacji do rozwiązania. Obecnie używam MATLAB's Optimization Toolbox (konkretnie fmincon()z algorytmem = 'sqp'), co jest dość skuteczne . Jednak większość mojego kodu znajduje się w języku Python i chciałbym również przeprowadzić optymalizację w języku Python. Czy istnieje solver NLP z powiązaniami Pythona, z którym …
4
Jakich wskazówek należy przestrzegać przy wyborze rzadkiego solvera systemu liniowego?
W aplikacjach pojawiają się rzadkie systemy liniowe. Do rozwiązania tych systemów jest wiele procedur do wyboru. Na najwyższym poziomie istnieje przełom między metodami bezpośrednimi (np. Rzadka eliminacja Gaussa lub rozkład Choleskiego, ze specjalnymi algorytmami porządkowania i metodami wielopłaszczyznowymi) i iteracyjnymi (np. GMRES, gradient (sprzężony)). Jak określić, czy zastosować metodę bezpośrednią, …
7
Jaki jest najszybszy sposób obliczenia największej wartości własnej ogólnej macierzy?
EDYCJA: Testuję, czy jakieś wartości własne mają wartość jednego lub więcej. Muszę znaleźć największą absolutną wartość własną dużej, rzadkiej, niesymetrycznej macierzy. Korzystam z eigen()funkcji R , która korzysta z algo QR z EISPACK lub LAPACK, aby znaleźć wszystkie wartości własne, a następnie używam, abs()aby uzyskać wartości bezwzględne. Jednak muszę to …
3
Odległość euklidesowa w oktawie
Chciałbym wiedzieć, czy istnieje szybki sposób na obliczenie odległości euklidesowej dwóch wektorów w oktawie. Wydaje się, że nie ma do tego żadnej specjalnej funkcji, więc czy powinienem po prostu użyć formuły z sqrt?
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. …
4
Jak zmienić kolejność zmiennych, aby uzyskać pasmową macierz minimalnej przepustowości?
Próbuję rozwiązać równanie Poissona 2D na podstawie różnic skończonych. W tym procesie otrzymuję rzadką macierz z tylko zmiennymi w każdym równaniu. Na przykład, jeśli zmienne byłyby , dyskretyzacja dałaby:555UUU Ui−1,j+Ui+1,j−4Ui,j+Ui,j−1+Ui,j+1=fi,jUi−1,j+Ui+1,j−4Ui,j+Ui,j−1+Ui,j+1=fi,jU_{i-1,j} + U_{i+1,j} -4U_{i,j} + U_{i,j-1} + U_{i,j+1} = f_{i,j} Wiem, że mogę rozwiązać ten system za pomocą metody iteracyjnej, ale …
1
Dlaczego SciPy eigsh () wytwarza błędne wartości własne w przypadku oscylatora harmonicznego?
Rozwijam jakiś większy kod do wykonywania obliczeń wartości własnych ogromnych rzadkich macierzy w kontekście fizyki obliczeniowej. Moje procedury sprawdzam na prostym oscylatorze harmonicznym w jednym wymiarze, ponieważ wartości własne są dobrze znane analitycznie. Robiąc to i porównując własne procedury z wbudowanymi rozwiązaniami SciPy, natknąłem się na osobliwość pokazaną na poniższym …
2
Oblicz wszystkie wartości własne bardzo dużej i bardzo rzadkiej macierzy przylegania
Mam dwa wykresy z prawie n ~ 100000 węzłów każdy. Na obu wykresach każdy węzeł jest podłączony dokładnie do 3 innych węzłów, więc macierz przylegania jest symetryczna i bardzo rzadka. Najtrudniejsze jest to, że potrzebuję wszystkich wartości własnych macierzy przylegania, ale nie wektorów własnych. Mówiąc dokładniej, będzie to raz w …
3
Czy algorytm Thomasa jest najszybszym sposobem rozwiązania symetrycznego, ukośnego, rzadkiego, tridiagonalnego układu liniowego
Zastanawiam się, czy algorytm Thomasa jest najszybszym (możliwym do udowodnienia?) Rozwiązaniem symetrycznego, dominującego po przekątnej, rzadkiego systemu tridiagonalnego pod względem złożoności algorytmicznej (nie szukając pakietów implementacyjnych takich jak LAPACK itp.). Wiem, że zarówno algorytm Thomasa, jak i multigrid mają złożoność , ale może stały współczynnik dla multigrid jest mniejszy? Nie …
3
Rzadki liniowy solver dla wielu prawych stron
Muszę rozwiązać ten sam rzadki układ liniowy (300 x 300 do 1000 x 1000) z wieloma prawymi bokami (300 do 1000). Oprócz tego pierwszego problemu chciałbym również rozwiązać różne systemy, ale z tymi samymi niezerowymi elementami (tylko różne wartości), to znaczy wiele rzadkich systemów o stałym wzorcu rzadkości. Moje macierze …
4
Czy jest jakiś pakiet podwójnej arytmetycznej macierzy rzadkich?
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 …
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.