Computational Science

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

6
Ograniczenia dotyczące
Przypuszczać minv e c (U)z zastrzeżeniem Uja , j≤ max { Uja , k, Uk , j} ,i , j , k = 1 , … , nminAvec(U)subject to Ui,j≤max{Ui,k,Uk,j},i,j,k=1,…,n\begin{align*} \min A &\mathrm{vec}(U) \\ &\text{subject to } U_{i,j} \leq \max\{U_{i,k}, U_{k,j}\}, \quad i,j,k = 1, \ldots, n \end{align*} gdzie jest …
1
Jak stosować warunki brzegowe, stosując metodę o skończonej objętości?
Zgodnie z moim poprzednim pytaniem próbuję zastosować warunki brzegowe do tej niejednorodnej siatki o skończonej objętości, Chciałbym zastosować warunek brzegowy typu Robin do lhs domeny ( , tak aby:x=xL)x=xL)x=x_L) σL=(dux+au)∣∣∣x=xLσL=(dux+au)|x=xL \sigma_L = \left( d u_x + a u \right) \bigg|_{x=x_L} gdzie jest wartością graniczną; a , d są współczynnikami zdefiniowanymi …
1
Alternatywy dla hdf5
Korzystam z HDF5 od lat, ale wraz ze wzrostem rozmiaru zbioru danych zaczynam doświadczać tych samych problemów wymienionych tutaj http://cyrille.rossant.net/moving-away-hdf5/ Czy możesz wskazać mi format taki jak HDF5 z - Dobra obsługa zapisu równoległego - Obsługa dostępu fragmentarycznego wielkich matryc Mój typowy przypadek użycia to macierz liczb całkowitych 100k x …
2
Opcje rozwiązywania systemów ODE na GPU?
Chciałbym wdrożyć systemy rozwiązywania ODE na GPU w „trywialnie równoległym” ustawieniu. Na przykład, wykonując analizę wrażliwości z 512 różnymi zestawami parametrów. Idealnie chciałbym wykonać rozwiązywanie ODE za pomocą inteligentnego adaptacyjnego solvera pomiaru czasu, takiego jak CVODE, zamiast stałego pomiaru czasu, takiego jak Forward Euler, ale działającego na GPU NVIDIA zamiast …
16 ode  gpu 
3
Konkursy programowania naukowego
Regularnie biorę udział w tak zwanych „konkursach programistycznych”, w których rozwiązujesz trudne problemy algorytmiczne za pomocą własnego kodu i umiejętności rozwiązywania problemów w ograniczonym czasie. Aby zapoznać się z przykładowymi przykładami tego, jak mogą one wyglądać, wyszukaj konkursy takie jak np. Google Code Jam lub ACM-ICPC. (Jeśli wiesz, jakie są …
1
Przydatność elementów o stabilności zależnej od siatki
Po wykonaniu około matematyki związane ze stabilnością elementów problemu 3D Stokesa I nieznacznie szoku sobie sprawę, że nie jest stabilna przez dowolną czworościenne siatki. Mówiąc dokładniej, jeśli masz element, w którym wszystkie węzły i trzy z czterech aspektów leżą na granicy domeny z warunkiem Dirichleta, otrzymujesz pojedynczą macierz. Jest to …
5
Dlaczego numeryczne rozwiązanie ODE odchodzi od niestabilnej równowagi?
Chciałbym zasymulować zachowanie systemu podwójnie wahadłowego. System jest manipulatorem robota o 2 stopniach swobody, który nie jest uruchamiany i dlatego będzie zachowywał się głównie jak podwójne wahadło pod wpływem grawitacji. Jedyną główną różnicą w przypadku podwójnego wahadła jest to, że składa się on z dwóch sztywnych ciał o właściwościach masy …
4
Wybieranie najbardziej rozproszonych punktów z zestawu punktów
Czy istnieje (skuteczny) algorytm do wybierania podzbioru punktów z zestawu punktów ( ) tak, aby „obejmowały” większość obszaru (we wszystkich możliwych podzbiorach rozmiaru )?M.M.MN.N.NM.&lt; NM.&lt;N.M < NM.M.M Zakładam, że punkty są w płaszczyźnie 2D. Naiwny algorytm jest prosty, ale zbyt skomplikowany pod względem złożoności czasowej: for each subset of N …
5
Algorytmy optymalizacji równoległej dla problemu z bardzo kosztowną funkcją celu
Optymalizuję funkcję 10-20 zmiennych. Zła wiadomość jest taka, że ​​ocena każdej funkcji jest kosztowna, około 30 minut obliczeń szeregowych. Dobrą wiadomością jest to, że mam do dyspozycji klaster z kilkudziesięcioma węzłami obliczeniowymi. Zatem pytanie: czy są dostępne algorytmy optymalizacji, które pozwolą mi efektywnie wykorzystać całą moc obliczeniową? Po jednej stronie …
7
Solidne obliczenie średniej z dwóch liczb zmiennoprzecinkowych?
Pozwól x, ybyć dwóch liczb zmiennoprzecinkowych. Jaki jest właściwy sposób na obliczenie ich średniej? Naiwny sposób (x+y)/2może doprowadzić do przepełnienia kiedy xi ysą zbyt duże. Myślę, że 0.5 * x + 0.5 * ymoże lepiej, ale wiąże się to z dwoma multiplikacjami (co może być nieefektywne) i nie jestem pewien, …
2
Czy jest jakiś sposób na „podwójne przygotowanie wstępne”
Pytanie: Załóżmy, że masz dwa różne (uwzględnione) warunki wstępne dla symetrycznej dodatniej macierzy : i gdzie odwrotność czynników wynosi łatwy do zastosowania.ZAZAAA ≈ BT.bZA≈bT.bA \approx B^TBA ≈ CT.do,ZA≈doT.do,A \approx C^TC,B , BT., C., C.T.b,bT.,do,doT.B, B^T, C, C^T Gdy jest to możliwe, aby wykorzystać informacje z zarówno i do budowania lepszego …
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.