Mam proste pytanie, które jest naprawdę trudne dla Google (oprócz kanonicznego Co każdy informatyk powinien wiedzieć o arytmetyki zmiennoprzecinkowej ). Kiedy należy używać funkcji takich jak log1plub expm1należy używać zamiast logi exp? Kiedy nie należy ich używać? Czym różnią się różne...
Czy ktoś może mi polecić dobry system zarządzania przepływem pracy (WMS), najlepiej w języku Python? Do tej pory korzystałem z GNU Make, ale wprowadza warstwę złożoności, której chcę uniknąć. Dobry WMS powinien mieć następujące funkcje: łatwa integracja z narzędziami wiersza poleceń i skryptami...
Inżynierowie często nalegają na stosowanie lokalnie konserwatywnych metod, takich jak skończona objętość, konserwatywna skończona różnica lub nieciągłe metody Galerkina do rozwiązywania PDE. Co może pójść nie tak przy użyciu metody, która nie jest lokalnie zachowawcza? Okej, więc lokalna ochrona...
Jako osoba posiadająca licencjat z fizyki byłem nieco zgorszony, kiedy zacząłem pracować z symulacjami molekularnymi. Szokiem było odkrycie, że nawet najbardziej szczegółowe i kosztowne obliczeniowo symulacje nie są w stanie odtworzyć ilościowo pełnego zachowania wody z pierwszych zasad. Wcześniej...
Problem Obecnie pracuję nad symulacją elementów skończonych Naviera Stokesa i chciałbym zbadać wpływ różnych parametrów. Niektóre parametry są określone w pliku wejściowym lub za pomocą opcji wiersza poleceń; inne parametry są dostarczane jako flagi w pliku Makefile, więc mój kod musi być...
Tło: Myślę, że mógłbym chcieć przenieść część kodu, który oblicza macierzowe produkty wykładnicze-wektorowe za pomocą metody podprzestrzeni Kryłowa od MATLAB-a do Pythona. (W szczególności funkcja expmvp Jitse Niesen , która wykorzystuje algorytm opisany w tym artykule ). Wiem jednak, że chyba, że...
Ostatnio spotkałem się z sformułowaniem meta-zjawiska : „ dwa są łatwe, trzy są trudne ” (sformułowane w ten sposób przez Federico Poloni), które można opisać następująco: Kiedy sformułowany jest pewien problem dla dwóch podmiotów, jest on stosunkowo łatwy do rozwiązania; jednak algorytm...
Journal of Computational Physics był w przeszłości ważnym rynkiem zbytu dla nauk obliczeniowych, który opublikowałem wcześniej. Z korzyścią dla tych (takich jak ja), którzy podpisali bojkot Elsevier , jakie czasopisma inne niż Elsevier byłyby odpowiednimi miejscami do publikowania artykułów, które...
det(A)≈0det(A)≈0\det(A) \approx 0 Czy odwrotność jest również prawdą? Czy źle uwarunkowana matryca ma prawie zerową determinantę? Oto coś, czego próbowałem w Octave: a = rand(4,4); det(a) %0.008 cond(a)%125 a(:,4) = 1*a(:,1) + 2*a(:,2) =
Chcesz poprawić ten post? Podaj szczegółowe odpowiedzi na to pytanie, w tym cytaty i wyjaśnienie, dlaczego Twoja odpowiedź jest poprawna. Odpowiedzi bez wystarczającej ilości szczegółów mogą być edytowane lub usuwane. Wersja 3.0 standardu MPI formalnie usunęła interfejs...
W statystyce i jej różnych zastosowaniach często obliczamy macierz kowariancji , która jest pozytywnie określona (w rozważanych przypadkach) i symetryczna, dla różnych zastosowań. Czasami potrzebujemy odwrotności tej macierzy do różnych obliczeń (na przykład formy kwadratowe z odwrotnością jako...
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...
To pytanie zostało przeniesione z przepełnienia stosu, ponieważ można na nie odpowiedzieć w programie Computational Science Stack Exchange. Migrował 7 lat temu . Czy są jakieś znane problemy / algorytmy w obliczeniach naukowych, których nie można przyspieszyć przez...
Wiem i użyłem f2py2e do owinięcia starego kodu Fortran 77, ale rozumiem, że nie działa z nowszym kodem Fortran 95. Zbadałem, co powinienem użyć, i natknąłem się na fwrap i G3 f2py, z których żadne nie wydaje się wyjaśniać ich obecnego stanu ani sposobu ich użycia (poza podstawowym...
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...
Twórcy oprogramowania mogą wybrać odpowiednią licencję zgodnie z celem (celami) pracy. Czy ktoś może podać jakieś rekomendacje / doświadczenia dotyczące licencji na oprogramowanie? Jakie są zalety / wady „rozdawania” całej zakodowanej pracy jako kodów open source? Jak radzić sobie z graczami...
Co może pójść nie tak, gdy zastosuje się wstępne metody Kryłowa z KSP ( pakiet solvera liniowego PETSc ) do rozwiązania rzadkiego układu liniowego, takiego jak te uzyskane przez dyskretyzację i linearyzację równań różniczkowych cząstkowych? Jakie kroki mogę podjąć, aby ustalić, co jest nie tak z...
Jednym z głównych problemów, z którymi musimy się zmierzyć w symulacjach molekularnych, jest obliczanie sił zależnych od odległości. Jeśli możemy ograniczyć funkcje siły i odległości do uzyskania równych mocy odległości separacji , możemy po prostu obliczyć kwadrat odległości r 2 = r ⋅ r i nie...
Nie znam zbyt dobrze wspólnych schematów dyskretyzacji dla PDE. Wiem, że Crank-Nicolson jest popularnym schematem dyskretyzacji równania dyfuzji. Czy jest to również dobry wybór na okres doradztwa? Interesuje mnie rozwiązanie równania Reakcja-Dyfuzja-Doradztwo
Czy są jakieś narzędzia zaprojektowane specjalnie do kompresji danych naukowych zmiennoprzecinkowych? Jeśli funkcja jest płynna, istnieje oczywiście duża korelacja między liczbami reprezentującymi tę funkcję, więc dane powinny być dobrze kompresowane. Binarne dane zmiennoprzecinkowe zip /...