Gleichungslöser linearer Gleichungssysteme |
05.08.2010, 12:03 | ola | Auf diesen Beitrag antworten » | |||||||
Gleichungslöser linearer Gleichungssysteme Ich bin auf der Suche nach "schnellen" und "nicht aufwändigen" (wenig Speichernutzung) Gleichungslösern für Probleme der Art b=Ax. Es handelt sich hierbei um ein sehr großes Problem der etwaigen Größe n=1.000.000. Mein Interesse liegt ersteinmal darin zu erfahren welche Löser oder auch Libraries es gibt und was die angeknüpften Bedingungen an die Matrix A sind für die jeweiligen Löser, bei welchen Matrizen sie also versagen und bei welchen sie anwendbar sind. Weiter hin ob die Matrix gut Konditioniert sein muss und welcher Rechenaufwand bzw auch Speicheraufwand besteht. Beispiel sind PARDISO, Lapack, Mehrgitter-Verfahren diese habe ich bisher gefunden und auch ein paar Eigenschaften. Ich höre gerne auch nochmal Komentare zu diesen Beispielen.... |
|||||||||
05.08.2010, 12:53 | tigerbine | Auf diesen Beitrag antworten » | |||||||
Ich kann dir leider nur 2 Begriffe hinwerfen: Krylovraum Methoden (https://lp.uni-goettingen.de/get/text/2020) und Vorkonditionierung mit BPX. Vielleicht findest du damit schon was für dich. |
|||||||||
25.08.2010, 00:45 | Formelmonster | Auf diesen Beitrag antworten » | |||||||
Hallo, Zu zweien deiner Beispiele:
Bei den Krylow-Unterraum-Methoden gibt es einen ganzen Haufen von Methoden. Welche man nimmt hängt von den Eigenschaften deiner Matrix ab. Ein paar dieser Methoden sind:
Diese beiden sind neben anderen in MATLAB enthalten. (CG in einer vorkonditionierten Variante mit dem Namen „pcg“). Sie lassen sich sehr einfach benutzen, du musst lediglich deine Matrix A in einer Sparse Matrix speichern und kannst das System dann bspw. mit
lösen. Welche Methode am besten klappt musst Du einfach ausprobieren Es gibt auch noch lineare Iterationsverfahren, die jedoch mit größer werdender Konditionszahl immer langsamer werden. Diese sind zwar langsam, dafür aber leicht zu implementieren. Das Jacobi-Verfahren funktioniert mit irreduzibel diagonaldominanten Matrizen. Das Gauß-Seidel-Verfahren und das SOR-Verfahren funktionieren zusätzlich mit lediglich positiv definiten Matrizen. Grüße und viel Erfolg! |
|
Verwandte Themen
Die Beliebtesten » |
|
Die Größten » |
Die Neuesten » |