Metoda, która tworzy sekwencję przybliżeń numerycznych, które są zbieżne (pod warunkiem spełnienia warunków technicznych) do rozwiązania problemu, zazwyczaj poprzez wielokrotne stosowanie pewnej procedury. Przykłady obejmują metodę Newtona do znajdowania pierwiastków i iterację Jacobiego do rozwiązywania macierzy-wektorów.
Według mojej wiedzy istnieją 4 sposoby rozwiązania układu równań liniowych (popraw mnie, jeśli jest więcej): Jeśli macierz systemowa jest kwadratową matrycą pełnego rzędu, można użyć reguły Cramera; Oblicz odwrotność lub pseudoinwersję macierzy systemowej; Użyj metod rozkładu macierzowego...
Co może pójść nie tak, gdy zastosuje się wstępne metody Kryłowa z KSP ( pakiet solvera liniowego PETSc ) do rozwiązania rzadkiego układu liniowego, takiego jak te uzyskane przez dyskretyzację i linearyzację równań różniczkowych cząstkowych? Jakie kroki mogę podjąć, aby ustalić, co jest nie tak z...
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...
O ile mi wiadomo, solwery wielosieciowe używają iteracyjnych wygładzaczy, takich jak Jacobi, Gauss-Seidel i SOR, aby tłumić błąd przy różnych częstotliwościach. Czy zamiast tego można zastosować metodę podprzestrzeni Kryłowa (jak gradient sprzężony, GMRES itp.)? Nie sądzę, że są one klasyfikowane...
Zastanawiam się, co stało się z wielomianowymi warunkami wstępnymi. Interesuję się nimi, ponieważ wydają się być względnie eleganckie z matematycznego punktu widzenia, ale o ile czytałem w ankietach na temat metod kryłowa, generalnie wypadają bardzo słabo jako warunki wstępne. Według słów Saada i...
Prowadzę zajęcia z analizy numerycznej i szukam motywacji do metody BFGS dla studentów z ograniczonym zapleczem / intuicją w optymalizacji! Chociaż nie mam czasu, aby rygorystycznie udowodnić, że wszystko się zbiega, staram się uzasadnić, dlaczego może pojawić się aktualizacja Hesji BFGS....
Według Wikipedii współczynnik konwergencji wyraża się jako konkretny stosunek norm wektorowych. Próbuję zrozumieć różnicę między szybkościami „liniowymi” i „kwadratowymi” w różnych punktach czasu (w zasadzie „na początku” iteracji i „na końcu”). Czy można stwierdzić,...
tło Rozwiązuję wariant równania Ornsteina-Zernike z teorii płynów. Abstrakcyjnie problem można przedstawić jako rozwiązanie problemu punktu stałego , gdzie A jest operatorem alegbraicznym integro, a c ( r ) jest funkcją rozwiązania (funkcja korelacji bezpośredniej OZ). Rozwiązuję za pomocą...
Załóżmy, że miałem problem z wartością graniczną: dud2udx2+dvdx=f in Ωd2udx2+dvdx=f in Ω\frac{d^2u}{dx^2} + \frac{dv}{dx}=f \text{ in } \Omega u=h in ∂Ωdudx+d2vdx2=g in Ωdudx+d2vdx2=g in Ω\frac{du}{dx} +\frac{d^2v}{dx^2} =g \text{ in } \Omega u=h in ∂Ωu=h in ∂Ωu=h \text{ in } \partial\Omega Moim...
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ą...
Załóżmy, że jest symetryczną, dodatnią określoną macierzą. A jest na tyle duże, że rozwiązanie A x = b jest drogie .A∈Rn×nA∈Rn×nA\in\mathbb{R}^{n\times n}AAAAx=bAx=bAx=b Czy istnieje iteracyjny algorytm znajdowania najmniejszej wartości własnej , który nie obejmuje odwracania A w każdej...
Rozważ dwa komputery z różnymi konfiguracjami sprzętu i oprogramowania. Podczas uruchamiania dokładnie tego samego seryjnego kodu Naviera-Stokesa na każdej platformie wykonanie x ity odpowiednio dla komputera 1 i 2 wymaga czasu xiy. W tym przypadku jest różnicą czasu iteracji między komputerem 1 a...
W informatyce często spotykamy duże układy liniowe, które musimy rozwiązać za pomocą niektórych (skutecznych) środków, np. Metod bezpośrednich lub iteracyjnych. Jeśli skupimy się na tym drugim, jak możemy ustalić, że iteracyjna metoda rozwiązywania dużych układów liniowych jest w praktyce...
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
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 }
Jak rozumiem, sukcesywne nad relaksacją działa poprzez wybranie parametru i użycie liniowej kombinacji (quasi) iteracji Gaussa-Seidela i wartości w poprzednim kroku czasu ... 0≤ω≤20≤ω≤20\leq\omega\leq2 uk+1=(ω)ugsk+1+(1−ω)ukuk+1=(ω)ugsk+1+(1−ω)uk{u}^{k+1} = (\omega){u_{gs}}^{k+1} + (1-\omega)u^{k}...
Mam układ liniowy z macierzą, której wartości własne są równomiernie rozmieszczone na okręgu jednostkowym w następujący sposób: Czy możliwe jest skuteczne rozwiązanie tego rodzaju systemu za pomocą metody iteracyjnej, może z pewnym warunkiem
Zanurzam się w fascynujący świat analizy elementów skończonych i chciałbym rozwiązać duży problem termomechaniczny (tylko mechaniczny efekt termiczny , brak sprzężenia zwrotnego).→→\rightarrow W przypadku problemu mechanicznego zrozumiałem już odpowiedź Geoffa , że muszę użyć iteracyjnego solvera...
Metoda Newtona rozwiązywania równań nieliniowych jest znana, że zbiega się kwadratowo, gdy domysły początkowe są „wystarczająco blisko” do rozwiązania. Co jest „wystarczająco blisko”? Czy istnieje literatura na temat struktury tego basenu
Załóżmy, że mam oryginalny duży, rzadki układ liniowy: . Teraz nie mam ponieważ A jest zbyt duże, aby uwzględnić czynnik lub jakikolwiek rozkład , ale zakładam, że mam rozwiązanie z rozwiązaniem iteracyjnym.ZAx0=b0ZAx0=b0A\textbf{x}_0=\textbf{b}_0ZA- 1ZA-1A^{-1}ZAZAAx0x0\textbf{x}_0 Teraz chcę...