Pytania otagowane jako numpy

5
Zezwól na macierz w miejscu w numpy
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 obliczeniowo rozwiązanie. W tej chwili ręcznie …
1
Jak skaluje się wydajność operacji tablicowych Python / Numpy wraz ze wzrostem wymiarów tablicy?
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 tablicy. Odkryłem, że zalety używania (niezbyt dobrych) …
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 …
3
Jak wyrazić to skomplikowane wyrażenie za pomocą numpy
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= …
14 python  numpy 
1
Czy przybliżony jakobian ze skończonymi różnicami może powodować niestabilność w metodzie Newtona?
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 możliwe!). Korzystając z opisów podanych w Ascher i Petzold 1998 , napisałem …
2
Złożoność inwersji macierzy w liczbach
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 wynika, że ​​skaluje się jako gdzie n jest liczbą rzędów, więc metoda wydaje się być …
4
Wydajne pod względem pamięci implementacje częściowych dekompozycji wartości osobliwych (SVD)
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 10^3AAA Jeśli po prostu wywołam svdprocedurę …
1
Całkowanie numeryczne krzywej modelowania dla nadprzewodników (Python)
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) = \frac{1}{eR_{\mathrm{n-n}}}\int_{-\infty}^{\infty}\frac{|E|}{[E^{2} - \Delta_{1}^{2}]^{1/2}}\frac{|E + eV|}{[(E + eV)^{2} - \Delta_{2}^{2}]^{1/2}}[f(E) - f(E + eV)]\,\mathrm{d}E \end{align} Wartości prądu ( lub w kodzie) oblicza się, oceniając tę ​​całkę dla danych napięć ( lub w kodzie).IIIIVVVv Próbowałem …
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.