Computational Science

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

2
Co z tym prostym oszacowaniem błędu dla liniowego PDE?
Niech ΩΩ\Omega być wypukła wielobocznie ograniczony Lipschitz domeny w R2R2\mathbb R^2 , niech f∈L2(Ω)f∈L2(Ω)f \in L^2(\Omega) . Δu=fΔu=f\Delta u = fΩΩ\Omegatraceu=0trace⁡u=0\operatorname{trace} u = 0∂Ω∂Ω\partial\OmegaH2H2H^2CCC∥u∥H2≤C∥f∥L2‖u‖H2≤C‖f‖L2\|u\|_{H^2} \leq C \|f\|_{L^2} Dla niektórych przybliżeń elementów skończonych , powiedzmy, z elementami węzłowymi na jednolitej siatce, mamy oszacowanie błęduuhuhu_h ∥u−uh∥H1≤Ch∥u∥H2‖u−uh‖H1≤Ch‖u‖H2\| u - u_h \|_{H^1} \leq C …

4
Nieliniowe najmniejszych kwadratów z ograniczeniami ramek
Jakie są zalecane sposoby wykonywania nieliniowych najmniejszych kwadratów, min ∑ e r rja( p )2)∑erri(p)2\sum err_i(p)^2 , z ograniczeniami ramek l ojot&lt; = pjot&lt; = godz jajotloj&lt;=pj&lt;=hijlo_j <= p_j <= hi_j ? Wydaje mi się (głupcy się spieszą), że można by uczynić kwadratowe ograniczenia kwadratowymi i zminimalizować ∑jae r rja( …

1
Czy można użyć schematu numerycznego do określenia poprawności problemu wartości początkowej lub brzegowej?
Wiem, że możemy zastosować techniki analizy matematycznej, aby udowodnić, czy IVP lub BVP ma rozwiązanie, jest unikalne i zależy w sposób ciągły od wartości granicznych / początkowych. W przypadku niektórych PDE, szczególnie nieliniowych pde, bardzo trudno jest, jeśli nie niemożliwe, udowodnić dobrą postawę. Czy istnieje jakaś technika numeryczna pozwalająca sprawdzić, …

3
Numeryczna integracja kompaktowo obsługiwanej funkcji na trójkącie
jak sugeruje tytuł, próbuję obliczyć całkę kompaktowo obsługiwanej funkcji (kwintyczny wielomian Wendlanda) na trójkącie. Zauważ, że środek funkcji znajduje się gdzieś w przestrzeni 3D. Integruję tę funkcję na dowolnym, ale małym trójkącie ( a r e a &lt; ( r a di U y / 4 )2)2)area&lt;(radius/4)22area < \frac{(radius/4)^2}{2} ). …
10 quadrature 

2
Śledź Isoline z drogiej funkcji 2D
Mam problem podobny do sformułowania do tego postu, z kilkoma znaczącymi różnicami: Jakie są proste metody adaptacyjnego próbkowania funkcji 2D? Jak w tym poście: Mam a ocena tej funkcji jest dość droga do obliczeniafa( x , y)fa(x,y)f(x,y) W przeciwieństwie do tego postu: Nie interesuje mnie wartość funkcji dokładnie wszędzie, ale …

1
Dlaczego musimy ponownie uruchomić solver CFD dla wyższej liczby Reynoldsa?
Zacząłem uczyć się OpenFOAM z samouczka Cavity, który jest dostępny na stronie internetowej . Podczas eksperymentowania z różnymi liczbami Reynoldsa w sekcji „2.1.8.2 Uruchamianie kodu” samouczek mówi, aby ponownie uruchomić solver, ponieważ „Rozsądne jest zwiększenie czasu rozwiązania”. Ale kiedy to zrobiłem, nie mogłem znaleźć żadnej różnicy między przepływem we wnęce …


