Computational Science

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

2
Zamieszanie na temat Quantum Monte Carlo
Moje pytanie dotyczy wydobycia obserwowalnych metod QMC, jak opisano w tym odnośniku . Rozumiem formalne wyprowadzenie różnych metod QMC, takich jak Path Integral Monte Carlo. Jednak pod koniec dnia wciąż nie rozumiem, jak skutecznie korzystać z tych technik. Podstawową ideą wyprowadzania metod MC kwantowych jest dyskrecja, za pomocą aproksymacji Trottera, …
4
Wydajne pod względem pamięci implementacje częściowych dekompozycji wartości osobliwych (SVD)
W celu zmniejszenia modelu chcę obliczyć lewe wektory osobliwe związane z - powiedzmy 20 - największymi wartościami osobliwymi macierzy , gdzie N ≈ 10 6 i k ≈ 10 3 . Niestety moja matryca A będzie gęsta bez żadnej struktury.A∈RN,kA∈RN,kA \in \mathbb R^{N,k}N≈106N≈106N\approx 10^6k≈103k≈103k\approx 10^3AAA Jeśli po prostu wywołam svdprocedurę …
1
Jakie są różnice między Parareal, PITA i PFASST?
Algorytmy Parareal, PITA i PFASST są technikami obejmującymi całą domenę , służącymi do równoległego rozwiązywania problemów zależnych od czasu w czasie. Jakie są główne zasady tych metod? Jakie są główne różnice między nimi? Czy mogę powiedzieć, że jedna opiera się na innej? W jaki sposób? Co z ich aplikacjami? Wiem, …
2
Nieciągły Galerkin / Poisson / Fenics
Próbuję rozwiązać równanie Poissona 2D za pomocą nieciągłej metody Galerkina (DG) i następującej dyskretyzacji (mam plik png, ale nie mogę go przesłać, przepraszam): Równanie: ∇ ⋅ ( κ ∇ T) + f= 0∇⋅(κ∇T)+f=0\nabla \cdot( \kappa \nabla T) + f = 0 Nowe równania: q= κ ∇ T∇ ⋅ q= - …
1
Cuda i metody numeryczne z domniemaną dyskretyzacją czasu
Chcę zaimportować kod, który rozwiązuje zbiór równań różniczkowych cząstkowych (PDE) metodą objętości skończonej w formie IMPLICIT (dla dyskretyzacji czasu). W wyniku tego powstaje tridiagonalny układ równań w kierunkach x, y, z, który jest obsługiwany przez schemat ADI / TDMA. Nie mogę znaleźć niczego, co dotyczy niejawnego rozwiązania PDE z CUDA. …
2
Kod Octree w Fortranie
Jestem nowy w informatyce naukowej. Szukam implementacji Fortran (najlepiej f90) Octree. Mój problem wymaga Octree, który dzieli moją domenę, dopóki nie będzie więcej niż niektóre cząsteczki N (lub źródła, w których znam wartość gęstości, którą można podłączyć do metody równania całkowego w celu rozwiązania równania takiego jak równanie Laplace'a) w …
4
Programowanie liniowe z ograniczeniami macierzowymi
Mam problem z optymalizacją, który wygląda następująco minJ,Bs.t.∑ij|Jij|MJ+BY=XminJ,B∑ij|Jij|s.t.MJ+BY=X \begin{array}{rl} \min_{J,B} & \sum_{ij} |J_{ij}|\\ \textrm{s.t.} & MJ + BY =X \end{array} Tutaj moje zmienne są macierzami JJJ i BBB , ale cały problem jest nadal programem liniowym; pozostałe zmienne są ustalone. Kiedy próbuję wprowadzić ten program do moich ulubionych narzędzi do …
2
Czy zasada maksimum / minimum równania ciepła jest utrzymywana przez dyskretyzację Cranka-Nicolsona?
Używam schematu różnic skończonych Cranka-Nicolsona do rozwiązania równania cieplnego 1D. Zastanawiam się, czy zasada maksimum / minimum równania ciepła (tj. Że maksimum / minimum występuje w stanie początkowym lub na granicach) również obowiązuje dla rozwiązania dyskretnego. Prawdopodobnie wynika to z faktu, że Crank-Nicolson jest stabilnym i zbieżnym schematem. Wygląda jednak …
2
Rozwiązywanie układu liniowego z argumentami macierzowymi
Wszyscy znamy wiele metod obliczeniowych do rozwiązania standardowego układu liniowego Ax=b.Ax=b. Ax=b. Jestem jednak ciekawy, czy istnieją jakieś „standardowe” metody obliczeniowe do rozwiązania bardziej ogólnego (skończonego wymiaru) układu liniowego formy LA=B,LA=B, LA=B, gdzie, powiedzmy, jest macierzą , jest macierzą , a jest operatorem liniowym przenoszącym macierze do macierzy , co …
5
Definicja przepływu nieściśliwego
Ponieważ wszyscy wiedzą, że przepływ nieściśliwy nie istnieje w rzeczywistości, jego założenie wprowadzono w celu uproszczenia równań rządzących. Nie możemy zastosować tego założenia bezpośrednio. Ogólnie liczba Macha (M <0,3 dla przepływu nieściśliwego), zmiana gęstości (zmiana gęstości zerowej) i rozbieżność prędkości (jest równa zeru dla przepływu nieściśliwego) są wspólnym kryterium definiującym …
3
wyprodukowane rozwiązania dla nieściśliwych Navier-Stokesa - jak znaleźć wolne od dywergencji pola prędkości?
W metodzie wytwarzanych rozwiązań (MMS) postuluje się dokładne rozwiązanie, zastępuje je w równaniach i oblicza odpowiedni termin źródłowy. Rozwiązanie służy następnie do weryfikacji kodu. W przypadku nieściśliwych równań Naviera-Stokesa MMS łatwo prowadzi do (niezerowego) terminu źródłowego w równaniu ciągłości. Ale nie wszystkie kody dopuszczają terminy źródłowe w równaniach ciągłości, więc …
1
w Matlabie, jakie są różnice między Linsolve a miliivide?
w matlabie zarówno linsolve , jak i miliiwide są używane do rozwiązywania układu równań liniowych we wszystkich określonych, zbyt określonych i nieokreślonych przypadkach. Czytając ich dokumenty, zastanawiałem się, jakie są między nimi różnice? Czy używają prawie tych samych algorytmów rozkładu macierzy i trójkątowania w trzech przypadkach? Jeśli A ma właściwoś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.