Computational Science

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

2
rozwiąż
Przenoszę istniejący kod z MATLAB do C ++ i mam liniowy system do rozwiązania (zamiast bardziej typowej formy A x = b )x A = bxZA=bxA=bA x = bZAx=bAx=b Matryca jest gęsta i ma ogólną formę, ale nie jest większa niż 1000 x 1000. Tak więc w MATLAB rozwiązaniem jest …
3
Rzadki liniowy solver dla wielu prawych stron
Muszę rozwiązać ten sam rzadki układ liniowy (300 x 300 do 1000 x 1000) z wieloma prawymi bokami (300 do 1000). Oprócz tego pierwszego problemu chciałbym również rozwiązać różne systemy, ale z tymi samymi niezerowymi elementami (tylko różne wartości), to znaczy wiele rzadkich systemów o stałym wzorcu rzadkości. Moje macierze …
5
Wielokrotnie rozwiązując
Korzystam z MATLAB, aby rozwiązać problem polegający na rozwiązywaniu za każdym razem, gdy b zmienia się z czasem. Obecnie realizuję to za pomocą MATLAB-a :Ax=bAx=b\mathbf{A} \mathbf{x}=\mathbf{b}bb\mathbf{b}mldivide x = A\b Mam elastyczność, aby wykonać tyle wstępnych obliczeń, ile potrzeba, więc zastanawiam się, czy istnieje szybsza i / lub dokładniejsza metoda niż …
1
Jak dokładnie działa * pełny * algorytm wielosiatkowy?
Więc rozumiem (a przynajmniej tak mi się wydaje), jak przebiega cykl V. Napisałem w Matlabie 1-D, rekurencyjną wersję cyklu V. Jednak kiedy uruchomiłem kod dla FMG, moje rozwiązanie nie było zbieżne. Wierzę, że mój problem polega na zrozumieniu faktycznej części FMG. To, co obecnie wiem, to: Tuż przed FMG interpolacji …
12 multigrid 
2
Strategie metody Newtona, gdy jakobian przy rozwiązaniu jest osobliwy
Próbuję rozwiązać następujący układ równań dla zmiennych i (wszystkie pozostałe są stałymi):x 2P., x1P.,x1P,x_1x2)x2)x_2 A ( 1 - P)2)- k1x1= 0A P.2)- k2)x2)= 0( 1 - P) ( r1+ x1)4L.1- P( r1+ x2))4L.2)= 0ZA(1-P.)2)-k1x1=0ZAP.2)-k2)x2)=0(1-P.)(r1+x1)4L.1-P.(r1+x2))4L.2)=0\frac{A(1-P)}{2}-k_1x_1=0 \\ \frac{AP}{2}-k_2x_2=0 \\ \frac{(1-P)(r_1+x_1)^4}{L_1}-\frac{P(r_1+x_2)^4}{L_2}=0 Widzę, że mogę przekształcić ten układ równań w pojedyncze równanie jednej zmiennej …
2
Jakie są dobre typy danych dla nieustrukturyzowanego kodu FVM CFD skoncentrowanego na komórkach?
Interesuje mnie rada dotycząca wydajnych struktur danych do przeglądania komórek w niestrukturalnej CFD o skończonej objętości opartej na komórkach. Jeden przykład, który napotkałem (w kodzie dolfyn cfd) wygląda następująco (pokażę odpowiedni segment) Mamy więc tablicę NFaces, w której przechowywana jest liczba twarzy dla każdej komórki. Następnie tablica CFace, która odwzorowuje …
1
Solwery PDE dla dyfuzji dryfu i modeli pokrewnych
Próbuję symulować podstawowe modele półprzewodnikowe do celów pedagogicznych - zaczynając od modelu dyfuzyjnego Drifta. Chociaż nie chcę używać gotowego symulatora półprzewodników - nauczę się innych (powszechnych, najnowszych lub niejasnych) modeli, ale chcę używać gotowego solvera PDE. Ale nawet w prostym przypadku 1D model dyfuzyjno-dyfuzyjny składa się z szeregu sprzężonych nieliniowych …
12 pde 
2
Integracja numeryczna - obsługa NaNs (C / Fortran)
Mam do czynienia z trudną całką, która wykazuje wartości NaN przy pewnych wartościach zbliżonych do zera, a w tej chwili radzę sobie z nimi dość brutalnie, używając instrukcji ISNAN, która ustawia intandand na zero, gdy to nastąpi. Próbowałem tego z biblioteką NMS w FORTRAN (procedura q1da - q1dax nie różni …
12 quadrature 
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
Jak zintegrować ekspresję wielomianową z 4-węzłowym elementem 3D?
Chcę zintegrować wyrażenie wielomianowe z elementem 4-węzłowym w 3D. Kilka książek na temat MES dotyczy przypadku, w którym integracja jest przeprowadzana na dowolnym płaskim 4-elementowym elemencie. Zwykłą procedurą w tym przypadku jest znalezienie macierzy Jacobiego i użycie jej wyznacznika do zmiany podstawy całkowania na znormalizowaną, w której mam prostsze granice …
4
Arbitralna dokładność Symulacja skalowalnej liny
Próbuję symulować obiekt linowy. Rozumiem, że sformułowanie to szereg cząstek połączonych sprężynami. Sprężyny te mają bardzo duże wartości k, tak że linia odkształca się, ale bardzo mało się rozciąga. Doszedłem do wniosku, że rozwiązanie tego w funkcji czasu nie jest możliwe w formie zamkniętej, ponieważ lina jest uogólnieniem wahadła (które …
12 simulation 
3
Heurystyczna kontrola stabilności numerycznej
Załóżmy, że mam funkcję o wartościach rzeczywistych niektórych zmiennych które chcę ocenić liczbowo. Zasadniczo wzór na może zawierać produkty, racjonalności, funkcje transcendentalne itp. I będzie musiał długo badać analityczną stabilność numeryczną. Albo będzie to przynajmniej czasochłonne, aby zrobić to w praktyce. Załóżmy, że nie mam krótszego odpowiednika ze stabilnością gwarantowaną. …
12 stability 
2
Wartość bezwzględna w ograniczeniach liniowych
Mam następujący problem z optymalizacją, w którym mam bezwzględną wartość w moich ograniczeniach: x∈Rnx∈Rn\mathbf{x} \in \mathbb{R}^nf0,f1,…,fmf0,f1,…,fm\mathbf{f}_0, \mathbf{f}_1, \ldots, \mathbf{f}_mnnnmins.t.fT0x|fT1x|≤|fT2x|≤…≤|fTmx|minf0Txs.t.|f1Tx|≤|f2Tx|≤…≤|fmTx|\begin{align} \min &\mathbf{f}_0^T \mathbf{x} \notag \\ \text{s.t.} &|\mathbf{f}_1^T \mathbf{x}| \leq |\mathbf{f}_2^T \mathbf{x}| \leq \ldots \leq |\mathbf{f}_m^T \mathbf{x}| \end{align} Wiem, że możliwa przestrzeń nie będzie wypukła i prawdopodobnie będę potrzebować MILP, aby rozwiązać …
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.