Computational Science

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

2
Jakie metody integracji czasowej powinniśmy stosować w przypadku hiperbolicznych PDE?
Jeśli zastosujemy metodę linii do dyskretyzacji (odrębna dyskretyza czasu i przestrzeni) hiperbolicznych PDE, uzyskamy po dyskretyzacji przestrzennej naszą ulubioną metodą numeryczną (fx. Finite Volume Method), czy w praktyce ma znaczenie to, który solver ODE zastosujemy do dyskretyzacji czasowej (TVD / SSP / etc)? Dodano kilka dodatkowych informacji: Problem z dokładnością …
3
Matematycznie dlaczego działa zbieranie macierzy masy / wektora obciążenia?
Wiem, że ludzie często zastępują spójne macierze masy bryłami diagonalnymi. W przeszłości zaimplementowałem również kod, w którym wektor obciążenia jest składany w skupieniu, a nie w sposób zgodny z MES. Ale nigdy nie zastanawiałem się, dlaczego wolno nam to robić. Jaka intuicja kryje się za zbijaniem, która pozwala zastosować ją …
2
Projektowanie wzorców w pisaniu oprogramowania numerycznego w C ++
Szukam zasobów na temat wzorców projektowych i zasad oprogramowania numerycznego, potencjalnie, ale niekoniecznie, ze szczególnym uwzględnieniem obiektowego podejścia do kodów numerycznych. Zdaję sobie sprawę z tekstu, Pisanie oprogramowania naukowego: Przewodnik po dobrym stylu , ale wciąż dotyczy to wielu rzeczy, takich jak problemy z liczbą zmiennoprzecinkową, jak i kiedy optymalizować …
3
Czy pliki HDF5 nadają się do kontroli wersji git?
Nie znam formatu pliku używanego w HDF5, ale zastanawiam się, czy pliki HDF5 nadają się do kontroli wersji za pomocą git (lub na przykład Mercurial lub Subversion)? Wydaje mi się, że mam na myśli: czy pliki HDF5 nadają się do różnicowania na podstawie linii, czy też git będzie musiał traktować …
1
Dlaczego przypięcie punktu do usunięcia pustego miejsca jest złe?
Równanie Poissona ze wszystkimi warunkami brzegowymi Neumanna ma pojedynczą stałą przestrzeń zerową. Podczas rozwiązywania metodą Kryłowa, pustą przestrzeń można usunąć, odejmując średnią z każdego rozwiązania iteracji lub przypinając wartość pojedynczego wierzchołka. Przypinanie pojedynczego wierzchołka ma tę zaletę, że jest proste, a także pozwala uniknąć dodatkowej globalnej redukcji na rzut. Jednak …
3
Dopasowywanie powierzchni niejawnych do zestawów punktów zorientowanych
Mam pytanie dotyczące dopasowania kwadratowego do zbioru punktów i odpowiadających normalnych (lub równoważnie stycznych). Dopasowanie kwadratowych powierzchni do danych punktowych jest dobrze zbadane. Niektóre prace są następujące: Bezpośrednie dopasowanie powierzchni kwadratowych o ograniczonym typie , James Andrews, Carlo H. Sequin Komputerowe projektowanie i zastosowania, 10 (a), 2013, bbb-ccc Algebraiczne dopasowanie …
1
Czy równanie doradcze o zmiennej prędkości może być zachowawcze?
Próbuję nieco lepiej zrozumieć równanie doradcze ze zmiennym współczynnikiem prędkości. W szczególności nie rozumiem, jak to równanie może być konserwatywne. Równanie adwekcja , ∂u∂t+∂∂x(vu)=0∂u∂t+∂∂x(vu)=0 \frac{\partial u}{\partial t} + \frac{\partial}{\partial x}(\boldsymbol{v}u) = 0 Zinterpretujmy u(x,t)u(x,t)u(x,t) jako koncentrację niektórych gatunków fizycznych ( cm−3cm−3cm^{-3} ) lub innej wielkości fizycznej, której nie można stworzyć …
2
Jaki jest najbardziej efektywny sposób pisania pętli „for” w Matlabie?
Czytałem, że jeśli na przykład mam podwójną forpętlę, która przebiega nad indeksami macierzy, umieszczenie indeksu pracy kolumny w zewnętrznej pętli jest bardziej wydajne. Na przykład: a=zeros(1000); for j=1:1000 for i=1:1000 a(i,j)=1; end end Jaki jest najskuteczniejszy sposób kodowania, jeśli mam trzy lub więcej forpętli? Na przykład: a=zeros(100,100,100); for j=1:100 for …
2
obserwacje punktowe vs. ciągłe w odwrotnym problemie PDE
Pracuję nad odwrotnym problemem dla mojego doktoratu. badania, które dla uproszczenia powiemy, określają wββ\beta L ( β) u ≡ - ∇ ⋅ ( k0miβ∇ u ) = fL(β)u≡−∇⋅(k0eβ∇u)=fL(\beta)u \equiv -\nabla\cdot(k_0e^\beta\nabla u) = f od kilku obserwacji ; k 0 jest stałą i F są znane. Zwykle jest to formułowane jako …
2
Algebraiczna multigrid: Dlaczego iloczyn interpolacji i ograniczenia nie prowadzi do czegoś z normą 1?
Obecnie pracuję z „A Multigrid Tutorial” Briggsa i in., Rozdział 8. Konstrukcja operatora interpolacji jest podana jako: Następnie konstrukcję operatora ograniczenia i operatora drobnej sieci podano jako: Załóżmy, że mamy trzy punkty siatki x0, x1, x2 ze środkowym jeden x1 jest w porządku, a pozostałe są gruboziarniste. Środkowy jest interpolowany …
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.