Chcę zmodyfikować gęstą kwadratową macierz przejścia w miejscu, zmieniając kolejność kilku jej wierszy i kolumn, używając biblioteki numpy Pythona. Matematycznie odpowiada to pomnożeniu macierzy przez macierz permutacji P i pomnożeniu jej przez P ^ -1 = P ^ T, ale nie jest to uzasadnione...
Jak skalują się tablice Python / Numpy wraz ze wzrostem wymiarów tablicy? Jest to oparte na niektórych zachowaniach, które zauważyłem podczas testowania kodu Python dla tego pytania: Jak wyrazić to skomplikowane wyrażenie za pomocą numpy Problem polegał głównie na indeksowaniu w celu zapełnienia...
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...
Chcę zaimplementować następujące wyrażenie w Pythonie: gdzie x i y są tablicami liczbowymi o rozmiarze n , a k to tablica liczbowa o wielkości n × n . Rozmiar n może wynosić do około 10000, a funkcja jest częścią wewnętrznej pętli, która będzie oceniana wiele razy, więc szybkość jest ważna.xja=...
Zaimplementowałem solver z Eulerem wstecznym w Pythonie 3 (używając numpy). Dla własnej wygody i jako ćwiczenie napisałem również małą funkcję, która oblicza przybliżoną różnicę skończoną różnicy gradientu, aby nie zawsze musiałem określać analitycznie jakobian (jeśli to w ogóle...
Rozwiązuję równania różniczkowe, które wymagają odwrócenia gęstych macierzy kwadratowych. Ta inwersja macierzy zużywa najwięcej mojego czasu obliczeniowego, więc zastanawiałem się, czy używam najszybszego dostępnego algorytmu. Mój obecny wybór to numpy.linalg.inv . Z moich danych liczbowych...
W celu zmniejszenia modelu chcę obliczyć lewe wektory osobliwe związane z - powiedzmy 20 - największymi wartościami osobliwymi macierzy , gdzie N ≈ 10 6 i k ≈ 10 3 . Niestety moja matryca A będzie gęsta bez żadnej struktury.A∈RN,kA∈RN,kA \in \mathbb R^{N,k}N≈106N≈106N\approx 10^6k≈103k≈103k\approx...
Jestem fizykiem, który próbuje modelować charakterystykę prądowo-napięciową złącza nadprzewodnik-nadprzewodnik. Równanie dla tego modelu to: I(V)=1eRn−n∫∞−∞|E|[E2−Δ21]1/2|E+eV|[(E+eV)2−Δ22]1/2[f(E)−f(E+eV)]dEI(V)=1eRn−n∫−∞∞|E|[E2−Δ12]1/2|E+eV|[(E+eV)2−Δ22]1/2[f(E)−f(E+eV)]dE\begin{align} I(V) =...