Computational Science

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

5
Zwiększenie trwałości archiwalnej kodu
Czy istnieje opublikowana lista najlepszych praktyk zapewniających długowieczność kodu, z myślą o odtwarzalnych wynikach naukowych? (np. open source, praktyki dokumentacji, wybór zależności, wybór języka, maszyny wirtualne itp.). Wiedz o wszelkich badaniach (lub ich brak, przykładach / anegdotach), które próbowały oszacować okres półtrwania typowego kodu naukowego lub innego oprogramowania (czy to …
3
W jakich okolicznościach integracja Monte Carlo jest lepsza niż quasi-Monte Carlo?
Proste pytanie: wykonać całkę wielowymiarową, biorąc pod uwagę, że zdecydowano, że jakaś metoda Monte Carlo jest odpowiednia, czy istnieje jakakolwiek zaleta, że ​​regularna integracja MC przy użyciu liczb pseudolosowych ma przewagę nad integracją quasi-Monte Carlo przy użyciu sekwencji quasirandom ? Jeśli tak, to jak rozpoznałbym sytuacje, w których ta korzyść …
2
Przewodniki po języku Python dotyczące programowania równoległego w pamięci współużytkowanej
Mam doświadczenie w kodowaniu OpenMP dla maszyn z pamięcią współużytkowaną (zarówno w C, jak i FORTRAN), aby wykonywać proste zadania, takie jak dodawanie macierzy, mnożenie itp. (Aby zobaczyć, jak konkuruje z LAPACK). Wiem, że OpenMP wystarcza do wykonywania prostych zadań bez konieczności przeglądania dokumentacji. Ostatnio przeszedłem na Python dla moich …
1
W jaki sposób w przypadku oprogramowania przesłanego do ACM TOMS umowa licencyjna na oprogramowanie ACM współdziała z innymi licencjami?
Czasopismo Association for Computing Machinery Transactions on Mathematical Software (ACM TOMS) publikuje wiele artykułów na temat algorytmów numerycznych obejmujących implementacje oprogramowania. Zgodnie z ich polityką redakcyjną , przesłanie artykułu z algorytmem zawiera kod źródłowy dla implementacji algorytmu opisanego we wspomnianym artykule. Ten kod źródłowy podlega Umowie o prawie autorskim i …
1
Rozwiązywanie ogromnego gęstego układu liniowego?
Czy jest jakaś nadzieja na skuteczne rozwiązanie następującego układu liniowego za pomocą metody iteracyjnej? A∈Rn×n,x∈Rn,b∈Rn, with n>106A∈Rn×n,x∈Rn,b∈Rn, with n>106A \in \mathbb{R}^{n \times n}, x \in \mathbb{R}^n, b \in \mathbb{R}^n \text{, with } n > 10^6 Ax=bAx=bAx=b z Δ - 6 6 1A=(Δ−K)A=(Δ−K) A=(\Delta - K) , gdzie jest bardzo rzadką …
2
Jak poprawić dokładność metody różnic skończonych w celu znalezienia systemu eigens pojedynczej liniowej ODE
Próbuję rozwiązać równanie typu: (−∂2∂x2−f(x))ψ(x)=λψ(x)(−∂2∂x2−f(x))ψ(x)=λψ(x) \left( -\tfrac{\partial^2}{\partial x^2} - f\left(x\right) \right) \psi(x) = \lambda \psi(x) Gdzie f(x)f(x)f(x) ma prosty biegun o wartości 000 , dla najmniejszych NNN wartości własnych i wektorów własnych. Warunki brzegowe to: ψ(0)=0ψ(0)=0\psi(0) = 0 i ψ(R)=0ψ(R)=0\psi(R)=0 , a ja patrzę tylko na funkcję ponad (0,R](0,R](0,R] . …
3
Biblioteki do rozwiązywania równania Lapunowa
Poniższe równanie macierzowe w dla danych macierzy i pojawia się w mojej pracy jako charakterystyka macierzy kowariancji. Dowiedziałem się, że to równanie jest znane, zwłaszcza w teorii ciągłego sterowania czasem, jako równanie Lapunowa i że istnieją różne dobrze znane algorytmy jego rozwiązywania, które wykorzystują szczególną naturę tego równania liniowego. Σ …
1
Jak wykryć krotność wartości własnych?
Załóżmy, że A jest ogólną macierzą rzadką i chcę obliczyć wartości własne. Nie wiem, jak wykryć wielokrotność wartości własnych. O ile mi wiadomo, w specjalnym przypadku znalezienia wielomianowych pierwiastków metodą macierzy towarzyszącej, możemy zastosować RRQR w celu wykrycia wielokrotności pierwiastków.
1
Wystawianie zerowej przestrzeni
Biorąc pod uwagę system gdzie , przeczytałem, że w przypadku gdy iteracja Jacobiego jest używana jako solver, metoda nie zbiegnie się, jeśli ma wartość niezerową składnik zerowa przestrzeni . Jak zatem można formalnie stwierdzić, że pod warunkiem, że ma niezerowy składnik obejmujący zerową przestrzeń , metoda Jacobiego jest niespójna? Zastanawiam …
2
Jak interpolować dane wielopunktowe do centrów komórkowych nieustrukturyzowanej siatki?
Mam zestawy danych pola wielopunktowego, każdy zestaw danych punktowych dotyczy pojedynczej komórki nieustrukturyzowanej siatki. Celem jest interpolacja danych do centrum komórki, bezpośrednio lub pośrednio, w jak najdokładniejszy sposób. Jeśli użyję interpolacji z odwrotnym ważeniem odległości, w przypadku, gdy odległość między źródłem a celem (centrum komórki) jest bardzo mała, może dojść …
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.