scipy.sparse.linalg.

bicg#

scipy.sparse.linalg.bicg(A, b, x0=None, *, rtol=1e-05, atol=0.0, maxiter=None, M=None, callback=None)[Quelle]#

Löse Ax = b mit der BIConjugate Gradient-Methode.

Parameter:

A{sparse array, ndarray, LinearOperator}: Die reelle oder komplexe N-by-N-Matrix des linearen Systems. Alternativ kann A ein linearer Operator sein, der Ax und A^T x erzeugen kann, z. B. mit scipy.sparse.linalg.LinearOperator.
bndarray: Rechte Seite des linearen Systems. Hat die Form (N,) oder (N,1).
x0ndarray: Startschätzung für die Lösung.
rtol, atolfloat, optional: Parameter für den Konvergenztest. Für Konvergenz sollte norm(b - A @ x) <= max(rtol*norm(b), atol) erfüllt sein. Der Standardwert ist atol=0. und rtol=1e-5.
maxiterinteger: Maximale Anzahl von Iterationen. Die Iteration wird nach maxiter-Schritten beendet, auch wenn die angegebene Toleranz nicht erreicht wurde.
M{sparse array, ndarray, LinearOperator}: Vorkonditionierer für A. Er sollte die Inverse von A approximieren (siehe Hinweise). Effektive Vorkonditionierung verbessert die Konvergenzrate dramatisch, was bedeutet, dass weniger Iterationen benötigt werden, um eine gegebene Fehlertoleranz zu erreichen.
callbackfunction: Vom Benutzer bereitgestellte Funktion, die nach jeder Iteration aufgerufen wird. Sie wird als callback(xk) aufgerufen, wobei xk der aktuelle Lösungsvektor ist.

Rückgabe:

xndarray

Die konvergierte Lösung.

infointeger

Liefert Konvergenzinformationen: 0 : erfolgreicher Abbruch >0 : Konvergenz zur Toleranz nicht erreicht, Anzahl der Iterationen <0 : Parameterzusammenbruch

Hinweise

Der Vorkonditionierer M sollte eine Matrix sein, so dass M @ A eine kleinere Konditionszahl als A hat, siehe [1] .

Referenzen

[1]

„Vorkonditionierer“, Wikipedia, https://de.wikipedia.org/wiki/Vorkonditionierer

[2]

„Bikonjugierter Gradientenverfahren“, Wikipedia, https://de.wikipedia.org/wiki/Bikonjugiertes_Gradientenverfahren

Beispiele

>>> import numpy as np
>>> from scipy.sparse import csc_array
>>> from scipy.sparse.linalg import bicg
>>> A = csc_array([[3, 2, 0], [1, -1, 0], [0, 5, 1.]])
>>> b = np.array([2., 4., -1.])
>>> x, exitCode = bicg(A, b, atol=1e-5)
>>> print(exitCode)  # 0 indicates successful convergence
0
>>> np.allclose(A.dot(x), b)
True