Pytania otagowane jako linear-solver

Nawiązanie do metod rozwiązywania liniowych układów równań.

2
Wydajny warunek wstępny dla Augmented Lagrangian
Chcę rozwiązać nieliniowy problem z nieliniowymi ograniczeniami równości i używam rozszerzonego Lagrangiana z terminem regularnej kary, który, jak wiadomo, psuje liczbę warunków moich zlinearyzowanych układów (przy każdej iteracji Newtona) . Im dłuższy okres kary, tym gorszy numer warunku. Czy ktoś znałby skuteczny sposób na pozbycie się tego złego uwarunkowania w …
2
Czy istnieją wstępne warunki wstępne dla metod wolnych od matrycy?
Metody Jakuba bez Newtona-Kryłowa (JFNK) i ogólnie metody Kryłowa mogą być bardzo przydatne, ponieważ nie wymagają jawnego przechowywania lub budowy matrycy, a jedynie wyniki produktów macierz-wektor. Jeśli faktycznie tworzysz rzadki system, istnieje dla ciebie wiele warunków wstępnych. Co jest dostępne dla prawdziwych metod bez macierzy? Googling pojawia się kilka odniesień …
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ą …
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 …
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 …
5
Najlepszy wybór solvera dla dużego rzadkiego symetrycznego (ale nie pozytywnie określonego) systemu
Obecnie pracuję nad rozwiązaniem bardzo dużych systemów symetrycznych (ale nie pozytywnie określonych), generowanych przez niektóre pewne algorytmy. Te macierze mają niezłą rzadkość blokową, którą można wykorzystać do rozwiązywania równoległego. Ale nie mogę zdecydować, czy powinienem zastosować podejście bezpośrednie (takie jak Multi-frontal) czy iteracyjne (wstępnie uwarunkowane GMRES lub MINRES). Wszystkie moje …
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 …
4
najszybsze rozwiązanie układu liniowego dla małych macierzy kwadratowych (10x10)
Jestem bardzo zainteresowany optymalizacją rozwiązania liniowego rozwiązywania problemów dla małych matryc (10x10), czasami nazywanych drobnymi matrycami. Czy istnieje na to gotowe rozwiązanie? Matryca może być przyjęta jako niejedna. Ten solver ma zostać wykonany ponad 1 000 000 razy w mikrosekundach na procesorze Intel. Mówię o poziomie optymalizacji stosowanym w grach …
4
Szybkie, wyraźne rozwiązanie dla
Szukam szybkiego (ośmielę się powiedzieć optymalnego?) Jawnego liniowego problemu 3x3, , . Ax=bAx=b\mathbf{A}\mathbf{x} = \mathbf{b}A∈R3×3,b∈R3A∈R3×3,b∈R3)\mathbf{A} \in \mathbf{R}^{3 \times 3}, \mathbf{b} \in \mathbf{R}^{3} Macierz jest ogólna, ale zbliżona do macierzy tożsamości z numerem warunku bliskim 1. Ponieważ to tak naprawdę pomiary czujnika z około 5 cyfrową precyzją, nie mam nic przeciwko …
3
Jakieś problemy, w których SOR jest szybszy niż Gaussa-Seidela?
Czy jest jakaś prosta zasada, by powiedzieć, czy warto zrobić SOR zamiast Gaussa-Seidela? (i możliwy sposób oszacowania parametru realxation )ωω\omega Mam na myśli po prostu patrząc na matrycę lub znajomość konkretnego problemu, jaki reprezentuje matryca? Czytałem odpowiedź na te pytania: Czy są jakieś heurystyki dla optymalizacji metody sukcesywnej nadmiernej relaksacji …
1
warunek wstępny dla metody bez macierzy rozwiązania Ax = b
Muszę rozwiązać Ax = b, ale zdaję sobie sprawę, że nawet jeśli jest on rzadki, przechowywanie współczynników macierzy mojego problemu zajmie zbyt dużo pamięci. Zastanawiam się teraz nad użyciem metody bez macierzy, ponieważ te same współczynniki pojawiają się w macierzy dużo czasu, więc mogłem użyć własnego prywatnego schematu przechowywania (i …
3
Zalecenia dotyczące lekkiego / bezinstalacyjnego gęstego solwera algebry liniowej opartego na C lub C ++
Większość mojego programowania to jednorazowe kody badawcze w C na własny użytek. Nigdy nie rozpowszechniałem żadnego kodu poza bliskimi współpracownikami. Opracowałem algorytm, który publikuję w czasopiśmie naukowym. Chcę podać kod źródłowy i być może kod wykonywalny w internetowym suplemencie do artykułu. Kolega poprosił mnie o uogólnienie algorytmu, który wymagał ode …
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.