Ostatnio napotkałem dziwny problem z FORTRAN95. Zainicjowałem zmienne X i Y w następujący sposób: X=1.0 Y=0.1 Później dodaję je razem i wydrukuję wynik: 1.10000000149012 Po zbadaniu zmiennych wydaje się, że 0,1 nie jest reprezentowane w podwójnej precyzji z pełną dokładnością. Czy jest jakiś...
Po kilku krokach aktualizacji bayesowskiej pozostaję z tylnym rozkładem postaci mieszaniny rozkładów normalnych,Oznacza to, że parametr \ theta jest pobierany z rozkładu, którego plik PDF jest podany jako ważona mieszanina normalnych plików PDF i nie jest sumą normalnych wartości RV. Chciałbym...
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...
Rozwiązuję wiele liczb PDE, ale zastosowana matematyka nie jest moją dziedziną. Nie wybrałem czasopism matematycznych, które powinienem przeczytać, aby nadążyć za najnowszymi osiągnięciami w tej dziedzinie. Jakie są dobre czasopisma do przeczytania, aby nadążyć za najnowszymi osiągnięciami w...
Mam ODE: u′=−1000u+sin(t)u′=−1000u+sin(t)u'=-1000u+sin(t) u(0)=−11000001u(0)=−11000001u(0)=-\frac{1}{1000001} Wiem, że ten konkretny ODE jest sztywny, analitycznie. Wiem również, że jeśli użyjemy jawnej (krokowej) metody krokowej (Euler, Runge-Kutta, Adams itp.), Metoda powinna zwrócić bardzo...
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...
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 }
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ść...
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...
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...
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...
Chcę wiedzieć, które z rozwiązują klasyczny liniowych (np Gauss-Seidel Jacobiego, SOR) gwarantowane są zbieżne do problemu , gdzie jest pozytywne pół definitywna i oczywiścieA x = bZAx=bAx=bZAZAAb ∈ i m ( A )b∈jam(ZA)b \in im(A) (Uwaga jest półokreślona i
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...
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...
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ć...
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 < ( r a di U y / 4...
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< = pjot< = godz jajotloj<=pj<=hijlo_j
Załóżmy, że wiedziałem, że losowa sekwencja liczb została wygenerowana przez liniowy generator kongruencjalny. To jest, xn + 1= ( a Xn+ c ) mod mxn+1=(aXn+c)modmx_{n+1}=(aX_n+c) \bmod m Jeśli jestem biorąc pod uwagę cały okres (lub przynajmniej znaczna przyległe podciąg z nim), w jaki sposób...
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=0traceu=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...
Zastanawiam się, jak znaleźć wartości własne macierzy rzadkich w danym przedziale [a, b] metodą iteracyjną. W moim osobistym rozumieniu bardziej oczywiste jest stosowanie metody podprzestrzeni Kryłowa w celu znalezienia ekstremalnych wartości własnych, a nie