Na moich kursach analizy numerycznej nauczyłem się analizować wydajność algorytmów, licząc liczbę wymaganych operacji zmiennoprzecinkowych (klap) w stosunku do wielkości problemu. Na przykład w tekście Trefethen & Bau na temat numerycznej algebry liniowej są nawet trójwymiarowe zdjęcia liczby flopów. Teraz modne jest stwierdzenie, że „flopy są bezpłatne”, ponieważ opóźnienie pamięci …
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 …
Jak można rozwiązać problem grawitacyjnego n-ciała równolegle numerycznie? Czy możliwy jest kompromis precyzji i złożoności? Jak precyzja wpływa na jakość modelu?
Mam serię punktów danych których spodziewam się (w przybliżeniu) podążania za funkcją która asymptota do linii na dużym . Zasadniczo zbliża się do zera jako x \ do \ infty i to samo można prawdopodobnie powiedzieć o wszystkich pochodnych f '(x) , f' '(x) itp. Ale nie wiem, jaka jest …
Szukam algorytmów do rysowania standardowych wykresów 2D dla funkcji, które mogą, ale nie muszą mieć osobliwości. Celem jest napisanie „Mini-CAS”, więc nie mam a priori wiedzy na temat rodzajów funkcji, które użytkownicy chcą wyświetlać. Ten problem jest bardzo stary, więc wyobrażam sobie, że w literaturze muszą być pewne standardowe algorytmy. …
Czy istnieje metoda do rozwiązania układów liniowych formy gdzie jest stałą macierzą SPD, a są dodatnimi macierzami diagonalnymi?O(n3+n2k)O(n3+n2k)O(n^3+n^2 k)kkk(Di+A)xi=bi(Di+A)xi=bi(D_i + A) x_i = b_iAAADiDiD_i Na przykład, jeżeli każdy jest skalarem, wystarczy obliczyć SVD . Jest to jednak podział na ogólne powodu braku przemienności.DiDiD_iAAADDD Aktualizacja : Jak dotąd odpowiedzi są „nie”. …
Nie mogę znaleźć literatury na temat algorytmów, która mogłaby zostać wykorzystana do rozwiązania problemu uogólnionego przypisania wiele do wielu (GAP), tj. Modeli, w których nie tylko można przypisać więcej zadań do jednego agenta, ale także wielu agentów przypisane do jednego zadania (punkty AP jeden do jednego i jeden do wielu …
Często piszę bardzo podobny kod dla jedno-, dwu- i trójwymiarowych wersji danej operacji / algorytmu. Utrzymanie wszystkich tych wersji może być nudne. Proste generowanie kodu działa dość dobrze, ale wydaje się, że istnieje lepszy sposób. Czy istnieje względnie prosty sposób, aby napisać operację raz i generalizować ją do wyższych lub …
Ostatnia odpowiedź wspomniała o zastosowaniu generatorów liczb losowych Fortuna lub Mersenne Twister ( RNG ) do zaszczepienia symulacji Monte Carlo . Nie słyszałem o Fortunie, więc spojrzałem na nią - wygląda na to, że jest przeznaczona głównie do użytku kryptograficznego. Obecnie używam Mersenne Twister w kodzie produkcyjnym do uruchomienia algorytmu …
W większości (wszystkich?) Implementacji szybkiej metody wielobiegunowej (FMM) do dekompozycji odpowiedniej domeny używa się oktetów. Teoretycznie oktany zapewniają proste wiązanie wolumetryczne, które jest przydatne do udowodnienia czasu działania O (n) FMM. Poza tym teoretycznym uzasadnieniem, czy istnieją korzyści z używania Octree w porównaniu do innych struktur danych drzewa lub trie? …
Załóżmy, że mam równanie macierzowe rekurencyjnie zdefiniowane jako A[n] = inverse([1 - b[n]A[n+1]]) * a[n] Następnie równanie dla A [1] wygląda podobnie do ułamka ciągłego, dla którego istnieje kilka wysoce wydajnych metod, które pozwalają uniknąć żmudnego ponownego obliczania (patrz „Przepisy numeryczne” dla niektórych przykładów). Zastanawiam się jednak, czy istnieją analogiczne …
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?
Opracowywanie algorytmów naukowych jest wysoce iteracyjnym procesem, który często wymaga zmiany wielu parametrów, które będę chciał zmieniać albo w ramach mojego projektu eksperymentalnego, albo w ramach ulepszenia wydajności algorytmu. Jakie strategie mogę zastosować, aby ustrukturyzować te parametry, aby móc je łatwo zmieniać między iteracjami i aby móc łatwo dodawać nowe?
Jestem z dziedziny fizyki akceleratora, szczególnie związanej z okrągłymi pierścieniami magazynowymidla synchrotronowych źródeł światła. Elektrony o wysokiej energii krążą wokół pierścienia, kierowane przez pola magnetyczne. Elektrony krążą miliardy razy i chce się przewidzieć stabilność. Można opisać ruch elektronów w jednym punkcie pierścienia w kategoriach przestrzeni fazowej (pozycja, przestrzeń pędu). Z …
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.