Computational Science

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

2
Reprezentowanie liczb Eisensteina bez liczb zmiennoprzecinkowych
Mam projekt, w którym muszę użyć pól kwadratowych. Dokładnie numery formularzaa + b- 3---√a+b−3a + b \sqrt{-3} z a , b ∈ Qa,b∈Qa,b \in \mathbb{Q}. Oto na przykład liczby pierwsze w liczbach całkowitych Eisensteina : Nie chcę używać szałwii. Chciałbym napisać własny typ danych do włączenia numpy. PARI byłby użyteczny …
2
Jakie nowe struktury danych są stosowane w adaptacyjnej MES?
Wiele adaptacyjnych bibliotek MES wykorzystuje bardziej zaawansowane struktury danych siatki do obsługi dodawania / usuwania węzłów, krawędzi, trójkątów, czworościanów itp. Na przykład biblioteka p4est używa struktur danych oktree do adaptacyjnego udoskonalania siatki; często nie ma oktetów używanych do obliczeń na siatce statycznej. Jakie zmiany po stronie algebry liniowej dla adaptacyjnego …
1
Wdrożenie metody Jacobi-Davidsona dla problemu wartości własnej sześciennej
Mam duży problem wartości własnej sześciennej: (A0+λA1+λ2A2+λ3A3)x=0.(A0+λA1+λ2A2+λ3A3)x=0.\left(\mathbf{A}_0 + \lambda\mathbf{A}_1 + \lambda^2\mathbf{A}_2 + \lambda^3\mathbf{A}_3\right)\mathbf{x} = 0. Mógłbym to rozwiązać, przechodząc na liniowy problem wartości własnych, ale spowodowałoby to, że układ byłby tak duży:3)2)323^2 ⎡⎣⎢-ZA0000ja000ja⎤⎦⎥⎡⎣⎢xyz⎤⎦⎥= λ⎡⎣⎢ZA1ja0ZA2)0jaZA3)00⎤⎦⎥⎡⎣⎢xyz⎤⎦⎥,[−A0000I000I][xyz]=λ[A1A2A3I000I0][xyz],\begin{bmatrix} -\mathbf{A}_0 & 0 & 0 \\ 0 & \mathbf{I} & 0 \\ 0 & 0 & …
2
Wysiłek obliczeniowy algorytmów
Rozważ ściśle wypukły, nieograniczony problem optymalizacjiO:=minx∈Rnf(x).O:=minx∈Rnf(x).\mathcal{O} := \min_{x \in \mathbb{R}^n} f(x).Niech oznacza jego unikalne minima, a x_0 będzie początkowym przybliżeniem do x_ \ text {opt}. Wywołamy wektor x an \ epsilon - zamknij rozwiązanie \ mathcal {O} if \ begin {equation} \ frac {|| x - x _ {\ text …
3
Jakieś problemy, w których SOR jest szybszy niż Gaussa-Seidela?
Czy jest jakaś prosta zasada, by powiedzieć, czy warto zrobić SOR zamiast Gaussa-Seidela? (i możliwy sposób oszacowania parametru realxation )ωω\omega Mam na myśli po prostu patrząc na matrycę lub znajomość konkretnego problemu, jaki reprezentuje matryca? Czytałem odpowiedź na te pytania: Czy są jakieś heurystyki dla optymalizacji metody sukcesywnej nadmiernej relaksacji …
3
Czy mój projekt może przetrwać bez Object Orientation?
Piszę mały pakiet MATLAB, który rozwiąże pewne problemy numeryczne klasy. Istnieją 3 etapy algorytmu, a użytkownik ma 5 opcji do wyboru dla każdego etapu. Zaimplementowałem cały problem za pomocą funkcji i 3 skrzynek przełączników (po jednym dla każdego etapu algorytmu). Działa dobrze, ale zastanawiam się nad tym, aby zrobić więcej …
9 matlab 
2
sugestia do zarządzania przebiegami symulacji?
Te pytania mogą być nieco nie na temat w comp-sci. jeśli to konieczne, proszę wskazać, gdzie to pasuje. Pytanie dotyczy tego, jak skutecznie zarządzać wszystkimi przebiegami symulacji. powiedzmy, na przykład, że symulacja wymaga ustalenia 2 parametrów, które należy zdefiniować w pewnym sugerowanym zakresie wartości. Aby znaleźć lepszy wynik uzyskany przez …
2
Osobliwy błąd przy rozwiązywaniu równania Poissona w metodzie nieskończonej objętości siatki (tylko 1D)
W ciągu ostatnich kilku dni próbowałem debugować ten błąd. Zastanawiałem się, czy ktoś ma porady, jak postępować. Rozwiązuję równanie Poissona dla rozkładu ładunku krokowego (powszechny problem w elektrostatyce / fizyce półprzewodników) na niejednorodnej siatce o skończonej objętości, w której nieznane są zdefiniowane w centrach komórek i strumienie na powierzchniach komórek. …
4
Uczenie się Fortran
Zastanawiałem się, jaki jest najlepszy sposób na naukę programowania w Fortran? Szukałem książek, przykłady, które wydawały mi się w porządku, to „Wyjaśnienie Fortran 95/2003” i „Fortran dla naukowców i inżynierów” . Nie wiem, co o tym myślicie? Zastanawiam się też, jakiego kompilatora powinienem użyć? Pracuję na komputerze z systemem Windows …
9 fortran 
1
Czy są jakieś zalety liczbowe w rozwiązywaniu macierzy symetrycznej w porównaniu do macierzy bez symetrii?
Stosuję metodę różnic skończonych do układu 3 równań sprzężonych. Dwa równania nie są sprzężone, jednak trzecie równanie łączy się z dwoma pozostałymi. Zauważyłem, że zmieniając porządek równań, powiedzmy od do że macierz współczynników staje się symetryczna.( x , y, z)(x,y,z)(x, y, z)( x , z, y)(x,z,y)(x, z, y) Czy jest …
3
FAS-multigrid wolniej niż liniowa korekcja defektów?
Wdrożyłem solver wielosiatkowy V-Cycle przy użyciu zarówno liniowej korekcji defektów (LDC), jak i pełnego schematu aproksymacji (FAS). Mój problem jest następujący: przy użyciu LDC reszta jest zmniejszana o współczynnik ~ 0,03 na cykl. Implementacja FAS również jest zbieżna z czynnikiem liniowym, ale współczynnik ten wynosi tylko ~ 0,58. Zatem FAS …
9 multigrid 
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.