In welchen packages ist bicgstab gut implementiert?

Wenn du dir nicht sicher bist, in welchem der anderen Foren du die Frage stellen sollst, dann bist du hier im Forum für allgemeine Fragen sicher richtig.
Antworten
Anfängerin
User
Beiträge: 5
Registriert: Mittwoch 4. Juli 2012, 17:35

Hallo,

ich möchte bicgstab verwenden um große lineare Gleichungssysteme zu lösen.
Dazu verwende ich bisher bicgstab aus scipy.sparse.linalg.isolve.iterative.
Gibt es da noch eine bessere Alternative?
Denn leider konvergiert das Verfahren nur bei sehr wenigen meiner Matrizen. Auch bei einigen nicht, die eigentlich sehr gut vorkonditioniert sind.
Man muss dazu sagen, dass ich den Vorkonditionierer der Verfahrens selber ausschalten will, um die Vorkonditionierung selbst mit ILU durchzuführen.(from scipy.sparse.linalg.dsolve.linsolve import spilu)

Viele Grüße und vielen Dank im Voraus
die Anfängerin
schaeffkoch
User
Beiträge: 38
Registriert: Dienstag 21. August 2012, 10:59

hallo,

die frage ist doch sehr speziell. eventuell bekommst du in http://dir.gmane.org/gmane.comp.python.scientific.user mehr hilfe.
Antworten