Pytania otagowane jako algorithms

Opis konkretnych kroków potrzebnych do jednoznacznego rozwiązania określonego problemu, wyrażony w formie abstrakcyjnej.

2
(jak) pisać symulacje, które działają szybciej?
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 …
1
Czy są jakieś wielopoziomowe implementacje ILU oparte na odwrotnych źródłach?
Jestem pod wielkim wrażeniem seryjnej wydajności wielopoziomowych odwracalnych warunków wstępnych ILU , szczególnie w przypadku heterogenicznego Helmholtza , ale jestem zaskoczony, że nie mogłem znaleźć żadnych implementacji typu open source. W szczególności ILUPACK udostępnia naukowo pliki binarne, ale nie wydaje się, aby publikowali swój kod źródłowy. Czy tak naprawdę nikt …
1
Jakie są względne zalety korzystania z Adamsa-Moultona w porównaniu z algorytmem Adamsa-Bashfortha?
Rozwiązuję układ dwóch sprzężonych PDE w dwóch wymiarach przestrzennych iw czasie obliczeniowo. Ponieważ oceny funkcji są drogie, chciałbym zastosować metodę wieloetapową (zainicjowaną przy użyciu Runge-Kutta 4-5). Metoda Adamsa-Bashfortha wykorzystująca pięć poprzednich ocen funkcji ma błąd globalny (jest to przypadek, gdy s = 5 w artykule w Wikipedii, do którego odwołuje …
5
Powtarzane obliczenia najbliższego sąsiada dla milionów punktów danych za wolno
Mam zestaw danych obejmujący miliony punktów danych w 3D. Aby wykonać obliczenia, muszę obliczyć sąsiada (wyszukiwanie zakresu) dla każdego punktu danych w promieniu, spróbować dopasować funkcję, obliczyć błąd dopasowania, powtórzyć to dla następnego punktu danych i tak dalej. Mój kod działa poprawnie, ale jego uruchomienie zajmuje bardzo dużo czasu, około …
2
Stabilny numerycznie sposób obliczania kątów między wektorami
Podczas stosowania klasycznej formuły kąta między dwoma wektorami: α=arccosv1⋅v2∥v1∥∥v2∥α=arccos⁡v1⋅v2‖v1‖‖v2‖\alpha = \arccos \frac{\mathbf{v_1} \cdot \mathbf{v_2}}{\|\mathbf{v_1}\| \|\mathbf{v_2}\|} stwierdzono, że dla bardzo małych / ostrych kątów występuje utrata precyzji, a wynik nie jest dokładny. Jak wyjaśniono w tej odpowiedzi Przepełnienie stosu , jednym rozwiązaniem jest użycie arcus tangens zamiast: α=arctan2(∥v1×v2∥,v1⋅v2)α=arctan⁡2(‖v1×v2‖,v1⋅v2)\alpha = \arctan2 \left(\|\mathbf{v_1} …
5
Jakie są zalety i wady związane z używaniem klas do enkapsulacji algorytmów numerycznych?
Wiele algorytmów wykorzystywanych w obliczeniach naukowych ma inną wewnętrzną strukturę niż algorytmy powszechnie uważane za mniej wymagające matematyki formy inżynierii oprogramowania. W szczególności poszczególne algorytmy matematyczne są zwykle bardzo złożone, często obejmują setki lub tysiące wierszy kodu, ale mimo to nie obejmują żadnego stanu (tj. Nie działają na złożoną strukturę …
13 algorithms 
1
Algorytmy dla dużych rzadkich macierzy całkowitych
Szukam biblioteki, która wykonuje operacje macierzowe na dużych macierzach rzadkich bez poświęcania stabilności numerycznej. Macierze będą miały wartości 1000+ na 1000+, a wartości macierzy będą zawierać się w przedziale od 0 do 1000. Będę wykonywać algorytm rachunku indeksu, więc będę generował (rzadkie) wektory rzędowe macierzy szeregowo. Gdy rozwijam każdy rząd, …
1
Wyliczenie wykresów pochodzących z teselacji Delaunaya w 3D
Czy istnieje algorytm, który wylicza wykresy odpowiadające pewnej teselacji punktów Delaunaya w 3D? Jeśli tak, to czy istnieje wydajna parametryzacja geometrii, która odpowiada dowolnemu „grafowi Delaunaya”? Staram się wyliczyć systematycznie wszystkie stabilne geometrie cząsteczek o określonym składzie bez żadnej wiedzy z zakresu wiązania itp. EDYCJA: Niech będzie zbiorem wykresów z …
1
Numeryczne metody odwracania transformat całkowych?
Próbuję odwrócić liczbowo następującą transformację całkowitą: F(y)=∫∞0yexp[−12(y2+x2)]I0(xy)f(x)d xF(y)=∫0∞yexp⁡[−12(y2+x2)]I0(xy)f(x)dxF(y) = \int_{0}^{\infty} y\exp{\left[-\frac{1}{2}(y^2 + x^2)\right]} I_0\left(xy\right)f(x)\;\mathrm{d}x Więc dla danego F(y)F(y)F(y) muszę aproksymować f( x )f(x)f(x) gdzie: fa( x )f(x)f(x) ifa(y)F(y)F(y) są rzeczywiste i dodatnie(są to ciągłe rozkłady prawdopodobieństwa) x , yx,yx,y są rzeczywiste i dodatnie(są wielkościami) W tej chwili mam bardzo niechlujną …
2
Obliczanie współczynnika Cholesky'ego
Twierdzenie o rozkładzie Cholesky'ego stwierdza zatem, że każda prawdziwa symetryczna macierz dodatnio określona ma rozkład Choleskiego gdzie jest dolną macierzą trójkątną.M = L L ⊤ L.MMMM=LL⊤M=LL⊤M= LL^\topL.LL Biorąc pod uwagę , wiemy już, że istnieją szybkie algorytmy do obliczania jego współczynnika Cholesky'ego .LM.MMLLL Załóżmy teraz, że otrzymałem prostokątną macierz , …
3
Szukam równoległej biblioteki dynamicznych grafów w C ++
Witaj społeczności scicomp, Pracowałem w obszarze algorytmów graficznych z wykorzystaniem frameworków takich jak NetworkX (Python), JUNG i YFiles (Java). Wchodzę teraz w obszar obliczeń równoległych i wysokowydajnych. W przypadku nowego projektu szukam biblioteki grafów C ++ z następującymi funkcjami: ma intuicyjny interfejs, który umożliwia tworzenie algorytmów obsługuje operacje dynamiczne: np. …
3
Algorytm równoległy dla eigensystem matrycy tridiagonalnej
Robię diagonalizację Lanczosa dużej rzadkiej macierzy (~ 2 miliony elementów). Prawie wszystkie kroki w algorytmie Lanzcosa są wykonywane równolegle na GPU, z wyjątkiem diagonalizacji macierzy Lanczosa w celu sprawdzenia zbieżności. W tym celu korzystałem z algorytmu TQLI z receptur numerycznych. Czy istnieją metody znalezienia układu macierzystego macierzy tridiagonalnej, które są …
2
Jak porównuje się koszt obliczeniowy operacji mpi_allgather z operacją gromadzenia / rozpraszania?
Pracuję nad problemem, który można zrównoleglić za pomocą pojedynczej operacji mpi_allgather lub jednej operacji mpi_scatter i jednej operacji mpi_gather. Te operacje są wywoływane w pętli while, więc mogą być wywoływane wiele razy. W implementacji ze schematem MPI_allgather zbieram wektor rozproszony dla wszystkich procesów w celu zduplikowanego rozwiązywania macierzy. W drugiej …
1
Posortuj chmurę punktów w odniesieniu do nieustrukturyzowanej siatki komórek sześciościennych
Pytanie Jak posortowałbyś chmurę punktów w odniesieniu do nieustrukturyzowanej siatki komórek sześciościennych? Każda komórka ma centrum i unikalną etykietę do jej reprezentowania. Zasadniczo istnieją dwa punkty chmurowe (pierwotna chmura punktów i chmura punktów centrów komórek), ale informacje o geometrii komórki (obwiednia) mogą być przydatne, nie jestem pewien. Wyniki Zadałem trochę …
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.