Computational Science

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

6
Jaki jest najlepszy sposób na przeprowadzenie powtarzalnych badań, jeśli potrzebujesz zastrzeżonych bibliotek?
Powtarzalne badania w dziedzinie obliczeń mają na celu udostępnienie kodu potrzebnego do wygenerowania wyników w dokumencie obliczeniowym innym badaczom, aby mogli uruchomić ten kod w celu odtworzenia wyników w tym artykule. Chciałbym, aby wszystkie moje badania były powtarzalne, ale wpadam w drobną przeszkodę: kilka artykułów, nad którymi pracuję, korzysta z …
6
Jak rozpocząć korzystanie z OpenFOAM dla CFD
Patrzę na wykorzystanie OpenFOAM do rozwiązywania podstawowych przepływów wewnętrznych w CFD. Jaki jest najlepszy sposób na rozpoczęcie pracy i czy ktoś mógłby mi wskazać dobre internetowe odnośniki, na które można by odpowiedzieć w przypadku jakichkolwiek pytań, które mogą pojawić się po nurkowaniu? Słyszałem, że na razie jest dość niszowa i …
1
Aktualizacja diagonalna symetrycznej dodatniej określonej macierzy
AAA jest symetryczną macierzą dodatnią z oznaczeniem dodatnim (SPD). jest rzadką macierzą diagonalną. jest duże ( > 10000), a liczba niezerowych w wynosi zwykle 100 ~ 1000.n×nn×nn \times nGGGnnnnnnGGG AAAZostał factorized w Cholesky'iego postaci jako .LDLTLDLTLDL^T Jak skutecznie zaktualizować i gdy staje się ?LLLDDDAAAA+GA+solA+G
3
Jakich wskazówek powinienem użyć, szukając dobrych metod przygotowania do określonego problemu?
W przypadku rozwiązania dużych układów liniowych metodami iteracyjnymi często interesujące jest wprowadzenie wstępnego kondycjonowania, np. Zamiast tego rozwiąż , gdzie jest tutaj stosowane do lewego wstępnego kondycjonowania układu . Zazwyczaj powinniśmy mieć ten i zapewnić podstawę (znacznie bardziej) wydajnego rozwiązania lub zmniejszenia zasobów obliczeniowych (np. Pamięci) w porównaniu z rozwiązaniem …
2
Co to jest pseudo stepping?
Czytając literaturę o rozwiązaniach PDE, spotkałem się dziś z terminem pseudo-stepping . Wydaje się, że jest to powszechny termin, jednak nie znalazłem dla niego dobrej definicji ani artykułu wprowadzającego. Stąd: Co to jest pseudo stepping i jak się go zwykle stosuje?
1
Dlaczego ósemki są używane do rozkładu przestrzeni wielobiegunowej?
W większości (wszystkich?) Implementacji szybkiej metody wielobiegunowej (FMM) do dekompozycji odpowiedniej domeny używa się oktetów. Teoretycznie oktany zapewniają proste wiązanie wolumetryczne, które jest przydatne do udowodnienia czasu działania O (n) FMM. Poza tym teoretycznym uzasadnieniem, czy istnieją korzyści z używania Octree w porównaniu do innych struktur danych drzewa lub trie? …
18 algorithms 
3
Rozwiązywanie nieograniczonych problemów z optymalizacją nieliniową na GPU
Próbuję rozwiązać niektóre nieliniowe problemy optymalizacji nieliniowej na GPU (CUDA). Funkcja celu jest gładką funkcją nieliniową, a jej gradient jest stosunkowo tani do obliczeń analitycznych, więc nie muszę się przejmować przybliżeniem numerycznym. Chcę rozwiązać ten problem głównie z opcjami matematycznymi fp32 (z różnych powodów), więc która nieliniowa metoda optymalizacji jest …
1
Jak zastosować falki do PDE?
Chciałbym dowiedzieć się, w jaki sposób można zastosować metody falkowe do PDE, ale niestety nie znam dobrego zasobu, aby dowiedzieć się na ten temat. Wydaje się, że wiele wstępów do falek koncentruje się na teorii interpolacji, np. Zestawianie sygnału przez superpozycję, najlepiej kilku falek. Czasami wspomina się o aplikacjach do …
18 pde  wavelet 
1
Katastrofalne anulowanie w logarytmie
Próbuję zaimplementować następującą funkcję w zmiennoprzecinkowym podwójnej precyzji z niskim błędem względnym : logsum(x,y)=log(exp(x)+exp(y))logsum(x,y)=log⁡(exp⁡(x)+exp⁡(y))\mathrm{logsum}(x,y) = \log(\exp(x) + \exp(y)) Jest to szeroko stosowane w aplikacjach statystycznych w celu dodania prawdopodobieństw lub gęstości prawdopodobieństwa, które są reprezentowane w przestrzeni dziennika. Oczywiście albo exp(x)exp⁡(x)\exp(x) albo exp(y)exp⁡(y)\exp(y) mogą łatwo przepełnić lub niedopełnić, co byłoby …
2
F2Py z alokowanymi i zakładanymi tablicami kształtów
Chciałbym używać f2pyz nowoczesnym Fortranem. W szczególności próbuję uzyskać następujący podstawowy przykład do działania. To jest najmniejszy użyteczny przykład, jaki mogłem wygenerować. ! alloc_test.f90 subroutine f(x, z) implicit none ! Argument Declarations ! real*8, intent(in) :: x(:) real*8, intent(out) :: z(:) ! Variable Declarations ! real*8, allocatable :: y(:) integer …
18 python  fortran 
4
Czy istnieje biblioteka ogólnego przeznaczenia do ustrukturyzowanego dostosowywania siatki?
Chcesz poprawić ten post? Podaj szczegółowe odpowiedzi na to pytanie, w tym cytaty i wyjaśnienie, dlaczego Twoja odpowiedź jest poprawna. Odpowiedzi bez wystarczającej ilości szczegółów mogą być edytowane lub usuwane. Adaptacyjne udoskonalanie siatki (AMR) jest powszechną techniką rozwiązywania problemu bardzo różnych skal przestrzennych w numerycznym rozwiązaniu PDE. Jakie biblioteki ogólnego …
5
Język programowania równoległego obliczeń naukowych?
Chcę opracować od podstaw równoległe oprogramowanie do obliczeń naukowych. Chcę się zastanowić, od którego języka zacząć. Program obejmuje odczytywanie / zapisywanie danych do plików txt i wykonywanie ciężkich obliczeń równolegle, z wieloma faktoryzacjami LU i wykorzystaniem rzadkich liniowych solverów. Rozwiązania, o których myślałem, to Fortran 2003/2008 z OpenMP lub co-array, …
4
Jakie metody mogą zapewnić, że wielkości fizyczne pozostaną dodatnie podczas symulacji PDE?
Wielkości fizyczne, takie jak ciśnienie, gęstość, energia, temperatura i stężenie powinny zawsze być dodatnie, ale metody numeryczne czasami obliczają wartości ujemne podczas procesu rozwiązywania. Nie jest to w porządku, ponieważ równania obliczą wartości złożone lub nieskończone (zwykle powoduje awarię kodu). Jakie metody numeryczne można zastosować, aby zagwarantować, że te ilości …
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.