4
Jak możemy ocenić wyniki uczniów na kursach informatyki?
Jako ktoś, kto musi prowadzić kursy informatyki, mam do czynienia z odwiecznym pytaniem: w jaki sposób oceniam zdolność uczniów do uczenia się przedmiotu zależnego od aplikacji trudnych do przetestowania za pomocą „standardowych” metod testowania ( egzaminy pisemne lub ustne)? Część kursu zależy od zrozumienia teorii i metod na poziomie abstrakcyjnym, …
10 education 

2
Czym różni się programowanie geometryczne od programowania wypukłego?
Czym różni się (uogólnione) programowanie geometryczne od ogólnego programowania wypukłego? Program geometryczny można przekształcić w program wypukły i zazwyczaj rozwiązuje się go metodą punktu wewnętrznego. Ale jaka jest przewaga nad bezpośrednim sformułowaniem problemu jako programu wypukłego i rozwiązaniem go metodą punktu wewnętrznego? Czy klasa programów geometrycznych stanowi jedynie podzbiór klasy …

2
Wybór solwera liniowego do obliczeń GPGPU (OpenCL)
Opracowałem już działające rozwiązanie metody elementów skończonych do rozwiązywania problemów związanych z przenoszeniem ciepła za pomocą GPU i OpenCL przy użyciu metody gradientu sprzężonego. Główną wadą tej metody jest duże zapotrzebowanie na pamięć. Co więcej, w przypadku kart graficznych pamięć jest często bardzo ograniczona. Widzę dwie opcje: Twórz poddomeny i …

3
Macierz wykładniczy prawdziwej macierzy asymetrycznej z Fortran 95 i LAPACK
Niedawno zadałem pytanie w tym samym stylu dla matryc skośno-hermitowskich. Zainspirowany sukcesem tego pytania i po kilku godzinach uderzenia głową o ścianę patrzę na wykładniczą macierz prawdziwych matryc asymetrycznych. Droga do znalezienia wartości własnych i wektorów własnych wydaje się dość skomplikowana i obawiam się, że się zgubiłem. Tło: Jakiś czas …

3
Dlaczego iteracyjne rozwiązywanie równań Hartree-Focka prowadzi do zbieżności?
W metodzie samodzielnego pola Hartree-Focka rozwiązywania niezależnego od czasu elektronicznego równania Schroedingera dążymy do zminimalizowania energii stanu podstawowego, , układu elektronów w polu zewnętrznym w odniesieniu do wyboru orbitali spinowych, { χ i } .E0E0E_{0}{χi}{χi}\{\chi_{i}\} Czynimy to poprzez iteracyjne rozwiązywanie 1 elektronów Hartree-Focka gdzie x i jest wirowanie / przestrzenne …


1
Jak złożyć i rozwiązać system macierzowy równolegle z wartościami generowanymi w różnych procesorach?
Rozwiązuję problem wieloskalowy za pomocą heterogenicznej metody wieloskalowej (HMM) . Zasadniczo moja konkretna procedura wykorzystuje następujący proces iteracyjny: Rozwiąż wiele lokalnych systemów macierzowych. Oblicz wartość zainteresowania z rozwiązań lokalnych systemów. Złóż globalny system macierzy na podstawie lokalnych „wartości zainteresowania” Rozwiąż globalny system macierzy Skorzystaj z rozwiązania globalnego systemu macierzy, aby …

2
Jaki jest najskuteczniejszy sposób obliczenia wektora własnego gęstej macierzy odpowiadającego wartości własnej największej wielkości?
Mam gęstą prawdziwą symetryczną macierz kwadratową. Wymiar wynosi około 1000 x 1000. Muszę obliczyć pierwszy główny składnik i zastanawiać się, jaki może być najlepszy algorytm. Wygląda na to, że MATLAB korzysta z algorytmów Arnoldi / Lanczos (dla eigs). Ale po przeczytaniu o nich nie jestem pewien, czy mają one jakąkolwiek …

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.