Computational Science

Pytania i odpowiedzi dla naukowców używających komputerów do rozwiązywania problemów naukowych

2
FeniCS: Wizualizacja elementów wysokiego rzędu
Właśnie zacząłem bawić się w FEniCS. Rozwiązuję Poissona z elementami trzeciego rzędu i chciałbym wizualizować wyniki. Kiedy jednak używam wykresu (u), wizualizacja jest po prostu liniową interpolacją wyników. To samo pojawia się, gdy przesyłam dane do VTK. W innym kodzie, z którym pracuję, napisałem program wyjściowy VTK, który upsamplowałby elementy …
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 …
4
Optymalna metoda ODE dla stałej liczby ocen RHS
W praktyce czas wykonywania rozwiązywania liczbowego IVP x˙(t)=f(t,x(t)) for t∈[t0,t1]x˙(t)=f(t,x(t)) for t∈[t0,t1] \dot{x}(t) = f(t, x(t)) \quad \text{ for } t \in [t_0, t_1] x(t0)=x0x(t0)=x0 x(t_0) = x_0 jest często zdominowane przez czas oceny prawej strony (RHS)fff . Załóżmy zatem, że wszystkie inne operacje są natychmiastowe (tj. Bez kosztów obliczeniowych). …
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 
5
Przykłady obliczeń PDE z wykorzystaniem równoległości w przestrzeni i czasie
W numerycznym rozwiązaniu początkowych wartości granicznych PDE bardzo często stosuje się równoległość w przestrzeni . O wiele rzadziej stosuje się jakąś formę paralelizmu w dyskretyzacji czasu , a paralelizm ten jest zwykle znacznie bardziej ograniczony. Jestem świadomy rosnącej liczby kodów i opublikowanych prac wykazujących równoległość czasową, ale żaden z nich …
6
Przybliżone spektrum dużej matrycy
Chcę obliczyć widmo ( wszystkie wartości własne) dużej rzadkiej macierzy (setki tysięcy wierszy). To jest trudne. Jestem gotów zadowolić się przybliżeniem. Czy istnieją do tego metody przybliżenia? Chociaż mam nadzieję na ogólną odpowiedź na to pytanie, byłbym również zadowolony z odpowiedzi w następującym konkretnym przypadku. Moja matryca jest znormalizowanym Laplacianem …
4
Jaki byłby najbardziej pomocny sposób pisania kodu dla artykułu, aby czytelnicy mogli wyraźnie dopasować wyniki do kodu, który je generuje?
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ć …
4
Warunki brzegowe dla równania doradczego dyskretyzowanego metodą różnic skończonych
Próbuję znaleźć zasoby, które pomogą wyjaśnić, jak wybrać warunki brzegowe podczas korzystania z metod różnic skończonych do rozwiązywania PDE. Książki i notatki, do których mam obecnie dostęp, mówią podobne rzeczy: Ogólne zasady rządzące stabilnością w obecności granic są zdecydowanie zbyt skomplikowane, aby można było wprowadzić tekst wprowadzający; wymagają skomplikowanych mechanizmów …
4
Ilustrujące przykłady metod mimetycznych różnic skończonych
Chociaż staram się znaleźć zwięzłe wyjaśnienie w Internecie, nie mogę pojąć pojęcia mimetycznej skończonej różnicy ani tego, w jaki sposób odnosi się ona do standardowych różnic skończonych. Naprawdę pomocne byłoby zobaczenie kilku prostych przykładów ich implementacji w klasycznych liniowych PDE (hiperbolicznych, eliptycznych i parabolicznych).
1
Algorytm Remeza
Algorytm Remeza jest dobrze znaną procedurą iteracyjną przybliżającą funkcję wielomianem w normie minimax. Ale, jak mówi o tym Nick Trefethen [1]: Większość tych [wdrożeń] sięga wielu lat wstecz, a właściwie większość z nich nie rozwiązuje ogólnego problemu najlepszego przybliżenia przedstawionego powyżej, ale warianty obejmujące zmienne dyskretne lub cyfrowe filtrowanie. W …
2
Jak przydatny jest PETSc do gęstych matryc?
Gdziekolwiek widziałem, tutorial / dokumenty PETSc itp. Mówią, że jest to przydatne w algebrze liniowej i zwykle określa, że ​​korzyści będą miały systemy rzadkie. Co z gęstymi matrycami? Jestem zaniepokojony rozwiązywania dla gęstej A .A x = bAx=bAx=bZAZAA W Fortranie napisałem własny kod dla CG i QMR. Podstawowa konstrukcja jest …
2
Jakie są zalety / wady metod punktów wewnętrznych w porównaniu z metodą simpleks do optymalizacji liniowej?
Jak rozumiem, ponieważ rozwiązanie programu liniowego zawsze występuje w wierzchołku jego wielościennego wykonalnego zestawu (jeśli istnieje rozwiązanie, a optymalna wartość funkcji celu jest ograniczona od dołu, zakładając problem minimalizacji), w jaki sposób można przeszukać wnętrze realnego regionu może być lepsze? Czy zbiega się szybciej? W jakich okolicznościach korzystniejsze byłoby zastosowanie …
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.