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?
Który z nich jest lepszy: FORTRAN lub Python? I chyba w obu przypadkach potrzebujesz Gnuplot, mam rację? Obecnie pracuję na komputerze z systemem Windows. Chciałbym go użyć do uzyskania numerycznych rozwiązań problemów fizyki, w tym symulacji Monte-Carlo, integracji i różnicowania numerycznego, dynamiki molekularnej itp. Widziałem kurs fizyki obliczeniowej, który wprowadza …
Ostatnio porównywałem różne nieliniowe solwery z scipy i byłem pod szczególnym wrażeniem przykładu Newtona-Kryłowa w książce kucharskiej Scipy, w którym rozwiązują równanie różniczkowe drugiego rzędu z nieliniowym terminem reakcji w około 20 wierszach kodu. Zmodyfikowałem przykładowy kod, aby rozwiązać nieliniowe równanie Poissona ( zwane również równaniem Poissona-Boltzmanna , patrz strona …
W swojej pracy naukowej chciałbym być niezależny od oprogramowania komercyjnego. Uzależnienie pakietów komercyjnych, takich jak Matlab i jego zestawy narzędzi, jest dla mnie niezadowalające, ponieważ nie wiem, czy będę mieć dostęp do Matlaba w przyszłości, a ponieważ nie podoba mi się ten język. Dlatego szukam alternatyw. Na szczęście jestem dość …
Zacząłem używać Pythona jako języka programowania do wykonywania wszystkich moich zadań w CFD. Mam bardzo małe doświadczenie w programowaniu. Pochodzę z inżynierii mechanicznej i studiuję na kierunku inżynieria lotnicza. czasami obliczeniowy aspekt CFD staje się bardziej nużący niż manipulowanie równaniami lub robienie matematyki. Jakie są ogólne wytyczne, które przyspieszają uruchamianie …
Moja grupa badawcza koncentruje się na dynamice molekularnej, która oczywiście może generować gigabajty danych w ramach jednej trajektorii, którą należy następnie przeanalizować. Kilka problemów, którymi się zajmujemy, wiąże się z korelacjami w zbiorze danych, co oznacza, że musimy śledzić duże ilości danych w pamięci i analizować je, zamiast stosować podejście …
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 …
Czy ktoś tutaj używał obliczeń naukowych o podwójnej precyzji z procesorami graficznymi nowej generacji (np. K20) za pośrednictwem Pythona? Wiem, że ta technologia szybko się rozwija, ale jaki jest obecnie najlepszy sposób? GPU nie wchodzi w zakres popularnych naukowych bibliotek Pythona numpy i scipy, a ja chciałem użyć theano, ale …
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= …
Piszę powtarzalny artykuł, a artykuł ma wyniki obliczeń generowane przez skrypt Pythona (podobny skrypt MATLAB generuje prawie identyczne wyniki). Wydaje mi się, że artykuł byłby łatwiejszy do zrozumienia dla czytelników, gdyby mogli dopasować obliczenia w artykule do obliczeń w kodzie. Praca proponuje abstrakcyjny formalizm, a przykłady w artykule mają uczynić …
Mam problem z implementacją funkcji numerycznie. Cierpi na tym fakt, że przy dużych wartościach wejściowych wynik jest bardzo dużą liczbą razy bardzo małą liczbą. Nie jestem pewien, czy katastrofalne anulowanie jest właściwym terminem, więc proszę mnie poprawić, jeśli tak jest. Dowody na coś nie tak: Jak mogę uniknąć oscylacji i …
Mam frustrację związaną ze sposobem, w jaki Matlab obsługuje integrację numeryczną vs. Scipy. Poniżej obserwuję następujące różnice w moim kodzie testowym: Wersja Matlaba działa średnio 24 razy szybciej niż mój odpowiednik w Pythonie! Wersja Matlaba jest w stanie obliczyć całkę bez ostrzeżeń, podczas gdy Python powraca nan+nanj Co mogę zrobić, …
Stosowałem różne metody do obliczania zarówno rangi macierzy, jak i rozwiązania układu równań macierzowych. Natknąłem się na funkcję linalg.svd. Porównując to do własnych wysiłków związanych z rozwiązaniem systemu z eliminacją Gaussa, wydaje się on zarówno szybszy, jak i bardziej precyzyjny. Próbuję zrozumieć, jak to jest możliwe. O ile mi wiadomo, …
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 przylegania jest również binarna, co sprawia, że myślę, że istnieje szybszy sposób, …
Muszę rozwiązać s.t.minx∥Ax−b∥22,∑ixi=1,xi≥0,∀i.minx‖Ax−b‖22,s.t.∑ixi=1,xi≥0,∀i.\begin{alignat}{1} & \min_{x}\|Ax - b\|^2_{2}, \\ \mathrm{s.t.} & \quad\sum_{i}x_{i} = 1, \\ & \quad x_{i} \geq 0, \quad \forall{i}. \end{alignat} Myślę , że to kwadratowy problem, który powinien być rozwiązany za pomocą CVXOPT , ale nie potrafię zrozumieć, jak 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.